Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,252,805 C→T 5.9% T580T (ACC→ACT dhaR → dhaKLM operon transcription activator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,252,8050CT 5.9% 238.4 / 4.0 68T580T (ACC→ACTdhaRdhaKLM operon transcription activator
Reads supporting (aligned to +/- strand):  ref base C (43/21);  new base T (2/2);  total (45/23)
Fisher's exact test for biased strand distribution p-value = 5.99e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01

TGCATGGCCAGGCAACGATTTTGAACTTTACAGCGTCATCGAGAATCTTGCTCTGAGTAGTGATAACGGGCGCATTCGCGTCAGTGATTTGCCGGAACATCTGTTTACCGAGCAGGCGACAGATGATGTCAGCGCCACCCGCCTTTCCACCAGTCTGTCATTTGCGGAAGTTGAAAAAGAGGCAATTATTAACGCAGCCCAGGTCACAGGCGGTCGCATTCAGGAAATGTCGGCTTTACTTGGGATCGGCCGCACTACGCTGTGGCGGAA  >  NC_000913/1252667‑1252936
                                                                                                                                          |                                                                                                                                   
tgCATGGCCAGGCAACGATTTTGAACTTTACAGCGTCATCGAGAATCTTGCTCTGAGTAGTGATAACGGGCGCATTCGCGTCAGTGATTTGCCGGAACATCTGTTTACCGAGCAGGCGACAGATGATGTCAGCGCCAccc                                                                                                                                    <  1:1544795/140‑1 (MQ=255)
 gCATGGCCAGGCAACGATTTTGAACTTTACAGCGTCATCGAGAATCTTGCTCTGAGTAGTGATAACGGGCGCATTCGCGTCAGTGATTTGCCGGAACATCTGTTTACCGAGCAGGCGACAGATGATGTCAGCGCCACCCg                                                                                                                                   >  1:765795/1‑140 (MQ=255)
  cATGGCCTGGCAACGATTTTGAACTTTACAGCGTCATCGAGAATCTTGCTCTGAGTAGTGATAACGGGCGCATTCGCGTCAGTGATTTGCCGGAACATCTGTTTACCGAGCAGGCGACAGATGATGTCAGCGCCACCCGc                                                                                                                                  <  2:1650770/140‑1 (MQ=255)
  cATGGCCAGGCAACGATTTTGAACTTTACAGCGTCATCGAGAATCTTGCTCTGAGTAGTGATAACGGGCGCATTCGCGTCAGTGATTTGCCGGAACATCTGTTTACCGAGCAGGCGACAGATGATGTCAGCGCCACCCGc                                                                                                                                  >  2:1287924/1‑140 (MQ=255)
  cATGGCCAGGCAACGATTTTGAACTTTACAGCGTCATCGAGAATCTTGCTCTGAGTAGTGATAACGGGCGCATTCGCGTCAGTGATTTGCCGGAACATCTGTTTACCGAGCAGGCGACAGATGATGTCAGCGCCACCCGc                                                                                                                                  <  2:260395/140‑1 (MQ=255)
  cATGGCCAGGCAACGATTTTGAACTTTACAGCGTCATCGAGAATCTTGCTCTGAGTAGTGATAACGGGCGCATTCGCGTCAGTGATTTGCCGGAACATCTGTTTACCGAGCAGGCGACAGATGATGTCAGCGCCACCCGc                                                                                                                                  >  2:161267/1‑140 (MQ=255)
  cATGGCCAGGCAACGATTTTGAACTTTACAGCGTCATCGAGAATCTTGCTCTGAGTAGTGATAACGGGCGCATTCGCGTCAGTGATTTGCCGGAACATCTGTTTACCGAGCAGGCGACAGATGATGTCAGCGCCACCCGc                                                                                                                                  <  1:1397547/140‑1 (MQ=255)
  cATGGCCAGGCAACGATTTTGAACTTTACAGCGTCATCGAGAATCTTGCTCTGAGTAGTGATAACGGGCGCATTCGCGTCAGTGATTTGCCGGAACATCTGTTTACCGAGCAGGCGACAGATGATGTCAGCGCCACCCGc                                                                                                                                  >  2:1288913/1‑140 (MQ=255)
  cATGGCCAGGCAACGATTTTGAACTTTACAGCGTCATCGAGAATCTTGCTCTGAGTAGTGATAACGGGCGCATTCGCGTCAGTGATTTGCCGGAACATCTGTTTACCGAGCAGGCGACAGATGATGTCAGCGCCACCCGc                                                                                                                                  >  2:1287974/1‑140 (MQ=255)
      gCCAGGCAACGATTTTGAACTTTACAGCGTCATCGAGAATCTTGCTCTGAGTAGTGATAACGGGCGCATTCGCGTCAGTGATTTGCCGGAACATCTGTTTACCGAGCAGGCGACAGATGATGTCAGCGCCACCCGCCttt                                                                                                                              <  1:332169/140‑1 (MQ=255)
          ggCAACGATTTTGAACTTTACAGCGTCATCGAGAATCTTGCTCTGAGTAGTGATAACGGGCGCATTCGCGTCAGTGATTTGCCGGAACATCTGTTTACCGAGCAGGCGACAGATGATGTCAGCGCCACCCGCCTTTccac                                                                                                                          >  1:1355965/1‑140 (MQ=255)
          ggCAACGATTTTGAACTTTACAGCGTCATCGAGAATCTTGCTCTGAGTAGTGATAACGGGCGCATTCGCGTCAGTGATTTGCCGGAACATCTGTTTACCGAGCAGGCGACAGATGATGTCAGCGCCACCCGCCTTTccac                                                                                                                          >  1:513134/1‑140 (MQ=255)
          ggCAACGATTTTGAACTTTACAGCGTCATCGAGAATCTTGCTCTGAGTAGTGATAACGGGCGCATTCGCGTCAGTGATTTGCCGGAACATCTGTTTACCGAGCAGGCGACAGATGATGTCAGCGCCACCCGCCTTTccac                                                                                                                          >  1:183821/1‑140 (MQ=255)
          ggCAACGATTTTGAACTTTACAGCGTCATCGAGAATCTTGCTCTGAGTAGTGATAACGGGCGCATTCGCGTCAGTGATTTGCCGGAACATCTGTTTACCGAGCAGGCGACAGATGATGTCAGCGCCACCCGCCTTTccac                                                                                                                          >  1:280694/1‑140 (MQ=255)
          ggCAACGATTTTGAACTTTACAGCGTCATCGAGAATCTTGCTCTGAGTAGTGATAACGGGCGCATTCGCGTCAGTGATTTGCCGGAACATCTGTTTACCGAGCAGGCGACAGATGATGTCAGCGCCACCCGCCTTTccac                                                                                                                          >  1:280779/1‑140 (MQ=255)
               cGATTTTGAACTTTACAGCGTCATCGAGAATCTTGCTCTGAGTAGTGATAACGGGCGCATTCGCGTCAGTGATTTGCCGGAACATCTGTTTACCGAGCAGGCGACAGATGATGTCAGCGCCACCCGCCTTTCCACCAGTc                                                                                                                     >  1:1688127/1‑140 (MQ=255)
                 aTTTTGAACTTTACAGCGTCATCGAGAATCTTGCTCTGAGTAGTGATAACGGGCGCATTCGCGTCAGTGATTTGCCGGAACATCTGTTTACCGAGCAGGCGACAGATGATGTCAGCGCCACCCGCCTTTCCACCAGTCTg                                                                                                                   <  2:1416982/140‑1 (MQ=255)
                 aTTTTGAACTTTACAGCGTCATCGAGAATCTTGCTCTGAGTAGTGATAACGGGCGCATTCGCGTCAGTGATTTGCCGGAACATCTGTTTACCGAGCAGGCGACAGATGATGTCAGCGCCACCCGCCTTTCCACCAGTCTg                                                                                                                   <  2:1007109/140‑1 (MQ=255)
                    ttGAACTTTACAGCGTCATCGAGAATCTTGCTCTGAGTAGTGATAACGGGCGCATTCGCGTCAGTGATTTGCCGGAACATCTGTTTACCGAGCAGGCGACAGATGATGTCAGCGCCACCCGCCTTTCCACCAGTCTGTCa                                                                                                                <  2:314017/140‑1 (MQ=255)
                    ttGAACTTTACAGCGTCATCGAGAATCTTGCTCTGAGTAGTGATAACGGGCGCATTCGCGTCAGTGATTTGCCGGAACATCTGTTTACCGAGCAGGCGACAGATGATGTCAGCGCCACCCGCCTTTCCACCAGTCTGTCa                                                                                                                <  2:1089376/140‑1 (MQ=255)
                     tGAACTTTACAGCGTCATCGAGAATCTTGCTCTGAGTAGTGATAACGGGCGCATTCGCGTCAGTGATTTGCCGGAACATCTGTTTACCGAGCAGGCGACAGATGATGTCAGCGCCACCCGCCTTTCCACCAGTCTGTCAt                                                                                                               >  1:1424082/1‑140 (MQ=255)
                     tGAACTTTACAGCGTCATCGAGAATCTTGCTCTGAGTAGTGATAACGGGCGCATTCGCGTCAGTGATTTGCCGGAACATCTGTTTACCGAGCAGGCGACAGATGATGTCAGCGCCACCCGCCTTTCCACCAGTCTGTCAt                                                                                                               >  2:940898/1‑140 (MQ=255)
                            tACAGCGTCATCGAGAATCTTGCTCTGAGTAGTGATAACGGGCGCATTCGCGTCAGTGATTTGCCGGAACATCTGTTTACCGAGCAGGCGACAGATGATGTCAGCGCCACCCGCCTTTCCACCAGTCTGTCATTTGCGGa                                                                                                        <  2:724711/140‑1 (MQ=255)
                                   tCATCGAGAATCTTGCTCTGAGTAGTGATAACGGGCGCATTCGCGTCAGTGATTTGCCGGAACATCTGTTTACCGAGCAGGCGACAGATGATGTCAGCGCCACCCGCCTTTCCACCAGTCTGTCATTTGCGGAAGTTGaa                                                                                                 >  2:13263/1‑140 (MQ=255)
                                      tCGAGAATCTTGCTCTGAGTAGTGATAACGGGCGCATTCGCGTCAGTGATTTGCCGGAACATCTGTTTACCGAGCAGGCGACAGATGATGTCAGCGCCACCCGCCTTTCCACCAGTCTGTCATTTGCGGAAGTTGaaaaa                                                                                              >  1:656086/1‑140 (MQ=255)
                                      tCGAGAATCTTGCTCTGAGTAGTGATAACGGGCGCATTCGCGTCAGTGATTTGCCGGAACATCTGTTTACCGAGCAGGCGACAGATGATGTCAGCGCCACCCGCCTTTCCACCAGTCTGTCATTTGCGGAAGTTGaaaaa                                                                                              >  1:1294/1‑140 (MQ=255)
                                      tCGAGAATCTTGCTCTGAGTAGTGATAACGGGCGCATTCGCGTCAGTGATTTGCCGGAACATCTGTTTACCGAGCAGGCGACAGATGATGTCAGCGCCACCCGCCTTTCCACCAGTCTGTCATTTGCGGAAGTTGaaaaa                                                                                              >  1:1116145/1‑140 (MQ=255)
                                           aaTCTTGCTCTGAGTAGTGATAACGGGCGCATTCGCGTCAGTGATTTGCCGGAACTTCTGTTTACCGAGCAGGCGACAGATGATGTCAGCGCCACCCGCCTTTCCACCAGTCTGTCATTTGCGGAAGTTGAAAAAGAgga                                                                                         >  2:600113/1‑139 (MQ=255)
                                            aTCTTGCTCTGAGTAGTGATAACGGGCGCATTCGCGTCAGTGATTTGCCGGAACATCTGTTTACCGAGCAGGCGACAGATGATGTCAGCGCCACCCGCCTTTCCACCAGTCTGTCATTTGCGGAAGTTGAAAAAGAGGCa                                                                                        <  2:816298/140‑1 (MQ=255)
                                            aTCTTGCTCTGAGTAGTGATAACGGGCGCATTCGCGTCAGTGATTTGCCGGAACATCTGTTTACCGAGCAGGCGACAGATGATGTCAGCGCCACCCGCCTTTCCACCAGTCTGTCATTTGCGGAAGTTGAAAAAGAGGCa                                                                                        <  1:944818/140‑1 (MQ=255)
                                            aTCGTGCTCTGAGTAGTGATAACGGGCGCCTTCGCGTCAGTGAATTGCCGGAACATCTGTTTTCCGAGCAGGCGACAGATGAGGTCAGCGCCACCCGCCTTTCCACCAGACTGACATTTGCTGAAGTTGAAAAAGAGGCa                                                                                        <  2:1661648/140‑1 (MQ=255)
                                             tCTTGCTCTGAGTAGTGATAACGGGCGCATTCGCGTCAGTGATTTGCCGGAACATCTGTTTACCGAGCAGGCGACAGATGATGTCAGCGCCACCCGCCTTTCCACCAGTCTGTCATTTGCGGAAGTTGAAAAAGAGGCaa                                                                                       >  2:1169293/1‑140 (MQ=255)
                                                 gCTCTGAGTAGTGATAACGGGCGCATTCGCGTCAGTGATTTGCCGGAACATCTGTTTACCGAGCAGGCGACAGATGATGTCAGCGCCACCCGCCTTTCCACCAGTCTGTCATTTGCGGAAGTTGAAAAAGAGGCActtat                                                                                   >  1:742396/1‑140 (MQ=255)
                                                      gAGTAGTGATAACGGGCGCATTCGCGTCAGTGATTTGCCGGAACATCTGTTTACCGAGCAGGCGACAGATGATGTCAGCGCCACCCGCCTTTCCACCAGTCTGTCATTTGCGGAAGTTGAAAAAGAGGCAATTATTAACg                                                                              <  2:894777/140‑1 (MQ=255)
                                                      gAGTAGTGATAACGGGCGCATTCGCGTCAGTGATTTGCCGGAACATCTGTTTACCGAGCAGGCGACAGATGATGTCAGCGCCACCCGCCTTTCCACCAGTCTGTCATTTGCGGAAGTTGAAAAAGAGGCAATTATTAACg                                                                              <  2:703595/140‑1 (MQ=255)
                                                      gAGTAGTGATAACGGGCGCATTCGCGTCAGTGATGTGCCGGAACATCTGTTTACCGAGCAGGCGACAGATGATGTCAGCGCCACCCGCCTTTCCACCAGTCTGTCATTTGCGGAAGTTGAAAAAGAGGCAATTATTAACg                                                                              <  2:371575/140‑1 (MQ=255)
                                                       agtagtGATAACGGGCGCATTCGCGTCAGTGATTTGCCGGAACATCTGTTTACCGAGCAGGCGACAGATGATGTCAGCGCCACCCGCCTTTCCACCAGTCTGTCATTTGCGGAAGTTGAAAAAGAGGCAATTATTAACGc                                                                             >  2:452195/1‑140 (MQ=255)
                                                       agtagtGATAACGGGCGCATTCGCGTCAGTGATTTGCCGGAACATCTGTTTACCGAGCAGGCGACAGATGATGTCAGCGCCACCCGCCTTTCCACCAGTCTGTCATTTGCGGAAGTTGAAAAAGAGGCAATTATTAACGc                                                                             >  2:661725/1‑140 (MQ=255)
                                                           gtGATAACGGGCGCATTCGCGTCAGTGATTTGCCGGAACATCTGTTTACCGAGCAGGCGACAGATGATGTCAGCGCCACCCGCCTTTCCACCAGTCTGTCATTTGCGGAAGTTGAAAAAGAGGCAATTATTAACGCAGcc                                                                         <  2:258852/140‑1 (MQ=255)
                                                           gtGATAACGGGCGCATTCGCGTCAGTGATTTGCCGGAACATCTGTTTACCGAGCAGGCGACAGATGATGTCAGCGCCACCCGCCTTTCCACCAGTCTGTCATTTGCGGAAGTTGAAAAAGAGGCAATTATTAACGCAGcc                                                                         <  1:499474/140‑1 (MQ=255)
                                                           gtGATAACGGGCGCATTCGCGTCAGTGATTTGCCGGAACATCTGTTTACCGAGCAGGCGACAGATGATGTCAGCGCCACCCGCCTTTCCACCAGTCTGTCATTTGCGGAAGTTGAAAAAGAGGCAATTATTAACGCAGcc                                                                         <  1:1262280/140‑1 (MQ=255)
                                                            tGATAACGGGCGCATTCGCGTCAGTGATTTGCCGGAACATCTGTTTACCGAGCAGGCGACAGATGATGTCAGCGCCACCCGCCTTTCCACCAGTCTGTCATTTGCGGAAGTTGAAAAAGAGGCAATTATTAACGCAGccc                                                                        >  1:850880/1‑140 (MQ=255)
                                                                aaCGGGCGCATTCGCGTCAGTGATTTGCCGGAACATCTGTTTACCGAGCAGGCGACAGATGATGTCAGCGCCACCCGCCTTTCCACCAGTCTGTCATTTGCGGAAGTTGAAAAAGAGGCAATTATTAACGCAGCCCAGGt                                                                    <  1:35682/140‑1 (MQ=255)
                                                                 actggCGCATTCGCGTCAGTGATTTGCCGGAACATCTGTTTACCGAGCAGGCGACAGATGATGTCAGCGCCACCCGCCTTTCCACCAGTCTGTCATTTGCGGAAGTTGAAAAAGAGGCAATTATTAACTCAGCCCAGGTc                                                                   >  2:865182/4‑140 (MQ=255)
                                                                 aCGGGCGCATTCGCGTCAGTGATTTGCCGGAACATCTGTTTACCGAGCAGGCGACAGATGATGTCAGCGCCACCCGCCTTTCCACCAGTCTGTCATTTGCGGAAGTTGAAAAAGAGGCAATTATTAACGCAGCCCAGTTc                                                                   >  2:493845/1‑140 (MQ=255)
                                                                 aCGGGCGCATTCGCGTCAGTGATTTGCCGGAACATCTGTTTACCGAGCAGGCGACAGATGATGTCAGCGCCACCCGCCTTTCCACCAGTCTGTCATTTGCGGAAGTTGAAAAAGAGGCAATTATTAACGCAGCCCAGGTc                                                                   >  2:59053/1‑140 (MQ=255)
                                                                 aCGGGCGCATTCGCGTCAGTGATTTGCCGGAACATCTGTTTACCGAGCAGGCGACAGATGATGTCAGCGCCACCCGCCTTTCCACCAGTCTGTCATTTGCGGAAGTTGAAAAAGAGGCAATTATTAACGCAGCCCAGGTc                                                                   >  2:828253/1‑140 (MQ=255)
                                                                 aCGGGCGCATTCGCGTCAGTGATTTGCCGGAACATCTGTTTACCGAGCAGGCGACAGATGATGTCAGCGCCACCCGCCTTTCCACCAGTCTGTCATTTGCGGAAGTTGAAAAAGAGGCAATTATTAACGCAGCCCAGGTc                                                                   >  1:1123522/1‑140 (MQ=255)
                                                                 aCGGGCGCATTCGCGTCAGTGATTTGCCGGAACATCTGTTTACCGAGCAGGCGACAGATGATGTCAGCGCCACCCGCCTTTCCACCAGTCTGTCATTTGCGGAAGTTGAAAAAGAGGCAATTATTAACGCAGCCCAGGTc                                                                   >  2:1059181/1‑140 (MQ=255)
                                                                 aCGGGCGCATTCGCGTCAGTGATTTGCCGGAACATCTGTTTACCGAGCAGGCGACAGATGATGTCAGCGCCACCCGCCTTTCCACCAGTCTGTCATTTGCGGAAGTTGAAAAAGAGGCAATTATTAACGCAGCCCAGGTc                                                                   >  1:487349/1‑140 (MQ=255)
                                                                 aCGGGCGCATTCGCGTCAGTGATTTGCCGGAACATCTGTTTACCGAGCAGGCGACAGATGATGTCAGCGCCACCCGCCTTTCCACCAGTCTGTCATTTGCGGAAGTTGAAAAAGAGGCAATTATTAACGCAGCCCAGGTc                                                                   >  2:1586166/1‑140 (MQ=255)
                                                                 aCGGGCGCATTCGCGTCAGTGATTTGCCGGAACATCTGTTTACCGAGCAGGCGACAGATGATGTCAGCGCCACCCGCCTTTCCACCAGTCTGTCATTTGCGGAAGTTGAAAAAGAGGCAATTATTAACGCAGCCCAGGTc                                                                   >  2:1419187/1‑140 (MQ=255)
                                                                      cgcATTCGCGTCAGTGATTTGCCGGAACATCTGTTTACCGAGCAGGCGACAGATGATGTCAGCGCCACCCGCCTTTCCACCAGTCTGTCATTTGCGGAAGTTGAAAAAGAGGCAATTATTAACGCAGCCCAGGTCACAgg                                                              >  1:476489/1‑140 (MQ=255)
                                                                      cgcATTCGCGTCAGTGATTTGCCGGAACATCTGTTTACCGAGCAGGCGACAGATGATGTCAGCGCCACCCGCCTTTCCACCAGTCTGTCATTTGCGGAAGTTGAAAAAGAGGCAATTATTAACGCAGCCCAGGTCACAgg                                                              >  2:497441/1‑140 (MQ=255)
                                                                      cgcATTCGCGTCAGTGATTTGCCGGAACATCTGTTTACCGAGCAGGCGACAGATGATGTCAGCGCCACCCGCCTTTCCACCAGTCTGTCATTTGCGGAAGTTGAAAAAGAGGCAATTATTAACGCAGCCCAGGTCACAgg                                                              >  2:1256462/1‑140 (MQ=255)
                                                                        ccttCGCGTCAGGGCTTTGCCGGAACATCTGTTTACCGAGCAGGCGACAGATGAGGTCAGCGCCACCCGCCTTTCCACCAGTCTGGCATTTGCGGACGTGGAACTAGAGGCACTTAGTAACGCAGCCCAGGACACAGGCg                                                            >  2:1691404/3‑140 (MQ=255)
                                                                                          gCCGGTACATCTGTTTACCGAGCAGGCGACAGTTGATGTCAGCGCCACCCGCCTTTCCACCAGTCTGTCATTTGCGGAAGTTGAAAAAGAGGCAATTATTAACGCAGCCCAGGTCACAGGCGGTCGCATTCAGGAAATGt                                          <  2:243065/140‑1 (MQ=255)
                                                                                              gAACATCTGTTTACCGAGCAGGCGACAGATGATGTCAGCGCCACCCGCCTTTCCACCAGTCTGTCATTTGCGGAAGTTGAAAAAGAGGCAATTATTAACGCAGCCCAGGTCACAGGCGGTCGCATTCAGGAAATGTCGGc                                      >  1:1655428/1‑140 (MQ=255)
                                                                                                         tACCGAGCAGGCGACAGATGATGTCAGCGCCACCCGCCTTTCCACCAGTCTGTCATTTGCGGAAGTTGAAAAAGAGGCAATTATTAACGCAGCCCAGGTCACAGGCGGTCGCATTCAGGAAATGTCGGCTTTACTTGGGa                           >  2:538973/1‑140 (MQ=255)
                                                                                                             gAGCAGGCGACAGATGATGTCAGCGCCACCCGCCTTTCCACCAGTCTGTCATTTGCGGAAGTTGAAAAAGAGGCAATTATTAACGCAGCCCAGGTCACAGGCGGTCGCATTCAGGAAATGTCGGCTTTACTTGGGATCgg                       <  1:1066414/140‑1 (MQ=255)
                                                                                                              aGCAGGCGACAGATGATGTCAGCGCCACCCGCCTTTCCACCAGTCTGTCATTTGCGGAAGTTGAAAAAGAGGCAATTATTAACGCAGCCCAGGTCACAGGCGGTCGCATTCAGGAAATGTCGGCTTTACTTGGGATCGGc                      >  1:639645/1‑140 (MQ=255)
                                                                                                              aGCAGGCGACAGATGATGTCAGCGCCACCCGCCTTTCCACCAGTCTGTCATTTGCGGAAGTTGAAAAAGAGGCAATTATTAACGCAGCCCAGGTCACAGGCGGTCGCATTCAGGAAATGTCGGCTTTACTTGGGATCGGc                      >  1:639653/1‑140 (MQ=255)
                                                                                                                 aGGCGACAGATGATGTCAGCGCCACCCGCCTTTCCACCAGTCTGTCATTTGCGGAAGTTGAAAAAGAGGCAATTATTAACGCAGCCCAGGTCACAGGCGGTCGCATTCAGGAAATGTCGGCTTTACTTGGGATCGGCCGc                   >  2:1600608/1‑140 (MQ=255)
                                                                                                                                  aGCGCCACTCGCCTTTCCACCAGTCTGTCATTTGCGGAAGTTGAAAAAGAGGCAATTATTAACGCAGCCCAGGTCACAGGCGGTCGCATTCAGGAAATGTCGGCTTTACTTGGGATCgg                       >  2:945718/1‑119 (MQ=255)
                                                                                                                                  aGCGCCACTCGCCTTTCCACCAGTCTGTCATTTGCGGAAGTTGAAAAAGAGGCAATTATTAACGCAGCCCAGGTCACAGGCGGTCGCATTCAGGAAATGTCGGCTTTACTTGGGATCgg                       >  2:948460/1‑119 (MQ=255)
                                                                                                                                  aGCGCCACTCGCCTTTCCACCAGTCTGTCATTTGCGGAAGTTGAAAAAGAGGCAATTATTAACGCAGCCCAGGTCACAGGCGGTCGCATTCAGGAAATGTCGGCTTTACTTGGGATCgg                       <  1:945718/119‑1 (MQ=255)
                                                                                                                                  aGCGCCACTCGCCTTTCCACCAGTCTGTCATTTGCGGAAGTTGAAAAAGAGGCAATTATTAACGCAGCCCAGGTCACAGGCGGTCGCATTCAGGAAATGTCGGCTTTACTTGGGATCgg                       <  1:948460/119‑1 (MQ=255)
                                                                                                                                  aGCGCCACCCGCCTTTCCACCAGTCTGTCATTTGCGGAAGTTGAAAAAGAGGCAATTATTAACGCAGCCCAGGTCACAGGCGGTCGCATTCAGGAAATGTCGGCTTTACTTGGGATCGGCCGCACTACGCTGTGGCGGaa  >  1:1600309/1‑140 (MQ=255)
                                                                                                                                  aGCGCCACCCGCCTTTCCACCAGTCTGTCATTTGCGGAAGTTGAAAAAGAGGCAATTATTAACGCAGCCCAGGTCACAGGCGGTCGCATTCAGGAAATGTCGGCTTTACTTGGGATCGGCCGCACTACGCTGTGGCGGaa  >  1:188394/1‑140 (MQ=255)
                                                                                                                                          |                                                                                                                                   
TGCATGGCCAGGCAACGATTTTGAACTTTACAGCGTCATCGAGAATCTTGCTCTGAGTAGTGATAACGGGCGCATTCGCGTCAGTGATTTGCCGGAACATCTGTTTACCGAGCAGGCGACAGATGATGTCAGCGCCACCCGCCTTTCCACCAGTCTGTCATTTGCGGAAGTTGAAAAAGAGGCAATTATTAACGCAGCCCAGGTCACAGGCGGTCGCATTCAGGAAATGTCGGCTTTACTTGGGATCGGCCGCACTACGCTGTGGCGGAA  >  NC_000913/1252667‑1252936

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