Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,288,385 G→A V208I (GTT→ATT)  nrfA → cytochrome c552 nitrite reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,288,3850GA94.3% 193.2 / 2.4 72V208I (GTT→ATT) nrfAcytochrome c552 nitrite reductase
Reads supporting (aligned to +/- strand):  ref base G (2/2);  new base A (34/32);  total (38/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.69e-01

AAGCCCGGATGTGGCGCGTCTGATCCAGAAAGACGGCGAAGATGGCTACTTCCACGGTAAATGGGCGCGCGGCGGTCCGGAAATCGTCAACAACTTAGGTTGTGCCGATTGCCATAACACCGCCTCTCCAGAGTTCGCCAAAGGCAAACCGGAGTTAACCCTTTCCCGTCCGTATGCGGCTCGCGCGATGGAAGCCATTGGTAAACCTTTTGAGAAAGCCGGACGTTTCGACCAGCAATCGATGGTTTGCGGTCAGTGCCATGTGGAGTATTACTTCGACGGCAAAAACAAAGCGGTTAAATTCCCGTGGGATGACGGCATGAAAGTCGAAAATATGGAGCAGTATTACGACAAAATTGCCTTCTCTGACTGGACTAACTCCCTGTCGAAAACGCCAATGCTGAAAGCGCAGCACCCGGAATATGAAACCTGGACAGCGGGCATTCACGGTAAAAACAACGTGACCTGTATCGACTGCCATATGCCAAAAG  >  NC_000913/4288141‑4288631
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                      
aaGCCCGGATGTGGCGCGTCTGATCCAGAAAGACGGCGAAGATGGCTACTTCCACGGTAAATGGGCGCGCGGCGGTCCGGAAATCGTCAACAACTTAGGTTGTGCCGATTGCCATAACACCGCCTCTCCAGAGTTCGCCAAAGGCAAACCGGAGTTAACCCTTTCCCGTCCGTATGCGGCTCGCGCGATGGAAGCCATTGGTAAACCTTTTGAGAAAGCCGGACGTTTCGACCAGCAATCGATGATTTGCg                                                                                                                                                                                                                                                  >  1:674571/1‑251 (MQ=255)
      ggATGTGGCGCGTCTGATCCAGAAAGACGGCGAAGATGGCTACTTCCACGGTAAATGGGCGCGCGGCGGTCCGGAAATCGTCAACAACTTAGGTTGTGCCGATTGCCATAACACCGCCTCTCCAGAGTTCGCCAAAGGCAAACCGGAGTTAACCCTTTCCCGTCCGTATGCGGCTCGCGCGATGGAAGCCATTGGTAAACCTTTTGAGAAAGCCGGACGTTTCGACCAGCAATCGATGATTTGCGGTCAGt                                                                                                                                                                                                                                            >  1:242902/1‑251 (MQ=255)
                  tCTGATCCAGAAAGACGGCGAAGATGGCTACTTCCACGGTAAATGGGCGCGCGGCGGTCCGGAAATCGTCAACAACTTAGGTTGTGCCGATTGCCATAACACCGCCTCTCCAGAGTTCGCCAAAGGCAAACCGGAGTTAACCCTTTCCCGTCCGTATGCGGCTCGCGCGATGGAAGCCATTGGTAACCCTTTTGAGAAAGCCGGACGTTTCGACCAGCCATCGATGATTTGCGGTCAGTGCCATGTGGAGt                                                                                                                                                                                                                                >  2:595071/1‑251 (MQ=255)
                     gaacaaaaaaGACGGCGAAGATGGCTACTTCCACGGTAAATGGGCCCGCGCAGGTCCGAAACTCGCAAACAACTTAGTTTGTGCCGATTGCCATACCACCGCCTCTCCAGAGTTCGCCAAAGGCAAACCGGAGTTAACCCTTTCCCGTCCGTATGCGGCTCGCGCGATGGAAGCCATTGGTAAACCTTTTGAGAAAGCCGGACGTTTCGACCAGCAATCGATGATTTGCGGTCAGTGCCATGTGGAGt                                                                                                                                                                                                                                <  1:638379/241‑1 (MQ=255)
                         cAGAAAGACGGCGAAGATGGCTACTTCCACGGTAAATGGGCGCGCGGCGGTCCGGAAATCGTCAACAACTTAGGTTGTGCCGATTGCCATAACACCGCCTCTCCAGAGTTCGCCAAAGGCAAACCGGAGTTAACCCTTTCCCGTCCGTATGCGGCTCGCGCGATGGAAGCCATTGGTAAACCTTTTGAGAAAGCCGGACGTTTCGACCAGCAATCGATGATTTGCGGTCAGTGCCATGAGGAGTATTACtt                                                                                                                                                                                                                         >  2:177148/1‑251 (MQ=255)
                             aaGAGGGCGAAGGGGGCTGCTTCCACGTAAAATGGGCGCGCGGCGGTCCGGAAACCGTCAACATGTTAGGTGGTCCCGATTGCTATAACACCGCCTCTCCAGAGTTCGCCAAAGGCAAACCGGAGTTAACCCTTTCCCGTCCGTATGCGGCTCGCGCGATGGAAGCCATTGGTAAACCTTTTGAGAAAGCCGGACGTTTCGACCAGCAATCGATGATTTTCGGTCAGTGCCATGTGGAGTATTACTTCGAc                                                                                                                                                                                                                     <  1:516109/251‑1 (MQ=255)
                             aaGACGACGAAGATGGCTACTTCCACGGTAAATGGGCGCGCTGCGGTCCGGAAATCGTCAACAACTAAGGTTGTGCCGATTGCCATAACACCGCCTCTCCAGAGGTCGCCAAAGGAAAACCGGAGTTAACCCTTTCCCGTCCGTATGCGGCTCGCGCGATGGAAGCCATTGGTACAACTTTTGAGACAGCCGGACGTTTCGAACAGCAATAGATGATTTGAGGTAAGTGGCATGGGGCGTATTAATTCggc                                                                                                                                                                                                                     >  2:325534/1‑249 (MQ=255)
                                      aagttGGCTACTCCCACGGTAAAGGGGCGCGCGGCGGTCCGAAAATCGTCAACATCTTAGGTTGTGCCGATTGCCATAACACCGCCTCTCCAGAGTTCGCCAAAGGCAAACCGGAGTTAACCCTTTCCCGTCCGTATGCGGCTCGCGCGATGGAAGCCATTGGTAAACCTTTTGAGAAAGCCGGACGTTTCGACCAGCAATCGATGATTTGCGGTCAGTGCCATGTGGAGTATTACTTCGACGGCAAAAAc                                                                                                                                                                                                            <  1:364110/247‑1 (MQ=255)
                                         agggcctcctCCACGGTAAATGGTCGCTCGGCGGACCGGACATCGTTACCATCTCAGGTTGTGCCGTGTGCCATAACTCCGCCTCTCCAGAGTTCGCCAAAGGCAAACCGGAGTTAACCCTTTCCCGTCCGTATGCGGCTCGCGCGATGGAAGCCATTGGTAATCCTTTTGAGAAAGCCGGACGTTTCGACCAGCAATCGATGATTTGCGGTCAGTGCCATGTGGAGTATTACTTCGACGGCAAAAACaaa                                                                                                                                                                                                         <  1:387129/242‑1 (MQ=255)
                                           ggCTACTTCCACGGAAAATGGGCGCGCGGCGGTCCGGAAATCGTCAACAACTTAGGTTGTGCCGATTGCCATAACACCGCCTCTCCAGAGTTCGCCAAAGGCAAACCGGAGTTAACCCTTTCCCGTCCGTATGCGGCTCGCGCGATGGAAGCCATTGGTAAACCTTTTGAGAAAGCCGGACGTTTCGACCAGCACTCGATGATTTGCGGTCAGTGCCATGTGGAGTATTACTTcaccgtcaataacaacgg                                                                                                                                                                                                       >  1:466865/1‑234 (MQ=255)
                                            gCTACTTCCACGGTAAATGGGCGCGCGGCGGTCCGGAAATCGTCAACAACTTAGGTTGTGCCGATTGCCATAACACCGCCTCTCCAGAGTTCGCCAAAGGCAAACCGGAGTTAACCCTTTCCCGTCCGTATGCGGCTCGCGCGATGGAAGCCATTGGTAAACCTTTTGAGAAAGCCGGACGTTTCGACCAGCAATCGATGATTTGCGGTCAGTGCCATGTGGAGTATTACTTCGACGGCAAAAACAAAGCg                                                                                                                                                                                                      >  1:431429/1‑251 (MQ=255)
                                                  gccgcggaaaaggggggcgcgGGGGCCCGGAAATCGACAACAACTGGGGTTGTGCAGATTGCCAAAACACCGCCTCTCCAGAGTTCGCCAAAGGCAAACCGGAGTTAACCCTTTCCCGTCCGTATGCGGCTCGCGCGATGGAAGCCATTGGTAAACCTTTTGAGAAAGCCGGACGTTTCGACCAGCAATCGATGATTTGCGGTCAGTGCCATGTGGAGTATTACTTCGACGGCAAAAACAAAGCGGTTaa                                                                                                                                                                                                 <  1:184264/234‑1 (MQ=255)
                                                  cccACGTTAAAGGGGCGCGCGGGGGTCCGGAAATCGCAAACAACTTAGGTTGGGCCGATTGCCAAAACACCGCCTCTCCAGAGTTCGCCAAAGGCAAACCGGAGTTAACCCTTTCCCGTCCGTATGCGGCTCGCGCGATGGAAGCCATTGGTAAACCTTTTGAGAAAGCCGGACGTTTCGACCAGCAATCGATGATTTGCGGTCAGTGCCATGTGGAGTATTACTTCGACGGCAAAAACAAAGCGGTTa                                                                                                                                                                                                  <  1:595071/248‑1 (MQ=255)
                                                           aaTGGGCGCGCGGCGGTCCGGAAATCGTCAACAACTTAGGTTGTGCCGATTGCCATAACACCGCCTCTCCAGAGTTCGCCAAAGGCAAACCGGAGTTAACCCTTTCCCGTCCGTATGCGGCTCGCGCGATGGAAGCCATTGGTAAACCTTTTGAGAAAGCCGGACGTTTCGACCAGCAATCGATGATTTGCGGTCAGTGCCATGTGGAGTATTACTTCGACGGCAAAAACAAAGCGGTTAAATTCCCGTgg                                                                                                                                                                                       >  1:43172/1‑251 (MQ=255)
                                                            aTGGGCGCGCGGCGGTCCGGAAATCGTCAACAACTTAGGTTGTGCCGATTGCCATAACACCGCCTCTCCAGAGTTCGCCAAAGGCAAACCGGAGTTAACCCTTTCCCGTCCGTATGCGGCTCGCGCGATGGAAGCCATTGGTAAACCTTTTGAGAAAGCCGGACGTTTCGACCAGCAATCGATGATTTGCGGTCAGTGCCATGTGGAGTATTACTTCGACGGCAAAAACAAAGCGGTTAAATTCCCGTggg                                                                                                                                                                                      >  1:462246/1‑251 (MQ=255)
                                                             tGGGCGCGCGGCGGTCCGGAAATCGTCAACAACTTAGGATGTGCCGATTGCCATAACACCGCCTCTCCAGAGTTCGCCAAAGGCAAACCGGAGTTAACCCTTTCCCGTCCGTATGCGGCTCGCGCGATGGAAGCCATTGGTAAACCTTTTGAGAAAGCCGGACGTTTCGTCCAGCAATCGATGTTTTGAGGTCAGTGCCATGGGGAGTATTACTTCGCCGGCAAAAACAAAGAGGTTAAATTCCCGTggg                                                                                                                                                                                      >  2:216382/1‑250 (MQ=255)
                                                                    ggggcggcCCGAAAATCGTAAACACTTTGGGTTGTGCCTATTGCCATAACACCGCCTCTCCAGAGTTCGCCAAAGGCAAACCGGAGTTAACCCTTTCCCGTCCGTATGCGGCTCGCGCGATGGAAGCCATTGGTAAACCTTTTGAGAAAGCCGGACGTTTCGACCAGCAATCGATGATTTGCGGTCAGTGCCATGTGGAGTATTACTTCGACGGCAAAAACAAAGCGGTTAAATTCCCGTGGGATGACGGc                                                                                                                                                                              <  1:312246/249‑1 (MQ=255)
                                                                    gcggcgGGCCGAAAATCGTCAACAACTGAGGTTGGGCCGTTTGCCATAACACCGCCTCTCCAGAGTTCGCCAAAGGCAAACCGGAGTTAACCCTTTCCCGTCCGTATGCGGCTCGCGCGATGGAAGCCATTGGTAAACCTTTTGAGAAAGCCGGACGTTTCGACCAGCAATCGATGATTTGCGGTCAGTGCCATGTGGAGTATTACTTCGACGGCAAAAACAAAGCGGTTAAATTCCCGTGGGATGACGGc                                                                                                                                                                              <  1:444760/251‑1 (MQ=255)
                                                                      ggcggTCCGGCAATCGTCAACAACTTAGGTTGTGCCGATTGCCATAACACCGCCTCTCCAGAGTTCGCCAAAGGCAAACCGGAGTTAACCCTTTCCCGTCCGTATGCGGCTCGCGCGATGGAAGCCATTGGTAAACCTTTTGAGAAAGCCGGACGTTTCGACCAGCAATCGATGATTTGCGGTCAGTGCCATGTGGAGTATTACTTCGACGGCAAAAACAAAGCGGTTAAATTCCCGTGGGATGACGGCAt                                                                                                                                                                            >  2:286258/1‑251 (MQ=255)
                                                                       gcgtcccGGAAAGCATCAACAACTTAGGTTGGGCCGGTTGCCATAACACCGCCTCTCCAGAGTTCGCCAAAGGCAAACCGGAGTTAACCCTTTCCCGTCCGTATGCGGCTCGCGCGATGGAAGCCATTGGTAAACCTTTTGAGAAAGCCGGACGTTTCGACCAGCAATCGATGATTTGCGGTCAGTGCCATGTGGAGTATTACTTCGACGGCAAAAACAAAGCGGTTAAATTCCCGTGGGATGACGGCATg                                                                                                                                                                           <  1:591251/246‑1 (MQ=255)
                                                                          gtccgaaaaccttaaacatctaagtttGTGCCGATTCCCAAACCACCGCCTCTCCAGAGTTCGCCAAAGGCAAACCGGAGTTAACCCTTTCCCGTCCGTATGCGGCTCGCGCGATGGAAGCCATTGGTAAACCTTTTGAGAAAGCCGGACGTTTCGACCAGCAATCGATGATTTGCGGTCAGTGCCATGTGGAGTATTACTTCGACGGCAAAAACAAAGCGGTTAAATTCCCGTGGGATGACGGCATGaaa                                                                                                                                                                        <  1:49700/226‑1 (MQ=255)
                                                                            taggAAATCGTCAACAACTTAGGATGTGCCGATTGCCATAACACCGCCTCTCCAGAGTTCGCCAAAGGCAAACCGGAGTTAACCCTTTCCCGTCCGTATGAGGCTCGCGCGATGGAAGCCATTGGAAAACCTTTTGAGACAGCCGGACGTTTCGCCCACCAATCGCTGATTTGCGGTCAGTGCAATGTGGAGAATTACTTCGCCGACAAAAACAAAGCGGTTAAATTCCCGTGGGATACCGGCATGAAAGt                                                                                                                                                                      >  2:228427/3‑251 (MQ=255)
                                                                              ggAAATCGTCAACAACTTAGGTTGTGCCGATTGCCATAACACCGCCTCTCCAGAGTTCGCCAAAGGCAAACCGGAGTTAACCCTTTCCCGTCCGTATGCGGCTCGCGCGATGGAAGCCATTGGTAAACCTTTTGAGAAAGCCGACCGTTTCGCCAACCAATCGATGATTTGCGGTCAGGGCCATGTGGAGAATTACTTGGACGGCAAAAACAAACCGGTTAAATTCCCGTGGAATGACGCCATGAAAgccg                                                                                                                                                                    >  1:581803/1‑248 (MQ=255)
                                                                              agAAATCGTCAACAACGTGGGTTGTGCCGATTTCCATCACACCGCCTCTCCAGAGTTCGCCAAAGGCAAACCGGAGTTAACCCTTTCCCGTCCGTATGCGGCTCGCGCGATGGAAGCCATTGGTAAACCTTTTGAGAAAGCCGGACGTTTCGACCAGCAAGCGATGATTTGCGGTCAGTGCCATGTGGAGTATTACTTCGACGGCAAAAACAAAGCGGTTAAATTCCCGTGGGATGACGGCATGAAAGTCg                                                                                                                                                                    <  1:617403/250‑1 (MQ=255)
                                                                                aaaTCGTCAACAACTTAGGTTGTGCCGATTGCCATAACACCGCCTCTCCAGAGTTCGCCAAAGGCAAACCGGAGTTAACCCTTTCCCGTCCGTATGCGGCTCGCGCGATGGAAGCCATTGGTAAACCTTTTGAGAAAGCCGGACGTTTCGACCAGCAATCGATGATTTGCGGTCAGTGCCATGTGGAGTATTACTTCGACGGCAAAAACAAAGCGGTTAAATTCCCGTGGGATGACGGCATGAAAGTCGa                                                                                                                                                                   >  1:61319/1‑250 (MQ=255)
                                                                                  aTCGTCAACAACTTAGGTTGTGCCGATTGCCATAACACCGCCTCTCCAGAGTTCGCCAAAGGCAAACCGGAGTTACCCCTTTCCCGTCCGTATGCGGCTCGCGCGATGGAAGCCATTGGTAAACCTTTTGAGAAAGCCGGACGTTTCGACCAGCAATCGATGATTTGCGGTCAGTGCCATGTGGAGTATTACTTCGACGGCAAAAACAAAGCGGTTAAATTCCCGTGGGATGACGGCATGAAAGTCGaaa                                                                                                                                                                 >  2:93451/1‑250 (MQ=255)
                                                                                         acattagtgtgtgtgcagtgtGCCAAAACCCCGCCTCTCCAGAGTTCGCCAAAGGCAAACCGGAGTTAACCCTTTCCCGTCCGTATGCGGCTCGCGCGATGGAAGCCATTGGTAAACCTTTTGAGAAAGCCGGACGTTTCGACCAGCAATCGATGATTTGCGGTCAGTGCCATGTGGAGTATTACTTCGACGGCAAAAACAAAGCGGTTAAATTCCCGTGGGATGACGGCATGAAAGTCGAAAATATGGAg                                                                                                                                                         <  2:431429/231‑1 (MQ=255)
                                                                                         acaacTTTGGTTGTGCCGATTGCCATAACACCGCCTCTCCAGAGTTCGCCAAAGGCAAACCGGAGTTAACCCTTTCCCGTCCGTATGCGGCTCGCGCGATGGAAGCCATTGGTAAACCTTTTGAGAAAGCCGGACGTTTCGACCAGCAATCGATGATTTGCGGTCAGTGCCATGTGGAGTATTACTTCGACGGCAAAAACAAAGCGGTTAAATTCCCGTGGGATGACGGCATGAAAGTCGAAAATATGGAg                                                                                                                                                         <  1:116825/251‑1 (MQ=255)
                                                                                                ggtttgggCCGATTGCCATAACACCGCCTCTCCAGAGTTCGCCAAAGGCAAACCGGAGTTAACCCTTTCCCGTCCGTATGCGGCTCGCGCGATGGAAGCCATTGGTAAACCTTTTGAGAAAGCCGGACGTTTCGACCAGCAATCGATGATTTGCGGTCAGTGCCATGTGGAGTATTACTTCGACGGCAAAAACAAAGCGGTTAAATTCCCGTGGGATGACGGCATGAAAGTCGAAAATATGGAGCAGTAtt                                                                                                                                                  <  1:661262/244‑1 (MQ=255)
                                                                                                 ggTTGGGCCGATTGCCCTAACACCGCCTCTCCAGAGTTCGCCAAAGGCAAACCGGAGTTAACCCTTTCCCGTCCGTATGCGGCTCGCGCGATGGAAGCCATTGGTAAACCTTTTGAGAAAGCCGGACGTTTCGACCAGCAATCGATGGTTTGCGGTCAGTGCCATGTGGAGTATTACTTCGACGGCAAAAACAAAGCGGTTAAATTCCCGTGGGATGACGGCATGAAAGTCGAAAATATGGAGCAGTATTa                                                                                                                                                 <  1:152606/251‑1 (MQ=255)
                                                                                                  gttgtgtcttttgcctaaccacCGCCTCTCCAGAGTTCGCCAAAGGCAAACCGGAGTTAACCCTTTCCCGTCCGTATGCGGCTCGCGCGATGGAAGCCATTGGTAAACCTTTTGAGAAAGCCGGACGTTTCGACCAGCAATCGATGATTTGCGGTCAGTGCCATGTGGAGTATTACTTCGACGGCAAAAACAAAGCGGTTAAATTCCCGTGGGATGACGGCATGAAAGTCGAAAATATGGAGCAGTATTAc                                                                                                                                                <  1:246433/232‑1 (MQ=255)
                                                                                                           aTTGCCATAACACCGCCTCTCCAGAGTTCGCCAAAGGCAAACCGGAGTTAACCCTTTCCCGTCCGTATGCGGCTCGCGCGATGGAAGCCATTGGTAAACCTTTTGAGAAAGCCGGACGTTTCGACCAGCAATCGATGATTTGCGGTCAGTGCCATGTGGAGTATTACTTCGACGGCAAAAACAAAGCGGTTAAATTCCCGTGGGATGACGGCATGAAAGTCGAAAATATGGAGCAGTATTACGACAAAAtt                                                                                                                                       <  1:14281/251‑1 (MQ=255)
                                                                                                                attccaccctccTCTCCAGATTTCGCCAAAGGCAAACCGGAGTTTACCCTTTCCCGTCCGTATGCGGCTCGCGCGATGGAAGCCATTGGTAAACCTTTTGAGAAAGCCGGACGTTTCGACCAGCAATCGATGATTTGCGGTCAGTGCCATGTGGAGTATTACTTCGACGGCAAAAACAAAGCGGTTAAATTCCCGTGGGATGACGGCATGAAAGTCGAAAATATGGAGCAGTATTACGACAAAATTGCCtt                                                                                                                                  <  2:528466/241‑1 (MQ=255)
                                                                                                                 aaaaccccGCCTCTCCAGAGTTCGCCTAAGGCCAACCGGAGTTAACCCTTTCCCGTCCGTATGCGGCTCGCGCGATGGAAGCCATTGGTAAACCTTTTGAGAAAGCCGGACGTTTCGACCAGCAATCGATGATTTGCGGTCAGTGCCATGTGGAGTATTACTTCGACGGCAAAAACAAAGCGGTTAAATTCCCGTGGGATGACGGCATGAAAGTCGAAAATATGGAGCAGTATTACGACAAAATTGCCTtc                                                                                                                                 <  2:530996/245‑1 (MQ=255)
                                                                                                                    ccacCGCCTCTCCAGAGTTCGCCAAAGGCAAACCGGAGTTAACCCTTTCCCGTCCGTATGCGGCTCGCGCGATGGAAGCCATTGGTAAACCTTTTGAGAAAGCCGGACGTTTCGACCAGCAATCGATGATTTGCGGTCAGTGCCATGTGGAGTATTACTTCGACGGCAAAAACAAAGCGGTTAAATTCCCGTGGGATGACGGCATGAAAGTCGAAAATATGGAGCAGTATTACGACAAAATTGCCTTctct                                                                                                                              <  1:39854/250‑1 (MQ=255)
                                                                                                                      cccGCCTGTCCAGAGTTCACCAAAGGATAACCGGAGTTAACCCTGTCCCGTCAGTATGCGGCTCGCGCGATGGAAGCCTTTTGTAAACCTTGTGATAAATCCTGACTTTTCGACCAGCAATCGATGATTTGCGGTCAGTTCCATGTGGAGTATTACTTCGACGGCAAAAACAAAGCGGTTAAATGCCCGTGGGATGACGCCATTAAAGTCGACAAGATGGAGCAGTATTACGACAAAATTGCCTTCTCTGa                                                                                                                            <  2:250054/250‑1 (MQ=255)
                                                                                                                       ccGCCTCTCCAGAGTTCGCCAAAGGCAAACCGGAGTTAACCCTTTCCCGTCCGTATGCGGCTCGCGCGATGGAAGCCATTGGTAAACCTTTTGAGAAAGCCGGACGTTTCGACCAGCAATCGATGATTTGCGGTCCGTGCCATGTGGAGAATTACTTCGACGGCAAAAACAAAGCGGTTAAATTCCCGTGGGATGACGGCATGAAAGTCGAAAATATGGGGAAGTATTACGACAAAATTGCCTTCTCTGAc                                                                                                                           >  2:499690/1‑251 (MQ=255)
                                                                                                                           ctctCCAAAGTTCGCCAAAGGCAAACCGGAGTTAACCCTTTCCCGTCCGTATGCGGCTCGCGCGATGGAAACCATTGGTAAACCTTTTGAGAAAGCCGGACGTTTCGACCAACAGTCGATGATTTGCGGTCAGTGCCATGAGGCGAATTACTTCGACGGCAACAACAAAGCGGTTAAATTCCCGTGGGCTGAAGGCATGAAAAGCGAAAGTATGGGACATGATTACGACAAACTTGCCTTCTCTGACTgg                                                                                                                        >  2:528892/1‑250 (MQ=255)
                                                                                                                                      tCGCCAAAGGCAAACCGGAGTTAACCCTTTCCCGTCCGTATGCGGCTCGCGCGATGGAAGCCATTGGTAAACCTTTTGAGAAAGCCGGACGTTTCGACCAGCAATCGATGATTTGCGGTCAGTGCCATGTGGAGTATTACTTCGACGGCAAAAACAAAGCGGTTAAATTCCCGTGGGATGACGGCATGAAAGTCGAAAATATGGAGCAGTATTACGACAAAATTGCCTTCTCTGACTGGACTAACTCCCt                                                                                                             <  1:114884/250‑1 (MQ=255)
                                                                                                                                        gCCAAAGGCAAACCGGAGTTAACCCTTTCCCGTCCGTATGCGGCTCGCGCGATGGAAGCCATTGGTAAACCTTTTGAGAAAGCCGGACGTTTCGACCAGCAATCGATGATTTGCGGTCAGTGCCATGTGGAGTATTACTTCGACGGCAAAAACAAAGCGGTTAAATTCCCGTGGGATGACGGCATGAAAGTCGAAAATATGGAGCAGTATTACGACAAAATTGCCTTCTCTGACTGGACTAACTCCCTGTc                                                                                                          <  1:286258/251‑1 (MQ=255)
                                                                                                                                            aaGGCAAACCGGAGTTAACCCTTTCCCGTCCGTATGCGGCTCGCGCGATGGAAGCCATTGGTAAACCTTTTGAGAAAGCCGGACGTTTCGACCAGCAATCGATGATTTGCGGTCAGTGCCATGTGGAGTATTACTTCGACGGCAAAAACAAAGCGGTTAAATTCCCGTGGGATGACGGCATGAAAGTCGAAAATATGGAGCAGTATTACGACAAAATTGCCTTCTCTGACTGGACTAACTCCCTGTCGaaa                                                                                                      <  1:582456/251‑1 (MQ=255)
                                                                                                                                            aaGGCAAACCGGAGTTAACCCTTTCCCGTCCGTATGCGGCTCGCGCGATGGAAGCCATTGGTAAACCTTTTGAGAAAGCCGGACGTTTCGACCAGCAATCGATGATTTGCGGTCAGTGCCATGTGGAGTATTACTTCGACGGCAAAAACAAAGCGGTTAAATTCCCGTGGGATGACGGCATGAAAGTCGAAAATATGGAGCAGTATTACGACAAAATTGCCTTCTCTGACTGGACTAACTCCCTGTCGaa                                                                                                       >  1:206330/1‑250 (MQ=255)
                                                                                                                                             aGGCAAACCGGAGTTAACCCTTTCCCGTCCGTATGCGGCTCGCGCGCTGGAAGCCATTGGTAAACCTTTTGAGAAAGCCGGACGTTTCGACCAGCAATCGATGATTTGCGGTCCGTGCCATGTGGAGTATTACTTCGACGGCAAAAACAAAGCGGTTAAATTCCCGTGGGATGACGGCATGAAAGTCGACAATATGGAGCAGTATTACGACAAAATTGCCTTCTCTGACTGGACTAACTCCCTGCCGaaa                                                                                                      >  2:109845/1‑250 (MQ=255)
                                                                                                                                              ggCAAACCGGAGTTAACACTTTCCCGTCCGTATGCGGCTCGCGCGATGGAAGACATTGGTAAACCTTTTGAGAAAGCCGGACGTTTCGACCAGCAATCGATGATTTGCGGTCAGTGCCATGTGGAGTATTACTTCGACGGCAAAAACAAAGCGGTTAAATTCCCGTGGGATGACGGCATGAAAGTCGAAAATATGGAGCAGTATTACGCCAAAATTGCCTTCTCTGACTGGACTAACTCCCTGTCGAAAAc                                                                                                    >  2:36345/1‑251 (MQ=255)
                                                                                                                                                       gAGTTAACCCTTTCCCGTCCGTATGCGGCTCGCGCGATGGAAGCCATTGGTAAACCTTTTGAGAAAGCCGGACGTTTCGACCAGCAATCGATGATTTGCGGTCAGTGCCATGTGGAGTATTACTTCGACGGCAAAAACAAAGCGGTTAAATTCCCGTGGGATGACGGCATGAAAGTCGAAACTATGGAGCAGTATTACGACACAATTGCCTTCTCTGACTGGACTAACCACCTGTCGAAAACTCAAATGCt                                                                                           >  2:228650/1‑251 (MQ=255)
                                                                                                                                                        aGTTAACCCTTTCCCGTCCGTATGCGGCTCGCGCGATGGAAGCCATTGGTAAACCTTTTGAGAAAGCCGGACGTTTCGACCAGCAATCGATGATTTGCGGTCAGTGCCATGTGGAGTATTACTTCGACGGCAAAAACAAAGCGGTTAAATTCCCGTGGGATGACGGCATGAAAGTCGAAAATATGGAGCAGTATTACGACAAAATTGCCTTCTCTGACTGGACTAACTCCCTGTCGAAAACGTCAATg                                                                                             <  1:207098/248‑1 (MQ=255)
                                                                                                                                                        aGTTAACCCTTTCCCGTCCGTATGCGGCTCGCGCGATGGAAGCCATTGGTAAACCTTTTGAGAAAGACGGACGTTTCGACCAGCAATCGATGATTTGCGGTCAGTGCCATGTGGAGTATTACTTCGACGGCAAAAACAAAGCGGTTAAATTCCCGTGGGATGACGGCATGAAAGTCGAAAATATGGAGCAGTATTACGACAAAATTGCCTTCTCTGACTGGACTAACTCCCTGTCGAAAACGCCAATg                                                                                             >  2:207098/1‑248 (MQ=255)
                                                                                                                                                         gTTAACCCTTTCCCGTCCGTATGCGGCTCGCGCGATGGAAGCCATTGGTAAACCTTTTGAGAAAGCCGGACGTTTCGACCAGCAATCGATGATTTGCGGTCAGTGCCATGTGGAGTATTACTTCGACGGCAAAAACAAAGCGGTTAAATTCCCGTGGGATGACGGCATGAAAGTCGAAAATATGGAGCAGTATTACGACAAAATTGCCTTCTCTGACTGGACTAACTCCCTGTCGAAAACGCCAATGCTg                                                                                          >  1:591371/1‑250 (MQ=255)
                                                                                                                                                            aaCCCTTTCCCGTCCGTATGCGGCTCGCGCGATGGAAGCCATTGGTAAACCTTTTGAGAAAGCCGGACGTTTCGACCAGCAATCGATGATTTGCGGTCAGTGCCATGTGGAGTATTACTTCGACGGCAAAAACAAAGCGGTTAAATTCCCGTGGGATGACGGCATGAAAGTCGAAAATATGGAGCAGTATTACGACAAAATTGCCTTCTCTGACTGGACTAACTCCCTGTCGAAAACGCCAATGCTGaaa                                                                                       <  1:22832/250‑1 (MQ=255)
                                                                                                                                                             aCCCTTTCCCGTCCGTATGCGGCTCGCGCGATGGAAGCCATTGGTAAACCTTTTGAGAAAGCCGGACGTTTCGACCAGCAATCGATGATTTGCGGTCAGTGCCATGTGGAGTATTACTTCGACGGCAAAAACAAAGCGGTTAAATTCCCGTGGGATGACGGCATGAAAGTCGAAAATATGGAGCAGTATTACGACAAAATTGCCTTCTCTGACTGGACTAACTCCCTGTCGAAAACGCCAATGCTGAAAgc                                                                                     <  1:36345/251‑1 (MQ=255)
                                                                                                                                                             aCCCTTTCCCGTCCGTATGCGGCTCGCGCGATGGAAGCCATTGGTAAACCTTTTGAGAAAGCCGGACGTTTCGACCAGCAATCGATGATTTGCGGTCAGTGCCATGTGGAGTATTACTTCGACGGCAAAAACAAAGCGGTTAAATTCCCGTGGGATGACGGCATGAAAGTCGAAAATATGGAGCAGTATTACGACAAAATTGCCTTCTCTGACTGGACTAACTCCCTGTCGAAAACGCCAATGCTGAAAgc                                                                                     <  1:624579/251‑1 (MQ=255)
                                                                                                                                                              cccTTTCCCGTCCGTATGCGGCTCGCGCGATGGAAGCCATTGGTAAACCTTTTGAGAAAGCCCGACGTTTCGACCAGCAATCGATGCTTTGCGGTCAGTGACCTGTGGAGTATTACTTCGAAGGCAAAAACAAAGCGGTTAAATTCCCGTGGGATGACGGCATGAAAGTCGAAAATATGGAGCAGTATTACGACAAAATTGCCTTCTCTGACTGGACTAACTCCCTGTCGAAAACGCCAATGCTGAAAgcg                                                                                    >  2:191147/1‑251 (MQ=255)
                                                                                                                                                                   ccccGTACGTATGCGGCGCGCGCGATGGAAGCCATTGTTAAACCCTTTGAGAAAGCCGGACGTTTCGACCAGCAATCGATGATGTGCGGTCAGTGCCATGTGGAGTATTACTTCGACGGCAAAAACAAAGCGGTTAAATTCCCGTGGTATGACGGCATGAAAGTCGAAAATATGGAGCAGTATGACGACAAAATTGCCTTCGCTGACTTGACTAACTCCCTGTCGAAAACTCCAATGCTGAAAGCtcagca                                                                               <  2:594676/250‑1 (MQ=255)
                                                                                                                                                                      cgtccgtATGCGGCTCGCGCGATGGAAGCCATTGGTAAACCTTTTGAGAAAGCCGGACGTTTCGACCAGCAATCGATGATTTGCGGTCAGTGCCATGTGGAGTATTACTTCGACGGCAAAAACAAAGCGGTTAAATTCCCGTGGGATGACGGCATGAAAGTCGAAAATATGGAGCAGTATTACGACAAAATTGCCTTCTCTGACTGGACTAACTCCCTGTCGAAAACGCCAATGCTGAAAGCGCAGCcccc                                                                            >  1:597361/1‑251 (MQ=255)
                                                                                                                                                                            tATGCGGCTCGCGCGATGGAAGCCATTGGTAAACCTTTTGAGAAAGCCGGACGTTTCGACCAGCAATCGATGATTTGCGGTCAGTGCCAt                                                                                                                                                                                                                                       >  2:171888/1‑90 (MQ=255)
                                                                                                                                                                            tATGCGGCTCGCGCGATGGAAGCCATTGGTAAACCTTTTGAGAAAGCCGGACGTTTCGACCAGCAATCGATGATTTGCGGTCAGTGCCAt                                                                                                                                                                                                                                       <  1:171888/90‑1 (MQ=255)
                                                                                                                                                                                       cgcgATGGAAGCCATTGGTAAACCTTTTGAGAAAGCCGGACGTTTCGACCAGCAATCGATGGTTTGCGGTCAGTGCCATGTGGAGTATTACTTCGACGGCAAAAACAAAGCGGTTAAATTCCCGTGGGATGACGGCATg                                                                                                                                                                           >  1:357629/1‑139 (MQ=255)
                                                                                                                                                                                       cgcgATGGAAGCCATTGGTAAACCTTTTGAGAAAGCCGGACGTTTCGACCAGCAATCGATGGTTTGCGGTCAGTGCCATGTGGAGTATTACTTCGACGGCAAAAACAAAGCGGTTAAATTCCCGTGGGATGACGGCATg                                                                                                                                                                           <  2:357629/139‑1 (MQ=255)
                                                                                                                                                                                          gATGGAAGCCATTGGTAAACCTTTTGAGAAAGCCGGACGTTTCGACCAGCAATCGATGATTTGCGGTCAGTGCCATGTGGAGTATTACTTCGACGGCAAAAACAAAGCGGTTAAATTCCCGTGGGATGACGGCATGAAAGTCGAAAATATGGAGCAGTATTACGACAAAATTGCCTTCTCTGACTGGACTAACTCCCTGTCGAAAACGCCAATGCTGAAAGCGCAGCACCCGGAATATGAAACCTGGAc                                                          >  1:146768/1‑249 (MQ=255)
                                                                                                                                                                                               aaGCCATTGGTAAACCTTTTGAGAAAGCCGGACGATTCGACCAGCAATCGATGATTTGCGGTCAGTGCCATGTGGAGTATTACTTCGACGGCAAAAACAAATCGGTTAAATTCGCGTGGGATGACGGCATGAAAGTCGAAAATATGGAGCAGTATTACGACAAAATTTCCTTCTCTGACTGGACTAACTCCCTGTCGAAAACTCCAATGCTGAAAGCGCAGAACCCAGAATATGGAACCTGGACAGCggg                                                    >  1:374062/1‑250 (MQ=255)
                                                                                                                                                                                                 gCCATTGGTAAACCTTTTGAGAAAGCCGGACGTTTCGACCAGCAATCGATGATTTGCGGTCAGTGCCATGTGGAGTATTACTTCGACGGCAAAAACAAAGCGGTTAAATTCCCGTGGGATGACGGCATGAAAGTCGCAAATATGGAGCAGTATTACGACAAAATTGCCTTCTCTGCCTGGACTAACTCCCTGTCGAAAACGCCAATGCTGAAAGCGCAGCACCCGGAATATGAAACCTGGACAGCGGGCAt                                                 >  1:403927/1‑251 (MQ=255)
                                                                                                                                                                                                      tGGTAAACCTTTTGAGAAAGCCGGACGTTTCGACCAGCAATCGATGATTTGCGGTCAGTGCCATGTGGAGTATTACTTCGACGGCAAAAACAAAGCGGTTAAATTCCCGTGGGATGACGGCATGAAAGTCGAAAATATGGAGCAGTATTACGACAAAATTGCCTTCTCTGACTGGACTAACTCCCTGTCGAAAACGCCAATGCTGAAAGCGCAGCACCCGGAATATGAAACCTGGACAGCGGGCATTCAc                                             <  1:93451/250‑1 (MQ=255)
                                                                                                                                                                                                                      aatgCCGTACGTTTCGACCAGCAATTGATGATTTGCGGGCATTGCCATGTGGAGTATTACTTCGACGGCAAAAACAAAGCGGTTAAATTCACGTGGGATGACGGCATGAAAGTCGAAAATATGGAGCAGTATTACGACAAAATTGCCTTCTCTGACTGGACTAACTCCCTGTCGAAAACGCCATTGCTGAAAGCGCAGCACCCGGAATATGAAACCTGGACAGCGGTCATTCACGGTAAAAACAACGTGaa                            <  2:265785/248‑2 (MQ=255)
                                                                                                                                                                                                                        aGCCGGACGTTTCGACCAGCAATCGATGATTTGCGGTCAGTGCCATGTGGAGTATTACTTCGACGGCAAAAACAAAGCGGTTAAATTCCCGTGGGATGACGGCATGAAAGTCGAAAATATGGAGCAGTATTACGACAAAATTGCCTTCTCTGACTGGACTAACTCCCTGTCGAAAACGCCAATGCTGAAAGCGCAGCACCCGGAATATGAAACCTGGACAGCGGGCATTCACGGTAAAAACAACGTGACCt                          >  1:417847/1‑251 (MQ=255)
                                                                                                                                                                                                                          ccGGACGTTTCGACCAGCAATCGATGATTTGCGGTCAGTGCCATGTGGAGTATTACTTCGACGGCAAAAACAAAGCGGTTAAATTCCCGTGGGATGACGGCATGAAAGTCGAAAATATGGAGCAGTATTACGACAAAATTGCCTTCTCTGACTGGACTAACTCCCTGTCGAAAACGCCAATGCTGAAAGCGCAGCACCCGGCATATGAAACCTGGACAGCGGGCATTCACGGTAAAAACAACGTGACCTGt                        >  1:28055/1‑251 (MQ=255)
                                                                                                                                                                                                                          ccGGACGTTTCGACCAGCAATAGATGATTTGCGGTCAGTGCCATGTGGAGTATTACTTCGACGGCAAAAACAAAGCGGTTAAATTCCCGTGGGATGACGGCATGAAAGTCGAAAATATGGAGCAGTATTACGACAAAATTGCCTTCTCTGACTGGACTAACTCCCTGTCGTAAACGCAAATGCTGAAAGCGCAGCACCCGGAATAAGAAACCTGGACAGCGGGCATGCACGGTAAAAACAACGGGACCTGt                        >  2:250444/1‑251 (MQ=255)
                                                                                                                                                                                                                                  ttCGACCAGCAATCGATGATTTGCGGTCAGTGCCATGTGGAGTATTACTTCGACGGCAAAAACAAAGCGGTTAAATTCCCGTGGGATGACGGCATGAAAGTCGAAAATATGGAGCAGTATTACGACAAAATTGCCTTCTCTGACTGGACTAACTCCCTGTCGAAAACGCCAATGCTGAAAGCGCAGCACCCGGAATATGAAACCTGGACAGCGGGCATTCACGGTAAAAACAACGTGACCTGTATCGACTg                >  1:195620/1‑251 (MQ=255)
                                                                                                                                                                                                                                   tCGACCAGCAATCGATGATTTGCGGTCAGTGCCATGTGGAGTATTACTTCGACGGCAAAAACAAAGCGGTTAAATTCCCGTGGGATGACGGCATGAAAGTCGAAAATATGGAGCAGTATTACGACAAAATTGCCTTCTCTGACTGGACTAACTCCCTGTCGAAAACGCCAATGCTGAAAGCGCAGCACCCGGAATATGAAACCTGGAc                                                          <  1:391011/208‑1 (MQ=255)
                                                                                                                                                                                                                                   tCGACCAGAAATCGATGATTTGCGGTCAGTGCCATGTGGAGTATTACTTCGACGGCAAAAACAAAGCGGTTAAATTCCCGTGGGATGACGGCATGAAAGTCGAAAATATGGAGCAGTATTACGACAACATTGCCTTCTCTGACTGGACTAACTCCCTGTCGAAAACGCCAATGCTGAAAGCGCAGCACCCGGAATATGAAACCTGGAc                                                          >  2:391011/1‑208 (MQ=255)
                                                                                                                                                                                                                                     gACCAGCAATCGATGGTTTGCGGTCAGTGCCATGTGGAGTATTACTTCGCCGGCAAAAACAAAGCGGTTAAATTCCCGTGGGATGACGGAATGAAAGTCGAAAATATGGAGCAGAATTACGACAAAATGGCCTTCTCTGACTGGACTAACGAACTGTCGCAAACGCCAATTCTGACAGCGAAGCACCCGGAATATGAAACCAGGACAGCGGGCATTCACGGGAAAAACAACGTGAccggaatggacggccc             >  2:131585/1‑237 (MQ=255)
                                                                                                                                                                                                                                      aCCAGCAATCGATGATTTGCGGTCAGTGCCATGTGGAGTATTACTTCGACGGCAAAAACAAAGCGGTTAAATTCCCGTGGGATGACGGCATGAAAGTCGAAAATATGGAGCAGTATTACGACAAACTTGCCTTCTCTGACTGGACTAACTCCCTGTCGAAAACGCCAATGCTGAAAGCGCAGCACCCGGCATATGAAACCTGGACAGCGGGCATTCACGGTAAAAACAACGTGACCTGTATCGACTGcaat            >  1:645524/1‑248 (MQ=255)
                                                                                                                                                                                                                                       ccAGCAATCGATGATTTGCGGTCAGTGCCATGTGGAGTATTACTTCGACGGCAAAAACAAAGCGGTTAAATTCCCGAGGGATGACGGCATGAAAGTCGAAAATATGGAGCAGTATTACGACAAAATTGCCTTCTCTGACTGGACTAACTCCCTGTCGAAAACGCCAATGCTGAAAGCGCAGCACCCGGAATATGAAACCTGGACAGCGGGCATTCACGGTAAAAACAACGTGACCTGTATCGACTGCCat            >  1:508248/1‑250 (MQ=255)
                                                                                                                                                                                                                                                gatgatTTGCGGTCAGTGCCATGTGGAGTATTACTTCGACGGCAAAAACAAAGCGGTTAAATTCCCGTGGGATGACGGCATGAAAGTCGAAAATATGGAGCAGTATTACGACAAAATTGCCTTCTCTGACTGGACTAACTCCCTGTCGAAAACGCCAATGCTGAAAGCGCAGCACCCGGAATATGAAACCTGGACAGCGGGCATTCACGGTAAAAACAACGTGACCTGTATCGACTGCCATATGCCAAAAg  <  1:325534/251‑1 (MQ=255)
                                                                                                                                                                                                                                                gatgatTTGCGGTCAGTGCAATGTGGAGTATTACTTCGACGGCAAAAACAAAGCGGTTAAATTCCCGGGGGGTGACGGCATGAAAGGCGAAAATATGGAGCAGTATTACGACAAAATTGCCTTCTCTGACTGGACTAACTCCATGTCGAAACCGCCAATGCTGAAAGCGCAGCACCCGGAATATGAAACCTGGACAGCGGGCATTCACGGAAAAAACAACGTGACCTGTATCGACTGCCATATGACaaagg  >  2:72431/1‑249 (MQ=255)
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                      
AAGCCCGGATGTGGCGCGTCTGATCCAGAAAGACGGCGAAGATGGCTACTTCCACGGTAAATGGGCGCGCGGCGGTCCGGAAATCGTCAACAACTTAGGTTGTGCCGATTGCCATAACACCGCCTCTCCAGAGTTCGCCAAAGGCAAACCGGAGTTAACCCTTTCCCGTCCGTATGCGGCTCGCGCGATGGAAGCCATTGGTAAACCTTTTGAGAAAGCCGGACGTTTCGACCAGCAATCGATGGTTTGCGGTCAGTGCCATGTGGAGTATTACTTCGACGGCAAAAACAAAGCGGTTAAATTCCCGTGGGATGACGGCATGAAAGTCGAAAATATGGAGCAGTATTACGACAAAATTGCCTTCTCTGACTGGACTAACTCCCTGTCGAAAACGCCAATGCTGAAAGCGCAGCACCCGGAATATGAAACCTGGACAGCGGGCATTCACGGTAAAAACAACGTGACCTGTATCGACTGCCATATGCCAAAAG  >  NC_000913/4288141‑4288631

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