Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,568,5120CA25.0% 136.2 / 19.3 48R3L (CGT→CTT) mdtMmultidrug efflux pump/bile salt:H(+) antiporter/Na(+):H(+) antiporter/K(+):H(+) antiporter
Reads supporting (aligned to +/- strand):  ref base C (4/32);  new base A (12/0);  total (16/32)
Fisher's exact test for biased strand distribution p-value = 2.61e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

AGACTGACATCGGCATTAAAATCACGTACCACATTAATGATCCCAGGCTGGATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAACTCCTTGTCCGGAGAACCCCAACGGGGAAAACGAAATTCAATACCGTCGGTAATGACGAAAGGCAGGTGAGAAGGGCAACGAGACTGCAAATGTACTGAATTTCATGTAATCAACCACA  >  NC_000913/4568374‑4568647
                                                                                                                                          |                                                                                                                                         
aGACTGACATCGGCATTAAAATCACGTACCACATTAATGATCCCAGGCTGGATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACg                                                                                                                                          >  1:765627/1‑140 (MQ=255)
   cTGACATCGGCATTAAAATCACGTACCACATTAATGATCCCAGGCTGGATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTgg                                                                                                                                       <  1:1158875/140‑1 (MQ=255)
      aCATCGGCATTAAAATCACGTACCACATTAATGATCCCAGGCTGGATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCAt                                                                                                                                    >  1:374453/1‑140 (MQ=255)
      aCATCGGCATTAAAATCACGTACCACATTAATGATCCCAGGCTGGATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGGAAAAAAAAGAGGCaa                                                                                                                                    >  2:1142894/1‑139 (MQ=255)
        aTCGGCATTAAAATCACGTACCACATTAATGATCCCAGGCTGGATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATaa                                                                                                                                  <  2:194741/140‑1 (MQ=255)
        aTCGGCATTAAAATCACGTACCACATTAATGATCCCAGGCTGGATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATaa                                                                                                                                  <  1:1142894/140‑1 (MQ=255)
           ggCATTAAAATCACGTACCACATTAATGATCCCAGGCTGGATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGGAAAAAAAACCGGCATaacaa                                                                                                                               >  1:1201775/1‑140 (MQ=255)
            gCATTAAAATCACGTACCACATTAATGATCCCAGGCTGGATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATaacaaa                                                                                                                              <  2:500825/140‑1 (MQ=255)
              aTTAAAATCACGTACCACATTAATGATCCCAGGCTGGATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAAAGTGGCATaacaaaca                                                                                                                            >  1:249297/1‑140 (MQ=255)
                tAAAATCACGTACCACATTAATGATCCCAGGCTGGATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAAAGTGGGATAACAaacaac                                                                                                                          >  1:791326/1‑140 (MQ=255)
                  aaaTCACGTACCACATTAATGATCCCAGGCTGGATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAACTc                                                                                                                        <  1:527819/140‑1 (MQ=255)
                   aaTCACGTACCACATTAATGATCCCAGGCTGGATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAACTcc                                                                                                                       <  2:374453/140‑1 (MQ=255)
                       aCGTACCACATTAATGATCCCAGGCTGGATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAACTCCTTGt                                                                                                                  <  2:1031143/140‑1 (MQ=255)
                       aCGTACCACATTAATGATCCCAGGCTGGATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGGGAAAAAAAATGGGATAAAAAAAAAAACCtttt                                                                                                                  >  2:1280154/1‑138 (MQ=255)
                          tACCACATTAATGATCCCAGGCTGGATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAACTCCTTGTCCg                                                                                                               <  2:512319/140‑1 (MQ=255)
                               caTTAATGATCCCAGGCTGGATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAACTCCTTGTCCGGAGaa                                                                                                         <  1:738471/140‑1 (MQ=255)
                                 ttAATGATCCCAGGCTGGATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAACTCCTTGTCCGGAGAAcc                                                                                                       >  1:136458/1‑140 (MQ=255)
                                 ttAATGATCCCAGGCTGGATCAGATCCGTCGACAGATACGCAGCAAAGTCACACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGGAAAAAAAACTGGCATAAAAAAAAAAGCCTTTTCCGGAGAAcc                                                                                                       >  1:1120963/1‑140 (MQ=255)
                                       aTCCCAGGCTGGATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAACTCCTTGTCCGGAGAACCCCAACg                                                                                                 >  1:423674/1‑140 (MQ=255)
                                       aTCCCAGGCTGGATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGGAAAAAAAAGTGGCATAACAAACAACTCCTTGTCCGGAGAACCCCAACg                                                                                                 >  1:1240514/1‑140 (MQ=255)
                                       aTCCCAGGCTGGATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGGAAAAAAAAGTGGCATAAAAAAAAACTCCCTTTCCGGGGAACCCCAACg                                                                                                 >  1:716357/1‑140 (MQ=255)
                                       aTCCCAGGCTGGATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGGAAAAAAAAAGGGCATAAAAAAAAAATCCTTGTCCGGAGAACCCCAACg                                                                                                 >  2:744640/1‑140 (MQ=255)
                                         cccAGGCTGGATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAAAGTGGCATAACAAACAACTCCTTGTCCGGAGAACCCCAACggg                                                                                               >  1:1202592/1‑140 (MQ=255)
                                          ccAGGCTGGATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGGAAAAAAAAGTGGCATAAAAAAAAACTCCTTGTCCGGAGAACCCCCACgggg                                                                                              >  2:781525/1‑140 (MQ=255)
                                          ccAGGCTGGATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGGAAAAAAAACTGGCATAACAAAAAAACCCATGTCCGGTGAAACCCAACgggg                                                                                              >  1:1121703/1‑140 (MQ=255)
                                           cAGGCTGGATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAACTCCTTGTCCGGAGAACCCCAACGGGGa                                                                                             <  2:765627/140‑1 (MQ=255)
                                            aGGCTGGATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGGAAAAAAAAATGGCATAACAAAAAAACCCCTGTCCCGGGAACCCCCACcgggga                                                                                            >  2:365975/1‑138 (MQ=255)
                                              gCTGGATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAAAGTGGCATAACAAACAACTCCTTGTCCGGAGAACCCCAACGGGGaaaa                                                                                          >  1:510275/1‑140 (MQ=255)
                                                tGGATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAACTCCTTGTCCGGAGAACCCCAACGGGGAAAACg                                                                                        <  2:1202592/140‑1 (MQ=255)
                                                 ggATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAACTCCTTGTCCGGAGAACCCCAACGGGGAAAACGa                                                                                       <  1:227893/140‑1 (MQ=255)
                                                 ggATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAAACGAGCATGAAAAAAAACACCCTTTCCGGGGAAACCCCAAAGGGAaaaaaa                                                                                        >  2:140711/1‑138 (MQ=255)
                                                  gATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGGAAAAAAACGTGGCATAACAAAAAACACCTTGTCCCGAGAACCCCCACGGGGAAAAAGaa                                                                                      >  2:770465/1‑140 (MQ=255)
                                                  gATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGGAAAAAAAAGTGGGATAAAAAAAAACTCCTTTTCCGGAGAACCCCCACGGGGGAAACGaa                                                                                      >  1:938593/1‑140 (MQ=255)
                                                  gATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGGAAAAAAAAGTGGCATAAAAAAAAAATCCCTTTCCGGGGGAACCCCACGGGGGAAAAGaa                                                                                      >  2:400115/1‑140 (MQ=255)
                                                   aTCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGGAAAAAAAAGTGGCATAAAAAAAAAATCCTTGTCCGGAGAAACCCAACGGGGGAAAAGaaa                                                                                     >  2:320187/1‑140 (MQ=255)
                                                   aTCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGGAAAAAAAAGGGGGATAAAAAAAAACCCCTTTTCCGGGGAAACCCAAAGGGGGAAAAGaaa                                                                                     >  1:1187522/1‑140 (MQ=255)
                                                      aGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAACTCCTTGTCCGGAGAACCCCAACGGGGAAAACGAAATTc                                                                                  <  2:249297/140‑1 (MQ=255)
                                                       gATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGGAAAAAAAAGGGGCATAACAAAAAACTCCTTGTCCGGAGAACCCCAACGGGGGAAACGAAAATCa                                                                                 >  1:1097224/1‑140 (MQ=255)
                                                        aTCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGTCATAACAAACAACACCGTTTTCCGGGAACCCCAACTGGGGAAACGGAAATCaa                                                                                >  2:445987/1‑140 (MQ=255)
                                                        aTCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAAAGGGGCATAAAAAAAAACCCCCTGTCCGGAGAACCCCAACCGGGGAAACAAAAATCaa                                                                                >  2:989824/1‑140 (MQ=255)
                                                         tCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAAAGGGGCATAAAAAACAACTCCTTGTCCGGAGAACCCCAACGGGGAAAACGAAATTCaaa                                                                               >  1:221037/1‑139 (MQ=255)
                                                            gTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGGAAAAAAAAATGGCATAAAAAAAAACTCCTTGTCCGGAGAAACCCCACAGGGGAAAAGAAAATCAATAcc                                                                            >  2:1245762/1‑140 (MQ=255)
                                                            gTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGGAAAAAAAAAGGGCGTAAAAAAAAACCCCTTGTCCGGGGGACCCCCACGGGGGAAAAAAAAATCaaaaac                                                                            >  2:953270/1‑136 (MQ=255)
                                                             tCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAACTCCTTGTCCGGAGAACCCCAACGGGGAAAACGAAATTCAATACCg                                                                           <  1:400115/140‑1 (MQ=255)
                                                             tCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGGAAAAAAAAGTGGCATAAAAAAAAACACCTTGTCCGGAGGACCCCCACGGGGGAAAAGAAAATCAATAccc                                                                           >  2:1089484/1‑139 (MQ=255)
                                                                    aTACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAACTCCTTGTCCTGAGAACCCCAACGGGGAAAACGAAATTCAATACCGTCGGTaa                                                                    <  2:722300/140‑1 (MQ=255)
                                                                    aTACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAACTCCTTGTCCGGAGAACCCCAACGGGGAAAACGAAATTCAATACCGTCGGTaa                                                                    <  2:1010006/140‑1 (MQ=255)
                                                                         cagcaAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAACTCCTTGTCCGGAGAACCCCAACGGGGAAAACGAAATTCAATACCGTCGGTAATGACg                                                               <  1:953270/140‑1 (MQ=255)
                                                                           gcaAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGGAAAAAAAAGTGGCATAACAAACAACTCCTTTTCCGGAGAACCCCAACGGGGAAAACGAAATTCAATACCGTCGGGAATGACGaa                                                             >  1:57589/1‑140 (MQ=255)
                                                                              aaGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAACTCCTTGTCCGGAGAACCCCAACGGGGAAAACGAAATTCAATACCGTCGGTAATGACGAAAgg                                                          <  1:920899/140‑1 (MQ=255)
                                                                              aaGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGGAAAAAAACGGGGCATAACAAAAAACTCCTTTTCCGGAGAACCCCAACGGGGGAAAAGAAATTAAAAACCGTCGGTAATGACGaaaag                                                          >  2:470610/1‑138 (MQ=255)
                                                                                    tACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAAAGAGGCATAAAAAAAAACTCCCTTTTCCGAGAACCCCCACGGGGGAAAAGAAAATCAAAAACGCCGGGAATGACGAAAAGCAgggg                                                    >  2:1283939/1‑138 (MQ=255)
                                                                                    tACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGGAAAAAAAAGTGGGATAAAAAAAAAATCCTTGTCCGGGGAAACCCAACGGGGAAAAAGAAAATCAAAACCCTCGGTAATGACGAAAAGCAgggg                                                    >  2:753977/1‑138 (MQ=255)
                                                                                     cagaaaTCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAACTCCTTGTCCGGAGAACCCCAACGGGGAAAACGAAATTCAATACCGTCGGTAATGACGAAAGGCAGGTga                                                   <  1:1212603/137‑1 (MQ=255)
                                                                                       aaaaTCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAACTCCTTGTCCGGAGAACCCCAACGGGGAAAACGAAATTCAATACCGTCGGTAATGACGAAAGGCAGGTgaga                                                 <  2:713342/140‑1 (MQ=255)
                                                                                           tCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGGAAAAAAAAGTGGCATAACAAAAAACTCCTTGTCCCGAGAACCCCAACGGGGGAAAAGAAAATCAAAACCGTCGGGAATGACGAAAAGCAGGGGAGAAggg                                             >  2:320287/1‑140 (MQ=255)
                                                                                              aCGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAAAATGGCATAACAAACAACTCCTTGTCCGGAGAAACCCAACGGGGGAAAAGAAAATCAATACCGTCGGTAAAGAAGAAAAGCAGGTGGGAAGGGCaa                                          >  2:777531/1‑140 (MQ=255)
                                                                                                 ccATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAACTCCTTGTCCGGAGAACCCCAACGGGGAAAACGAAATTCAATACCGTCGGTAATGACGAAAGGCAGGTGAGAAGGGCAACga                                       <  1:140711/140‑1 (MQ=255)
                                                                                                  cATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAACTCCTTGTCCGGAGAACCCCAACGGGGAAAACGAAATTCAATACCGTCGGTAATGACGAAAGGCAGGTGAGAAGGGCAACgag                                      <  2:1120963/140‑1 (MQ=255)
                                                                                                  cATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAACTCCTTGTCCGGAGAACCCCAACGGGGAAAACGAAATTCAATACCGTCGGTAATGACGAAAGGCAGGTGAGAAGGGCAACgag                                      <  2:136458/140‑1 (MQ=255)
                                                                                                   aTCGGGAAAAACAGCGTGGCGGCATGGCGGGGAAAAAAACGTGGCATAAAAAAAAAATCCCTGTCCGGAGAAACCCAACGGGGGAAAAGAAAATCAATACCCTCGGGAATGACGAAAAGCAGGGGAGAAAGGGaaagccg                                     >  2:1212603/1‑135 (MQ=255)
                                                                                                    tCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAACTCCTTGTCCGGAGAACCCCAACGGGGAAAACGAAATTCAATACCGTCGGTAATGACGAAAGGCAGGTGAGAAGGGCAACGAGAc                                    <  2:510275/140‑1 (MQ=255)
                                                                                                               aGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAACTCCTTGTCCGGAGAACCCCAACGGGGAAAACGAAATTCAATACCGTCGGTAATGACGAAAGGCAGGTGAGAAGGGCAACGAGACTGCAAATGTAc                         <  1:1089484/140‑1 (MQ=255)
                                                                                                                gCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAACTCCTTGTCCGGAGAACCCCAACGGGGAAAACGAAATTCAATGCCGTCGGTAATGACGAAAGGCAGGTGAGAAGGGCAACGAGACTGCAAATGTACt                        <  1:320287/140‑1 (MQ=255)
                                                                                                                  gTGGCGGCATGGCGGGGAAAAAAAAGTGGCATAACAAACAACTCCTTGTCCGGAGAACCCCAACGGGGGAAACGAAATTCAATACCGTCGGGAATGACGAAAAGCAGGGGAGAAGGGCAACGAGACTGCAAATGTACTGa                      >  2:376958/1‑140 (MQ=255)
                                                                                                                          aTGGCGGGTAAAAAAACGTGGCATAACAAACAACTCCTTGTCCGGAGAACCCCAACGGGGAAAACGAAATTCAATACCGTCGGTAATGACGAAAGGCAGGTGAGAAGGGCAACGAGACTGCAAATGTACTGAATTTCATg              <  2:57589/140‑1 (MQ=255)
                                                                                                                          aTGGCGGGTAAAAAAACGTGGCATAACAAACAACTCCTTGTCCGGAGAACCCCAACGGGGAAAACGAAATTCAATACCGTCGGTAATGACGAAAGGCAGGTGAGAAGGGCAACGAGACTGCAAATGTACTGAATTTCATg              <  1:1245762/140‑1 (MQ=255)
                                                                                                                            ggCGGGTAAAAAAACGTGGCATAACAAACAACTCCTTGTCCGGAGAACCCCAACGGGGAAAACGAAATTCAATACCGTCGGTAATGACGAAAGGCAGGTGAGAAGGGCAACGAGACTGCAAATGTACTGAATTTCATGTa            <  1:1280154/140‑1 (MQ=255)
                                                                                                                            ggCGGGTAAAAAAACGTGGCATAACAAACAACTCCTTGTCCGGAGAACCCCAACGGGGAAAACGAAATTCAATACCGTCGGTAATGACGAAAGGCAGGTGAGAAGGGCAACGAGACTGCAAATGTACTGAATTTCATGTa            <  1:13973/140‑1 (MQ=255)
                                                                                                                            ggCGGGTAAAAAAACGTGGCATAACAAACAACTCCTTGTCCGGAGAACCCCAACGGGGAAAACGAAATTCAATACCGTCGGTAATGACGAAAGGCAGGTGAGAAGGGCAACGAGACTGCAAATGTACTGAATTTCATGTa            <  1:744640/140‑1 (MQ=255)
                                                                                                                                ggTAAAAAAACGTGGCATAACAAACAACTCCTTGTCCGGAGAACCCCAACGGGGAAAACGAAATTCAATACCGTCGGTAATGACGAAAGGCAGGTGAGAAGGGCAACGAGACTGCAAATGTACTGAATTTCATGTAATCa        <  1:1199743/140‑1 (MQ=255)
                                                                                                                                      aaaaCGTGGCATAACAAACAACTCCTTGTCCGGAGAACCCCAACGGGGAAAACGAAATTCAATACCGTCGGTAATGACGAAAGGCAGGTGAGAAGGGCAACGAGACTGCAAATGTACTGAATTTCATGTAATCAACcaca  <  2:716357/140‑1 (MQ=255)
                                                                                                                                          |                                                                                                                                         
AGACTGACATCGGCATTAAAATCACGTACCACATTAATGATCCCAGGCTGGATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAACTCCTTGTCCGGAGAACCCCAACGGGGAAAACGAAATTCAATACCGTCGGTAATGACGAAAGGCAGGTGAGAAGGGCAACGAGACTGCAAATGTACTGAATTTCATGTAATCAACCACA  >  NC_000913/4568374‑4568647

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