Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,939,836 A→C F259L (TTT→TTG hsrA ← putative multidrug or homocysteine efflux system

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,939,8360AC98.4% 211.2 / ‑3.2 63F259L (TTT→TTGhsrAputative multidrug or homocysteine efflux system
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base C (34/28);  total (34/29)
Fisher's exact test for biased strand distribution p-value = 4.60e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.55e-01

TCACCGTGATCCCCACTAACGTATGGCGATAGCCCAGACGACGTAAGACTTGGGTAACCATCGATTTTGCAATAATGGAACCTAACGCTGTCGGTGCCATCATACAGCCAGCAATAAACGCCTGATAACCAAATCCTACCTGTAACATCAATGGCATAAGGAACGGTACACAGCCGGTCCCCAGACGGGTTGCAATATTGCCTACGATACCGATCGAGAAAGTGCGGGTTTTAAATAAATCTAATGAAATTAATGGGTTTGGCGTGCGTCGTGCATGGAGAATATAGAGAAGCAGTAACCCGATGCTGGTGACAATTACCGTCAAGGCAATCCAGCTGGCGACAATCTTTTCCCCGAATAGCTCTATTCCGCTTGAGAAGAGAACAAGGCTGAGGCCAAACAGCAAAAAGCCAGTGATATCGAATCTGCGTCGTGCGGTGGTGAAATTGGGC  >  NC_000913/3939604‑3940055
                                                                                                                                                                                                                                        |                                                                                                                                                                                                                           
tCACCGTGATCCCCACTAACGTATGGCGATAGCCCAGACGACGTAAGACTTGGGTAACCATCGATTTTGCAATAATGGAACCTAACGCTGTCGGTGCCATCATACAGCCAGCAATAAACGCCTGATAACCAAATCCTACCTGTAACATCAATGGCATAAGGAACGGTACACAGCCGGTCCCCAGACGGGTTGCAATATTGCCTACGATACCGATCGAGAAAGTGCGGGTTTtcaat                                                                                                                                                                                                                          <  2:345126/236‑1 (MQ=255)
  aCCGTGATCCCCACTAACGTATGGCGATAGCCCAGACGACGTAAGACTTGGGTAACCATCGATTTTGCAATAATGGAACCTAACGCTGTCGGTGCCATCATACAGCCAGCAATAAACGCCTGATAACCAAATCCTACCTGTAACATCAATGGCATAAGGAACGGTACACAGCCGGTCCCCAGACGGGTTGCAATATTGCCTACGATACCGATCGAGAAAGTGCGGGTTTtcaataa                                                                                                                                                                                                                        >  1:526414/1‑236 (MQ=255)
      tGATCCCCACTAACGTATGGCGATAGCCCAGACGACGTAAGACTTGGGTAACCATCGATTTTGCAATAATGGAACCTAACGCTGTCGGTGCCATCATACAGCCAGCAATAAACGCCTGATAACCAAATCCTACCTGTATCATCAATGGCATAAGGAACGGTACACAGCCGGTCCCCAGACGGGTTGCAATATTGCCTACGATACCGATCGAGAAAGTGCGGGTTTTCAATAAATCt                                                                                                                                                                                                                    >  2:240077/1‑236 (MQ=255)
        aTCCCCACTAACGTATGGCGATAGCCCAGACGACGTAAGACTTGGGTAACCATCGATTTTGCAATAATGGAACCTAACGCTGTCGGTGCCATCATACAGCCAGCAATAAACGCCTGATAACCAAATCCTACCTGTAACATCAATGGCATAAGGAACGGTACACAGCCGGTCCCCAGACGGGTTGCAATATTGCCTACGATACCGATCGGGAAAGTGCGGGGTTTCAATAAATCTaa                                                                                                                                                                                                                  >  1:646538/1‑236 (MQ=255)
          cccaaatAAGGTATGGCAATAGCCCAGACGTAGTAAGACTTGGGTAACCATCGATTTTGCAATAATGGAACCTAACGCTGTCGGTGCCATCATACAGCCAGCAATAAACGCCTGATAACCAAATCCTACCTGTAACATCAATGGCATAAGGAACGGTACACAGCCGGTCCCCAGACGGGTTGCAATATTGCCTACGATACCGATCGAGAAAGTGCGGGTTTTCAATAAATCTAATg                                                                                                                                                                                                                <  2:419340/230‑1 (MQ=255)
                 aaCGTATGGCGATAGCCCAGACGACGTAAGACTTGGGTAACCATCGATTTTGCAATAATGGAACCTAACGCTGTCGGTGCCATCATACAGCCAGCAATAAACGCCTGATAACCAAATCCTACCTGTAACATCAATGGCATAAGGAACGGTACACAGCCGGTCCCCAGACGGGTTGCAATATTGCCTACGATACCGATCGAGAAAGTGCGGGTTTTCAATAAATCTAATGAAATTaa                                                                                                                                                                                                         <  1:112981/236‑1 (MQ=255)
                   cGTATGGCGATAGCCCAGACGACGTAAGACTTGGGTAACCATCGATTTTGCAATATTGGAACCTAACGCTGTCGGTGCCATCATACAGCCAGCAATAAACGCCTGATAACCAAATCCTACCTGTAACATCAATGGCATAAGGAACGGTACACAGCCGGTCCCCAGACGGGTTGCAATATTGCCTACGATAACGATCGAGAAAGTGCGGGTTTTCAATAAATCTAATGAAATTAATg                                                                                                                                                                                                       >  2:716311/1‑236 (MQ=255)
                      tgggCGATAGCCCAGACGACGTAAGGCTTGGGTAACCATCGATTTTGCAATAATGGAACCTAACGCTGTCGGTGCCATCATACAGCCAGCAATAAACGCCTGATAACCAAATCCTACCTGTAACATCAATGGCATAAGGAACGGTACACAGCCGGTCCCCAGACGGGTTGCAATATTGCCTACGATACCGATCGAGAAAGTGCGGGTTTTCAATAAATCTAATGAAATTAATGGGt                                                                                                                                                                                                    <  2:17495/234‑1 (MQ=255)
                      taggCAATAGCCCAGAAGACGTAAGTCTTGGGTAACCATCGATTTTGCAATAATGGAACCTAACGCTGTCGGTGCCATCATACAGCCAGCAATAAACGCCTGATAACCAAATCCTACCTGTAACATCAATGGCATAAGGAACGGTACACAGCCGGTCCCCAGACGGGTTGCAATATTGCCTACGATACCGATCGAGAAAGTGCGGGTTTTCAATAAATCTAATGAAATTAATGGGt                                                                                                                                                                                                    <  2:646538/234‑1 (MQ=255)
                      aTGGCGATAGCCCGGACAACGTAGGACTTGGGTAACCATCGATTTTGCAATAATGGAACCTAACGCTGTCGGTGCCATCATACAGCCAGCAATAAACGCCTGATAACCAAATCCTACCTGTAACATCAATGGCATAAGGAACGGTACACAGCCGGTCCCCAGACGGGTTGCAATATTGCCTACGATACCGATCGAGAAAGTGCGGGTTTTCAATAAATCTAATGAAATTAATGGGt                                                                                                                                                                                                    <  2:722316/236‑1 (MQ=255)
                                cccAGACGACGTAAGACTTGGGTAACCATCGATTTTGCAATAATGGAACCTAACGCTGTCGGTGCCATCATACAGCCAGCAATAAACGCCTGATAACCAAATCCTACCTGTAACATCAATGGCATAAGGAACGGTACACAGCCGGTCCCCAGACGGGTTGCAATATTGCCTACGGTACCGATCGAGAAAGTGCGGGTTTTCAATAAATCTAATGAAATTAATGGGTTTGgcgt                                                                                                                                                                                             >  2:588818/1‑233 (MQ=255)
                                    gacgacGTAAGACTTGGGTAACCATCGATTTTGCAATAATGGAACCTAACGCTGTCGGTGCCATCATACAGCCAGCAATAAACGCCTGATAACCAAATCCTACCTGTAACATCAATGGCATAAGGAACGGTACACAGCCGGTCCCCAGACGGGTTGCAATATTGCCTACGATACCGATCGAGAAAGTGCGGGTTTTAAATAAATCTAATGAAATTAATGGGTTTGgcgtgcg                                                                                                                                                                                          <  2:34143/232‑1 (MQ=255)
                                         cgTAAGACTTGGGTAACCATCGATTTTGCAATAATGGAACCTAACGCTGTCGGTGCCATCATACAGCCAGCAATAAACGCCTGATAACCAAATCCTACCTGTAACATCAATGGCATAAGGAACGGTACACAGCCGGTCCCCAGACGGGTTGCAATATTGCCTACGATACCGATCGAGAAAGTGCGGGTTTTCAATAAATCTAATGAAATTAATGGGTTTGGCGTGCGTCGTGCATg                                                                                                                                                                                 >  1:474453/1‑236 (MQ=255)
                                         cgTAAGACTTGGGTAACCATCGATTTTGCAATAATGGAACCTAACGCTGTCGGTGCCATCATACAGCCAGCAATAAACGCCTGATAACCAAATCCTACCTGTAACATCAATGGCATAAGGAACGGTACACAGCCGGTCCCCAGACGGGTTGCAATATTGCCTACGATACCGATCGAGAAAGGCCGGGTTTTCAATAAATCTAATGAAATTAATTGGTTTGGCGTGCGCCGTGCATg                                                                                                                                                                                 >  1:712633/1‑236 (MQ=255)
                                                   gggTAACCATCGATTTTGCAATAATGGAACCTAACGCTGTCGGTGCCATCATACAGCCAGCAATAAACGCCTGATAACCAAATCCTACCTGTAACATCAATGGCATAAGGAACGGTACACAGCCGGTCCCCAGACGGGTTGCAATATTGCCTACGATACCGATCGAGAAAGTGCGGGTTTTCAATAAATCTAATGAAATTAATGGGTTTGGCGTGCGTCGTGCATg                                                                                                                                                                                 <  2:176847/226‑1 (MQ=255)
                                                   gggTAACCATCGATTTTGCAATAATGGAACCTAACGCTGTCGGTGCCATCATACAGCCAGCAATAAACGCCTGATAACCAAATCCTACCTGTAACATCAATGGCATAAGGAACGGTACACAGCCGGTCCCCAGACGGGTTGCAATATTGCCTACGATACCGATCGAGAAAGTGCGGGTTTTCAATAAATCTAATGAAATTAATGGGTTTGGCGTGCGTCGTGCATg                                                                                                                                                                                 >  1:176847/1‑226 (MQ=255)
                                                   gggTAACCATCGATTTTGCAATAATGGAACCTAACGCTGTCGGTGCCATCATACAGCCAGCAATAAACGCCTGATAACCAAATCCTACCTGTAACATCAATGGCATAAGGAACGGTACACAGCCGGTCCCCAGACGGGTTGCAATATTGCCTACGATACCGATCGAGAAAGTGCGGGTTTTCAATAAATCTAATGAAATTAATGGGTTTGGCGTGCGTCGTGCATGGAGAATATag                                                                                                                                                                       >  1:100249/1‑236 (MQ=255)
                                                      tAACCATCGATTTTGCAATAATGGAACCTAACGCTGTCGGTGCCATCATACAGCCAGCAATAAACGCCTGATAACCAAATCCTACCTGTAACATCAATGGCATAAGGAACGGTACACAGCCGGTCCCCAGACGGGTTGCAATATTGCCTACGATACCGATCGAGAAAGTGCGGGTTTTCAATAAATCTAATGAAATTAATGGGTTTGGCGTGCGTCGTGCATGGAGAATATAgaga                                                                                                                                                                    <  2:526414/236‑1 (MQ=255)
                                                                   tGCAATAATGGAACCTAACGCTGTCGGTGCCATCATACAGCCAGCAATAAACGCCTGATAACCAAATCCTACCTGTAACATCAATGGCATAAGGAACGGCACACAGCCGGTCCCCAGACGGGTTGCAATATTGCCTACGACACCGATCGAGGAAGTCCGGGTTTTCAATAAATCTAATGAAATTAAAGGGATTGACGGTCGTCGTGCATTGCAAATATCTAGAATCatgaaacc                                                                                                                                                         >  2:355669/1‑227 (MQ=255)
                                                                          atGGAACCTAACGCTGTCGGTGCCATCATACAGCCAGCAATAAACGCCTGATAACCAAATCCTACCTGTAACATCAATGGCATAAGGAACGGTACACAGCCGGTCCCCAGACGGGTTGCAATATTGCCTACGATACCGATCGAGAAAGTGCGGGTTTTCAATAAATCTAATGAAATTAATGGGTTTGGCGTGCGTCGTGCATGGAGAATATAGAGAAGCAGTAACCCGATGCTGGt                                                                                                                                                <  2:484952/236‑1 (MQ=255)
                                                                                          tCGGTGCCATCATACAGCCAGCAATAAACGCCTGATAACCAAATCCTACCTGTAACATCAATGGCATAAGGAACGGTACACAGCCGGTCCCCAGACGGGTTGCAATATTGCCTACGATACCGATCGAGAAAGTGCGGGTTTTCAATAAATCTAATGAAATTAATGGGTTTGGCGTGCGTCGTGCATGGAGAATATAGAGAAGCAGTAACCCGATGCTGGTGACAATTACCGTCAAg                                                                                                                                >  2:522296/1‑236 (MQ=255)
                                                                                                  atcatACAGCCAGCAATAAACGCCTGATAACCAAATCCTACCTGTAACATCAATGGCATAAGGAACGGTACACAGCCGGTCCCCAGACGGGTTGCAATATTGCCTACGATACCGATCGAGAAAGTGCGGGTTTTCAATAAATCTAATg                                                                                                                                                                                                                <  1:247135/148‑1 (MQ=255)
                                                                                                  atcatACAGCCAGCAATAAACGCCTGATAACCAAATCCTACCTGTAACATCAATGGCATAAGGAACGGTACACAGCCGGTCCCCAGACGGGTTGCAATATTGCCTACGATACCGATCGAGAAAGTGCGGGTTTTCAATAAATCTAATg                                                                                                                                                                                                                >  2:247135/1‑148 (MQ=255)
                                                                                                             agcAATAAACGCCTGATAACCAAATCCTACCTGTAACATCAATGGCATAAGGAACGGTACACAGCCGGTCCCCAGACGGGTTGCAATATTGCCTACGATACCGATCGAGAAAGTGCGGGTTTTCAATAAATCTAATGAAATTAATGGGTTTGGCGTGCGTCGTGCATGGAGAATATAgaga                                                                                                                                                                    >  1:277966/1‑181 (MQ=255)
                                                                                                             agcAATAAACGCCTGATAACCAAATCCTACCTGTAACATCAATGGCATAAGGAACGGTACACAGCCGGTCCCCAGACGGGTTGCAATATTGCCTACGATACCGATCGAGAAAGTGCGGGTTTTCAATAAATCTAATGAAATTAATGGGTTTGGCGTGCGTCGTGCATGGAGAATATAgaga                                                                                                                                                                    <  2:277966/181‑1 (MQ=255)
                                                                                                             agcAATAAACGCCTGATAACCAAATCCTACCTGTAACATCAATGGCATAAGGAACGGTACACAGCCGGTCCCCAGACGGGTTGCAATATTGCCTACGATACCGATCGAGAAAGTGCGGGTTTTCAATAAATCTAATGAAATTAATGGGTTTGGCGTGCGTCGTGCATGGAGAATATAGAGAAGCAGTAACCCGATGCTGGTGACAATTACCGTCAAGGCAATCCAGCTGGCGACaa                                                                                                             >  1:496770/1‑236 (MQ=255)
                                                                                                             agcAATAAACGCCTGATAACCAAATCCTACCTGTAACATCAATGGCATAAGGAACGGTACACAGCCGGTCCCCAGACGGGTTGCAATATTGCCTACGATACCGATCGAGAAAGTGCGGGTTTTCAATAAATCTAATGAAATTAATGGGTTTGGCGTGCGTCGTGCATGGAGAATATAGAGAAGCAGTAACCCGATGCTGGTGACAATTACCGTCAAGGCAATCCAGCTGGCGACaa                                                                                                             >  1:709284/1‑236 (MQ=255)
                                                                                                             agcAATAAACGCCTGATAACCAAATCCTACCTGTAACATCAATGGCATAAGGAACGGTACACAGCCGGTCCCCAGACGGGTTGCAATATTGCCTACGATACCGATCGAGAAAGTGCGGGTTTTCAATAAATCTAACGAAATTAATGGGTTTGGCGTGCGTCGTGCATGGAGAATATAGAGAAGCAGTAACCCGATGCTGGTGACAATTACCGTCAAGGCAATCCAGCTGGCGCCaa                                                                                                             >  1:75874/1‑236 (MQ=255)
                                                                                                                     aCGCCTGATAACCAAATCCTACCTGTAACATCAATGGCATAAGGAACGGTACACAGCCGGTCCCCAGACGGGTTGCAATATTGCCTACGATACCGATCGAGAAAGTGCGGGTTTTCAATAAATCTAATGAAATTAATGGGTTTGgcgtgcgt                                                                                                                                                                                         <  2:709144/152‑1 (MQ=255)
                                                                                                                     aCGCCTGATAACCAAATCCTACCTGTAACATCAATGGCATAAGGAACGGTACACAGCCGGTCCCCAGACGGGTTGCAATATTGCCTACGATACCGATCGAGAAAGTGCGGGTTTTCAATAAATCTAATGAAATTAATGGGTTTGgcgtgcgt                                                                                                                                                                                         >  1:709144/1‑152 (MQ=255)
                                                                                                                         cTGATAACCAAATCCTACCTGTAACATCAATGGCATAAGGAACGGTACACAGCCGGTCCCCAGACGGGTTGCAATATTGCCTACGATACCGATCGAGAAAGTGCGGGTTTTCAATAAATCTAATGAAATTAATGGGTTTGGCGTGCGTCGTGCATGGAGAATATAGAGAAGCAGTAACCCGATGCTGGTGACAATTACCGTCAAGGCAATCCAGCTGGGCACAAACTTTTcccct                                                                                                  >  2:304573/1‑234 (MQ=255)
                                                                                                                         cTGATAACCAAATCCTACCTGTAACATCAATGGCATAAGGAACGGTACACAGCCGGTCCCCAGACGGGTTGCAATATTGCCTACGATACCGATCGAGAAAGTGCGGGTTTTCAATAAATCTAATGAAATTAATGGGTTTGGCGTGCGTCGTGCATGGAGAATATAGAGAAGCAGTAACCCGATGCTGGTGACAATTACCGTCAAGGCAATCCAGCTGGCGACAATCTTTTCCCCg                                                                                                  >  1:443908/1‑235 (MQ=255)
                                                                                                                          taaTAACCAAATCCTACCTGTAACATCAATGGCATAAGGAACGGTACACAGCCGGTCCCCAGACGGGTTGCAATATTGCCTACGATACCGATCGAGAAAGTGCGGGTTTTCAATAAATCTAATGAAATTAATGGGTTTGGCGTGCGTCGTGCATGGAGAATATAGAGAAGCAGTAACCCGATGCTGGTGACAATTACCGTCAAGGCAATCCAGCTGGCGACAATCTTTTCCCCGaa                                                                                                <  2:296045/234‑1 (MQ=255)
                                                                                                                            aTAACCAAATCCTACCTGTAACATCAATGGCATAAGGAACGGTACACAGCCGGTCCCCAGACGGGTTGCAATATTGCCTACGATACCGATCGAGAAAGTGCGGGTTTTCAATAAATCTAATGAAATTAATGGGTTTGGCGTGCGTCGTGCATGGAGAATATAGAGAAGCAGTAACCCGATGCTGGTGCCAATTACCTTAAAGGCAATACAGCTGGCGACAAACTTTTCCCCgcaac                                                                                              >  1:639096/1‑232 (MQ=255)
                                                                                                                             tAACCAAATCCTACCTGTAACATCAATGGCATAAGGAACGGTACACAGCCGGTCCCCAGACGGGTTGCAATATTGCCTACGATACCGATCGAGAAAGTGCGGGTTTTCAATAAATCTAATGAAATTAATGGGTTTGGCGTGCGTCGTGCATGGAGAATATAGAGAAGCAGTAACCCGATGCTGGTGACAATTACCGTCAAGGCAATCCAGCTGGCGACAATCTTTTCCCCGAATAg                                                                                             <  2:100249/236‑1 (MQ=255)
                                                                                                                             tAACCAAATCCATCCTTGAAGATCAATGGCATAAGGAACGGTACACAGCCGGTCCCCAGACGGGTTGCAATATTGCCTACGATACCGATCGAGAAAGTGCGGGTTTTCAATAAATCTAATGAAATTAATGGGTTTGGCGTGCGTCGTGCATGGAGAATATAGAGAAGCAGTAACCCGATGCTGGTGACAATTACCGTCAAGGCAATCCAGCTGGCGACAATCTTTTCCCCGAATAg                                                                                             <  2:474453/236‑1 (MQ=255)
                                                                                                                                           cTGTAACATCAATGGCATAAGGAACGGTACACAGCCGGTCCCCAGACGGGTTGCAATATTGCCTACGATACCGATCGAGAAAGTGCGGGTTTTCAATAAATCTAATGAAATTAATGGGt                                                                                                                                                                                                    >  2:571071/1‑119 (MQ=255)
                                                                                                                                           cTGTAACATCAATGGCATAAGGAACGGTACACAGCCGGTCCCCAGACGGGTTGCAATATTGCCTACGATACCGATCGAGAAAGTGCGGGTTTTCAATAAATCTAATGAAATTAATGGGt                                                                                                                                                                                                    <  1:571071/119‑1 (MQ=255)
                                                                                                                                               aaCATCAATGGCATAAGGAACGGTACACAGCCGGTCCCCAGACGGGTTGCAATATTGCCTACGATACCGATCGAGAAAGTGCGGGTTTTCAATAAATCTAATGAAATTAATGGGTTTGGCGTGCGTCGTGCATGGAGAATATAGAGAAGCAGTAACCCGATGCTGGTGACAATTACCGTCAAGGCAATCCAGCTGGCGACAATCTTTTCCCCGAATAGCTCTATTCCGCTTgagaa                                                                           <  2:496770/236‑1 (MQ=255)
                                                                                                                                                      aTGGCATAAGGAACGGTACACAGCCGGTCCCCAGACGGGTTGCAATATTGCCTACGATACCGATCGAGAAAGTGCGGGTTTTCAATAAATCTAATGAAATTAATGGGTTTGGCGTGCGTCGTGCATGGAGAATATAGAGAAGCAGTAACCCGATGCTGGTGACAATTAGCGTCAAGGCAATCCAGCTGGCGACAATCTTTTCCCCGAATAGCTCTATTCCGCTTGAGAAGCGAAc                                                                     >  1:359865/1‑235 (MQ=255)
                                                                                                                                                          cATAAGGAACGGTACACAGCCGGTCCCCAGACGGGTTGCAATATTGCCTACGATACCGATCGAGAAAGTGCGGGTTTTCAATAAATCTAATGAAATTAATGGGTTTGGCGTGCGTCGTGCATGGAGAATATAGAGAAGCAGTAACCCGATGCTGGTGACAATTACCGTCAAGGCAATCCAGCTGGCGACAATCTTTTCCCCGAATAGCTCTATTCCGCTTGAGAAGAGACCAAGGc                                                                >  1:227977/1‑236 (MQ=255)
                                                                                                                                                            tAAGGAACGGTACACAGCCGGTCCCCAGACGGGTTGCAATATTGCCTACGATACCGATCGAGAAAGTGCGGGTTTTCAATAAATCTAATGAAATTAATGGGTTTGGCGTGCGTCGTGCATGGAGAATATAGAGAAGCAGTAACCCGATGCTGGTGACAATTACCGTCAAGGCAATCCAGCTGGCGACAATCTTTTCCCCGAATAGCTCTATTCCGCTTGAGAAGAGAACAAGGCt                                                               >  1:406066/1‑235 (MQ=255)
                                                                                                                                                             aaGGAACGGTACACAGCCGGTCCCCAGACGGGTTGCAATATTGCCTACGATACCGATCGAGAAAGTGCGGGTTTTCAATAAATCTAATGAAATTAATGGGTTTGGCGTGCGTCGTGCATGGAGAATATAGAGAAGCAGTAACCCGATGCTGGTGACAATTACCGTCAAGGCAATCCAGCTGGCGACAATCTTTTCCCCGAATAGCTCTATTCCGCTTGAGAAGAGAACAAGGCTg                                                              <  1:647440/235‑1 (MQ=255)
                                                                                                                                                              aGGAACGGTACACAGCCGGTCCCCAGACGGGTTGCAATATTGCCTACGATACCGATCGAGAAAGTGCGGGTTTTCAATAAATCTAATGAAATTAATGGGTTTGGCGTGCGTCGTGCATGGAGAATATAGAGAAGCAGTAACCCGATGCTGGTGACAATTACCGTCAAGGCAATCCAGCTGGCGACAATCTTTTCCCCGAATAGCTCTATTCCGCTTGAGAAGAGAACAAGGCTGAg                                                            >  2:294385/1‑236 (MQ=255)
                                                                                                                                                                    ggTACACAGCCGGTCCCCAGACGGGTTGCAATATTGCCTACGATACCGATCGAGAAAGTGCGGGTTTTCAATAAATCTAATGAAATTAATGGGTTTGGCGTGCGTCGTGCATGGAGAATATAGAGAAGCAGTAACCCGATGCTGGTGACAATTACCGTCAAGGCAATCCAGCTGGCGACAATCTTTTCCCCGAATAGCTCTATTCCGCTTGAGAAGAGAACAAGGCTGAGGCCaa                                                       >  2:249700/1‑235 (MQ=255)
                                                                                                                                                                    ggTACACAGCCGGTCCCCAGACGGGTTGCAATATTGCCTACGATACCGATCGAGAAAGTGCGGGTTTTCAATAAATCTAATGAAATTAATGGGTTTGGCGTGCGTCGTGCATGGAGAATATAGAGAAGCAGTAACCCGATGCTGGTGACAATTACCGTCAAGGCAATCCAGCTGGCGACAATCTTTTCCCCGAATAGCTCTATTCCGCTTGAGAAGAGAACAAGGCTGAGGCCa                                                        <  2:406066/234‑1 (MQ=255)
                                                                                                                                                                     gTACACAGCCGGTCCACAGACGGGTTGCAATATTGCATACGATACCGATCGAGAAAGTGCGGGTTTTCAATAAATCTAATGAAATTAATGGGTTTGGCGTGCGTCGTGCATGGAGAATATAGAGAAGCAGTAACCCGATGCTGGTGACAATTACCGTCAAGGCAATCCAGCTGGCGACAATCTTTTCCCCGAATAGCTCTATTCCGCTTGAGAAGAGAACAAGGCTGAGGCCaa                                                       <  2:722524/234‑1 (MQ=255)
                                                                                                                                                                          caGCCGGTCCCCAGACGGGTTGCAATATTGCCTACGATACCGATCGAGAAAGTGCGGGTTTTCAATAAATCTAATGAAATTAATGGGTTTGGCGTGCGTCGTGCATGGAGAATATAGAGAAGCAGTAACCCGATGCTGGTGACAATTACCGTCAAGGCAATCCAGCTGGCGACAATCTTATCCCCGAATAGCTCTATTCCGCTTGAGAAGAGAACAAGGCTGAGGCCAAACAGCa                                                 <  2:709284/235‑1 (MQ=255)
                                                                                                                                                                                gTCCCCAGACGGGTTGCAATATTGCCTACGATACCGATCGAGAAAGTGCGGGTTTTCAATAAATCTAATGAAATTaa                                                                                                                                                                                                         <  2:201863/77‑1 (MQ=255)
                                                                                                                                                                                gTCCCCAGACGGGTTGCAATATTGCCTACGATACCGATCGAGAAAGTGCGGGTTTTCAATAAATCTAATGAAATTaa                                                                                                                                                                                                         >  1:201863/1‑77 (MQ=255)
                                                                                                                                                                                gTCCCCAGACGGGTTGCAATATTGCCTACGATACCGATCGAGAAAGTGCGGGTTTTCAATAAATCTAATGAAATTAATGGGTTTGGCGTGCGTCGTGCATGGAGAATATAGAGAAGCAGTAACCCGATGCTGGTGACAATTACCGTCAAGGCAATCCAGCTGGCGACAATCTTTTCCCCGAATAGCTCTATTCCGCTTgagaa                                                                           >  1:720522/1‑203 (MQ=255)
                                                                                                                                                                                gTCCCCAGACGGGTTGCAATATTGCCTACGATACCGATCGAGAAAGTGCGGGTTTTCAATAAATCTAATGAAATTAATGGGTTTGGCGTGCGTCGTGCATGGAGAATATAGAGAAGCAGTAACCCGATGCTGGTGACAATTACCGTCAAGGCAATCCAGCTGGCGACAATCTTTTCCCCGAATAGCTCTATTCCGCTTgagaa                                                                           <  2:720522/203‑1 (MQ=255)
                                                                                                                                                                                gTCCCCAGACGGGTTGCAATATTGCCTACGATACCGATCGAGAAAGTGCGGGTTTTCAATAAATCTAATGAAATTAATGGGTTTGGCGTGCGTCGTGCATGGAGAATATAGAGAAGCAGTAACCCGATGCTGGTGACAATTACCGTCAAGGCAATCCAGCTGGCGACAATCTTTTCCCCGAATAGCTCTATTCCGCTTGAGAAGAGAACAAGGCTGAGGCCAAACAGCAAAAAGcc                                          >  1:109758/1‑236 (MQ=255)
                                                                                                                                                                                    ccAGACGGGTTGCAATATTGCCTACGATACCGATCGAGAAAGTGCGGGTTTTCAATAAATCTAACGAAATTAATGGGTTTGGCGTGCGTCGTGCATGGAGAATATAGAGAAGCAGTAACCCGATGCTGGTGACAATTACCGTCAAGGCAATCCAGCTGGCGACAATCTTTTCCCCGAATAGCTCTATTCCGCTTGAGAAGAGAACAAGGCTGAGGCCAAACAGCAAAAAGCCAGt                                       <  2:75874/235‑1 (MQ=255)
                                                                                                                                                                                     cAGACGGGTTGCAATATTGCCTACGATACCGATCGAGAAAGTGCGGGTTTTCAATAAATCTAATGAAATTAATGGGTTTGGCGTGCGTCGTGCATGGAGAATATAGAGAAGCAGTAACCCGATGCTGGTGACAATTACCGTCAAGGCAATCCAGCTGGCGACAATCTTTTCCCCGAATAGCTCTATTCCGCTTGAGAAGAGAACAAGGCTGAGGCCAAACAGCAAAAAGCCAGTg                                      >  2:169579/1‑235 (MQ=255)
                                                                                                                                                                                           ggTTGCAATATTGCCTACGATACCGATCGAGAAAGTGCGGGTTTTCAATAAATCTAATGAAATTAATGGGTTTGGCGTGCGTCGTGCATGGAGAATATAGAGAAGCAGTAACCCGATGCTGGTGACAATTACCGTCAAGGCAATCCAGCTGGCGACAATCTTTTCCCCGAATAGCTCTATTCCGCTTGAGAAGAGAACAAGGCTGAGGCCAAACAGCAAAAAGCCAGTGATATc                                 >  1:49063/1‑234 (MQ=255)
                                                                                                                                                                                                      tGCCTACGATACCGATCGAGAAAGTGCGGGTTTTCAATAAATCTAATGAAATTAATGGGTTTGGCGTGCGTCGTGCATGGAGAATATAGAGAAGCAGTAACCCGATGCTGGTGACAATTACCGTCAAGGCAATCCAGCTGGCGACAATCTTTTCCCCGAATAGCTCTATTCCGCTTGAGAAGAGAACAAGGCTGAGGCCaaa                                                      >  2:44457/1‑202 (MQ=255)
                                                                                                                                                                                                      tGCCTACGATACCGATCGAGAAAGTGCGGGTTTTCAATAAATCTAATGAAATTAATGGGTTTGGCGTGCGTCGTGCATGGAGAATATAGAGAAGCAGTAACCCGATGCTGGTGACAATTACCGTCAAGGCAATCCAGCTGGCGACAATCTTTTCCCCGAATAGCTCTATTCCGCTTGAGAAGAGAACAAGGCTGAGGCCaaa                                                      <  1:44457/202‑1 (MQ=255)
                                                                                                                                                                                                      tGCCTACGATACCGATCGAGAAAGTGCGGGTTTTCAATAAATCTAATGAAATTAATGGGTTTGGCGTGCGTCGTGCATGGAGAATATAGAGAAGCAGTAACCCGATGCTGGTGACAATTACCGTCAAGGCAATCCAGCTGGCGACAATCTTTTCCCCGAATAGCTCTATTCCGCTTGAGAAGAGAACAAGGCTGAGGCCAAACAGCAAAAAGCCACTGATATCGAAACTGcg                        >  2:217143/1‑232 (MQ=255)
                                                                                                                                                                                                       gCCTACGATACCGATCGAGAAAGTGCGGGTTTTCAATAAATCTAATGAAATTAATGGGTTTGGCGTGCGTCGTGCATGGAGAATATAGAGAAGCAGTAACCCGATGCTGGTGACAATTACCGTCAAGGCAATCCAGCTGGCGACAATCTTTTCCCCGAATAGCTCTATTCCGCTTGAGAAGAGAACAAGGCTGAGGCCAAACAGCAAAAAGCCAGTGATATCGAATCTGCGTCGTg                   <  2:49063/236‑1 (MQ=255)
                                                                                                                                                                                                          tACGATACCGATCGAGAAAGTGCGGGTTTTCAATAAATCTAATGAAATTAATGGGTTTGGCGTGCGTCGTGCATGGAGAATATAGAGAAGCAGTAACCCGATGCTGGTGACAATTACCGTCAAGGCAATCCAGCTGGCGACAATCTTTTCCCCGAATAGCTCTATTCCGCTTGAGAAGAGAACAAGGCTGAGGCCAAACAGCAAAAAGCCAGTGATATCGAATCTGCGTCGTGCgg                <  1:522296/236‑1 (MQ=255)
                                                                                                                                                                                                          tACGATACCGATCGAGAAAGTGCGGGTTTTCAATAAATCTAATGAAATTAATGGGTTTGGCGTGCGTCGTGCATGGAGAATATAGAGAAGCAGTAACCCGATGCTGGTGACAATTACCGTCAAGGCAATCCAGCTGGCGACAATCTTTTCCCCGAATAGCTCTATTCCGCTTGAGAAGAGAACAAGGCTGAGGCCAAACAGCAAAAAGCCAGTGATATCGAATCTGCGTCGTGCgg                <  1:240077/236‑1 (MQ=255)
                                                                                                                                                                                                                    aTCGAGAAAGTGCGGGTTTTCAATAAATCTAATGAAATTAATGGGTTTGGCGTGCGTCGTGCATGGAGAATATAGAGAAGCAGTAACCCGATGCTGGTGACAATTACCGTCAAGGCAATCCAGCTGGCGACAATCTTTTCCCCGAATAGCTCTATTCCGCTTGTGAAGAGAACAAGGCTGAGGCCAAACAGCCAAAAGCCAGTGATATCGAATCTGCGTCGCGCGGTGGCGAAAtt      >  1:746465/1‑236 (MQ=255)
                                                                                                                                                                                                                    aTCGAGAAAGTGCGGGTTTTCAATAAATCTAATGAAATTAATGGGTTTGGCGTGCGTCGTGCATGGAGAATATAGAGAAGCAGTAACCCGATGCTGGTGACAATTACCGTCAAGGCAATCCAGCTGGCGACAATCTTTTCCCCGAATAGCTCTATTCCGCTTGAGAAGAGAACAAGGCTGAGGCCAAACAGCAAAAAGCCAGTGATATCGAATCTGCGTCGTGCGGTGGTGAAAt       >  1:425448/1‑235 (MQ=255)
                                                                                                                                                                                                                        agaAAGTGCGGGTTTTCAATAAATCTAATGAAATTAATGGGTTTGGCGTGCGTCGTGCATGGAGAATATAGAGAAGCAGTAACCCGATGCTGGTGACAATTACCGTCAAGGCAATCCAGCTGGCGACAATCTTTTCCCCGAATAGCTCTATTCCGCTTGAGAAGAGAACAAGGCTGAGGCCAAACAGCAAAAAGCCAGTGATATCGAATCTGCGTCGTGCGGTGGTGAAATTGGGc  <  1:716311/236‑1 (MQ=255)
                                                                                                                                                                                                                                        |                                                                                                                                                                                                                           
TCACCGTGATCCCCACTAACGTATGGCGATAGCCCAGACGACGTAAGACTTGGGTAACCATCGATTTTGCAATAATGGAACCTAACGCTGTCGGTGCCATCATACAGCCAGCAATAAACGCCTGATAACCAAATCCTACCTGTAACATCAATGGCATAAGGAACGGTACACAGCCGGTCCCCAGACGGGTTGCAATATTGCCTACGATACCGATCGAGAAAGTGCGGGTTTTAAATAAATCTAATGAAATTAATGGGTTTGGCGTGCGTCGTGCATGGAGAATATAGAGAAGCAGTAACCCGATGCTGGTGACAATTACCGTCAAGGCAATCCAGCTGGCGACAATCTTTTCCCCGAATAGCTCTATTCCGCTTGAGAAGAGAACAAGGCTGAGGCCAAACAGCAAAAAGCCAGTGATATCGAATCTGCGTCGTGCGGTGGTGAAATTGGGC  >  NC_000913/3939604‑3940055

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