Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 884,868 T→A 6.6% L399* (TTA→TAA)  mdfA → multidrug efflux system protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913884,8680TA 6.6% 334.9 / 11.8 121L399* (TTA→TAA) mdfAmultidrug efflux system protein
Reads supporting (aligned to +/- strand):  ref base T (53/60);  new base A (6/2);  total (59/62)
Fisher's exact test for biased strand distribution p-value = 1.57e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.05e-01

Only 100 of 123 total aligned reads displayed.

TCTGCCGCGATGGGAATGCTGCAAATGCTGATCTTTACCGTTGGTATTGAAATCAGCAAACATGCCTGGCTGAACGGGGGCAACGGACTGTTTAATCTCTTCAACCTTGTCAACGGAATTTTGTGGCTGTCGCTGATGGTTATCTTTTTAAAAGATAAACAGATGGGAAATTCTCACGAAGGGTAAAAAAATGCCTGACTGCTTTGTGCGATCAGGCATTCTCGAATTAATGGTGATGGTCGTCAATCTGGTGTTCGATAACCATCCCTTCACCTACGCTGGCAAGATGGCGAAC  >  NC_000913/884720‑885014
                                                                                                                                                    |                                                                                                                                                   
tCTGACTCGATTGTAATGCTTCAAATGCTGATCTTTACCGTTGGTATTGAAATCAGCAAAAATGCCTGGATTAACTGGGTCAACGGACTGTATAATCTCTTCAACCTTGTCAACGGAATTTTGTGGCTGTCGCTGATGGTTATCttttt                                                                                                                                                     <  2:939098/149‑1 (MQ=255)
 cTGCCGCGATGGGAATGCTGCAAATGCTGATCTTTACCGTTGGTATTGAAATCAGCAAACATGCCTGGCTGAACGGGGGCAACGGACTGTTTAATCTCTTCAACCTTGTCAACTGAATTTTGTGGATGTAGCTGATGGTTATCTTTTTa                                                                                                                                                    >  2:967972/1‑149 (MQ=255)
 cTGCCGCGATGGGAATGCTGCAAATGCTGATCTTTACCGTTGGTATTGAAATCAGCAAACATGCCTGGCTGAACGGGGGCAACGGACTGTTTAATCTCTTCAACCTTGTCAACGGAATTTTGTGGCTGTCGCTGATGGTTATCTTTTTa                                                                                                                                                    >  2:1045304/1‑149 (MQ=255)
 cTGCCGCGATGGGAATGCTGCAAATGCTGATCTTTACCGTTGGTATTGAAATCAGCAAACATGCCTGGCTGAACGGGGGCAACGGACTGTTTAATCTCTTCAACCTTGTCAACGGAATTTTGTGGCTGTCGCTGATGGTTATCTTTTTa                                                                                                                                                    >  1:1552937/1‑149 (MQ=255)
   gCCGCGATGGGAATGCTGCAAATGCTGATCTTTACCGTTGGTATTGAAATCAGCAAACATGCCTGGCTGAACGGGGGCAACGGACTGTTTAATCTCTTCAACCTTGTCAACGGAATTTTGTGGCTGTCGCTGATGGTTATCTTTTTaaa                                                                                                                                                  <  2:138110/149‑1 (MQ=255)
   gCCGCGATGGGAATGCTGCAAATGCTGATCTTTACCGTTGGTATTGAAATCAGCAAACATGCCTGGCTGAACGGGGGCAACGGACTGTTTAATCTCTTCAACCTTGTCAACGGAATTTTGTGGCTGTCGCTGATGGTTATCTTTTTaaa                                                                                                                                                  <  2:1201229/149‑1 (MQ=255)
    ccGCGATGGGAATGCTGCAAATGCTGATCTTTACCGTTGGTATTGAAATCAGCAAACATGCCTGGCTGAACGGGGGCAACGGACTGTTTAATCTCTTCAACCTTGTCAACGGAATTTTGTGGCTGTCGCTGATGGTTATCTTTTTaaaa                                                                                                                                                 >  1:1028787/1‑149 (MQ=255)
        gATGGGAATGCTGCAAATGCTGATCTTTACCGTTGGTATTGAAATCAGCAAACATGCCTGGCTGAACGGGGGCAACGGACTGTTTAATCTCTTCAACCTTGTCAACGGAATTTTGTGGCTGTCGCTGATGGTTATCTTTTTAAAAGATa                                                                                                                                             <  1:1544551/149‑1 (MQ=255)
        gATGGGAATGCTGCAAATGCTGATCTTTACCGTTGGTATTGAAATCAGCAAACATGCCTGGCTGAACGGGGGCAACGGACTGTTTAATCTCTTCAACCTTGTCAACGGAATTTTGTGGCTGTCGCTGATGGTTATCTTTTTAAAAGATa                                                                                                                                             >  2:318648/1‑149 (MQ=255)
         aTGGGAATGCTGCAAATGCTGATCTTTACCGTTGGTATTGAAATCAGCAAACATGCCTGGCTGAACGGGGGCAACGGACTGTTTAATCTCTTCAACCTTGTCAACGGAATTTTGTGGCTGTCGCTGATGGTTATCTTTTAAAAAGATaa                                                                                                                                            <  2:821519/149‑1 (MQ=255)
          tGGGAATGCTGCAAATGCTGATCTTTACCGTTGGTATTGAAATCAGCAAACATGCCTGGCTGAACGGGGGCAACGGACTGTTTAATCTCTTCAACCTTGTCAACGGAATTTTGTGGCTGTCGCTGATGGTTATCTTTTTAAAAGATaaa                                                                                                                                           <  1:1560739/149‑1 (MQ=255)
          tGGGAATGCTGCAAATGCTGATCTTTACCGTTGGTATTGAAATCAGCAAACATGCCTGGCTGAACGGGGGCAACGGACTGTTTAATCTCTTCAACCTTGTCAACGGAATTTTGTGGCTGTCGCTGATGGTTATCTTTTAAAAAGATaac                                                                                                                                           <  2:1472108/149‑2 (MQ=255)
           gggAATGCTGCAAATGCTGATCTTTACCGTTGGTATTGAAATCAGCAAACATGCCTGGCTGAACGGCGGCAACGGACTGTTAAATCTCTTCAACCTTGTCAACGGAATTTTGTGGCTGTCGCTGATGGTTATCTTTTTAAAAGATAAAc                                                                                                                                          >  2:1443948/1‑149 (MQ=255)
             gAATGCTGCAAATGCTGATCTTTACCGTTGGTATTGAAATCAGCAAACATGCCTGGCTGAACGGGGGCAACGGACTGTTTAATCTCTTCAACCTTGTCAACGGAATTTTGTGGCTGTCGCTGATGGTTATCTTTTTAAAAGATAAACAg                                                                                                                                        <  1:1045304/149‑1 (MQ=255)
              aaTGCTGCAAATGCTGATCTTTACCGTTGGTATTGAAATCAGCAAACATGCCTGGCTGAACGGGGGCAACGGACTGTTTAATCTCTTCAACCTTGTCAACGGAATTTTGTGGCTGTCGCTGATGGTTATCTTTTTAAAAGATAAACAGa                                                                                                                                       >  2:699859/1‑149 (MQ=255)
                   tgcAAATGCTGATCTTTACCGTTGGTATTGAAATCAGCAAACATGCCTGGCTGAACGGGGGCAACGGACTGTTTAATCTCTTCAACCTTGTCAACGGAATTTTGTGGCTGTCGCTGATGGTTATCTTTTTAAAAGATAAACAGATGGGa                                                                                                                                  >  1:1428098/1‑149 (MQ=255)
                     cAAATGCTGATCTTTACCGTTGGTATTGAAATCAGCAAACATGCCTGGCTGAACGGGGGCAACGGACTGTTTAATCTCTTCAACCTTGTCAACGGAATTTTGTGGCTGTCGCTGATGGTTATCTTTTTAAAAGATAAACAGATGGGaaa                                                                                                                                >  2:1333856/1‑149 (MQ=255)
                     cAAATGCTGATCTTTACCGTGGGTATTGAAATCAGCAAACATGCCTGGCTGAACGGGGGCAACGGACTGTTTAATCTCTTCAACCTTGTCAACGGAATTTTGTGGCTGTCGCTGATGGTTATCTTTTTAAAAGATAAACAGATGGGaaa                                                                                                                                <  2:634110/149‑1 (MQ=255)
                         tgatGATCTTTACCGTTGGTAATGAAATCATCAAACATGCCTTGCTGAACGGGGGCAACGGACTGTTTAATCTCTTCAACCTTGTCAACGGAATTTTGTGGCTGTCGCTGATGGTTATCTTTTTAAAAGATAAACAGATGGGAAATtct                                                                                                                            <  2:806211/146‑1 (MQ=255)
                           gtGATCATTACAGTTGGTATTGAACTCAGCAAACATGCCTGGCTGAACGGTGGCAACTGACTGTTTAATCTCTTCAACCTTGTCAACGGAATTTTGTGGCTGTCGCTGATGGTTATCTTTTTAAAAGATAAACAGAGGGGAAATTCTCa                                                                                                                          <  1:251707/148‑1 (MQ=255)
                             gATCTTTACCGTTGGTATTGAAATCAGCAAACATGCCTGGCTGAACGGGGGCAACGGACTGTTTAATCTCTTCAACCTTGTCAACGGAATTTTGTGGCTGTCGCTGATGGTTATCTTTTTAAAAGATAAACAGATGGGAAATTCTCACg                                                                                                                        >  1:316282/1‑149 (MQ=255)
                                 tttAACGTTGGTATTGAAATCAGCAAACATGCCTCGCTGAACGGGGGCAACGGACTGTTTAATCTCTTCAACCTTGTCAACGGAATTTAGAGGCTGTCGCTGATGGTTATCTTTTAAAAAGATTAACAAATGGGAAATTCTCACGAAgg                                                                                                                    >  1:308387/1‑149 (MQ=255)
                                  ttACCGTTGGTATTGAAATCAGCAAACATGCCTGGCTGAACGGGGGCAACGGACTGTTTAATCTCTTCAACCTTGTCAACGGAATTTTGTGGCTGTCGCTGATGGTTATCTTTTTAAAAGATAAACAGATGGGAAATTCTCACGAAggg                                                                                                                   >  2:1209547/1‑149 (MQ=255)
                                        ttGGTATTGAAATCAGCAAACATGCCTGGCTGAACGGGGGCAACGGACTGTTTAATCTCTTCAACCTTGTCAACGGAATTTTGTGGCTGTCGCTGATGGTTATCTTTTTAAAAGATAAACAGATGGGAAATTCTCACGAAGGGTaaaaa                                                                                                            >  1:659337/1‑149 (MQ=255)
                                        ttGGTATTGAAATCAGCAAACATACCTGGCTGAACGGGGGCAACGGACTGTTTAATCTCTTCAACCTTGTCAACGGAATTTTGTGGATGTCGCTGATGGTTATCTTTTTAAAAGATAAACAGATGGGAAATTCTCACGAAGGTTaaaaa                                                                                                            >  2:1036681/1‑149 (MQ=255)
                                             tttGAAATCAGCAAACATGCCTGGCTGAACGGGGGCAACGGACTGTTTAATCTCTTCAACCTTGTCAACGGAATTTTGTGGCTGTCGCTGATGGTTATCTTTTTAAAAGATAAACAGATGGGAAATTCTCACGAAGGGTAAAAAAATGc                                                                                                       <  2:1337615/148‑1 (MQ=255)
                                              ttGAAATCAGCAAACATGCCTGGCTGAACGGGGGCAACGGACTGTTTATTCTCTTCAACCTTGTCAACGGAATTTTGTGGCTGTCGCTGATGGTTATCTTTTTAAAAGATAAACAGATGGGAAATTCTCACGAAGGGTAAAAAAATGcc                                                                                                      >  1:102960/1‑149 (MQ=255)
                                              ttGAAATCAGCAAACATGCCTGGCTGAACGGGGGCAACGGACTGTTTAATCTCTTCAACCTTGTCAACGGAATTTTGTGGCTGTCGCTGATGGTTATCTTTTTAAAAGATAAACAGATGGGAAATTCTCACGAAGGGTAAAAAAATGcc                                                                                                      >  1:254835/1‑149 (MQ=255)
                                              ttGAAATCAGCAAACATGCCTGGCTGAACGGGGGCAACGGACTGTTTAATCTCTTCAACCTTGTCAACGGAATTTTGTGGCTGTCGCTGATGGTTATCTTTTTAAAAGATAAACAGATGGGAAATTCTCACGAAGGGTAAAAAAATGcc                                                                                                      >  2:641987/1‑149 (MQ=255)
                                              ttGAAATCAGCAAACATGCCTGGCTGAACGGGGGCAACGGACTGTTTAATCTCTTCAACCTTGTCAACGGAATTTTGTGGCTGTCGCTGATGGTTATCTTTTTAAAAGATAAACAGATGGGAAATTCTCACGAAGGGTAAAAAAATGcc                                                                                                      <  1:1249071/149‑1 (MQ=255)
                                                 aaaTCAGCAAACATGCCTGGCTGAACGGGGGCAACGGACTGTTTAATCTCTTCAACCTTGTCAACGGAATTTTGTGGCTGTCGCTGATGGTTATCTTTTTAAAAGATAAACAGATGGGAAATTCTCACGATGGGTAAAAAAATCCATGa                                                                                                   >  1:746247/1‑149 (MQ=255)
                                                 aaaTCAGCAAACATGCATGGCTGTACGGGGGCAACGGACTGTTTTATCTCTTCAACCTTGTCAACGGAATTTTGTGGCTGTCGCTGATGGTTATCTTTTTAAAAGATAAACAGATGGGAAATTCTCACGAAGGGTAAAAAAATGCCTGa                                                                                                   <  1:1020871/149‑1 (MQ=255)
                                                      aGCAAACATGCCTGGCTGAACGGGTGCAACGGACTGTTTAATCTCTTCAACCTTGTCAACGGAATTTTGTGGCTGTCGCTGATGGTTATATTTTTAAAAGATAAACAGATGGTAAATTCTCACGAAGGGTAAAAAACTGCCTGACTGCt                                                                                              >  2:73819/1‑149 (MQ=255)
                                                      aGCAAACATGCCTGGCTGAACGGGGGCAACGGACTGTTTAATCTCTTCAACCTTGTCAACGGAATTTTGTGGCTGTCGCTGATGGTTATCTTTTTAAAAGATAAACAGATGGGAAATTCTCGCGAAGGGTAAAAAAATGCCTGACTGCt                                                                                              <  1:565048/149‑1 (MQ=255)
                                                       gCAAACATGCCTGGCTGAACGGGGGCAACGGACTGTTTAATCTCTTCAGCCTTGTCAACGGAATTTTGTGGCTGTCGCTGATGGTTATCTTTTTAAAAGATAAACAGATGGGAAATTCTCACGAAGGGTAAAAAAATGCCTGACTGCtt                                                                                             >  1:1259163/1‑149 (MQ=255)
                                                       gCAAACATGCCTGGCTGAACGGGGGCAACGGACTGTTTAATCTCTTCAACCTTGTCAACGGAATTTTGTGGCTGTCGCTGATGGTTATCTTTTTAAAAGATAAACAGATGGGAAATTCTCACGAAGGGTAAAAAAATGCCTGACTGCtt                                                                                             >  2:223030/1‑149 (MQ=255)
                                                       gCAAACATGCCTGGCTGAACGGGGGCAACGGACTGTTTAATCTCTTCAACCTTGTCAACGGAATTTTGTGGCTGTCGCTGATGGTTATCTTTTTAAAAGATAAACAGATGGGAAATTCTCACGAAGGGTAAAAAAATGCCTGACTGCtt                                                                                             >  2:793930/1‑149 (MQ=255)
                                                        cAAACATGCCTGGCTGAACGGGGGCAACGGACTGTTTAATCTCTTCACCCTTGTCAACGGAATTTTGTGGCTGTCGCTGATGGTTATCTTTTTAAAAGATAAACAGATGGGAAATTCTCACGAAGGGTAAAAAAAATGCCTGACTGCtt                                                                                             >  2:293109/1‑149 (MQ=255)
                                                        cAAACATGCCTGGCTGAACGGGGGCAACGGACTGTTTAATCTCTTCAACCTTGTCAACGGAATTTTGTGGCTGTCGCTGATGGTTATCTTTTTAAAAGATAAACAGATGGGAAATTCTCACGAAGGGTAAAAAAATGCCTGACTGCttt                                                                                            >  1:1413174/1‑149 (MQ=255)
                                                        cAAACATGCCTGGCTGAACGGGGGCAACGGACTGTTTAATCTCTTCAACCTTGTCAACGGAATTTTGTGGCTGTCGCTGATGGTTATCTTTTTAAAAGATAAACAGATGGGAAATTCTCACGAAGGGTAAAAAAATGCCTGACTGCttt                                                                                            >  2:815327/1‑149 (MQ=255)
                                                           aCATGCCTGGCTGAACGGGGGCAACGGACTGTTTAATCTCTTCAACCTTGTCAACGGAATTTTGTGGCTGTCGCTGATGGTTATCTTTTTAAAAGATAAACAGATGGGAAATTCTCACGAAGGGTAAAAAAATGCCTGACTGCTTtgtg                                                                                         >  2:1383377/1‑149 (MQ=255)
                                                                 cTGGCTGAACGGGGGCAACGGACTGTTTAATCTCTTCAACCTTGTCAACGGAATTTTGTGGCTGTCGCTGATGGTTATCTTTTTAAAAGATAAACAGATGGGAAATTCTCACGAAGGGTAAAAAAATGCCTGACTGCTTTGTGCGATCa                                                                                   <  1:1279090/149‑1 (MQ=255)
                                                                 cTGGCTGAACGGGGGCAACGGACTGTTTAATCTCTTCAACCTTGTCAACGGAATTTTGTGGCTGTCGCTGATGGTTATCTTTTTAAAAGATAAACAGATGGGAAATTCTCACGAAGGGTAAAAAAATGCCTGACTGCTTTGTGCGATCa                                                                                   >  1:330125/1‑149 (MQ=255)
                                                                  tGGCTGAACGGGGGCAACGGACTGTTTAATCTCTTCAACCTTGTCAACGGAATTTTGTGGCTGTCGCTGATGGTTATCTTTTTAAAAGATAAACAGATGGGAAATTCTCACGAAGGGTAAAAAAATGCCTGACTGCTTTGTGCGATCAg                                                                                  >  2:247324/1‑149 (MQ=255)
                                                                  tGGCTGAACGGGGGCAACGGACTGTTTAATCTCTTCAACCTTGTCAACGGAATTTTGTGGCTGTCGCTGATGGTTATCTTTTTAAAAGATAAACAGATGGGAAATTCTCACGAAGGGTAAAAAAATGCCTGACTGCTTTGTGCGATCAg                                                                                  <  1:491752/149‑1 (MQ=255)
                                                                   ggCTGAACGGGGGCAACGGACTGTTTAATCTCTTCAACCTTGTCAACGGAATTTTGTGGCTGTCGCTGATGGTTATCTTTTTAAAAGATAAACAGATGGGAAATTCTCACGAAGGGTAAAAAAATGCCTGACTGCTTTGTGCGATCAgg                                                                                 <  2:886332/149‑1 (MQ=255)
                                                                    gCTGAACGGGGGCAACGGACTGTTTAATCTCTTCAACCTTGTCAACGGAATTTTGTGGCTGTCGCTGATGGTTATCTTTTTAAAAGATAAACAGATGGGAAATTCTCACGAAGGGTAAAAAAATGCCTGACTGCTTTGTGCGATCAGGc                                                                                <  1:1217009/149‑1 (MQ=255)
                                                                    gCTGAACGGGGGCAACGGACTGTATAATCTCTTCAACCTTGTCAACGGAATTCTGTGGCTGTCGCTGATGGTTATCTTTTTAAAAGATAAACAGATGGGAAATTCTCACGAAGGGTAAAAAAATGCCTGACTGCTTTGTGCGATCAGGc                                                                                <  1:715771/149‑1 (MQ=255)
                                                                      tGAACGGGGGCAACGGACTGTTTAATCTCTTCAACCTTGTCAACGGAATTTTGTGGCGGTCGCTGATGGTTTTCTTTTTAAAATATAAACAGATGGGAAATTATCACGAAGGGTAAAAAAATGCCTGCCTGCTTTGTGCGATCAGGCAt                                                                              >  2:54093/1‑149 (MQ=255)
                                                                          cGGGGGCAACGGACTGTTTAATCTCTTCAACCTTGTCAACGGAATTTTGTGGCTGTCGCTGATGGTTATCTTTTTAAAAGATAAACAGATGGGAAATTCTCACGAAGGGTAAAAAAATGCCTGACTGCTTTGTGCGATCAGGCATtctc                                                                          <  1:1193811/149‑1 (MQ=255)
                                                                            ggggCAACGGACTGTTTAATCTCTTCAACCTTGTCAACGGAATTTTGTGGCTGTCGCTGATGGTTATCTTTTTAAAAGATAAACAGATGGGAAATTCTCACGAAGGGTAAAAAAATGCCTGACTGCTTTGTGCGATCAGGCATTCTCGa                                                                        <  1:483406/149‑1 (MQ=255)
                                                                             gttcAACTGACTGTTTAATCTATTATACATTTTCAACTGAATTTTGTTGCTGTCTCTGATGGTTATCTTTTTAAAAGATAATCAGATGGGAAATTCTCACTAAGGGTAAAAAAATGCCTGACTGCTTTGTGCGATCAGGCATTCTCGaa                                                                       <  1:1229619/146‑1 (MQ=255)
                                                                             gggCAACGGACTGTTTAATCTCTTCAACCTTGTCAACGGAATTTTGTGGCTGTCGCTGATGGTTATCTTTTTAAAAGATAAACAGATGGGAAATTCTCACGAAGGGTAAAAAAATGCCTGACTGCTTTGTGCGATCAGGCATTCTCGaa                                                                       <  2:706321/149‑1 (MQ=255)
                                                                              ggCAACGGACTGTTTAATCTCTTCCACCTTGTCAACGGAATTTTGTGGCTGTCGCTGATGGTTATCTTTTTAAAAGATAAACAGATGGGAAATTCTCACGAAGGGTAAAAAAATGCCTGACTGCTTTGTGCGATCAGGCATTCTCGAAt                                                                      >  2:1441883/1‑149 (MQ=255)
                                                                                 aaCGGACTGTTTAATCTCTTCAACCTTGTCAACGGAATTTTGTGGCTGTCGCTGATGGTTATCTTTTTAAAAGATAAACAGATGGGAAATTCTCACGAAGGGTAAAAAAATGCCTGACTGCTTTGTGCGATCAGGCATTCTCGAATTaa                                                                   <  1:793930/149‑1 (MQ=255)
                                                                                  aCGGACTGTGTAATCTCTTCAACCTTGTCAACGGAATTTTGTGGCTGTCGCTGATGGTTATCTTTTTAAAAGATAAACAGATGGGAAATTCTCACGAAGGGTAAAAAAATGCCTGACTGCTTTGTGCGATCAGGCATTCTCGAATTAAt                                                                  <  2:529822/149‑1 (MQ=255)
                                                                                   cGGACTGTTTAATCTCTTCAACCTTGTCAACGGAATTTTGTGGCTGTCGCTGATGGTTATCTTTTTAAAAGATAAACAGATGGGAAATTCTCACGAAGGGTAAAAAAATGCCTGACTGCTTTGTGCGATCAGGCATTCTCGAATTAATg                                                                 >  1:835687/1‑149 (MQ=255)
                                                                                   cGGACTGTTTAATCTCTTCAACCTTGTCAACGGAATTTTGTGGCTGTCGCTGATGGTTATCTTTTTAAAAGATAAACAGATGGGAAATTCTCACGAAGGGTAAAAAAATGCCTGACTGCTTTGTGCGATCAGGCATTCTCGAATTAATg                                                                 <  1:883775/149‑1 (MQ=255)
                                                                                      aCTGTTTAATCTCTTCAACCTTGTCAACGGAATTTTGTGGCTGTCGCTGATGGTTATCTTTTTAAAAGATAAACAGATGGGAAATTCTCACGAAGGGTAAAAAAATGCCTGACTGCTTTGTGCGATCAGGCATTCTCGAATTAATGGTg                                                              >  2:49791/1‑149 (MQ=255)
                                                                                      aCTGTTTAATATCTTCAACCTTTTCAACGGAATTTTGTGGCTGTCGCTGATGGTTATCTTTTTAAAAGATAAACAGATGGGAAATTCTCACTAAGGGTAAAAAAATGCCTGACTGCTTTGTGCGATCAGGCATTCTCGAATTAATGGTg                                                              <  1:967972/149‑1 (MQ=255)
                                                                                       cTGTTTAATCTCTTCAACCTTGTCAACGGAATTTTGTGGCTGTCGCTGATGGTTATCTTTTTAAAAGATAAACAGATGGGAAATTCTCACGAAGGGTAAAAAAATGCCTGACTGCTTTGTGCGATCAGGCATTCTCGAATTAATGGTGa                                                             <  2:1159548/149‑1 (MQ=255)
                                                                                          tttAATCTCTTCAACCTTGTCAACGGAATTTTGTGGCTGTCGCTGATGGTTATCTTTTTAAAAGATAAACAGATGGGAAATTCTCACGAAGGGTAAAAAAATGCCTGACTGCTTTGTGCGATCAGGCATTCTCGAATTAATGGTGATgg                                                          <  1:1396739/149‑1 (MQ=255)
                                                                                          tttAATCTCTTCAACCTTGTCAACGGAATTTTGTGGCTGTCGCTGATGGTTATCTTTTTAAAAGATAAACAGATGGGAAATTCTCACGAAGGGTAAAAAAATGCCTGACTGCTTTGTGCGATCAGGCATTCTCGAATTAATGGTGATgg                                                          >  2:848185/1‑149 (MQ=255)
                                                                                            tAATCTCTTCAACCTTGTCAACGGAATTTTGTGGCTGTCGCTGATGGTTATCTTTTTAAAAGATAAACAGATGGGAAATTCTCACGAAGGGTAAAAAAATGCCTGACTGCTTTGTGCGATCAGGCATTCTCGAATTAATGGTGATGgtc                                                        <  1:460130/149‑1 (MQ=255)
                                                                                               tctcTTCAACCTTGTCAACGGAATTTTGTGGCTGTCGCTGATGGTTATCTTTTTAAAAGATAAACAGATGGGAAATTCTCACGAAGGGTAAAAAAATGCCTGACTGCTTTGTGCGATCAGGCATTCTCGAATTAATGGTGATGgtcgtc                                                     >  2:598189/1‑149 (MQ=255)
                                                                                               tctcTTCAACCTTGTCAACGGAATTTTGTGGCTGTCGCTGATGGTTATCTTTTTAAAAGATAAACAGATGGGAAATTCTCACGAAGGGTAAAAAAATGCCTGACTGCTTTGTGCGATCAGGCATTCTCGAATTAATGGTGATGgtcgtc                                                     >  2:773201/1‑149 (MQ=255)
                                                                                               tctcTTCAACCTTGTCAACGGAATTTTGTGGCTGTCGCTGATGGTTATCTTTTTAAAAGATAAACAGATGGGAAATTCTCACGAAGGGTAAAAAAATGCCTGACTGCTTTGTGCGATCAGGCATTCTCGAATTAATGGTGATGgtcgtc                                                     >  2:1356631/1‑149 (MQ=255)
                                                                                                 tctTCAACCTTGTCAACGGAATTTTGTGGCTGTCGCTGATGGTTATCTTTTTAAAAGATAAACAGATGGGAAATTCTCACGAAGGGTAAAAAAATGCCTGACTGCTTTGTGCGATCAGGCATTCTCGAATTAATGGTGATGGTCGTCaa                                                   <  1:1351103/149‑1 (MQ=255)
                                                                                                   ttCAACCTTGTCAACGGAATTTTGTGGCTGTCGCTGATGGTTATCTTTTTAAAAGATAAACAGATGGGAAATTCTCACGAAGGGTAAAAAAATGCCTGACTGCTTTGTGCGATCAGGCATTCTCGAATTAATGGTGATGTTCGTCAATc                                                 >  2:406245/1‑149 (MQ=255)
                                                                                                   ttCAACCTTGTCAACGGAATTTTGTGGCTGTCGCTGATGGTTATCTTTTTAAAAGATAAACAGATGGGAAATTCTCACGAAGGGTAAAAAAATGCCTGACTGCTTTGTGCGATCAGGCATTCTCGAATTAATGGTGATGGTCGTCAATc                                                 <  1:820948/149‑1 (MQ=255)
                                                                                                     cAACCTTGTCAACGGAATTTTGTGGCTGTCGCTGATGGTTATCTTTTTAAAAGATAAACAGATGGGAAATTCTCACGAAGGGTAAAAAAATGCCTGACTGCTTTGTGCGATCAGGCATTCTCGAATTAATGGTGATGGTCGTCAATCTg                                               <  1:752103/149‑1 (MQ=255)
                                                                                                     cAACCTTGTCAACGGAATTTTGTGGCTGTCGCTGATGGTTATCTTTTTAAAAGATAAACAGATGGGAAATTCTCACGAAGGGTAAAAAAATGCCTGACTGCTTTGTGCGATCAGGCATTCTCGAATTAATGGTGATGGTCGTCAATCTg                                               <  2:316282/149‑1 (MQ=255)
                                                                                                      aaCCTTGTCAACGGAATTTTGTGGCTGTCGCTGATGGTTATCTTTTTAAAAGATAAACAGATGGGAAATTCTCACGAAGGGTAAAAAAATGCCTGACTGCTTTGTGCGATCAGGCATTCTCGAATTAATGGTGATGGTCGTCAATCTgg                                              <  1:1333856/149‑1 (MQ=255)
                                                                                                      aaCCTTGTCAACGGAATTTTGTGGCTGTCGCTGATGGTTATCTTTTTAAAAGATAAACAGATGGGAAATTCTCACGAAGGGTAAAAAAATGCCTGACTGCTTTGTGCGATCAGGCATTCTCGAATTAATGGTGATCGTCGTCAATCTgg                                              <  2:1037170/149‑1 (MQ=255)
                                                                                                       aCCTTGTCAACGGAATTTTGTGGCTGTCGCTGATGGTTATCTTTTTAAAAGATAAACAGATGGGAAATTCTCACGAAGGGTAAAAAAATGCCTGACTGCTTTGTGCGATCAGGCATTCTCGAATTAATGGTGATGGTCGTCAATCTGgt                                             >  1:817714/1‑149 (MQ=255)
                                                                                                        ccTTGTCAACGGAATTTTGTGGCTGTCGCTGATGGTTATCTTTTTAAAAGATAAACAGATGGGAAATTCTCACGAAGGGTAAAAAAATGCCTGACTGCTTTGTGCGATCAGGCATTCTCGAATTAATGGTGATGGTCGTCAATCTGgtg                                            >  2:1480946/1‑149 (MQ=255)
                                                                                                         cTTGTCAACGGAATTTTGTGGCTGTCGCTGATGGTTATCTTTTTAAAAGATAAACAGATGGGAAATTCTCACGAAGGGTAAAAAAATGCCTGACTGCTTTGTGCGATCAGGCATTCTCGAATTAATGGTGATGGTCGTCAATCTGgtgt                                           <  1:1354318/149‑1 (MQ=255)
                                                                                                          ttGTCAACGGAATTTTGTGGCTGTCGCTGATGGTTATCTTTTTAAAAGATAAACAGATGGGAAATTCTCACGAAGGGTAAAAAAATGCCTGACTGCTTTGTGCGATCAGGCATTCTCGAATTAATGGTGATGGTCGTCAATCTGGTGtt                                          <  1:815327/149‑1 (MQ=255)
                                                                                                           tGTCAACGGAATTTTGTGGCTGTCGCTGATGGTTATCTTTTTAAAAGATAAACAGATGGGAAATTCTCACGAAGGGTAAAAAAATGCCTGACTGCTTTGGGCGATCAGGCATTCTCGAATAAATGGTGATGGTCGTCAATCTGGTGTTc                                         <  2:276853/149‑1 (MQ=255)
                                                                                                              cAACGGAATTTTGTGGCTGTCGCTGATGGTTATCTTTTTAAAAGATAAACAGATGGGAAATTCTCACGAAGGGTAAAAAAATGCCTGACTGCTTTGTGCGATCAGGCATTCTCGAATTAATGGTGATGGTCGTCAATCTGGTGTTCGAt                                      <  2:1428098/149‑1 (MQ=255)
                                                                                                                aCGGAATTTTGTGGCTGTCGCTGATGGTTATCTTTTTAAAAGATAAACAGATGGGAAATTCTCACGAAGGGTAAAAAAATGCCTGACTGCTTTGTGCGATCAGGCATTCTCGAATTAATGGTGATGGTCGTCAATCTGGTGTTCGATaa                                    >  1:1557097/1‑149 (MQ=255)
                                                                                                                 cGGAATTTTGTGGCTGTCGCTGATGGTTATCTTTTTAAAAGATAAACAGATGGGAAATTCTCACGAAGGGTAAAAAAATGCCTGACTGCTTTGTGCGATCAGGCATTCTCGAATTAATGGTGATGGTCGTCAATCTGGTGTTCGATAAc                                   >  1:519498/1‑149 (MQ=255)
                                                                                                                    aaTTTTGTGGCTGTCGCTGATGGTTATCTTTTTAAAAGATAAACAGATGGGAAATTCTCACGAAGGGTAAAAAAATGCCTGACTGCTTTGTGCGATCAGGCATTCTCGAATTAATGGTGATGGTCGTCAATCTGGTGTTCGATAACCAt                                <  2:538861/149‑1 (MQ=255)
                                                                                                                     aTTTTGTGGCTGTCGCTGATGGTTATCTTTTTAAAAGATAAACAGATGGGAAATTCTCACGAAGGGTAAAAAAAATGCCTGACTGCTTTGTGCGATCAGGCATTCTCGAATTAATGGTGATGGTCGTCAATCTGGTGTTCGATAACCAt                                <  2:1543481/149‑1 (MQ=255)
                                                                                                                      ttttGTGGCTGTCGCTGATGGTTATCTTTTTAAAAGATAAACAGATGGGAAATTCTCACGAAGGGTAAAAAAATGCCTGACTGCTTTGTGCGATCAGGCATTCTCGAATTAATGGTGATGGTCGTCAATCTGGTGTTCGATAACCTTcc                              <  1:690769/149‑1 (MQ=255)
                                                                                                                          gtgGCTGTCGCTGATGGTTATCTTTTTAAAAGATAAACAGATGGGAAATTCTCACGAAGGGTAAAAATATGCCTGACTGCTTTGTGCGATCAGGCATTCTCGAATTAATGGTGATGGTCGTCAATCTGGTGTTCGATAACCATCCCTTc                          <  2:102960/149‑1 (MQ=255)
                                                                                                                            ggatGTCGCTGATGGTTATCTTTTTAAAAGATAAACAGATGGGAAATTCTCACGAATGTTAAAAAAATGCCTGACTGCTTTGTGCGATCAGGCATTCTCGAATTAATGGTGATGGTCGTCAATCTGGTGTTCGATAACCATCCCTTCAc                        <  1:1443948/146‑1 (MQ=255)
                                                                                                                            ggCTGTCGCTGATGTTTATCTTTTTAAAAGATAAACAGATGGGAAATTCTCACGAAGGGTAAAAAAATGCCTGACTGCTTTGTGCGATCAGGCATTCTCGAATTAATGGGGATGGTCGTCAATCTGGTGTTCGATAACCATCCCTTCAc                        <  2:693884/149‑1 (MQ=255)
                                                                                                                            ggCTGTCGCTGATGGTTATCTTTTTAAAAGATAAACAGATGGGAAATTCTCACGAAGGGTAAAAAAATGCCTGACTGCTTTGTGCGATCAGGCATTCTCGAATTAATGGTGATGGTCGTCAATCTGGTGTTCGATAACCATCCCTTCAc                        <  1:1349539/149‑1 (MQ=255)
                                                                                                                                    catctgtttATCTTTTAAAAAGATAAACAGATGGGAAATTCTCACGAAGGGTAAAAAAATGCCTGACTGCTTTGTGCGATCAGGCATTCTCGAATTAATGGTGATGGTCGTCAATCTGGTGTTCGATAACCATCCCTTCACCTACGCTg                >  2:593925/8‑149 (MQ=255)
                                                                                                                                    catctgtttATCTTTTAAAAAGATAAACAGATGGGAAATTCTCACGAAGGGTAAAAAAATGCCTGACTGCTTTGTGCGATCAGGCATTCTCGAATTAATGGTGATGGTCGTCAATCTGGTGTTCGATAACCATCCCTTCACCTACGCTg                >  1:942432/8‑149 (MQ=255)
                                                                                                                                     tGATGGTTATCTTTTTAAAAGATAAACAGATGGGAAATTCTCACGAAGGGTAAAAAAATGCCTGACTGCTTTGTGCGATCAGGCATTCTCGAATTAATGGTGATGGTCGTCAATCTGGTGTTCGATAACCATCCCTTCACCTACGCTgg               >  2:1366040/1‑149 (MQ=255)
                                                                                                                                          tttATCTTTTAAAAAGATAAACAGATGGGAAATTCTCACGAAGGGTAAAAAAATGCCTGACTGCTTTGTGCGATCAGGCATTCTCGAATTAATGGTGATGGTCGTCAATCTGGTGTTCGATAACCATCCCTTCACCTACGCTGGCAAGa          >  2:1593284/2‑149 (MQ=255)
                                                                                                                                             aTCTTTTTAAAAGATAAACAGATGGGAAATTCTCACGAAGGGTAAAAAAATGCCTGACTGCTTTGTGCGATCAGGCATTCTCGAATTAATGGTGATGGTCGTCAATCTGGTGTTCGATAACCATCCCTTCACCTACGCTGGCAAGATgg       >  1:490540/1‑149 (MQ=255)
                                                                                                                                             aTCTTTTAAAAAGATAAACAGATGGGAAATTCTCACGAAGGGTAAAAAAATGCCTGACTGCTTTGTGCTATCAGGCATTCTCGAATTAATGGTGATGGTCGTCAATCTGGTGTTCGATAACCATCCCTTCACCTACGCTGGCAAGATgg       >  1:573665/1‑149 (MQ=255)
                                                                                                                                             aTCTTTTAAAAAGATAAACAGATGGGAAATTCTCACGAAGGGTAAAAAAATGCCTGACTGCTTTGTGCAATCAGGCATTCTCGAATTAATGGTGATGGTCGTCAATCTGGTGTTCGATAACCATCCCTTCACCTACGCTGGCAAGATgg       >  1:568074/1‑149 (MQ=255)
                                                                                                                                              tCTTTTTAAAAGATAAACAGATGGGAAATTCTCACGAAGGGTAAAAAAATGCCTGACTGCTTTGTGCGATCAGGCATTCTCGAATTAATGGTGATGGTCGTCAATCTGGTGTTCGATAACCATCCCTTCACCTACGCTGGCAAGATGGc      <  2:330125/149‑1 (MQ=255)
                                                                                                                                              tCTTTTTAAAAGATAAACAGATGGGAAATTCTCACGAAGGGTAAAAAAATGCCTGACTGCTTTGTGCGATCAGGCATTCTCGAATTAATGGTGATGGTCGTCAATCTGGTGTTCGATAACCATCCCTTCACCTACGCTGGCAAGATGGc      <  2:1557097/149‑1 (MQ=255)
                                                                                                                                              tCTTTTTAAAAGATAAACAGATGGGAAATTCTCACGAAGGGTAAAAAAATGCCTGACTGCTTTGTGCGATCAGGCATTCTCGAATTAATGGGGATGGTCGTCAATCTGGTGTTCGATAACCATCCCTTCACCTACGCTGGCAAGATGGc      <  2:254835/149‑1 (MQ=255)
                                                                                                                                                  tttAAAAGATAAACAGATGGGAAATTCTCACGAAGGGTAAAAAAATGCCTGACTGCTTTGTGCGATCAGGCATTCTCGAATTAATGGTGATGGTCGTCAATCTGGTGTTCGATAACCATCCCTTCACCTACGCTGGCAAGATGGCGAAc  <  2:318193/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
TCTGCCGCGATGGGAATGCTGCAAATGCTGATCTTTACCGTTGGTATTGAAATCAGCAAACATGCCTGGCTGAACGGGGGCAACGGACTGTTTAATCTCTTCAACCTTGTCAACGGAATTTTGTGGCTGTCGCTGATGGTTATCTTTTTAAAAGATAAACAGATGGGAAATTCTCACGAAGGGTAAAAAAATGCCTGACTGCTTTGTGCGATCAGGCATTCTCGAATTAATGGTGATGGTCGTCAATCTGGTGTTCGATAACCATCCCTTCACCTACGCTGGCAAGATGGCGAAC  >  NC_000913/884720‑885014

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