Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,469,199 T→A 5.0% intergenic (‑50/‑42) insC1 ← / → ydbA IS2 repressor TnpA/pseudogene, autotransporter homolog; interrupted by IS2 and IS30

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,469,1990TA 5.0% 352.0 / 6.9 99intergenic (‑50/‑42)insC1/ydbAIS2 repressor TnpA/pseudogene, autotransporter homolog; interrupted by IS2 and IS30
Reads supporting (aligned to +/- strand):  ref base T (45/49);  new base A (3/2);  total (48/51)
Fisher's exact test for biased strand distribution p-value = 6.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.78e-01

Only 100 of 101 total aligned reads displayed.

TTTTCCTGTGTGGTACGCCGTCTGCGTTTCTCCGGCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAAACCCCAGAACTTACTTATGCTGATTCCGGTGCGAAAATTGTTAATAAAGGTACTTGTAGATTCAATTGGTCAACGCAACAGTTATGTGAAAACATGGGGTTG  >  NC_000913/1469097‑1469343
                                                                                                      |                                                                                                                                                
ttttCCTGTGTGGTACGCCGTCTGCGTTTCTCCGGCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAAAccc                                                                                                    >  1:691122/1‑149 (MQ=31)
ttttCCTGTGTGGTACGCCGTCTGCGTTTCTCCGGCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAAAccc                                                                                                    >  3:172238/1‑149 (MQ=31)
ttttCCTGTGTGGTACGCCGTCTGCGTTTCTCCGGCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAAAccc                                                                                                    >  1:441311/1‑149 (MQ=31)
  ttCCTGTGTGGTACGCCGTCTGCGTTTCTCCGGCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAAACCCCa                                                                                                  <  1:154414/149‑1 (MQ=32)
  ttCCTGTGTGGTACGCCGTCTGCGTTTCTCCGGCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAAACCCCa                                                                                                  <  1:62003/149‑1 (MQ=32)
  ttCCTGTGTGGTACGCCGTCTGCGTTTCTCCGGCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAAACCCCa                                                                                                  <  4:647917/149‑1 (MQ=32)
     cTGTGTGGTACGCCGTCTGCGTTTCTCCGGCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAAACCCCAGaa                                                                                               >  2:654397/1‑149 (MQ=35)
     cTGTGTGGTACGCCGTCTGCGTTTCTCCGGCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAAACCCCAGaa                                                                                               <  1:374257/149‑1 (MQ=35)
     cTGTGTGGTACGCCGTCTGCGTTTCTCCGGCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAAACCCCAGaa                                                                                               <  2:368290/149‑1 (MQ=35)
      tgtgtgGTACGCCGTCTGCGTTTCTCCGGCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAAACCCCAGAac                                                                                              >  2:292391/1‑149 (MQ=35)
      tgtgtgGTACGCCGTCTGCGTTTCTCCGGCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAAACCCCAGAac                                                                                              <  1:921054/149‑1 (MQ=35)
      tgtgtgGTACGCCGTCTGCGTTTCTCCGGCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAAACCCCAGAac                                                                                              <  3:564346/149‑1 (MQ=35)
      tgtgtgGTACGCCGTCTGCGTTTCTCCGGCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAAACCCCAGAac                                                                                              >  4:767950/1‑149 (MQ=35)
         gtgGTACGCCGTCTGCGTTTCTCCGGCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAAACCCCAGAactta                                                                                           <  1:574696/149‑1 (MQ=38)
          tgGTACGCCGTCTGCGTTTCTCCGGCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAAACCCCAGAacttac                                                                                          <  1:249110/149‑1 (MQ=39)
              aCGCCGTCTGCGTTTCTCCGGCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAAACCCCAGAACTTACTTAt                                                                                      <  3:925061/149‑1 (MQ=255)
                 ccGTCTGCGTTTCTCCGGCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAAACCCCAGAACTTACTTATGCt                                                                                   <  4:213550/149‑1 (MQ=255)
                 ccGTCTGCGTTTCTCCGGCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAAACCCCAGAACTTACTTATGCt                                                                                   >  1:109892/1‑149 (MQ=255)
                   gTCTGCGTTTCTCCGGCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAACCCCCAGAACTTACTTATGCTGa                                                                                 <  1:568940/149‑1 (MQ=255)
                   gTCTGCGTTTCTCCGGCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAAACCCCAGAACTTACTTATGCTGa                                                                                 <  4:544234/149‑1 (MQ=255)
                    tCTGCGTTTCTCCGGCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAAACCCCAGAACTTACTTATGCTGAt                                                                                <  2:100904/149‑1 (MQ=255)
                      tGCGTTTCTCCGGCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAACACTATCACTTATTTAAGTCTTATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAAACCCCAGAACTTACTTATGCTGATTc                                                                              <  4:788484/149‑1 (MQ=255)
                         gTTTCTCCGGCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCTAGTCTAAACCCCAGAACTTACTTATGCTGATTCCgg                                                                           >  3:883079/1‑149 (MQ=255)
                          tttCTCCGGCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAAACCCCAGAACTTACTTATGCTGATTCCGGt                                                                          <  4:617576/149‑1 (MQ=255)
                                  gCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGTTATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAAACCCCAGAACTTACTTATGCTGATTCCGGTGCGACAAt                                                                  >  2:185684/1‑149 (MQ=255)
                                    ccTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAAACCCCa                                                                                                  >  1:646711/1‑115 (MQ=33)
                                    ccTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAAACCCCa                                                                                                  <  2:646711/115‑1 (MQ=33)
                                    ccTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAAACCCCAGAACTTACTTATGCTGATTCCGGTGCGAAAATTg                                                                >  3:86645/1‑149 (MQ=255)
                                    ccTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAAACCCCAGAACTTACTTATGCTGATTCCGGTGCGAAAATTg                                                                >  1:343397/1‑149 (MQ=255)
                                     cTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAAACCCCAGAACTTACTTATGCTGATTCCGGTGCGAAAATTGt                                                               <  2:607870/149‑1 (MQ=255)
                                     cTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAAACCCCAGAACTTACTTATGCTGATTCCGGTGCGAAAATTGt                                                               >  1:893705/1‑149 (MQ=255)
                                     cTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAAACCCCAGAACTTACTTATGCTGATTCCGGTGCGAAAATTGt                                                               <  3:362596/149‑1 (MQ=255)
                                      tAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAAACCCCAGAACTTACTTATGCTGATTCCGGTGCGAAAATTGtt                                                              <  1:651061/149‑1 (MQ=255)
                                      tAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAAACCCCAGAACTTACTTATGCTGATTCCGGTGCGAAAATTGtt                                                              >  2:934837/1‑149 (MQ=255)
                                      tAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAAACCCCAGAACTTACTTATGCTGATTCCGGTGCGAAAATTGtt                                                              <  2:181662/149‑1 (MQ=255)
                                      tAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAAACCCCAGAACTTACTTATGCTGATTCCGGTGCGAAAATTGtt                                                              >  1:502634/1‑149 (MQ=255)
                                      tAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAAACCCCAGAACTTACTTATGCTGATTCCGGTGCGAAAATTGt                                                               >  2:62231/1‑148 (MQ=255)
                                      tAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAAACCCCAGAACTTACTTATGCTGATTCCGGTGCGAAAATTGt                                                               <  1:62231/148‑1 (MQ=255)
                                        aGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAAACCCCAGAACTTACTTATGCTGATTCCGGTGCGAAAATTGTtaa                                                            >  1:321537/1‑149 (MQ=255)
                                         gACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAAACCCCAGAACTTACTTATGCTGATTCCGGTGCGAAAATTGTtaat                                                           <  4:90412/149‑1 (MQ=255)
                                            atcaatcaTCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAAACCCCAGAACTTACTTATGCTGATTCCGGTGCGAAAATTGTTAATaaa                                                        >  3:94755/1‑149 (MQ=255)
                                            atcaatcaTCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAAACCCCAGAACTTACTTATGCTGATTCCGGTGCGAAAATTGTTAATaaa                                                        <  1:540116/149‑1 (MQ=255)
                                            atcaatcaTCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAAACCCCAGAACTTACTTATGCTGATTCCGGTGCGAAAATTGTTAATaaa                                                        <  4:29752/149‑1 (MQ=255)
                                            atcaagcaTCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAAACCCCAGAACTTACTTATGCTGATTCCGGTGCGAAAATTGTTAATaaa                                                        <  4:192563/149‑1 (MQ=255)
                                             tcaatcaTCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAAACCCCAGAACTTACTTATGCTGATTCCGGTGCGAAAATTGTTAATAAAg                                                       <  2:314580/149‑1 (MQ=255)
                                             tcaatcaTCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAAACCCCAGAACTTACTTATGCTGATTCCGGTGCGAAAATTGTTAATAAAg                                                       <  4:264760/149‑1 (MQ=255)
                                                   atcTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAAACCCCAGAACTTACTTATGCTGATTCCGGTGCGAAAATTGTTAATAAAGGTACtt                                                 >  3:243820/1‑149 (MQ=255)
                                                   atcTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAAACCCCAGAACTTACTTATGCTGATTCCGGTGCGAAAATTGTTAATAAAGGTACtt                                                 <  1:877009/149‑1 (MQ=255)
                                                   atcTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAAACCCCAGAACTTACTTATGCTGATTCCGGTGCGAAAATTGTTAATAAAGGTACtt                                                 >  2:83256/1‑149 (MQ=255)
                                                      tGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAAACCCCAGAACTTACTTATGCTGATTCCGGTGCGAAAATTGTTAATAAAGGTACTTGTa                                              >  3:711321/1‑149 (MQ=255)
                                                      tGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAAACCCCAGAACTTACTTATGCTGATTCCGGTGCGAAAATTGTTAATAAAGGTACTTGTa                                              >  4:652045/1‑149 (MQ=255)
                                                        ttCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAAACCCCAGAACTTACTTATGCTGATTCCGGTGCGAAAATTGTTAATAAAGGTACTTGTAGa                                            <  3:767950/149‑1 (MQ=255)
                                                        ttCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAAACCCCAGAACTTACTTATGCTGATTCCGGTGCGAAAATTGTTAATAAAGGTACTTGTAGa                                            >  1:582492/1‑149 (MQ=255)
                                                        ttCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAAACCCCAGAACTTACTTATGCTGATTCCGGTGCGAAAATTGTTAATAAAGGTACTTGTAGa                                            <  1:498958/149‑1 (MQ=255)
                                                          ctcCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAAACCCCAGAACTTACTTATGCTGATTCCGGTGCGAAAATTGTTAATAAAGGTACTTGTAGAtt                                          <  3:701789/149‑1 (MQ=255)
                                                          ctcCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAAACCCCAGAACTTACTTATGCTGATTCCGGTGCGAAAATTGTTAATAAAGGTACTTGTAGAtt                                          <  3:630596/149‑1 (MQ=255)
                                                           tcCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAAACCCTAGAACTTACTTATGCTGATTCCGGTGCGAAAATTGTTAATAAAGGTACTTGTAGATTc                                         >  2:924003/1‑149 (MQ=255)
                                                           tcCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAAACCCCAGAACTTACTTATGCTGATTCCGGTGCGAAAATTGTTAATAAAGGTACTTGTAGATTc                                         <  4:265469/149‑1 (MQ=255)
                                                             cAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAAACCCCAGAACTTACTTATGCTGATTCCGGTGCGAAAATTGTTAATAAAGGTACTTGTAGATTCaa                                       >  4:487544/1‑149 (MQ=255)
                                                             cAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAAACCCCAGAACTTACTTATGCTGATTCCGGTGCGAAAATTGTTAATAAAGGTACTTGTAGATTCaa                                       <  1:344890/149‑1 (MQ=255)
                                                             cAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAAACCCCAGAACTTACTTATGCTGATTCCGGTGCGAAAATTGTTAATAAAGGTACTTGTAGATTCaa                                       <  3:487544/149‑1 (MQ=255)
                                                               aTGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAAACCCCAGAACTTACTTATGCTGATTCCGGTGCGAAAATTGTTAATAAAGGTACTTGTAGATTCAAtt                                     <  1:363/149‑1 (MQ=255)
                                                                tGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAAACCCCAGAACTTACTTATGCTGATTCCGGTGCGAAAATTGTTAATAAAGGTACTTGTAGATTCAATTg                                    <  2:691122/149‑1 (MQ=255)
                                                                tGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAAACCCCAGAACTTACTTATGCTGATTCCGGTGCGAAAATTGTTAATAAAGGTACTTGTAGATTCAATTg                                    <  3:426628/149‑1 (MQ=255)
                                                                tGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAAACCCCAGAACTTACTTATGCTGATTCCGGTGCGAAAATTGTTAATAAAGGTACTTGTAGATTCAATTg                                    <  1:242291/149‑1 (MQ=255)
                                                                tGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAAACCCCAGAACTTACTTATGCTGATTCCGGTGCGAAAATTGTTAATAAAGGTACTTGTAGATTCAATTg                                    <  1:654397/149‑1 (MQ=255)
                                                                tGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAAACCCCAGAACTTACTTATGCTGATTCCGGTGCGAAAATTGTTAATAAAGGTACTTGTAGATTCAATTg                                    <  2:489610/149‑1 (MQ=255)
                                                                       tCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAAACCCCAGAACTTACTTATGCTGATTCCGGTGCGAAAATTGTTAATAAAGGTACTTGTAGATTCAATTGGTCAACg                             >  1:98024/1‑149 (MQ=255)
                                                                       tCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAAACCCCAGAACTTACTTATGCTGATTCCGGTGCGAAAATTGTTAATAAAGGTACTTGTAGATTCAATTGGTCAACg                             >  4:119424/1‑149 (MQ=255)
                                                                         tataaaCTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAAACCCCAGAACTTACTTATGCTGATTCCGGTGCGAAAATGGTTAATAAAGGTACTTGTAGATTCAATTGGTCAACGCa                           >  4:303744/4‑149 (MQ=255)
                                                                         tAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAAACCCCAGAACTTACTTATGCTGATTCCGGTGCGAAAATTGTTAATAAAGGTACTTGTAGATTCAATTGGTCAACGCa                           >  3:758576/1‑149 (MQ=255)
                                                                         tAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAAACCCCAGAACTTACTTATGCTGATTCCGGTGCGAAAATTGTTAATAAAGGTACTTGTAGATTCAATTGGTCAACGCa                           >  1:315456/1‑149 (MQ=255)
                                                                         tAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAAACCCCAGAACTTACTTATGCTGATTCCGGTGCGAAAATTGTTAATAAAGGTACTTGTAGATTCAATTGGTCAACGCa                           >  3:819189/1‑149 (MQ=255)
                                                                         tAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAAACCCCAGAACTTACTTATGCTGATTCCGGTGCGAAAATTGTTAATAAAGGTACTTGTAGATTCAATTGGTCAACGCa                           >  1:928966/1‑149 (MQ=255)
                                                                           aaaaCTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAAACCCCAGAACTTACTTATGCTGATTCCGGTGCGAAAATTGTTAATAAAGGTACTTGTAGATTCAATTGGTCAACGCAAc                         >  3:870553/1‑149 (MQ=255)
                                                                            aaaCTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAAACCCCAGAACTTACTTATGCTGATTCCGGTGCGAAAATTGTTAATAAAGGTACTTGTAGATTCAATTGGTCAACGCAACa                        >  2:908742/1‑149 (MQ=255)
                                                                              aCTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAAACCCCAGAACTTACTTATGCTGATTCCGGTGCGAAAATTGTTAATAAAGGTACTTGTAGATTCAATTGGTCAACGCAACAGt                      >  1:817798/1‑149 (MQ=255)
                                                                                  gTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAAACCCCAGAACTTACTTATGCTGATTCCGGTGCGAAAATTGTTAATAAAGGTACTTGTAGATTCAATTGGTCAACGCAACAGTTAtg                  >  4:723077/1‑149 (MQ=255)
                                                                                   tATTAAGACTATCACTTATTTATGTGATATTGGTTGTCTGGAGATTAAGGGGGCCAGTCTAAACCCCAGAACTTACTTATGCTGATTCCGTTGCGAAAATTGTTAATAAAGGAACATGTAGATTCAATTGGTCAACGCAACAGTTAtgt                 <  3:552440/149‑1 (MQ=37)
                                                                                   tATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAAACCCCAGAACTTACTTATGCTGATTCCGGTGCGAAAATTGTTAATAAAGGTACTTGTAGATTCAATTGGTCAACGCAACAGTTAtgt                 <  4:151724/149‑1 (MQ=255)
                                                                                     ttAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAAACCCCAGAACTTACTTATGCTGATTCCGGTGCGAAAATTGTTAATAAAGGTACTTGTAGATTCAATTGGTCAAc                              <  1:123705/134‑1 (MQ=255)
                                                                                     ttAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAAACCCCAGAACTTACTTATGCTGATTCCGGTGCGAAAATTGTTAATAAAGGTACTTGTAGATTCAATTGGTCAAc                              >  2:123705/1‑134 (MQ=255)
                                                                                     ttAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAAACCCCAGAACTTACTTATGCTGATTCCGGTGCGAAAATTGTTAATAAAGGTACTTGTAGATTCAATTGGTCAACGCAACAGTTATGTGa               >  2:450724/1‑149 (MQ=255)
                                                                                      tAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAAACCCCAGAACTTACTTATGCTGATTCCGGTGCGAAAATTGTTAATAAAGGTACTTGTAGATTCAATTGGTCAACGCAACAGTTATGTGaa              >  3:794165/1‑149 (MQ=255)
                                                                                      aaccaatATCACTTAAATAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAAACCCCAGAACTTACTTATGCTGATTCCGGTGCGAAAATTGTTAATAAAGGTACTTGTAGATTc                                         <  2:350557/116‑1 (MQ=255)
                                                                                      aaccaatATCACTTAAATAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAAACCCCAGAACTTACTTATGCTGATTCCGGTGCGAAAATTGTTAATAAAGGTACTTGTAGATTc                                         >  1:350557/7‑122 (MQ=37)
                                                                                       aaGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAAACCCCAGAACTTACTTATGCTGATTCCGGTGCGAAAATTGTTAATAAAGGTACTTGTAGATTCAATTGGTCAACGCAACAGTTATGTGaaa             >  4:437786/1‑149 (MQ=255)
                                                                                        ccaatATCACTTAAATAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAAACCCCAGAACTTACTTATGCTGATTCCGGTGCGAAAATTGTTAATAAAGGTACTTGTAGATTCAATTGGTCAACGCAACAGTTATGTGaaaa            >  1:131946/5‑149 (MQ=255)
                                                                                        ccaatATCACTTAAATAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAAACCCCAGAACTTACTTATGCTGATTCCGGTGCGAAAATTGTTAATAAAGGTACTTGTAGATTCAATTGGTCAACGCAACAGTTATGTGaaaa            >  3:734289/5‑149 (MQ=255)
                                                                                         gACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCTAGTCTAAACCCCAGAACTTACTTATGCTGATTCCGGTGCGAAAATTGTTAATAAAGGTACTTGTAGATTCAATTGGTCAACGCAACAGTTATGTGAAAAc           <  4:883079/149‑1 (MQ=255)
                                                                                          aCTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAAACCCCAGAACTTACTTATGCTGATTCCGGTGCGAAAATTGTTAATAAAGGTACTTGTAGATTCAATTGGTCAACGCAACAGTTATGTGAAAACa          <  2:98024/149‑1 (MQ=255)
                                                                                            tATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAAACCCCAGAACTTACTTATGCTGATTCCGGTGCGAAAATTGTTAATAAAGGTACTTGTAGATTCAATTGGTCAACGCAACAGTTATGTGAAAACATg        >  1:258388/1‑149 (MQ=255)
                                                                                              tCACTTAAATAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAAACCCCAGAACTTACTTATGCTGATTCCGGTGCGAAAATTGTTAATAAAGGTACTTGTAGATTCAATTGGTCAACGCAACAGTTATGTGAAAACATggg      <  4:734289/149‑1 (MQ=255)
                                                                                                aCTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAAACCCCAGAACTTACTTATGCTGATTCCGGTGCGAAAATTGTTAATAAAGGTACTTGTAGATTCAATTGGTCAACGCAACAGTTATGTGAAAACATGGGGt    >  2:913372/1‑149 (MQ=255)
                                                                                                 cTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAAACCCCAGAACTTACTTATGc                                                                                    >  1:778124/1‑68 (MQ=255)
                                                                                                 cTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAAACCCCAGAACTTACTTATGc                                                                                    <  2:778124/68‑1 (MQ=255)
                                                                                                  ttATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAAACCCCAGAACTTACTTATGCTGATTCCGGTGCGAAAATTGTTAATAAAGGTACTTGTAGATTCAATTGGTCAACGCAACAGTTATGTGAAAACATGGGGTTg  >  1:80444/1‑149 (MQ=255)
                                                                                                  ttATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAAACCCCAGAACTTACTTATGCTGATTCCGGTGCGAAAATTGTTAATAAAGGTACTTGTAGATTCAATTGGTCAACGCAACAGTTATGTGAAAACATGGGGTTg  >  1:260115/1‑149 (MQ=255)
                                                                                                     tttAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAAACCCCAGAACTTACTTATGCTGATTCCGGTGCGAAAATTGTTAATAAAGGt                                                     <  3:661181/95‑1 (MQ=255)
                                                                                                     tttAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAAACCCCAGAACTTACTTATGCTGATTCCGGTGCGAAAATTGTTAATAAAGGt                                                     >  4:661181/1‑95 (MQ=255)
                                                                                                      |                                                                                                                                                
TTTTCCTGTGTGGTACGCCGTCTGCGTTTCTCCGGCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAAACCCCAGAACTTACTTATGCTGATTCCGGTGCGAAAATTGTTAATAAAGGTACTTGTAGATTCAATTGGTCAACGCAACAGTTATGTGAAAACATGGGGTTG  >  NC_000913/1469097‑1469343

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