Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,282,906 C→A 5.0% intergenic (+105/+34) rplY → / ← yejK 50S ribosomal subunit protein L25/spermidine nucleoid‑associated protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,282,9060CA 5.0% 332.7 / 3.1 96intergenic (+105/+34)rplY/yejK50S ribosomal subunit protein L25/spermidine nucleoid‑associated protein
Reads supporting (aligned to +/- strand):  ref base C (42/49);  new base A (3/2);  total (45/51)
Fisher's exact test for biased strand distribution p-value = 6.63e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.90e-01

ACAAACCGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTCCGGCGCTGCAATTGGTCGCGCAAATTCGGCGGCGTGCCTTTAATGGTCAAGGTATCGGTTGCTGGATCCCAGAAGATGCGTTCACCTAACAGCA  >  NC_000913/2282761‑2283052
                                                                                                                                                  |                                                                                                                                                  
aCAAACCGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCAcg                                                                                                                                                 >  4:886753/1‑149 (MQ=255)
  aaaCCGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCAcgcc                                                                                                                                               >  3:722727/1‑149 (MQ=255)
          gCTGCAGCACATCGACTTCGTTCGCGCTTAATAGCTGAATAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTgg                                                                                                                                       <  4:323955/149‑1 (MQ=255)
            tGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAGAACCCCGCTCCG‑TCGGGGTTTTTTGTATCTGCAGATTATGCCTGACTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCa                                                                                                                                     >  4:908569/1‑149 (MQ=255)
            tGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATGCGGCGTGAACGCCTTATCTGGc                                                                                                                                      <  1:759738/148‑1 (MQ=255)
                gcaCATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAATCCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTc                                                                                                                                  <  2:161251/148‑1 (MQ=255)
                gcaCATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTc                                                                                                                                  >  1:497318/1‑148 (MQ=255)
                gcaCATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCa                                                                                                                                 >  4:459203/1‑149 (MQ=255)
                gcaCATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCa                                                                                                                                 <  3:372801/149‑1 (MQ=255)
                 cacaGCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAg                                                                                                                                >  4:850906/1‑149 (MQ=255)
                   caTCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGaa                                                                                                                              <  4:486734/149‑1 (MQ=255)
                   caTCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGaa                                                                                                                              <  3:846358/149‑1 (MQ=255)
                   caTCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGa                                                                                                                               <  2:379857/148‑1 (MQ=255)
                   caTCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACACCGCATCTGGCATTCAGaa                                                                                                                              <  3:789505/149‑1 (MQ=255)
                    aTCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGaag                                                                                                                             >  3:85748/1‑148 (MQ=255)
                    aTCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTAACGCCGCATCTGGCATTCAGaa                                                                                                                              >  1:247758/1‑148 (MQ=255)
                     tCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACCCCCCATCTGGCATTCAGAac                                                                                                                             >  2:876324/1‑148 (MQ=255)
                       gACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAacac                                                                                                                           >  1:688152/1‑148 (MQ=255)
                               tcgCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGa                                                                                                                   <  1:783946/148‑1 (MQ=255)
                               tcgCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGAt                                                                                                                  >  3:495258/1‑149 (MQ=255)
                                   gcTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTaa                                                                                                               <  1:938731/148‑1 (MQ=255)
                                      tAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTgc                                                                                                           <  4:948533/149‑1 (MQ=255)
                                      tAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTg                                                                                                            >  2:921664/1‑148 (MQ=255)
                                      tAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTg                                                                                                            >  2:591828/1‑148 (MQ=255)
                                       aaTTGCTGAATAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTgcc                                                                                                          <  3:403931/149‑1 (MQ=255)
                                       aaTTGCTGAATAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTgc                                                                                                           <  1:310494/148‑1 (MQ=255)
                                          tGCTGAATAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTgccgc                                                                                                        >  1:434124/1‑148 (MQ=255)
                                           gCTGAATAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCa                                                                                                      >  3:754982/1‑149 (MQ=255)
                                             tGAATAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGa                                                                                                    <  3:831049/149‑1 (MQ=255)
                                              gAATAAGTTGTAAAAAACCCCGCTCCG‑GCTGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGATCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCTGCCAGa                                                                                                    <  2:50551/148‑1 (MQ=255)
                                              gAATAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGa                                                                                                    <  2:594443/148‑1 (MQ=255)
                                               aaTAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGAc                                                                                                   >  2:819976/1‑148 (MQ=255)
                                               aaTAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACg                                                                                                  >  4:151289/1‑149 (MQ=255)
                                                aTAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACCCAAAACCGATTAATTGCCGCCAGACg                                                                                                  <  2:818825/148‑1 (MQ=255)
                                                aTAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACg                                                                                                  <  2:141435/148‑1 (MQ=255)
                                                aTAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACg                                                                                                  >  1:170446/1‑148 (MQ=255)
                                                aTAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGt                                                                                                 >  4:79431/1‑149 (MQ=255)
                                                 tAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGt                                                                                                 <  1:287378/148‑1 (MQ=255)
                                                 tAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTc                                                                                                >  4:634311/1‑149 (MQ=255)
                                                 tAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTc                                                                                                >  3:844528/1‑149 (MQ=255)
                                                     ttGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATGCGGCGTGAACGCCTTATCTGGcc                                                                                                                                     <  4:163327/108‑2 (MQ=38)
                                                     ttGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATGCGGCGTGAACGCCTTATCTGGcc                                                                                                                                     >  3:163327/1‑107 (MQ=38)
                                                      tGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATGCGGCGTGAACGCCTTATCTGGc                                                                                                                                      >  1:926866/1‑106 (MQ=38)
                                                      tGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTCCGgc                                                                                            >  2:663899/1‑148 (MQ=255)
                                                        tAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTCCGGCGCt                                                                                         >  4:917664/1‑149 (MQ=255)
                                                         aaaaaaCCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTCCGGCGCTg                                                                                        >  4:409824/1‑149 (MQ=255)
                                                         aaaaaaCCCCGC‑CGGAGCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTCCGGCGCTg                                                                                        >  3:357310/1‑149 (MQ=39)
                                                          aaaaaTCCAGCTCCG‑GCGGGGTTTTTTAGATGTGCAGATTATGCCTGTTTAGGGTATTCCTATTTTTTACAGGCCAGATAAGGTGTTCACGCCGCATATGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTCCGGCGCTGc                                                                                       >  4:237438/1‑149 (MQ=255)
                                                            aaaCCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTCCGGCGCTGCaa                                                                                     >  4:763246/1‑149 (MQ=255)
                                                                 ccGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTCCGGCGCTGCAATTgg                                                                                 >  1:671594/1‑148 (MQ=255)
                                                                  cGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTCCGGCGCTGCAATTGGt                                                                                >  1:694450/1‑148 (MQ=255)
                                                                          gcggGTTTTTTTGTGTCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCCCGCCGCACCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCATACGCCCGGCGCTCCCATTGCTCTTGCaa                                                                         <  2:781766/148‑1 (MQ=255)
                                                                           cggGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTCCGGCGCTGCAATTGGTCGCGCAAAt                                                                       >  3:475044/1‑149 (MQ=255)
                                                                               gTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTCCGGCGCTGCAATTGGTCGCGCAAATTcg                                                                    <  2:33549/148‑1 (MQ=255)
                                                                                ttttttGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTCCGGCGCTGCAATTGGTCGCGCAAATTcggc                                                                  <  3:371091/149‑1 (MQ=255)
                                                                                ttttttGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTCCGGCGCTGCAATTGGTCGCGCAAATTcgg                                                                   <  1:265815/148‑1 (MQ=255)
                                                                                  ttttGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTCCGGCGCTGCAATTGGTCGCGCAAATTcggcgg                                                                >  4:617870/1‑149 (MQ=255)
                                                                                   tttGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTCCGGCGCTGCAATTGGTCGCGCAAATTCggcggc                                                               <  4:517026/149‑1 (MQ=255)
                                                                                       tATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTCCGGCGCTGCAATTGGTCGCGCAAATTCGGCGGCGTg                                                            <  1:314494/148‑1 (MQ=255)
                                                                                         tCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTaa                                                                                                               <  1:209281/95‑1 (MQ=255)
                                                                                         tCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTaa                                                                                                               >  2:209281/1‑95 (MQ=255)
                                                                                         tCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTCCGGCGCTGCAATTGGTCGCGCAAATTCGGCGGCGTGCCt                                                         <  3:409824/149‑1 (MQ=255)
                                                                                            gCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTCCGGCGCTGCAATTGGTCGCGCAAATTCGGCGGCGTGCCttt                                                       >  2:404433/1‑148 (MQ=255)
                                                                                             cAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTCCGGCGCTGCAATTGGTCGCGCAAATTCGGCGGCGTGCCTTTaa                                                     <  4:475044/149‑1 (MQ=255)
                                                                                             cAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTCCGGCGCTGCAATTGGTCGCGCAAATTCGGCGGCGTGCCTTTa                                                      <  1:921664/148‑1 (MQ=255)
                                                                                              aGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTCCGGCGCTGCAATTGGTCGCGCAAATTCGGCGGCGTGCCTTTAAt                                                    <  3:766131/149‑1 (MQ=255)
                                                                                              aGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTCCGGCGCTGCAATTGGTCGCGCAAATTCGGCGGCGTGCCTTTAAt                                                    <  3:908569/149‑1 (MQ=255)
                                                                                                 ttATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTCCGGCGCTGCAATTGGTCGCGCAAATTCGGCGGCGTGCCTTTAATGGt                                                 <  3:886753/149‑1 (MQ=255)
                                                                                                    tGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTCCGGCGCTGCAATTGGTCGCGCAAATTCGGCGGCGTGCCTTTAATGGTCaa                                              <  3:151289/149‑1 (MQ=255)
                                                                                                    tGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTCCGGCGCTGCAATTGGTCGCGCAAATTCGGCGGCGTGCCTTTAATGGTCa                                               <  1:630072/148‑1 (MQ=255)
                                                                                                     gCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTCCGGCGCTGCAATTGGTCGCGCAAATTCGGCGGCGTGCCTTTAATGGTCaa                                              <  1:644289/148‑1 (MQ=255)
                                                                                                     gCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCCATTAATTGCCGCCAGACGTCCGGCGCTGCAATTGGTCGCGCAAATTCGGCGGCGTGCCTTTAACGGTCAAg                                             <  4:236296/149‑1 (MQ=255)
                                                                                                          aTTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTCCGGCGCTGCAATTGGTCGCGCAAATTCGGCGGCGTGCCTTTAATGGTCAAGGTAt                                         <  1:819976/148‑1 (MQ=255)
                                                                                                           ttACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTCCGGCGCTGCAATTGGTCGCGCAAATTCGGCGGCGTGCCTTTAATGGTCAAGGTATCg                                       <  4:357310/149‑1 (MQ=255)
                                                                                                            tACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTCCGGCGCTGCAATTGGTCGCGCAAATTCGGCGGCGTGCCTTTAATGGTCAAGGTATCg                                       >  2:484845/1‑148 (MQ=255)
                                                                                                                 tATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTCCGGCGCTGCAATTGGTCGCGCAAATTCGGCGGCGTGCCTTTAATGGTCAAGGTATCGGTTGCt                                 <  3:617870/149‑1 (MQ=255)
                                                                                                                   ttGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTCCGGCGCTGCAATTGGTCGCGCAAATTCGGCGGCGTGCCTTTAATGGTCAAGGTATCGGTTGCTg                                >  1:217369/1‑148 (MQ=255)
                                                                                                                    tGCTATTTTTTGCATGCCAAAGAAGGCGTTCACGCGGCATCTGTCATTCATTACACAAAACCGATTAATTGCCGCCAGACGACCGGCGCTGCAATTAATCGCGCAAATTCGGCGGCGTGCCTTTAATGGTCAAGGTATCGGTTGCTgg                               >  2:224051/1‑148 (MQ=255)
                                                                                                                     gCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTCCGGCGCTGCAATTGGTCGCGCAAATTCGGCGGCGTGCCTTTAATGGTCAAGGTATCGGTTGCTGGa                              <  2:434124/148‑1 (MQ=255)
                                                                                                                       tATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTCCGGCGCTGCAATTGGTCGCGCAAATTCGGCGGCGTGCCTTTAATGGTCAAGGTATCGGTTGCTGGATcc                           >  3:754197/1‑149 (MQ=255)
                                                                                                                       tATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTCCGGCGCTGCAATTGGTCGCGCAAATTCGGCGGCGTGCCTTTAATGGTCAAGGTATCGGTTGCTGGATc                            >  2:368456/1‑148 (MQ=255)
                                                                                                                              tGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTCCGGCGCTGCAATTGGTCGCGCAAATTCGGCGGCGTGCCTTTAATGGTCAAGGTATCGGTTGCTGGATCCCagaag                     <  1:368456/148‑1 (MQ=255)
                                                                                                                              tGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTCCGGCGCTGCAATTGGTCGCGCAAATTCGGCGGCGTGCCTTTAATGGTCAAGGTATCGGTTGCTGGATCCCagaag                     <  2:763564/148‑1 (MQ=255)
                                                                                                                               gCAGGCCAGATAAGGCGTTCACGACGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTCCGGCGCTGCAATTGGTCGCGCAAATTCGGCGGCGTGCCTTTAATGGTCAAGGTATCGGTTGCTGGATCCCagaaga                    <  2:130169/148‑1 (MQ=255)
                                                                                                                                 aGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTCCGGCGCTGCAATTGGTCGCGCAAATTCGGCGGCGTGCCTTTAATGGTCAAGGTATCGGTTGCTGGATCCCAGAAGATGc                 <  4:779968/149‑1 (MQ=255)
                                                                                                                                 aGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTCCGGCGCTGCAATTGGTCGCGCAAATTCGGCGGCGTGCCTTTAATGGTCAAGGTATCGGTTGCTGGATCCCAGAAGATGc                 >  4:118580/1‑149 (MQ=255)
                                                                                                                                  ggCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTCCGGCGCTGCAATTGGTCGCGCAAATTCGGCGGCGTGCCTTTAATGGTCAAGGTATCGGTTGCTGGATCCCAGAAGATGc                 >  2:262605/1‑148 (MQ=255)
                                                                                                                                  ggCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTCCGGCGCTGCAATTGGTCGCGCAAATTCGGCGGCGTGCCTTTAATGGTCAAGGTATCGGTTGCTGGATCCCAGAAGATGCg                >  4:696339/1‑149 (MQ=255)
                                                                                                                                   gCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTCCGGCGCTGCAATTGGTCGCGCAAATTCGGCGGCGTGCCTTTAATGGTCAAGGTATCGGTTGCTGGATCCCAGAAGATGCg                <  1:663899/148‑1 (MQ=255)
                                                                                                                                      aGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTCCGGCGCTGCAATTGGTCGCGCAAATTCGGCGGCGTGCCTTTAATGGTCAAGGTATCGGTTGCTGGATCCCAGAAGATGCGTTc             <  2:492624/148‑1 (MQ=255)
                                                                                                                                       gATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTCCGGCGCTGCAATTGGTCGCGCAAATTCGGCGGCGTGCCTTTAATGGTCAAGGTATCGGTTGCTGGATCCCAGAAGATGCGTTCa            <  2:688152/148‑1 (MQ=255)
                                                                                                                                           aGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTCCGGCGCTGCAATTGGTCGCGCAAATTCGGCGGCGTGCCTTTAATGGTCAAGGTATCGGTTGCTGGATCCCAGAAGATGCGTTCACCTa        <  2:694450/148‑1 (MQ=255)
                                                                                                                                            ggCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTCCGGCGCTGCAATTGGTCGCGCAAATTCGGCGGCGTGCCTTTAATGGTCAAGGTATCGGTTGCTGGATCCCAGAAGATGCGTTCACCTaa       <  1:606638/148‑1 (MQ=255)
                                                                                                                                             gCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTCCGGCGCTGCAATTGGTCGCGCAAATTCGGCGGCGTGCCTTTAATGGTCAAGGTATCGGTTGCTGGATCCCAGAAGATGCGTTCACCTAAc      <  1:404433/148‑1 (MQ=255)
                                                                                                                                               gTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTCCGGCGCTGCAATTGGTCGCGCAAATTCGGCGGCGTGCCTTTAATGGTCAAGGTATCGGTTGCTGGATCCCAGAAGATGCGTTCACCTAACAgc   >  4:282977/1‑149 (MQ=255)
                                                                                                                                                 tCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTCCGGCGCTGCAATTGGTCGCGCAAATTCGGCGGCGTGCCTTTAATGGTCAAGGTATCGGTTGCAGGATCCCAGAAGATGCGTTCACGTAACAgca  >  2:611528/1‑148 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
ACAAACCGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTCCGGCGCTGCAATTGGTCGCGCAAATTCGGCGGCGTGCCTTTAATGGTCAAGGTATCGGTTGCTGGATCCCAGAAGATGCGTTCACCTAACAGCA  >  NC_000913/2282761‑2283052

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑