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

ACAAACCGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCGGCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTCCGGCGCTGCAATTGGTCGCGCAAATTCGGCGGCGTGCCTTTAATGGTCAAGGTATCGGTTGCTGGATCCCAGAAGATGCGTTCACCTAACAGCA  >  NC_000913/2282761‑2283052
                                                                                                                                                  |                                                                                                                                                  
aCAAACCGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCGGCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCAcg                                                                                                                                                 >  4:886753/1‑149 (MQ=255)
  aaaCCGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCGGCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCAcgcc                                                                                                                                               >  3:722727/1‑149 (MQ=255)
          gCTGCAGCACATCGACTTCGTTCGCGCTTAATAGCTGAATAAGTTGTAAAAAACCCCGCTCCGGCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTgg                                                                                                                                       <  4:323955/149‑1 (MQ=255)
            tGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAGAACCCCGCTCCGTCGGGGTTTTTTGTATCTGCAGATTATGCCTGACTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCa                                                                                                                                     >  4:908569/1‑149 (MQ=255)
            tGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCGGCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATGCGGCGTGAACGCCTTATCTGGc                                                                                                                                      <  1:759738/148‑1 (MQ=255)
                gcaCATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAATCCCCGCTCCGGCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTc                                                                                                                                  <  2:161251/148‑1 (MQ=255)
                gcaCATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCGGCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTc                                                                                                                                  >  1:497318/1‑148 (MQ=255)
                gcaCATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCGGCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCa                                                                                                                                 >  4:459203/1‑149 (MQ=255)
                gcaCATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCGGCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCa                                                                                                                                 <  3:372801/149‑1 (MQ=255)
                 cacaGCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCGGCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAg                                                                                                                                >  4:850906/1‑149 (MQ=255)
                   caTCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCGGCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGaa                                                                                                                              <  4:486734/149‑1 (MQ=255)
                   caTCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCGGCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGaa                                                                                                                              <  3:846358/149‑1 (MQ=255)
                   caTCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCGGCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGa                                                                                                                               <  2:379857/148‑1 (MQ=255)
                   caTCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCGGCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACACCGCATCTGGCATTCAGaa                                                                                                                              <  3:789505/149‑1 (MQ=255)
                    aTCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCGGCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGaag                                                                                                                             >  3:85748/1‑148 (MQ=255)
                    aTCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCGGCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTAACGCCGCATCTGGCATTCAGaa                                                                                                                              >  1:247758/1‑148 (MQ=255)
                     tCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCGGCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACCCCCCATCTGGCATTCAGAac                                                                                                                             >  2:876324/1‑148 (MQ=255)
                       gACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCGGCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAacac                                                                                                                           >  1:688152/1‑148 (MQ=255)
                               tcgCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCGGCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGa                                                                                                                   <  1:783946/148‑1 (MQ=255)
                               tcgCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCGGCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGAt                                                                                                                  >  3:495258/1‑149 (MQ=255)
                                   gcTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCGGCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTaa                                                                                                               <  1:938731/148‑1 (MQ=255)
                                      tAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCGGCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTgc                                                                                                           <  4:948533/149‑1 (MQ=255)
                                      tAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCGGCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTg                                                                                                            >  2:921664/1‑148 (MQ=255)
                                      tAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCGGCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTg                                                                                                            >  2:591828/1‑148 (MQ=255)
                                       aaTTGCTGAATAAGTTGTAAAAAACCCCGCTCCGGCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTgcc                                                                                                          <  3:403931/149‑1 (MQ=255)
                                       aaTTGCTGAATAAGTTGTAAAAAACCCCGCTCCGGCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTgc                                                                                                           <  1:310494/148‑1 (MQ=255)
                                          tGCTGAATAAGTTGTAAAAAACCCCGCTCCGGCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTgccgc                                                                                                        >  1:434124/1‑148 (MQ=255)
                                           gCTGAATAAGTTGTAAAAAACCCCGCTCCGGCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCa                                                                                                      >  3:754982/1‑149 (MQ=255)
                                             tGAATAAGTTGTAAAAAACCCCGCTCCGGCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGa                                                                                                    <  3:831049/149‑1 (MQ=255)
                                              gAATAAGTTGTAAAAAACCCCGCTCCGGCTGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGATCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCTGCCAGa                                                                                                    <  2:50551/148‑1 (MQ=255)
                                              gAATAAGTTGTAAAAAACCCCGCTCCGGCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGa                                                                                                    <  2:594443/148‑1 (MQ=255)
                                               aaTAAGTTGTAAAAAACCCCGCTCCGGCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGAc                                                                                                   >  2:819976/1‑148 (MQ=255)
                                               aaTAAGTTGTAAAAAACCCCGCTCCGGCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACg                                                                                                  >  4:151289/1‑149 (MQ=255)
                                                aTAAGTTGTAAAAAACCCCGCTCCGGCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACCCAAAACCGATTAATTGCCGCCAGACg                                                                                                  <  2:818825/148‑1 (MQ=255)
                                                aTAAGTTGTAAAAAACCCCGCTCCGGCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACg                                                                                                  <  2:141435/148‑1 (MQ=255)
                                                aTAAGTTGTAAAAAACCCCGCTCCGGCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACg                                                                                                  >  1:170446/1‑148 (MQ=255)
                                                aTAAGTTGTAAAAAACCCCGCTCCGGCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGt                                                                                                 >  4:79431/1‑149 (MQ=255)
                                                 tAAGTTGTAAAAAACCCCGCTCCGGCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGt                                                                                                 <  1:287378/148‑1 (MQ=255)
                                                 tAAGTTGTAAAAAACCCCGCTCCGGCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTc                                                                                                >  4:634311/1‑149 (MQ=255)
                                                 tAAGTTGTAAAAAACCCCGCTCCGGCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTc                                                                                                >  3:844528/1‑149 (MQ=255)
                                                     ttGTAAAAAACCCCGCTCCGGCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATGCGGCGTGAACGCCTTATCTGGcc                                                                                                                                     <  4:163327/108‑2 (MQ=38)
                                                     ttGTAAAAAACCCCGCTCCGGCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATGCGGCGTGAACGCCTTATCTGGcc                                                                                                                                     >  3:163327/1‑107 (MQ=38)
                                                      tGTAAAAAACCCCGCTCCGGCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATGCGGCGTGAACGCCTTATCTGGc                                                                                                                                      >  1:926866/1‑106 (MQ=38)
                                                      tGTAAAAAACCCCGCTCCGGCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTCCGgc                                                                                            >  2:663899/1‑148 (MQ=255)
                                                        tAAAAAACCCCGCTCCGGCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTCCGGCGCt                                                                                         >  4:917664/1‑149 (MQ=255)
                                                         aaaaaaCCCCGCTCCGGCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTCCGGCGCTg                                                                                        >  4:409824/1‑149 (MQ=255)
                                                         aaaaaaCCCCGCCGGAGCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTCCGGCGCTg                                                                                        >  3:357310/1‑149 (MQ=39)
                                                          aaaaaTCCAGCTCCGGCGGGGTTTTTTAGATGTGCAGATTATGCCTGTTTAGGGTATTCCTATTTTTTACAGGCCAGATAAGGTGTTCACGCCGCATATGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTCCGGCGCTGc                                                                                       >  4:237438/1‑149 (MQ=255)
                                                            aaaCCCCGCTCCGGCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTCCGGCGCTGCaa                                                                                     >  4:763246/1‑149 (MQ=255)
                                                                 ccGCTCCGGCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTCCGGCGCTGCAATTgg                                                                                 >  1:671594/1‑148 (MQ=255)
                                                                  cGCTCCGGCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTCCGGCGCTGCAATTGGt                                                                                >  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)
                                                                                                                                                  |                                                                                                                                                  
ACAAACCGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCGGCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTCCGGCGCTGCAATTGGTCGCGCAAATTCGGCGGCGTGCCTTTAATGGTCAAGGTATCGGTTGCTGGATCCCAGAAGATGCGTTCACCTAACAGCA  >  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: