Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,939,662 C→T M317I (ATG→ATA hsrA ← putative transporter HsrA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,939,6620CT98.6% 241.3 / NA 70M317I (ATG→ATAhsrAputative transporter HsrA
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (31/38);  minor base G (0/1);  total (31/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.85e-01

GCGTGACCGCCAGAACACTGTTACCGCTGCTGGCATTGTCATCGGTCAGATCGGCAAGTGTGATGGTATTCATCGCGGTAAACTGCGTCGACATAGCCATCCCTAATATAAACAACGGCAAGATCAGCATCCATATAGCCATTGCTGGTGATTGCAAAGAGAACTGAGCGATCATTAGCCCAATAATCACCGTGATCCCCACTAACGTATGGCGATAGCCCAGACGACGTAAGACTTGGGTAACCATCGATTTTGCAATAATGGAACCTAACGCTGTCGGTGCCATCATACAGCCAGCAATAAACGCCTGATAACCAAATCCTACCTGTAACATCAATGGCATAAGGAACGGTACACAGCCGGTCCCCAGACGGGTTGCAATATTGCCTACGATACCGATCGAGAAAGTGCGGGTTTTAAATAAATCTAATGAAATTAATGGGTTTGGCGTGCGTCGTGCATGGAGAATATAGAGAAGCAGTAACC  >  NC_000913/3939418‑3939903
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                 
ccGTGACCGCCAGAACTATCTTACCGCTGCTGGCATTGTTAACAGTGAAATCGGCAAGTGGGATGGTCTTCATCGCGGTAAACTGCGTCGACATAGCCATCCCTAATATAAACAACGGCAAGATCAGCATCCATATAGCCATTGCTGGTGATTGCAAAGAGAACTGAGCGATCATTAGCCCAATAATCACCGTGATCCCCACTATCGTATGGCGATAGCCCAGACGACGTAAGACTTGGGTAACTATCGa                                                                                                                                                                                                                                              <  2:380337/249‑1 (MQ=255)
 cGTGACCGCCAGAACACTGTTACCGCTGCTGGCCTTGTCATCGGTCAGATCGGCAAGTGTGATGGTATTCATCGCGGTAAACTGCGTCGACATAGCCATCCCTAATATAAACAACGGCAAGATCAGCATCCATATAGCCATTGCTGGTGATTGCAAAGAGAACTGAGCGATCATTAGCCCAATAATCACCGTGATCCCCACTAACGTATGGCGATAGCCCAGACGACGTAAGACTTGGGTAACTATCGAtt                                                                                                                                                                                                                                            >  2:84990/1‑251 (MQ=255)
  gTGACCGCCAGAACACTGTTACCGCTGCTGGCATTGTCATCGGTCAGATCGGCAAGTGTGATGGTATTCATCGCGGTAAACTGCGTCGACATAGCCATCCCTAATATAAACAACGGCAAGATCAGCATCCATATAGCCATTGCTGGTGATTGCAAAGAGAACTGAGCGATCATTAGCCCAATAATCACCGTGATCCCCACTAACGTATGGCGATAGCCCAGACGACGTAAGACTTGGGTAACTATCGAttt                                                                                                                                                                                                                                           >  1:97051/1‑251 (MQ=255)
     aCCGCCAGAACACTGTTACCGCTGCTGGCATTGTCATCGGTCAGATCGGCAAGTGTGATGGTATTCATCGCGGTAAACTGCGTCGACATAGCCATCCCTAATATAAACAACGGCAAGATCAGCATCCATATAGCCATTGCTGGTGATTGCAAAGAGAACTGAGCGATCATTAGCCCAATAATCACCGTGATCCCCACTAACGTATGGCGATAGCCCAGACGACGTAAGACTTGGGTAACTATCGATTTTGc                                                                                                                                                                                                                                        <  1:93173/251‑1 (MQ=255)
        gCCAGAACACTGTTACCGCTGCTGGCATTGTCATCGGTCAGATCGGCAAGTGTGATGGTATTCATCGCGGTAAACTGCGTCGACATAGCCATCCCTAATATAAACAACGGCAAGATCAGCATCCATATAGCCATTGCTGGTGATTGCAAAGAGAACTGAGCGATCATTAGCCCAATAATCACCGTGATCCCCACTAACGTATGGCGATAGCCCAGACGACGTAAGACTTGGGTAACTATCGATTTTGCaat                                                                                                                                                                                                                                     >  2:167592/1‑251 (MQ=255)
        gCCAGAACACTGTTACCGCTGCTGGCATTGTCATCGGTCAGATCGGCAAGTGTGATGGTATTCATCGCGGTAAACTGCGTCGACATAGCCATCCCTAATATAAACAACGGCAAGATCAGCATCCATATAGCCATTGCTGGTGATTGCAAAGAGAACTGAGCGATCATTAGCCCAATAATCACCGTGATCCCCACTAACGTATGGCGATAGCCCAGACGACGTAAGAATTGGGTAACTATCGATTTTGCaa                                                                                                                                                                                                                                      >  2:228577/1‑250 (MQ=255)
                   gTTACCGCTGCTGGCATTGTCATCGGTCAGATCGGCAAGTGTGATGGTATTCATCGCGGTAAACTGCGTCGACATAGCCATCCCTAATATAAACAACGGCAAGATCAGCATCCATATAGCCATTGCTGGTGATTGCAAAGAGAACTGAGCGATCATTAGCCCAATAATCACCGTGATCCCCACTAACGTATGGCGATAGCCCAGACGACGTAAGACTTGGGTAACTATCGATTTTGCAATAATGGAACCTa                                                                                                                                                                                                                          <  2:357929/251‑1 (MQ=255)
                    ttACCGCTGCTGGCATTGTCATCGGTCAGATCGGCAAGTGTGATGGTATTCATCGCGGTAAACTGCGTCGACATAGCCATCCCTAATATAAACAACGGCAAGATCAGCATCCATATAGCCATTGCTGGTGATTGCAAAGAGAACTGAGCGATCATTAGCCCAATAATCACCGTGATCCCCACTAACGTATGGCGATAGCCCAGACGACGTAAGACTTGGGTAACTATCGATTTTGCAATAATGGAAACTa                                                                                                                                                                                                                          >  2:327024/1‑250 (MQ=255)
                     tACCGCTGCTGGCATTGTCATCGGTCAGATCGGCAAGTGTGATGGTATTCATCGCGGTAAACTGCGTCGACATAGCCATCCCTAATATAAACAACGGCAAGATCAGCATCCATATAGCCATTGCTGGTGATTGCAAAGAGAACTGAGCGATCATTAGCCCAATAATCACCGTGATCCCCACTAACGTATGGCGATAGCCCAGACGACGTAAGACTTGGGTAACTATCGATTTTGCAATAATGGAACCTAAc                                                                                                                                                                                                                        >  2:184345/1‑251 (MQ=255)
                      aCCGCTGCTGGCATTGTCATCGGTCAGATCGGCAAGTGTGATGGTATTCATCGCGGTAAACTGCGTCGACATAGCCATCCCTAATATAAACAACGGCAAGATCAGCATCCATATAGCCATTGCTGGTGATTGCAAAGAGAACTGAGCGATCATTAGCCCAATAATCACCGTGATCCCCACTAACGGATGGCGATAGCCCAGACGACGTAAGACTTGGGTAACTATCGATTTTGCAATAATGGAACCTAAc                                                                                                                                                                                                                        <  2:211689/250‑1 (MQ=255)
                               ggCATTGTCATCGGTCAGATCGGCAAGTGTGATGGTATTCATCGCGGTAAACTGCGTCGACATAGCCATCCCTAATATAAACAACGGCAAGATCAGCATCCATATAGCCATTGCTGGTGATTGCAAAGAGAACTGAGCGAGCATTAGCCCAATAATCACCGTGATCCCCACTAACGTAGGGCGATAGCCCAGACGACGTAAGACTTGGGTGACGATCGATTTTGCAATAATGGAACCTAACGCTGTCGGt                                                                                                                                                                                                               <  1:259356/250‑1 (MQ=255)
                               ggCATCGTCATCGGTCAGACCGGCAAGTGTGATGGTATTCATCGCGGTAAACTGCGTCGACATAGCAATCCCTACTATAAACAACGGCAAGATCAGCATCCATATAGCCATTGCTGGTGATTGCAAAGAGAACTGAGCGATCATTAGCCCAATAATCACCGTGATCCCCACTAACGTATGGTGATAGCCCAGACGACGTCAGACTTGGGTAACTATCGATTTTGCAATAATGGAACCTAACGCTGTCgggg                                                                                                                                                                                                              >  2:181119/1‑249 (MQ=255)
                                gCATTGTCATCGGTCAGATCGGCAAGTGTTATGGTATTCAGCGCGGTAAACTGCGTCGACATAGCCAGCCCTAATAGAAACAACGGCAAGATCAGCATCCATATAGCCATTGCTGGTGATTGCATAGAGAACTGAGGGATCTTTAGCCGAATAGTGGCCGTGAGCCCCATTAACGTATGGCGATAGCCTAGACGACGTAATTCTTGGGTAAGTATCGATTTTGCAATAATGGAACCTAACGCTGTCGGTGc                                                                                                                                                                                                             <  1:173483/251‑1 (MQ=255)
                                 cATTGTCATCGGTCAGATCGGCAAGTGTGATGGTATTCATCGCGGTAAACTGCGTCGACATAGCCATCCCTAATATAAACAACGGCAAGATCAGCATCCATATAGCCATTGCTGGTGATTGCAAAGAGAACTGAGCGATCATTAGCCCAATAATCACCGTGATCCCCACTAACGTATGGCGATAGCCCAGACGACGTAAGACTTGGGTAACTATCGATTTTGCAATAATGGAACCTAACGCTGTCGGTGcc                                                                                                                                                                                                            <  1:69307/251‑1 (MQ=255)
                                    tGTCATCGGTCAGATCGGCAAGTGTGATGGTATTCATCGCGGTAAACTGCGCCGACATAGCCATCCCTACTATAAACAACGGCAAGATCAGCATCCATATAGCCATTGCTGGTGATTGCAAAGAGAACTGAGCGATCATTAGCCCAATAATCACCGTGATCCCCACTAACGTATGGCGATAGCCCAGACGACGTAAGACTTGGGAAACTATCGATTTTGCAATAATGGACCCTAACGCTGTCGGTGCcatc                                                                                                                                                                                                         >  2:19993/1‑251 (MQ=255)
                                     gTCATCGGTCAGATCGGCAAGTGTTATGGTATTCATCGCGGTAAACTGCGTCGACATAGCCAGCCCTAATATAAACAACGGCAAGATCAGCATCCATATAGCCATTGCTGGTGATTGCAAAGAGAACTGAGCGATCATTAGCCCAATAAGCACCGTGATCCCCACTAACGTATGGCGATAGCGCAGACGACGGAAGACGTGGGTAACTATCGAGGTGTCAATAATGGAAGCTAACGCGGTCTGTGCcatca                                                                                                                                                                                                        <  2:401187/251‑1 (MQ=255)
                                        ctCGGTCAGATCGGCAAGTGTGATGGTATTCATCGCGGTAAACTGCGTCGACATAGCCATCCCTAATATAAACAACGGCAAGATCAGCATCCATATAGCCATTGCTGGTGATTGCAAAGAGAACTGAGCGATCATTAGCCCAATAATCACCGTGATCCCCACTAACGTATGGCGAGAGCCCAGACGACGTAAGACTTGGGTAACTATCGATTTTGCAATAATGGAACCTAACGCTGTCGGTGCCATCATa                                                                                                                                                                                                      >  2:92806/2‑250 (MQ=255)
                                                   cGGCAAGTGTGATGGTATTCATCGCGGTAAACTGCGTCGACATAGCCATCCCTAATATAAACAACGGCAAGATCAGCATCCATATAGCCATTGCTGGTGATTGCAAAGAGAACTGAGCGATCATTAGCCCAATAATCACCGTGATCCCCACTAACGTATGGCGATAGCCCAGACGACGTAAGACTTGGGTAACTATCGATTTTGCAATAATGGAACCTAACGCTGTCGGTGCCATCATACAGCCAGCAAt                                                                                                                                                                                           <  1:327819/250‑1 (MQ=255)
                                                      cAAGTGTGATGGTATTCATCGCGGTAAACTGCGTCGACATAGCCATCCCTAATATAAACAACGGCAAGATCAGCATCCATATAGCCATTGCTGGTGATTGCAAAGAGAACTGAGCGATCATTAGCCCAATAATCACCGTGATCCCCACTAACGTATGGCGATAGCCCAGACGACGTAAGACTTGGGTAACTATCGATTTTGCAATAATGGAACCTAACGCTGTCGGTGCCATCATACAGCCAGCAATAAAc                                                                                                                                                                                       <  1:112988/251‑1 (MQ=255)
                                                             gATGGTATTCATCGCGGTAAACTGCGTCGACATAGCCATCCCTAATATAAACAACGGCAAGATCAGCATCCATATAGCCATTGCTGGTGATTGCAAAGAGAACTGAGCGATCATTAGCCCAATAATCACCGTGATCCCCACTAACGTATGTCGATAGCCCAGACGACGTGAGACTTGGGTAACTAGCGAGTTTGCAATAATGGAATCTATCGCTTTCGGTGCCATCATACAGCCAGCAATAAACGCCTGAt                                                                                                                                                                                <  1:196080/251‑1 (MQ=255)
                                                             gATGGTATTCATCGCGGTAAACTGCGTCGACATAGCCATCACTAATATAACCAACGGCAAGATCAGCATCCATATAGCCATTGCTGGTGATTGCAAAGAGAACTGAGCGATCATTAGCCCAATAATCACCGTGATCCCCACTAACGTATGGCGATAGCCCAGACGACGTAAGACTTGGGTAACTATCGATTTTGCAATAATGGAACCTAACGCTGTCGGTGCCATCATACAGCCAGCAATAAACGCCTgtt                                                                                                                                                                                >  1:412121/1‑249 (MQ=255)
                                                                       aTCGCGGTAAACTGCGTCGACATAGCCATCCCTAATATAAACCACGGCAAGATCAGCATCCATATAGCCATTGCTGGTGATTGCAAAGAGAACTGAGCGATCCTTAGCCCAATAATCACCGTGATCCCCACTAACGTATGGCGATAGCCCAGACGACGTAAGACTTGGGTAACTATCGATTTTGCAATAATGGAACCTAACGCTGTCGGTGCCATCATACAGCCAGCAATAAACGCCTGATAACCAAATcc                                                                                                                                                                      >  2:202061/1‑251 (MQ=255)
                                                                                 aCTGCGTCGACATAGCCATCCCTAATATAAACAACGGCAAGATCAGCATCCATATAGCCATTGCTGGTGATTGCAAAGAGAACTGAGCGATCATTAGCCCAATAATCACCGTGATCCCCACTAACGTATGGCGATAGCCCAGACGACGTAAGACTTGGGTAACTATCGATTTTGCAATAATGGAACCTAACGCTGTCGGTGCCATCATACAGCCAGCAATAAACGCCTGATAACCAAATCCTACCTGTAAc                                                                                                                                                            >  2:145691/1‑251 (MQ=255)
                                                                                 aCTGCGACGACATAGACATCCCCCATATAACCAAAGGCAAGCTCAAGATCCATATAGACATTGCTGGAGATAGTAAAGAGAACTGAGGGATCATTAGTTCAATAATCACCGTGCTACCACCTACCGTATGGCGCTAGCCCAGACGTCGTCAGACTTGGGTAACTATCGATTTTGCACTAATGGAACCTAACGCTGACGGTGCCATCATACAGCCAGCAATAAACGCCTGCTAACCACATCATACCTGTAAc                                                                                                                                                            >  2:383555/1‑251 (MQ=255)
                                                                                  cTGCGTCGACATAGCCATCCCTAATATAAACAACGGCAAGATCAGCATCCATATAGCCATTGCTGGTGATTGCAAAGAGAACTGAGCGATCATTAGCCCAATAATCACCGTGATCCCCACTAACGTATGGCGATAGCCCAGACGACGTAAGACTTGGGTAACTATCGATTTTGCAATAATGGAACCTAACGCTGTCGGTGCCATCATACAGCCAGCAATAAACGCCTGATAACCAAATCCTACCTGtcac                                                                                                                                                            >  2:118608/1‑247 (MQ=255)
                                                                                                           ataAACAACGGCAAGATCAGCATCCATATAGCCATTGCTGGTGATTGCAAAGAGAACTGAGCGATCATTAGCCCAATAATCACCGTGATCCCCACTAACGTATGGCGATAGCCCAGACGACGTAAGACTTGGGTAACTATCGATTTTGCAATAATGGAACCTAACGCTGTCGGTGCCATCATACAGCCAGCAATAAACGCCTGATAACCAAATCCTACCTGTAACATCAATGGCATAAGGAACGGTacac                                                                                                                                   >  1:352605/1‑250 (MQ=255)
                                                                                                            taAACAACGGCAAGATCAGCATCCATATAGCCATTGCTGGTGATTGCAAAGAGAACTGAGCGATCATTAGCCCAATAATCACCGTGATCCCCACTAACGTATGGCGATAGCCCAGACGACGTAAGACTTGGGTAACTATCGATTTTGCAATAATGGAACCTAACGCTGTCGGTGCCATCATACAGCCAGCAATAAACGCCTGATAACCAAATCCTACCTGt                                                                                                                                                               <  1:347805/221‑1 (MQ=255)
                                                                                                            taAACAACGGCAAGATCAGCATCCATATAGCCATTGATGGTGATTGCAACGAGAACTGAGCGATCATTAGCCCAATAATCACCGTGATCCCCACTAACGTATGGCGATAGCCCAGACGACGTAAGACTTGGGTAACTATCGATTTTGCAATAATGGAACCTAACGCTGTCGGTGCCATCATACAGCCAGCAATAAACGCCTGATAACCAAATCCTACCTGt                                                                                                                                                               >  2:347805/1‑221 (MQ=255)
                                                                                                              aacaacGGCAAGATCAGCATCCATATAGCCATTGCTGGTGATTGCAAAGAGAACTGAGCGATCATTAGCCCAATAATCACCGTGATCCCCACTAACGTATGGCGATAGCCCAGACGACGTAAGACTTGGGTAACTATCGATTTTGCAATAATGGAACCTAACGCTGTCGGTGCCATCATACAGCCAGCAATAAACGCCTGATAACCAAATCCTACCTGTAACATCAATGGCATAAGGAACGGTACACAGcc                                                                                                                               >  2:171454/1‑251 (MQ=255)
                                                                                                                 aacGGCAAGATCAGCATCCATATAGCCATTGCTGGTGATTGCAAAGAGAACTGAGCGATCATTAGCCCAATAATCACCGTGATCCCCACTAACGTATGGCGATAGCCCAGACGACGTAAGACTTGGGTAACTATCGATTTTGCAATAATGGAACCTAACGCTGTCGGTGCCATCATACAGCCAGCAATAAACGCCTGATAACCAAATCCTACCTGTAACATCAATGGCATAAGGAACGGTACACAGCCGGt                                                                                                                            <  1:352502/251‑1 (MQ=255)
                                                                                                                    ggCAAGATCAGCATCCATATAGCCATTGCTGGTGATTGCAAAGAGAACTGAGCGATCATTAGCCCAATAATCACCGTGATCCCCACTAACGTATGGCTATAGCCCAGACGACGTAAGACTTGGGTAACTATCGATTTTGCAATAATGGAACCTAACGCTGTc                                                                                                                                                                                                                  >  2:114832/1‑162 (MQ=255)
                                                                                                                    ggCAAGATCAGCATCCATATAGCCATTGCTGGTGATTGCAAAGAGAACTGAGCGATCATTAGCCCAATAATCACCGTGATCCCCACTAACGTATGGCTATAGCCCAGACGACGTAAGACTTGGGTAACTATCGATTTTGCAATAATGGAACCTAACGCTGTc                                                                                                                                                                                                                  <  1:114832/162‑1 (MQ=255)
                                                                                                                      cAAGATCAGCATCCATATAGCCATTGCTGGTGATTGCAAAGAGAACTGAGCGATCATTAGCCCAATAATCACCGTGATCCCCACTAACGTATGGCGATAGCCCAGACGACGTAAGACTTGGGTAACTATCGATTTTGCAATAATGGAACCTAACGCTGTCGGTGCCATCATACAGCCAGCAATAAACGCCTGATAACCAAATCCTACCTGTAACATCAATGGCATAAGGAACGGTACACAGCCGGTCCCCa                                                                                                                       <  1:66084/251‑1 (MQ=255)
                                                                                                                       aaGATCAGCATCCATATAGCCATTGCTGGTGATTGCAAAGAGAACTGAGCGATCATTAGCCCAATAATCACCGTGATCCCCACTAACGTATGGCGATAGCCCAGACGACGTAAGACTTGGGTAACTATCGATTTTGCAATAATGGAACCTAACGCTGTCGGTGCCATCATACAGCCAGCAATAAACGCCTGATAACCAAATCCTACCTGTAACATCAATGGCATAAGGAACGGTACACAGCCGGTcccc                                                                                                                        >  1:83501/1‑249 (MQ=255)
                                                                                                                         gATCAGCATCCATATAGCCATTGCTGGTGATTGCAAAGAGAACTGAGCGATCATTAGCCCAATAATCACCGTGATCCCCACTAACGTATGGCGATAGCCCAGACGACGTAAGACTTGGGTAACTATCGATTTTGCAATAATGGAACCTAACGCTGTCGGTGCCATCATACAGCCAGCAATAAACGCCTGATAACCAAATCCTACCTGTAACATCAATGGCATAAGGAACGGTACACAGCCGGTCCCCAGAc                                                                                                                    <  1:19993/251‑1 (MQ=255)
                                                                                                                             aGCATCCATATAGCCATTGCTGGTGATTGCAAAGAGAACTGAGCGATCATTAGCCCAATAATCACCGTGATCCCAACTAACGTATGGCGATAGCCCAGACGACGTAAGACTTGGGTAACTATCGATTTTGCAATAATGGAACCTAACGCTGTCGGTGCCATCATACAGCCAGCAATAAACGCCTGATAACCAAATCCTACCTGTAACATCAATGGCATAAGGAACGGTACACAGCCGGTCCCCAGACGGGt                                                                                                                >  2:302079/1‑251 (MQ=255)
                                                                                                                              gCATCCATATAGCCATTGCTGGTGATTGCAAAGAGAACTGAGCGATCATTAGCCCAATAAGCACCGTGATCCCCACTAACGTATGGCGATAGCCCAGACGACGTAAGACTTGGGTAACTATCGATTTTGCAATAATGGAACCTAACGCTGTCGGTGCCATCATACAGCCAGCAATAAACGCCTGATACCCAAATCCTACCTGTAACATCAATGGCATAAGGAACGGTACACAGCCGGTCCCCAGACGGGtt                                                                                                               <  2:306884/251‑1 (MQ=255)
                                                                                                                                 tCCATATAGCCATTGCTGGTGATTGCAAAGAGAACTGAGCGATCATTAGCCCAATAATCACCGTGATCCCCACTAACGTATGGCGATAGCCCAGACGACGTAAGACTTGGGTAACTATCGATTTTGCAATAATGGAACCTAACGCTGTCGGTGCCATCATACAGCCAGCAATAAACGCCTGATAACCAAATCCTACCTGTAACATCAATGGCATAAGGAACGGTACACAGCCGGTCCCCAGACGGGTTGc                                                                                                             >  2:315669/1‑250 (MQ=255)
                                                                                                                                      ataGCCATTGCTGGTGATTGCAAAGAGAACTGAGCGATCATTAGCCCAATAATCACCGTGATCCCCACTAACGTATGGCGATAGCCCAGACGACGTAAGACTTGGGTAACTATCGATTTTGCAATAATGGAACCTAACGCTGTCGGTGCCATCATACAGCCAGCAATAAACGCCTGATAACCAAATCCTACCTGTAACATCAATGGCATAAGGAACGGTACACAGCCGGTCCCCAGACGGGTTGCAATAtt                                                                                                       >  1:128865/1‑251 (MQ=255)
                                                                                                                                            aTTGCTGGTGATTGCAAAGAGAACTGAGCGATCATTAGCCCAATAATCACCGTGATCCCCACTAACGTATGGCGATAGCCCAGACGACGTAAGACTTGGGTAACTATCGATTTTGCAATAATGGAACCTAACGCTGTCGGTGCCATCATGCAGCCAGCAATAAACGCCTGATAACCAAATCCTACCTGTAACATCAATGGCAGAATGAACGGTACACAGCCGGTCCCCAGACGGGTTGCAATATTGCCTAc                                                                                                 <  1:397426/251‑1 (MQ=255)
                                                                                                                                                  ggTGATTGCAAAGAGAACTGAGCGATCATTAGCCCAATAATCACCGTGATCCCCACTAACGTATGGCGATAGCCCAGACGACGTAAGACTTGGGTAACTATCGATTTTGCAATAATGGAACCTAACGCTGTCGGTGCCATCATACAGCCAGCAATAAACGCCTGATAACCAAATCCTACCTGTAACATCAATGGCATAAGGAACGGTACACAGCCGGTCCCCAGACGGGTTGCAATATTGCCTACGATAc                                                                                            <  1:228577/250‑1 (MQ=255)
                                                                                                                                                  ggTGATTGAAAAGAGAACTGAGAGATCATTAGCCCAATAAACACCGTGATCCCCACTAACGTATGGCGATAGCCCAGACGACGTAAGACTGGGGTAACTATCGAGTTTGCAATAATGGAACCTAACGCTGTCGGGGCCATCATACAGCCAGCAATAAACGCCTGATAACCAAATCCTACCTGTAACATCAATGGCATAAGGAACGGTACACAGCCGGTCCCCAGACGGGTGGCAATATTGCCTACGATAcc                                                                                           <  2:3201/251‑1 (MQ=255)
                                                                                                                                                   gTGATTGCAAAGAGAACTGAGCGATCATTAGCCCAATAATCACCGTGATCCCCCCTAACGTATGGCGACAGCCCAGACGAAGTAAGCCTTGGGTAACTATCGATTTTGCAATAATGGAACCTAACGATGTCGGAGCCATCATACAGCCAGCAATAAACGCCTGATAACCAAATCCTACCTGTAACATCAATGGCATAAGGAACGGTACACAGCCGGTCCCCAGACGGGTTGCAATATTGCCTACGATACCg                                                                                          >  1:202855/1‑251 (MQ=255)
                                                                                                                                                      aTTGCAAAGAGAACTGAGCGATCATTAGCCCAATAATCACCGTGATCCCCACTAACGTATGGCGATAGCCCAGACGACGTAAGACTTGGGTAACTATCGATTTTGCAATAATGGAACCTAACGCTGTCGGTGCCATCATACAGCCAGCAATAAACGCCTGATAACCAAATCCTACCTGTAACATCAATGGCATAAGGAACGGTACACAGCCGGTCCCCAGACGGGTTGCAATATTGCCTACGATACCGATc                                                                                       <  1:21921/251‑1 (MQ=255)
                                                                                                                                                      aTTGCAAAAAGACATGAGCGATCATTAGCCCAATAATCACCGTTCTCCCCACTAACGTATGGCTATAGCCCAGACGACGTAAGACTTGGGTAACTATCGATTTTGCAATAATTGAACCTACCGCTGTCGGTGCCACCATACAGCCAGCAATAAACGCCTGTTAACCAAATCCTACCGGTAACATCAATGGCATAAGGAACGGTACACAGCCGGTCCCCAGTCGTGTGGCAATATTGCCTACGATACCGATc                                                                                       <  2:163957/251‑1 (MQ=255)
                                                                                                                                                        tGCAAAGAGAACTGAGCGATCATTAGCCCAATAATCACCGTGATCCCCACTAACGTATGGCGATAGCCCAGACGACGTAAGACTTGGGTAACTATCGATTTTGCAATAATGGAACCTAACGCTGTCGGTGCCATCATACAGCCAGCAATAAACGCCTGATAACCAAATCCTACCTGTAACATCAATGGCATAAGGAACGGTACACAGCCGGTCCCCAGACGGGTTGCAATATTGCCTACGATACCGATCg                                                                                      <  1:145691/250‑1 (MQ=255)
                                                                                                                                                                     gAGCGATCATTAGCCCAATAATCACCGTGATCCCCACTAACGTATGGCGATAGCCCAGACGACGTAAGACTTGGGTAACTATCGATTTTGCAATAATGGAACCTAACGCTGTCGGTGCCATCATACAGCCAGCAATAAACGCCTTATAACCAAATCCTACCTGTAACATCAATGGCATAAGGAACGGTACACAGCCGGTCCCCAGTCGGGTTGCAATATTGCCTACGATACCGATCGAGAAAGTGCggg                                                                          <  1:34052/249‑1 (MQ=255)
                                                                                                                                                                     gAGCGATCATTAGCCCAATAATCACCGTGATCCCCACTAACGTATGGCGATAGCCCAGACGACGTAAGACTTGGGTAACTATCGATTTTGCAATAATGGAACCTAACGCTGTCGGTGCCATCATACAGCCAGCAATAAACGCCTGATAACCAAATCCTACCTGTAACATCAATGGCATAAGGAACGGTACACAGCCGGTCCCCAGACGGGTTGCAATATTGCCTACGATACCGATCGAGAAAGTGCGGGt                                                                         <  1:171454/250‑1 (MQ=255)
                                                                                                                                                                         gATCATTAGCCCAATAATCACCGTGATCCCCACTAACGTATGGCGATAGCCCAGACGACGTAAGACTTGCGTAACTATCGATTTTGCAATAATGGAACCTAACGCTGTCGGTGCCATCATACAGCCAGCAATAAACGCCTGATAACCAAATCCTACCTGTAACATCAATGGCATAAGGAACGGTACACAGCCGGTCCCCAGACGGGTTGCAATATTGCCTACGATACCGATCGAGAAAGTGCGGGTTTtaa                                                                    <  1:84990/251‑1 (MQ=255)
                                                                                                                                                                               tAGCCCAATAATCACCGTGATCCCCACTAACGTATGGCGATAGCCCAGACGACGTAAGACTTGGGTAACTATCGATTTTGCAATAATGGAACCTAACGCTGTCGGTGCCATCATACAGCCAGCAATAAACGCCTGATAACCAAATCCTACCTGTAACATCAATGGCATAAGGAACGGTACACAGCCGGTCCCCAGACGGGTTGCAATATTGCCTACGATACCGATCGAGAAAGTGCGGGTTTTaaataaat                                                              >  1:132176/1‑251 (MQ=255)
                                                                                                                                                                                aGCCCAATAATCCCCGTGATCCCCACTAACGTATGGCGATAGCCCAGACGCCGTAAGACTTGGGTATCTATCGATTTTGCAATAATGGAACCTAACGCTGTCTGTGCCACTATACAGCCAGCAATAAACGCCTGAAATCCAGATCCTACCGTTATCATCAATGGCATAAGGAACTGTCCACAGCCGGTCCCCAGGCGTGTTTCAATAGTGCCTACGAGACCGATCTAGAAAGTGCGGGCGTTAAATAAATc                                                             <  2:261113/251‑1 (MQ=255)
                                                                                                                                                                                        aatCACCGTGATCCCCACTAACGTATGGCGATAGCCCAGACGACGTAAGACTTGGGTAACTATCGATTTTGCAATAATGGAACCTAACGCTGTCGGTGCCATCATACAGCCAGCAATAAACGCCTGATAACCAAATCCTACCTGTAACATCAATGGCATAAGGAACGGTACACAGCCGGTCCCCAGACGGGTTGCAATATTGCCTACGATACCGATCGAGAAAGAGAGGGTTTTAAATAAATCTAATGaa                                                      >  2:140480/1‑250 (MQ=255)
                                                                                                                                                                                            aCCGTGATCCCCACTAACGTATGGCGATAGCCCAGACGACGTAAGACTTGGGTAACTATCGATTTTGCAATAATGGAACCTAACGCTGTCGGTGCCATCATACAGCCAGCAATAAACGCCTGATAACCAAATCCTACCTGTAACATCAATGGCATAAGGAACGGTACACAGCCGGTCCCCAGACGGGTTGCAATATTGCCTACGATACCGATCGAGAAAGTGCGGGTTTTAAATAAATCTAATGAAATTa                                                  >  1:361410/1‑250 (MQ=255)
                                                                                                                                                                                              cGTGATCCCCACTAACGTATGGCGATAGCCCAGACGACGTAAGACTTGGGTAACTATCGATTTTGCAATAATGGAACCTAACGCTGTCGGTGCCATCATACAGCCAGCAATAAACGCCTGATAACCAAATCCTACCTGTAACATCAATGGCATAAGGAACGGTACACAGCCGGTCCCCAGACGGGTTGCAATATTGCCTACGATACCGATCGAGAAAGTGCGGGTTTTAAATAAATCTAATGAAATTAATg                                               <  1:302079/251‑1 (MQ=255)
                                                                                                                                                                                               gTGATCCCCACTAACGTATGGCGATAGCCCAGACGACGTAAGACTTGGGTAACTATCGATTTTGCAATAATGGAACCTAACGCTGTCGGTGCCATCATACAGCCAGCAATAAACGCCTGATAACCAAATCCTACCTGTAACATCAATGGCATAAGGAACGGTACACAGCCGGTCCCCAGACGGGTTGCAATATTGCCTACGATACCGATCGAGAAAGTGCGGGTTTTAAATAAATCGAATGAAATTAATgg                                              <  1:126622/251‑1 (MQ=255)
                                                                                                                                                                                               gTGATCCACACTAACGTATGGCGATAGCCCAGACGACGTAAGACTTGGGTAACTATCGATTTTGCAATACTGGAACCTAACGCTGTCGGTGCCATCATACAGCCAGCAATAAACGCCTGATAACCAAATCCTACCTGTAACATCAATGGCATAAGGAACGGTACACAGCCGGTCCCCAGACGGGTTGCAATATTGCCTACGATACCGATCGAGAAAGTGTGGGTTTTAAATAAATCTAATGAAATTAATgg                                              >  1:63871/1‑251 (MQ=255)
                                                                                                                                                                                                    ccccACTAACGTATGGCGATAGCCCAGACGACGTAAGACTTGGGTAACTATCGATTTTGCAATAATGGAACCTAACGCTGTCGGTGCCATCATACAGCCAGCAATAAACGCCTGATAACCAAATCCTACCTGTAACATCAATGGCATAAGGAACGGTACACAGCCGTTCCCCAGACGGGTTGCAATATGGCCTACGATACCGATCGAGAAAGTGCGGGTTTTAAATAAATCTAATGAAATTAATGGGTTTg                                         <  1:202061/251‑1 (MQ=255)
                                                                                                                                                                                                               taaggCGATAGCCCAGACGACGTAAGACTTGGGTAACTATCGATTTTGCAAGAATGGAACCTAACGCTGTCGGTGCCATCATACAGCCAGCAATAAACGCCTGATAACCAAATCCTACCGGTAACATCAATGGCATAAGGAACGGTACACAGCCGGTCCCCAGACGGGTTGCAATATTGCCTACGATACCGATCGAGAAAGTGCGGGTTTTAAATAAATCTAATGAAAGTAATGGGTTTGGCGTGcgtcg                               <  2:391953/247‑1 (MQ=255)
                                                                                                                                                                                                                aTGGCGATAGCCCCGACGACTTAAGACTTGGGTAACTATCGATTTTGCAATAATGAAACCTAACGCTGTCGGTGCCATCATACAGCCAGCAATAAACGCCTGATAACCAAATCCTACCTGTAACATCAATGGCATAAGGAACGGTACACAGCCGGTCCCCAGACGGGTTGCAATATTGCCTACGATACCGATCGAGAAAGTGCGGGTTTTAAATAAATCTAATGAAATTAATGTGTTTGGCGTGcgtcgt                              <  2:426455/250‑1 (MQ=255)
                                                                                                                                                                                                                     gATAGCCTAGACGACGTAAGACCTGGGAAACTATCGATTTTGCAATAATGGTACCTAAAGATGTGGGTGCCATCATACACCCAGCAATAAACGCCTGATAACCCAATCCTACCTGTAACATCAATGGCATAAGGCACGGCACAAAGCCGGTCCCCAGACCGGTAGCAATATAGCCTACGCTACCGAACGAGAAAGCGGGCGATTTAAATAAATCTCATGAAATTAATGGGTTTGGCGAGCGTCGTGCATgg                        >  2:77741/1‑251 (MQ=255)
                                                                                                                                                                                                                      aTAGCCCAGACGACGTAAGACTTGGGTAACTATCGATTTTGCACTAATGGAACCTAACGCTGTCGGTGCCATCATACAGCCAGCAATAAACGCCTGATAACCAAATCCTACCTGTAACATCAATGGCATAAGGAACGGTACACAGCCGGTCCCCAGACGGGTTGCAATATTGCCTACGATACCGATCGAGAAAGTGCGGGTTTTAAATAAATCTAATGAAATTAATGGGTTTGGCGTGCGTCGTGCATggg                       >  1:366155/1‑250 (MQ=255)
                                                                                                                                                                                                                            cAGACGACGTAAGACTTGGGTAACTATCGATTTTGCAATAATGGAACCTAACGCTGTCGGTGCCATCATACAGCCAGCAATAAACGCCTGATAACCAAATCCTACCTGTAACATCAATGGCATAAGGAACGGTACACAGCCGGTCCCCAGACGGGTTGCAATATTGCCTACGATACCGATCGAGAAAGTGCGGGTTTTAAATAAATCTAATGAAATTAATGGGTTTGGCGTGCGTCGTGCATGGAGAatat                 <  1:327024/251‑1 (MQ=255)
                                                                                                                                                                                                                                 gacgTAAGACTTGGGTAACTATCGATTTTGCAATAATGGAACCTAACGCTGTCGGTGCCATCATACAGCCAGCAATAAACGCCTGATAACCAAATCCTACCTGTAACATCAATGGCATAAGGAACGGTACACAGCCGGTCCCCAGACGGGTTGCAATATTGCCTACGATACCGATCGAGAAAGTGCGGGTTTTAAATAAATCTAATGAAATTAATGGGTTTGGCGTGCGTCGTGCATGGAGAATATag               >  2:71466/1‑248 (MQ=255)
                                                                                                                                                                                                                                 gacgTAAGACTTGGGTAACTATCGATTTTGCAATAATGGAACCTAACGCTGTCGGTGCCATCATACAGCCAGCAATAAACGCCTGATAACCAAATCCTACCTGTAACATCAATGGCATAAGGAACGGTACACAGCCGGTCCCCAGACGGGTTGCAATATTGCCTACGATACCGATCGAGAAAGTGCGGGTTTTAAATAAATCTAATGAAATTAATGGGTTTGGCGTGCGTCGTGCATGGAGAATATAgaga            <  2:83501/251‑1 (MQ=255)
                                                                                                                                                                                                                                  acgTAAGACTTGGGTAACTATCGATTTTGCAATAATGGAACCAAACGCTGTCGGTGCCATCATACAGCCAGCAATAAACGCCTAATAACCGAATACTACCTGTACCATCAATGGCATAAGGAACGGTACACAGCCGGTCCCCAGACGGGTTGCAATATTACCTACGCTACCGATCGAGAAAGTGCGGGTTTTAAATAAATCTAATGAAATTAATGGGTTTGGCGTGCGTCGTGCATGGCGAATATAgaga            >  1:171916/1‑250 (MQ=255)
                                                                                                                                                                                                                                   ggTAAGACTTGGGTAACTATCGATGTTACAATAATGGAACCTAACGCTGTCGGTGCCATCATACAGCCAGCAATAAACGCCTGATAACCAAATCCTACCTGTAACATCAATGGCATAAGGAACGGTACACAGCCGGTCCCCAGACGGGTTGCAATATTGCCTACGATACCGATCGAGAAAGTGCGGGTTTTAAATAAATCTAATGAAATTAATGGGTTTGGCGTGCGTCGTGCATGGAGAATATAGAGAAg          <  1:63107/250‑1 (MQ=255)
                                                                                                                                                                                                                                      aaGAGTTGGGTAACTATCTATTTTGCAATAATGGAACCTAACGCTGGCGGTGGCATCATATAGCCAGCAATAAGCGCCGGTTAACCAAATCCTACCTGTAACATCGATGTCATAAGGGCCGGTACACAGGCGGTTCCCGGACGGGGGGCAATGTTGGCTACGATATCGATCGAGAAGGTGCGGGTTGTTAATAAATCTTATGGAGTTATTTGATTTGGCGTGCGTCGTGCATGGAGAATATAGAGAAGCa        <  1:257338/250‑1 (MQ=255)
                                                                                                                                                                                                                                        gccTTGGGTAACTTTCCACTTTGCAATAATGGCACCTAACGCTGTTGGTGCCCTCATACAGCCAGCAACAACCCCCTGATATCCAACTCCTCCCTGTACCACCAATGGCATAAGTCCCGGTACACAGCCGGTTCCCAGACGGGTTGCAATACTGCCTACTATACCGATCCAGAAAGTGCGGTTTTCAAATAAACCTAATGAAATTAATGGGTTTGGCGTGCGTCGTGCATGGAGAATATAGAGAAGCAGTa     <  1:41535/249‑1 (MQ=255)
                                                                                                                                                                                                                                            tGGGTAACTATCGATTTTGCAATAATGGAACCTAACGCTGTCGGTGCCATCATACAGCCAGCAATAAACGCCTGATAACCAAATCCTACCTGTAACATCAATGGCATAAGGAACGGTACACAGCCGGTCCCCAGACGGGTTGCAATATTGCCTACGATACCGATCGAGAAAGTGCGGGTTTTAAATAAATCTAATGAAATTAATGGGTTTGGCGTGCGTCTTGCATGGAGAATATAGAGAAGCAGTAAcc  <  1:118608/250‑1 (MQ=255)
                                                                                                                                                                                                                                            tGGGTAACTATCGATTTTGCAATAATGGAACCTAACGCTGTCGGTGCCATCATACAGCCAGCAATAAACGCCTGATAACCAAATCCTACCTGTAACATCAATGGCATAAGGAACGGTACACAGCCGGTCCCCAGACGGGTTGCAATATTGCCTACGATACCGATCGAGAAAGTGCGGGTTTTAAATAAATCTAATGAAATTAATGGGTTTGGCGTGCGTCGTGCATGGAGAATATAGAGAAGCAGTAAc   <  1:184345/249‑1 (MQ=255)
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                 
GCGTGACCGCCAGAACACTGTTACCGCTGCTGGCATTGTCATCGGTCAGATCGGCAAGTGTGATGGTATTCATCGCGGTAAACTGCGTCGACATAGCCATCCCTAATATAAACAACGGCAAGATCAGCATCCATATAGCCATTGCTGGTGATTGCAAAGAGAACTGAGCGATCATTAGCCCAATAATCACCGTGATCCCCACTAACGTATGGCGATAGCCCAGACGACGTAAGACTTGGGTAACCATCGATTTTGCAATAATGGAACCTAACGCTGTCGGTGCCATCATACAGCCAGCAATAAACGCCTGATAACCAAATCCTACCTGTAACATCAATGGCATAAGGAACGGTACACAGCCGGTCCCCAGACGGGTTGCAATATTGCCTACGATACCGATCGAGAAAGTGCGGGTTTTAAATAAATCTAATGAAATTAATGGGTTTGGCGTGCGTCGTGCATGGAGAATATAGAGAAGCAGTAACC  >  NC_000913/3939418‑3939903

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: