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.

TCTGCCGCGATGGGAATGCTGCAAATGCTGATCTTTACCGTTGGTATTGAAATCAGCAAACATGCCTGGCTGAACGGGGGCAACGGACTGTTTAATCTCTTCAACCTTGTCAACGGAATTTTGTGGCTGTCGCTGATGGTTATCTTTTTAAAAGATAAACAGATGGGAAATTCTCACGAAGGGT‑AAAAAAATGCCTGACTGCTTTGTGCGATCAGGCATTCTCGAATTAATGGTGATGGTCGTCAATCTGGTGTTCGATAACCATCCCTTCACCTACGCTGGCAAGATGGCGAAC  >  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)
                                        ttGGTATTGAAATCAGCAAACATGCCTGGCTGAACGGGGGCAACGGACTGTTTAATCTCTTCAACCTTGTCAACGGAATTTTGTGGCTGTCGCTGATGGTTATCTTTTTAAAAGATAAACAGATGGGAAATTCTCACGAAGGGT‑aaaaa                                                                                                            >  1:659337/1‑149 (MQ=255)
                                        ttGGTATTGAAATCAGCAAACATACCTGGCTGAACGGGGGCAACGGACTGTTTAATCTCTTCAACCTTGTCAACGGAATTTTGTGGATGTCGCTGATGGTTATCTTTTTAAAAGATAAACAGATGGGAAATTCTCACGAAGGTT‑aaaaa                                                                                                            >  2:1036681/1‑149 (MQ=255)
                                             tttGAAATCAGCAAACATGCCTGGCTGAACGGGGGCAACGGACTGTTTAATCTCTTCAACCTTGTCAACGGAATTTTGTGGCTGTCGCTGATGGTTATCTTTTTAAAAGATAAACAGATGGGAAATTCTCACGAAGGGT‑AAAAAAATGc                                                                                                       <  2:1337615/148‑1 (MQ=255)
                                              ttGAAATCAGCAAACATGCCTGGCTGAACGGGGGCAACGGACTGTTTATTCTCTTCAACCTTGTCAACGGAATTTTGTGGCTGTCGCTGATGGTTATCTTTTTAAAAGATAAACAGATGGGAAATTCTCACGAAGGGT‑AAAAAAATGcc                                                                                                      >  1:102960/1‑149 (MQ=255)
                                              ttGAAATCAGCAAACATGCCTGGCTGAACGGGGGCAACGGACTGTTTAATCTCTTCAACCTTGTCAACGGAATTTTGTGGCTGTCGCTGATGGTTATCTTTTTAAAAGATAAACAGATGGGAAATTCTCACGAAGGGT‑AAAAAAATGcc                                                                                                      >  1:254835/1‑149 (MQ=255)
                                              ttGAAATCAGCAAACATGCCTGGCTGAACGGGGGCAACGGACTGTTTAATCTCTTCAACCTTGTCAACGGAATTTTGTGGCTGTCGCTGATGGTTATCTTTTTAAAAGATAAACAGATGGGAAATTCTCACGAAGGGT‑AAAAAAATGcc                                                                                                      >  2:641987/1‑149 (MQ=255)
                                              ttGAAATCAGCAAACATGCCTGGCTGAACGGGGGCAACGGACTGTTTAATCTCTTCAACCTTGTCAACGGAATTTTGTGGCTGTCGCTGATGGTTATCTTTTTAAAAGATAAACAGATGGGAAATTCTCACGAAGGGT‑AAAAAAATGcc                                                                                                      <  1:1249071/149‑1 (MQ=255)
                                                 aaaTCAGCAAACATGCCTGGCTGAACGGGGGCAACGGACTGTTTAATCTCTTCAACCTTGTCAACGGAATTTTGTGGCTGTCGCTGATGGTTATCTTTTTAAAAGATAAACAGATGGGAAATTCTCACGATGGGT‑AAAAAAATCCATGa                                                                                                   >  1:746247/1‑149 (MQ=255)
                                                 aaaTCAGCAAACATGCATGGCTGTACGGGGGCAACGGACTGTTTTATCTCTTCAACCTTGTCAACGGAATTTTGTGGCTGTCGCTGATGGTTATCTTTTTAAAAGATAAACAGATGGGAAATTCTCACGAAGGGT‑AAAAAAATGCCTGa                                                                                                   <  1:1020871/149‑1 (MQ=255)
                                                      aGCAAACATGCCTGGCTGAACGGGTGCAACGGACTGTTTAATCTCTTCAACCTTGTCAACGGAATTTTGTGGCTGTCGCTGATGGTTATATTTTTAAAAGATAAACAGATGGTAAATTCTCACGAAGGGT‑AAAAAACTGCCTGACTGCt                                                                                              >  2:73819/1‑149 (MQ=255)
                                                      aGCAAACATGCCTGGCTGAACGGGGGCAACGGACTGTTTAATCTCTTCAACCTTGTCAACGGAATTTTGTGGCTGTCGCTGATGGTTATCTTTTTAAAAGATAAACAGATGGGAAATTCTCGCGAAGGGT‑AAAAAAATGCCTGACTGCt                                                                                              <  1:565048/149‑1 (MQ=255)
                                                       gCAAACATGCCTGGCTGAACGGGGGCAACGGACTGTTTAATCTCTTCAGCCTTGTCAACGGAATTTTGTGGCTGTCGCTGATGGTTATCTTTTTAAAAGATAAACAGATGGGAAATTCTCACGAAGGGT‑AAAAAAATGCCTGACTGCtt                                                                                             >  1:1259163/1‑149 (MQ=255)
                                                       gCAAACATGCCTGGCTGAACGGGGGCAACGGACTGTTTAATCTCTTCAACCTTGTCAACGGAATTTTGTGGCTGTCGCTGATGGTTATCTTTTTAAAAGATAAACAGATGGGAAATTCTCACGAAGGGT‑AAAAAAATGCCTGACTGCtt                                                                                             >  2:223030/1‑149 (MQ=255)
                                                       gCAAACATGCCTGGCTGAACGGGGGCAACGGACTGTTTAATCTCTTCAACCTTGTCAACGGAATTTTGTGGCTGTCGCTGATGGTTATCTTTTTAAAAGATAAACAGATGGGAAATTCTCACGAAGGGT‑AAAAAAATGCCTGACTGCtt                                                                                             >  2:793930/1‑149 (MQ=255)
                                                        cAAACATGCCTGGCTGAACGGGGGCAACGGACTGTTTAATCTCTTCACCCTTGTCAACGGAATTTTGTGGCTGTCGCTGATGGTTATCTTTTTAAAAGATAAACAGATGGGAAATTCTCACGAAGGGTAAAAAAAATGCCTGACTGCtt                                                                                             >  2:293109/1‑149 (MQ=255)
                                                        cAAACATGCCTGGCTGAACGGGGGCAACGGACTGTTTAATCTCTTCAACCTTGTCAACGGAATTTTGTGGCTGTCGCTGATGGTTATCTTTTTAAAAGATAAACAGATGGGAAATTCTCACGAAGGGT‑AAAAAAATGCCTGACTGCttt                                                                                            >  1:1413174/1‑149 (MQ=255)
                                                        cAAACATGCCTGGCTGAACGGGGGCAACGGACTGTTTAATCTCTTCAACCTTGTCAACGGAATTTTGTGGCTGTCGCTGATGGTTATCTTTTTAAAAGATAAACAGATGGGAAATTCTCACGAAGGGT‑AAAAAAATGCCTGACTGCttt                                                                                            >  2:815327/1‑149 (MQ=255)
                                                           aCATGCCTGGCTGAACGGGGGCAACGGACTGTTTAATCTCTTCAACCTTGTCAACGGAATTTTGTGGCTGTCGCTGATGGTTATCTTTTTAAAAGATAAACAGATGGGAAATTCTCACGAAGGGT‑AAAAAAATGCCTGACTGCTTtgtg                                                                                         >  2:1383377/1‑149 (MQ=255)
                                                                 cTGGCTGAACGGGGGCAACGGACTGTTTAATCTCTTCAACCTTGTCAACGGAATTTTGTGGCTGTCGCTGATGGTTATCTTTTTAAAAGATAAACAGATGGGAAATTCTCACGAAGGGT‑AAAAAAATGCCTGACTGCTTTGTGCGATCa                                                                                   <  1:1279090/149‑1 (MQ=255)
                                                                 cTGGCTGAACGGGGGCAACGGACTGTTTAATCTCTTCAACCTTGTCAACGGAATTTTGTGGCTGTCGCTGATGGTTATCTTTTTAAAAGATAAACAGATGGGAAATTCTCACGAAGGGT‑AAAAAAATGCCTGACTGCTTTGTGCGATCa                                                                                   >  1:330125/1‑149 (MQ=255)
                                                                  tGGCTGAACGGGGGCAACGGACTGTTTAATCTCTTCAACCTTGTCAACGGAATTTTGTGGCTGTCGCTGATGGTTATCTTTTTAAAAGATAAACAGATGGGAAATTCTCACGAAGGGT‑AAAAAAATGCCTGACTGCTTTGTGCGATCAg                                                                                  >  2:247324/1‑149 (MQ=255)
                                                                  tGGCTGAACGGGGGCAACGGACTGTTTAATCTCTTCAACCTTGTCAACGGAATTTTGTGGCTGTCGCTGATGGTTATCTTTTTAAAAGATAAACAGATGGGAAATTCTCACGAAGGGT‑AAAAAAATGCCTGACTGCTTTGTGCGATCAg                                                                                  <  1:491752/149‑1 (MQ=255)
                                                                   ggCTGAACGGGGGCAACGGACTGTTTAATCTCTTCAACCTTGTCAACGGAATTTTGTGGCTGTCGCTGATGGTTATCTTTTTAAAAGATAAACAGATGGGAAATTCTCACGAAGGGT‑AAAAAAATGCCTGACTGCTTTGTGCGATCAgg                                                                                 <  2:886332/149‑1 (MQ=255)
                                                                    gCTGAACGGGGGCAACGGACTGTTTAATCTCTTCAACCTTGTCAACGGAATTTTGTGGCTGTCGCTGATGGTTATCTTTTTAAAAGATAAACAGATGGGAAATTCTCACGAAGGGT‑AAAAAAATGCCTGACTGCTTTGTGCGATCAGGc                                                                                <  1:1217009/149‑1 (MQ=255)
                                                                    gCTGAACGGGGGCAACGGACTGTATAATCTCTTCAACCTTGTCAACGGAATTCTGTGGCTGTCGCTGATGGTTATCTTTTTAAAAGATAAACAGATGGGAAATTCTCACGAAGGGT‑AAAAAAATGCCTGACTGCTTTGTGCGATCAGGc                                                                                <  1:715771/149‑1 (MQ=255)
                                                                      tGAACGGGGGCAACGGACTGTTTAATCTCTTCAACCTTGTCAACGGAATTTTGTGGCGGTCGCTGATGGTTTTCTTTTTAAAATATAAACAGATGGGAAATTATCACGAAGGGT‑AAAAAAATGCCTGCCTGCTTTGTGCGATCAGGCAt                                                                              >  2:54093/1‑149 (MQ=255)
                                                                          cGGGGGCAACGGACTGTTTAATCTCTTCAACCTTGTCAACGGAATTTTGTGGCTGTCGCTGATGGTTATCTTTTTAAAAGATAAACAGATGGGAAATTCTCACGAAGGGT‑AAAAAAATGCCTGACTGCTTTGTGCGATCAGGCATtctc                                                                          <  1:1193811/149‑1 (MQ=255)
                                                                            ggggCAACGGACTGTTTAATCTCTTCAACCTTGTCAACGGAATTTTGTGGCTGTCGCTGATGGTTATCTTTTTAAAAGATAAACAGATGGGAAATTCTCACGAAGGGT‑AAAAAAATGCCTGACTGCTTTGTGCGATCAGGCATTCTCGa                                                                        <  1:483406/149‑1 (MQ=255)
                                                                             gttcAACTGACTGTTTAATCTATTATACATTTTCAACTGAATTTTGTTGCTGTCTCTGATGGTTATCTTTTTAAAAGATAATCAGATGGGAAATTCTCACTAAGGGT‑AAAAAAATGCCTGACTGCTTTGTGCGATCAGGCATTCTCGaa                                                                       <  1:1229619/146‑1 (MQ=255)
                                                                             gggCAACGGACTGTTTAATCTCTTCAACCTTGTCAACGGAATTTTGTGGCTGTCGCTGATGGTTATCTTTTTAAAAGATAAACAGATGGGAAATTCTCACGAAGGGT‑AAAAAAATGCCTGACTGCTTTGTGCGATCAGGCATTCTCGaa                                                                       <  2:706321/149‑1 (MQ=255)
                                                                              ggCAACGGACTGTTTAATCTCTTCCACCTTGTCAACGGAATTTTGTGGCTGTCGCTGATGGTTATCTTTTTAAAAGATAAACAGATGGGAAATTCTCACGAAGGGT‑AAAAAAATGCCTGACTGCTTTGTGCGATCAGGCATTCTCGAAt                                                                      >  2:1441883/1‑149 (MQ=255)
                                                                                 aaCGGACTGTTTAATCTCTTCAACCTTGTCAACGGAATTTTGTGGCTGTCGCTGATGGTTATCTTTTTAAAAGATAAACAGATGGGAAATTCTCACGAAGGGT‑AAAAAAATGCCTGACTGCTTTGTGCGATCAGGCATTCTCGAATTaa                                                                   <  1:793930/149‑1 (MQ=255)
                                                                                  aCGGACTGTGTAATCTCTTCAACCTTGTCAACGGAATTTTGTGGCTGTCGCTGATGGTTATCTTTTTAAAAGATAAACAGATGGGAAATTCTCACGAAGGGT‑AAAAAAATGCCTGACTGCTTTGTGCGATCAGGCATTCTCGAATTAAt                                                                  <  2:529822/149‑1 (MQ=255)
                                                                                   cGGACTGTTTAATCTCTTCAACCTTGTCAACGGAATTTTGTGGCTGTCGCTGATGGTTATCTTTTTAAAAGATAAACAGATGGGAAATTCTCACGAAGGGT‑AAAAAAATGCCTGACTGCTTTGTGCGATCAGGCATTCTCGAATTAATg                                                                 >  1:835687/1‑149 (MQ=255)
                                                                                   cGGACTGTTTAATCTCTTCAACCTTGTCAACGGAATTTTGTGGCTGTCGCTGATGGTTATCTTTTTAAAAGATAAACAGATGGGAAATTCTCACGAAGGGT‑AAAAAAATGCCTGACTGCTTTGTGCGATCAGGCATTCTCGAATTAATg                                                                 <  1:883775/149‑1 (MQ=255)
                                                                                      aCTGTTTAATCTCTTCAACCTTGTCAACGGAATTTTGTGGCTGTCGCTGATGGTTATCTTTTTAAAAGATAAACAGATGGGAAATTCTCACGAAGGGT‑AAAAAAATGCCTGACTGCTTTGTGCGATCAGGCATTCTCGAATTAATGGTg                                                              >  2:49791/1‑149 (MQ=255)
                                                                                      aCTGTTTAATATCTTCAACCTTTTCAACGGAATTTTGTGGCTGTCGCTGATGGTTATCTTTTTAAAAGATAAACAGATGGGAAATTCTCACTAAGGGT‑AAAAAAATGCCTGACTGCTTTGTGCGATCAGGCATTCTCGAATTAATGGTg                                                              <  1:967972/149‑1 (MQ=255)
                                                                                       cTGTTTAATCTCTTCAACCTTGTCAACGGAATTTTGTGGCTGTCGCTGATGGTTATCTTTTTAAAAGATAAACAGATGGGAAATTCTCACGAAGGGT‑AAAAAAATGCCTGACTGCTTTGTGCGATCAGGCATTCTCGAATTAATGGTGa                                                             <  2:1159548/149‑1 (MQ=255)
                                                                                          tttAATCTCTTCAACCTTGTCAACGGAATTTTGTGGCTGTCGCTGATGGTTATCTTTTTAAAAGATAAACAGATGGGAAATTCTCACGAAGGGT‑AAAAAAATGCCTGACTGCTTTGTGCGATCAGGCATTCTCGAATTAATGGTGATgg                                                          <  1:1396739/149‑1 (MQ=255)
                                                                                          tttAATCTCTTCAACCTTGTCAACGGAATTTTGTGGCTGTCGCTGATGGTTATCTTTTTAAAAGATAAACAGATGGGAAATTCTCACGAAGGGT‑AAAAAAATGCCTGACTGCTTTGTGCGATCAGGCATTCTCGAATTAATGGTGATgg                                                          >  2:848185/1‑149 (MQ=255)
                                                                                            tAATCTCTTCAACCTTGTCAACGGAATTTTGTGGCTGTCGCTGATGGTTATCTTTTTAAAAGATAAACAGATGGGAAATTCTCACGAAGGGT‑AAAAAAATGCCTGACTGCTTTGTGCGATCAGGCATTCTCGAATTAATGGTGATGgtc                                                        <  1:460130/149‑1 (MQ=255)
                                                                                               tctcTTCAACCTTGTCAACGGAATTTTGTGGCTGTCGCTGATGGTTATCTTTTTAAAAGATAAACAGATGGGAAATTCTCACGAAGGGT‑AAAAAAATGCCTGACTGCTTTGTGCGATCAGGCATTCTCGAATTAATGGTGATGgtcgtc                                                     >  2:598189/1‑149 (MQ=255)
                                                                                               tctcTTCAACCTTGTCAACGGAATTTTGTGGCTGTCGCTGATGGTTATCTTTTTAAAAGATAAACAGATGGGAAATTCTCACGAAGGGT‑AAAAAAATGCCTGACTGCTTTGTGCGATCAGGCATTCTCGAATTAATGGTGATGgtcgtc                                                     >  2:773201/1‑149 (MQ=255)
                                                                                               tctcTTCAACCTTGTCAACGGAATTTTGTGGCTGTCGCTGATGGTTATCTTTTTAAAAGATAAACAGATGGGAAATTCTCACGAAGGGT‑AAAAAAATGCCTGACTGCTTTGTGCGATCAGGCATTCTCGAATTAATGGTGATGgtcgtc                                                     >  2:1356631/1‑149 (MQ=255)
                                                                                                 tctTCAACCTTGTCAACGGAATTTTGTGGCTGTCGCTGATGGTTATCTTTTTAAAAGATAAACAGATGGGAAATTCTCACGAAGGGT‑AAAAAAATGCCTGACTGCTTTGTGCGATCAGGCATTCTCGAATTAATGGTGATGGTCGTCaa                                                   <  1:1351103/149‑1 (MQ=255)
                                                                                                   ttCAACCTTGTCAACGGAATTTTGTGGCTGTCGCTGATGGTTATCTTTTTAAAAGATAAACAGATGGGAAATTCTCACGAAGGGT‑AAAAAAATGCCTGACTGCTTTGTGCGATCAGGCATTCTCGAATTAATGGTGATGTTCGTCAATc                                                 >  2:406245/1‑149 (MQ=255)
                                                                                                   ttCAACCTTGTCAACGGAATTTTGTGGCTGTCGCTGATGGTTATCTTTTTAAAAGATAAACAGATGGGAAATTCTCACGAAGGGT‑AAAAAAATGCCTGACTGCTTTGTGCGATCAGGCATTCTCGAATTAATGGTGATGGTCGTCAATc                                                 <  1:820948/149‑1 (MQ=255)
                                                                                                     cAACCTTGTCAACGGAATTTTGTGGCTGTCGCTGATGGTTATCTTTTTAAAAGATAAACAGATGGGAAATTCTCACGAAGGGT‑AAAAAAATGCCTGACTGCTTTGTGCGATCAGGCATTCTCGAATTAATGGTGATGGTCGTCAATCTg                                               <  1:752103/149‑1 (MQ=255)
                                                                                                     cAACCTTGTCAACGGAATTTTGTGGCTGTCGCTGATGGTTATCTTTTTAAAAGATAAACAGATGGGAAATTCTCACGAAGGGT‑AAAAAAATGCCTGACTGCTTTGTGCGATCAGGCATTCTCGAATTAATGGTGATGGTCGTCAATCTg                                               <  2:316282/149‑1 (MQ=255)
                                                                                                      aaCCTTGTCAACGGAATTTTGTGGCTGTCGCTGATGGTTATCTTTTTAAAAGATAAACAGATGGGAAATTCTCACGAAGGGT‑AAAAAAATGCCTGACTGCTTTGTGCGATCAGGCATTCTCGAATTAATGGTGATGGTCGTCAATCTgg                                              <  1:1333856/149‑1 (MQ=255)
                                                                                                      aaCCTTGTCAACGGAATTTTGTGGCTGTCGCTGATGGTTATCTTTTTAAAAGATAAACAGATGGGAAATTCTCACGAAGGGT‑AAAAAAATGCCTGACTGCTTTGTGCGATCAGGCATTCTCGAATTAATGGTGATCGTCGTCAATCTgg                                              <  2:1037170/149‑1 (MQ=255)
                                                                                                       aCCTTGTCAACGGAATTTTGTGGCTGTCGCTGATGGTTATCTTTTTAAAAGATAAACAGATGGGAAATTCTCACGAAGGGT‑AAAAAAATGCCTGACTGCTTTGTGCGATCAGGCATTCTCGAATTAATGGTGATGGTCGTCAATCTGgt                                             >  1:817714/1‑149 (MQ=255)
                                                                                                        ccTTGTCAACGGAATTTTGTGGCTGTCGCTGATGGTTATCTTTTTAAAAGATAAACAGATGGGAAATTCTCACGAAGGGT‑AAAAAAATGCCTGACTGCTTTGTGCGATCAGGCATTCTCGAATTAATGGTGATGGTCGTCAATCTGgtg                                            >  2:1480946/1‑149 (MQ=255)
                                                                                                         cTTGTCAACGGAATTTTGTGGCTGTCGCTGATGGTTATCTTTTTAAAAGATAAACAGATGGGAAATTCTCACGAAGGGT‑AAAAAAATGCCTGACTGCTTTGTGCGATCAGGCATTCTCGAATTAATGGTGATGGTCGTCAATCTGgtgt                                           <  1:1354318/149‑1 (MQ=255)
                                                                                                          ttGTCAACGGAATTTTGTGGCTGTCGCTGATGGTTATCTTTTTAAAAGATAAACAGATGGGAAATTCTCACGAAGGGT‑AAAAAAATGCCTGACTGCTTTGTGCGATCAGGCATTCTCGAATTAATGGTGATGGTCGTCAATCTGGTGtt                                          <  1:815327/149‑1 (MQ=255)
                                                                                                           tGTCAACGGAATTTTGTGGCTGTCGCTGATGGTTATCTTTTTAAAAGATAAACAGATGGGAAATTCTCACGAAGGGT‑AAAAAAATGCCTGACTGCTTTGGGCGATCAGGCATTCTCGAATAAATGGTGATGGTCGTCAATCTGGTGTTc                                         <  2:276853/149‑1 (MQ=255)
                                                                                                              cAACGGAATTTTGTGGCTGTCGCTGATGGTTATCTTTTTAAAAGATAAACAGATGGGAAATTCTCACGAAGGGT‑AAAAAAATGCCTGACTGCTTTGTGCGATCAGGCATTCTCGAATTAATGGTGATGGTCGTCAATCTGGTGTTCGAt                                      <  2:1428098/149‑1 (MQ=255)
                                                                                                                aCGGAATTTTGTGGCTGTCGCTGATGGTTATCTTTTTAAAAGATAAACAGATGGGAAATTCTCACGAAGGGT‑AAAAAAATGCCTGACTGCTTTGTGCGATCAGGCATTCTCGAATTAATGGTGATGGTCGTCAATCTGGTGTTCGATaa                                    >  1:1557097/1‑149 (MQ=255)
                                                                                                                 cGGAATTTTGTGGCTGTCGCTGATGGTTATCTTTTTAAAAGATAAACAGATGGGAAATTCTCACGAAGGGT‑AAAAAAATGCCTGACTGCTTTGTGCGATCAGGCATTCTCGAATTAATGGTGATGGTCGTCAATCTGGTGTTCGATAAc                                   >  1:519498/1‑149 (MQ=255)
                                                                                                                    aaTTTTGTGGCTGTCGCTGATGGTTATCTTTTTAAAAGATAAACAGATGGGAAATTCTCACGAAGGGT‑AAAAAAATGCCTGACTGCTTTGTGCGATCAGGCATTCTCGAATTAATGGTGATGGTCGTCAATCTGGTGTTCGATAACCAt                                <  2:538861/149‑1 (MQ=255)
                                                                                                                     aTTTTGTGGCTGTCGCTGATGGTTATCTTTTTAAAAGATAAACAGATGGGAAATTCTCACGAAGGGTAAAAAAAATGCCTGACTGCTTTGTGCGATCAGGCATTCTCGAATTAATGGTGATGGTCGTCAATCTGGTGTTCGATAACCAt                                <  2:1543481/149‑1 (MQ=255)
                                                                                                                      ttttGTGGCTGTCGCTGATGGTTATCTTTTTAAAAGATAAACAGATGGGAAATTCTCACGAAGGGT‑AAAAAAATGCCTGACTGCTTTGTGCGATCAGGCATTCTCGAATTAATGGTGATGGTCGTCAATCTGGTGTTCGATAACCTTcc                              <  1:690769/149‑1 (MQ=255)
                                                                                                                          gtgGCTGTCGCTGATGGTTATCTTTTTAAAAGATAAACAGATGGGAAATTCTCACGAAGGGT‑AAAAATATGCCTGACTGCTTTGTGCGATCAGGCATTCTCGAATTAATGGTGATGGTCGTCAATCTGGTGTTCGATAACCATCCCTTc                          <  2:102960/149‑1 (MQ=255)
                                                                                                                            ggatGTCGCTGATGGTTATCTTTTTAAAAGATAAACAGATGGGAAATTCTCACGAATGTT‑AAAAAAATGCCTGACTGCTTTGTGCGATCAGGCATTCTCGAATTAATGGTGATGGTCGTCAATCTGGTGTTCGATAACCATCCCTTCAc                        <  1:1443948/146‑1 (MQ=255)
                                                                                                                            ggCTGTCGCTGATGTTTATCTTTTTAAAAGATAAACAGATGGGAAATTCTCACGAAGGGT‑AAAAAAATGCCTGACTGCTTTGTGCGATCAGGCATTCTCGAATTAATGGGGATGGTCGTCAATCTGGTGTTCGATAACCATCCCTTCAc                        <  2:693884/149‑1 (MQ=255)
                                                                                                                            ggCTGTCGCTGATGGTTATCTTTTTAAAAGATAAACAGATGGGAAATTCTCACGAAGGGT‑AAAAAAATGCCTGACTGCTTTGTGCGATCAGGCATTCTCGAATTAATGGTGATGGTCGTCAATCTGGTGTTCGATAACCATCCCTTCAc                        <  1:1349539/149‑1 (MQ=255)
                                                                                                                                    catctgtttATCTTTTAAAAAGATAAACAGATGGGAAATTCTCACGAAGGGT‑AAAAAAATGCCTGACTGCTTTGTGCGATCAGGCATTCTCGAATTAATGGTGATGGTCGTCAATCTGGTGTTCGATAACCATCCCTTCACCTACGCTg                >  2:593925/8‑149 (MQ=255)
                                                                                                                                    catctgtttATCTTTTAAAAAGATAAACAGATGGGAAATTCTCACGAAGGGT‑AAAAAAATGCCTGACTGCTTTGTGCGATCAGGCATTCTCGAATTAATGGTGATGGTCGTCAATCTGGTGTTCGATAACCATCCCTTCACCTACGCTg                >  1:942432/8‑149 (MQ=255)
                                                                                                                                     tGATGGTTATCTTTTTAAAAGATAAACAGATGGGAAATTCTCACGAAGGGT‑AAAAAAATGCCTGACTGCTTTGTGCGATCAGGCATTCTCGAATTAATGGTGATGGTCGTCAATCTGGTGTTCGATAACCATCCCTTCACCTACGCTgg               >  2:1366040/1‑149 (MQ=255)
                                                                                                                                          tttATCTTTTAAAAAGATAAACAGATGGGAAATTCTCACGAAGGGT‑AAAAAAATGCCTGACTGCTTTGTGCGATCAGGCATTCTCGAATTAATGGTGATGGTCGTCAATCTGGTGTTCGATAACCATCCCTTCACCTACGCTGGCAAGa          >  2:1593284/2‑149 (MQ=255)
                                                                                                                                             aTCTTTTTAAAAGATAAACAGATGGGAAATTCTCACGAAGGGT‑AAAAAAATGCCTGACTGCTTTGTGCGATCAGGCATTCTCGAATTAATGGTGATGGTCGTCAATCTGGTGTTCGATAACCATCCCTTCACCTACGCTGGCAAGATgg       >  1:490540/1‑149 (MQ=255)
                                                                                                                                             aTCTTTTAAAAAGATAAACAGATGGGAAATTCTCACGAAGGGT‑AAAAAAATGCCTGACTGCTTTGTGCTATCAGGCATTCTCGAATTAATGGTGATGGTCGTCAATCTGGTGTTCGATAACCATCCCTTCACCTACGCTGGCAAGATgg       >  1:573665/1‑149 (MQ=255)
                                                                                                                                             aTCTTTTAAAAAGATAAACAGATGGGAAATTCTCACGAAGGGT‑AAAAAAATGCCTGACTGCTTTGTGCAATCAGGCATTCTCGAATTAATGGTGATGGTCGTCAATCTGGTGTTCGATAACCATCCCTTCACCTACGCTGGCAAGATgg       >  1:568074/1‑149 (MQ=255)
                                                                                                                                              tCTTTTTAAAAGATAAACAGATGGGAAATTCTCACGAAGGGT‑AAAAAAATGCCTGACTGCTTTGTGCGATCAGGCATTCTCGAATTAATGGTGATGGTCGTCAATCTGGTGTTCGATAACCATCCCTTCACCTACGCTGGCAAGATGGc      <  2:330125/149‑1 (MQ=255)
                                                                                                                                              tCTTTTTAAAAGATAAACAGATGGGAAATTCTCACGAAGGGT‑AAAAAAATGCCTGACTGCTTTGTGCGATCAGGCATTCTCGAATTAATGGTGATGGTCGTCAATCTGGTGTTCGATAACCATCCCTTCACCTACGCTGGCAAGATGGc      <  2:1557097/149‑1 (MQ=255)
                                                                                                                                              tCTTTTTAAAAGATAAACAGATGGGAAATTCTCACGAAGGGT‑AAAAAAATGCCTGACTGCTTTGTGCGATCAGGCATTCTCGAATTAATGGGGATGGTCGTCAATCTGGTGTTCGATAACCATCCCTTCACCTACGCTGGCAAGATGGc      <  2:254835/149‑1 (MQ=255)
                                                                                                                                                  tttAAAAGATAAACAGATGGGAAATTCTCACGAAGGGT‑AAAAAAATGCCTGACTGCTTTGTGCGATCAGGCATTCTCGAATTAATGGTGATGGTCGTCAATCTGGTGTTCGATAACCATCCCTTCACCTACGCTGGCAAGATGGCGAAc  <  2:318193/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
TCTGCCGCGATGGGAATGCTGCAAATGCTGATCTTTACCGTTGGTATTGAAATCAGCAAACATGCCTGGCTGAACGGGGGCAACGGACTGTTTAATCTCTTCAACCTTGTCAACGGAATTTTGTGGCTGTCGCTGATGGTTATCTTTTTAAAAGATAAACAGATGGGAAATTCTCACGAAGGGT‑AAAAAAATGCCTGACTGCTTTGTGCGATCAGGCATTCTCGAATTAATGGTGATGGTCGTCAATCTGGTGTTCGATAACCATCCCTTCACCTACGCTGGCAAGATGGCGAAC  >  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: ‑