Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,568,5040TG33.9% 84.9 / 15.7 62T6P (ACC→CCC) mdtMmultidrug efflux pump/bile salt:H(+) antiporter/Na(+):H(+) antiporter/K(+):H(+) antiporter
Reads supporting (aligned to +/- strand):  ref base T (11/28);  new base G (20/0);  total (34/28)
Fisher's exact test for biased strand distribution p-value = 3.14e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

GCCAGACTGACATCGGCATTAAAATCACGTACCACATTAATGATCCCAGGCTGGATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAACTCCTTGTCCGGAGAACCCCAACGGGGAAAACGAAATTCAATACCGTCGGTAATGACGAAAGGCAGGTGAGAAGGGCAACGAGACTGCAAATGTACTGAATTTCATGTAATCA  >  NC_000913/4568371‑4568641
                                                                                                                                      |                                                                                                                                          
gCCAGACTGACATCGGCATTAAAATCACGTACCACATTAATGATCCCAGGCTGGATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTaaaaaa                                                                                                                                      >  1:713757/1‑140 (MQ=255)
 ccAGACTGACATCGGCATTAAAATCACGTACCACATTAATGATCCCAGGCTGGATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTaaaaaaa                                                                                                                                     <  2:852728/140‑1 (MQ=255)
      cTGACATCGGCATTAAAATCACGTACCACATTAATGATCCCAGGCTGGATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTgg                                                                                                                                <  1:549995/140‑1 (MQ=255)
        gACATCGGCATTAAAATCACGTACCACATTAATGATCCCAGGCTGGATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGcc                                                                                                                              >  1:323041/1‑139 (MQ=255)
          cATCGGCATTAAAATCACGTACCACATTAATGATCCCAGGCTGGATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATa                                                                                                                            <  1:110014/140‑1 (MQ=255)
          cATCGGCATTAAAATCACGTACCACATTAATGATCCCAGGCTGGATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATa                                                                                                                            <  1:396396/140‑1 (MQ=255)
           aTCGGCATTAAAATCACGTACCACATTAATGATCCCAGGCTGGATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATaa                                                                                                                           <  2:547150/140‑1 (MQ=255)
               gCATTAAAATCACGTACCACATTAATGATCCCAGGCTGGATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGGAAAAAAAAGTGGCATaacaaa                                                                                                                       >  1:111138/1‑140 (MQ=255)
               gCATTAAAATCACGTACCACATTAATGATCCCAGGCTGGATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGGAAAAAAAAATGGCATaaaaaa                                                                                                                       >  2:706303/1‑140 (MQ=255)
                  ttAAAATCACGTACCACATTAATGATCCCAGGCTGGATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAacaaacaa                                                                                                                    <  1:706303/140‑1 (MQ=255)
                              aCCACATTAATGATCCCAGGCTGGATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAACTCCTTGTCCgg                                                                                                       <  2:898553/140‑1 (MQ=255)
                                 acaTTAATGATCCCAGGCTGGATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGGAAAAAAACAGGGGATAACAAACAAACCCTCTTCCGgagg                                                                                                    >  2:902819/1‑139 (MQ=255)
                                  caTTAATGATCCCAGGCTGGATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGCAAAAAAACGTGGCATAACAAACAACCCCATTTCCGGAGaa                                                                                                   >  1:832775/1‑140 (MQ=255)
                                      aaTGATCCCAGGCTGGATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAAAGTGGGATAAAAAAAAACTCCTTTTCCGGGGAAaccc                                                                                               >  1:872781/1‑140 (MQ=255)
                                        tGATCCCAGGCTGGATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGGAAAAAAAAGTGGCATAAAAAAAAAATCCTTGTCCGGGGAACCCCaa                                                                                             >  2:734918/1‑140 (MQ=255)
                                          aTCCCAGGCTGGATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAAAAGGGCATAACAAAAAACTCCTTTTCCGGAGAACCCCAACg                                                                                           >  2:162788/1‑140 (MQ=255)
                                          aTCCCAGGCTGGATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGGGAAAAAAACCGGCATAACAAAAAAATCCTTTTCCGGAGAACCCCAACg                                                                                           >  2:610421/1‑140 (MQ=255)
                                          aTCCCAGGCTGGATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGGAAAAAAACGGGGCATAACAAACAAATCCTTGTCCCGAGGAcccccaaa                                                                                           >  2:176160/1‑136 (MQ=255)
                                            cccAGGCTGGATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGGGAAAAAAACACGGACAACAAACAACACCTTGTCCGGGACAACCCCAAggg                                                                                          >  2:99993/1‑140 (MQ=255)
                                            cccAGGCTGGATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAACTCCTTGTCCGGAGAACCCCAACggg                                                                                         >  1:719851/1‑140 (MQ=255)
                                            cccAGGCTGGATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAAAGTGGCATAACAAAAAACTCCCTGTCCGGAGAACCCCAACggg                                                                                         >  2:340083/1‑140 (MQ=255)
                                            cccAGGCTGGATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGGAAAAAAAAAGGGCATAAAAAAAAACCCCCTGTCCGGAGGACCCCAACagg                                                                                         >  2:294324/1‑140 (MQ=255)
                                            cccAGGCTGGATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGAAAAAAAACGTGGCATAAAAAACAACCCCCTGTCCGGAGAACCCCAAAggg                                                                                         >  2:84403/1‑140 (MQ=255)
                                             ccAGGCTGGATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGGATAACAAACAAACCCTTGTCCGGGGAACCCCAACcggg                                                                                        >  1:684938/1‑140 (MQ=255)
                                               aGGCTGGATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACATGGCCTAAAAAAAAAAACCTTTTATTGAGGAACCCAACggggga                                                                                      >  2:942419/1‑138 (MQ=255)
                                                 gCTGGATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGGAAAAAAAAGTGGCATAACAAACAACTCCTTGTCCGGAGAACCCCAACGGGGgaaa                                                                                    >  1:706391/1‑140 (MQ=255)
                                                   tGGATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAACTCCTTGTCCGGAGAACCCCAACGGGGAAAACg                                                                                  <  2:719851/140‑1 (MQ=255)
                                                     gATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGGAAAAAAAAGTGGCCAAAAAAACAAAACCTTTTCCGGAGAACCCCCACGGGGAAAACCaa                                                                                >  1:442370/1‑140 (MQ=255)
                                                     gATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGGAAAAAAAAGAGGCATAACAAAAAAATCCCTTTCCGGGGAACCCCCACGGGGGAAAAGaa                                                                                >  2:689056/1‑140 (MQ=255)
                                                     gATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGGAAAAAAAAATGGGATAAAAAACAAACCCTTTTCCGGGGAACCCCAACGGGGGAAACGaa                                                                                >  2:437619/1‑140 (MQ=255)
                                                      aTCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGGAAAAAAAAGTGGCATAAAAAACAACTCCTTGTCCGGGGAACCCCAACGGGGGAAAAGaaa                                                                               >  2:1529/1‑140 (MQ=255)
                                                           aTCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGGTAAAAAAAGGGGCATTACAAACAACCCCTTGTCCGGAGAACCCCAACGGGGGAAAAGAAATTCaa                                                                          >  1:608688/1‑140 (MQ=255)
                                                           aTCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGGAAAAAAAAGTGGCATAACAAAAAAATCCTTTTCCGGAGAACCCCCACGGGGAAAAAGAAAATCaa                                                                          >  2:16575/1‑140 (MQ=255)
                                                            tCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGGAAAAAAAAGAGGCATAACAAAAAACTCCTTGTCCGGAGAACCCCAACGGGGGAAAAGAAAATCaaa                                                                         >  2:482408/1‑139 (MQ=255)
                                                               gTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAACTCCTTGTCCGGAGAACCCCAACGGGGAAAACGAAATTCAATAcc                                                                      >  1:872107/1‑140 (MQ=255)
                                                               gTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGGAAAAAAAACGGGGATAACAAACAACACCCTGTCCGGGGAACCCCAAGGGGGGAAACGAAAATCAAAAcc                                                                      >  2:401214/1‑140 (MQ=255)
                                                                     aGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGGAAAAAAACGTGGCATAAAAAACAACTCCTTGTCCGGAGAACCCCAACGGGGAAAACGAAAATCAATACCCTCggg                                                                >  2:297056/1‑139 (MQ=255)
                                                                      gATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAACTCCTTGTCCGGAGAACCCCAACGGGGAAAACGAAATTCAATACCGTCGGTa                                                               <  1:482408/140‑1 (MQ=255)
                                                                      gATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAACTCCTTGTCCGGAGAACCCCAACGGGGAAAACGAAATTCAATACCGTCGGTa                                                               <  1:16575/140‑1 (MQ=255)
                                                                       aTACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAACTCCTTGTCCGGAGAACCCCAACGGGGAAAACGAAATTCAATACCGTCGGTaa                                                              <  1:99993/140‑1 (MQ=255)
                                                                       aTACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGGAAAAAAAAGGGGCATAAAAAAAAAACCCCTGTCCGGAGAAACCCAACGGGGGAAAAGAAAATCAAAACCGTCGGGaa                                                              >  2:335875/1‑140 (MQ=255)
                                                                        tACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAACTCCTTGTCCGGAGAACCCCAACGGGGAAAACGAAATTCAATACCGTCGGTAAt                                                             <  2:924510/140‑1 (MQ=255)
                                                                            cagcaAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAACTCCTTGTCCGGAGAACCCCAACGGGGAAAACGAAATTCAATACCGTCGGTAATGACg                                                         <  1:385805/140‑1 (MQ=255)
                                                                            cagcaAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAACTCCTTGTCCGGAGAACCCCAACGGGGAAAACGAAATTCAATACCGTCGGTAATGACg                                                         <  2:713757/140‑1 (MQ=255)
                                                                                 aaGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAACTCCTTGTCCGGAGAACCCCAACGGGGAAAACGAAATTCAATACCGTCGGTAATGACGAAAgg                                                    <  2:476327/140‑1 (MQ=255)
                                                                                  aGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAAAGGGGCATAAAAAAAAACACCTTTTCCGGGGAACCCCAACCGGGGAAACGAAAATTAATACCGTCGGGAATGACGaaaagg                                                   >  2:559561/1‑137 (MQ=255)
                                                                                  aGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGCAAAAAAAAAGGGCCTAACAAACAACTCCCTGGCCGGAGAAACCCAACCGGGAAAAAGAAAATTAATACCCACTGGAATGACGGAAAGc                                                   >  2:815339/1‑140 (MQ=255)
                                                                                        aCAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAACTCCTTGTCCGGAGAACCCCAACGGGGAAAACGAAATTCAATACCGTCGGTAATGACGAAAGGCAGGTga                                             <  1:937828/140‑1 (MQ=255)
                                                                                          aaaaTCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAACTCCTTGTCCGGAGAACCCCAACGGGGAAAACGAAATTCAATACCGTCGGTAATGACGAAAGGCAGGTgaga                                           <  2:523770/140‑1 (MQ=255)
                                                                                              tCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAAAGTGGCATAAAAAACAACTCCTTGTCCGGAGAACCCCAACGGGGAAAACGAAAATCAATACCGTCGGTAATGACGAAAGGCAGGTGAGAAggg                                       >  1:488198/1‑140 (MQ=255)
                                                                                               cAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAACTCCTTGTCCGGAGAACCCCAACGGGGAAAACGAAATTCAATACCGTCGGTAATGACGAAAGGCAGGTGAGAAGGGc                                      <  1:734918/140‑1 (MQ=255)
                                                                                                 aCGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGGAAAAAAAAGTGGCATAACAAAAAACTCCTTGTCCGGAGAACCCCAACGGGGGAAAAGAAAATCAATACCGTCGGGAATGACGAAAGGCAGGGGAGAAGGGGaa                                    >  2:192671/1‑140 (MQ=255)
                                                                                                 aCGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGGAAAAAAAAGAGGCATAACAAACAACTCCTTGTCCGGAGAACCCCCACGGGGGAAAAGAAAATCAAAAACGTCGGGAATGACGAAAAGCAGGGGAGAAAGGGaa                                    >  2:343710/1‑140 (MQ=255)
                                                                                                   gCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAACTCCTTGTCCGGAGAACCCCAACGGGGAAAACGAAATTCAATACCGTCGGTAATGACGAAAGGCAGGTGAGAAGGGCAACg                                  <  1:689056/140‑1 (MQ=255)
                                                                                                    ccATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAACTCCTTGTCCGGAGAACCCCAACGGGGAAAACGAAATTCAATACCGTCGGTAATGACGAAAGGCAGGTGAGAAGGGCAACga                                 <  1:181007/140‑1 (MQ=255)
                                                                                                    ccATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAACTCCTTGTCCGGAGAACCCCAACGGGGAAAACGAAATTCAATACCGTCGGTAATGACGAAAGGCAGGTGAGAAGGGCAACga                                 <  2:1890/140‑1 (MQ=255)
                                                                                                       tCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAACTCCTTGTCCGGAGAACCCCAACGGGGAAAACGAAATTCAATACCGTCGGTAATGACGAAAGGCAGGTGAGAAGGGCAACGAGAc                              <  1:401214/140‑1 (MQ=255)
                                                                                                       tCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAACTCCTTGTCCGGAGAACCCCAACGGGGAAAACGAAATTCAATACCGTCGGTAATGACGAAAGGCAGGTGAGAAGGGCAACGAGAc                              <  1:340083/140‑1 (MQ=255)
                                                                                                       tCGGGAAAAACAGCGTGGCGGCATGGCGGGGAAAAGAACGGGGCATAAAAAAAAAATCCTTGTCCGGGGAACCCCCACGGGGGAAAAGAAATTTAATACCCTCGGGAATGACGAAAAGCAGGGGAGAAGGGGAACgggaa                              >  2:820306/1‑139 (MQ=255)
                                                                                                            aaaaaCAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAACTCCTTGTCCGGAGAACCCCAACGGGGAAAACGAAATTCAATACCGTCGGTAATGACGAAAGGCAGGTGAGAAGGGCAACGAGACTGCaa                         <  2:174032/140‑1 (MQ=255)
                                                                                                                      tGGCGGCATGGCGGGTAAAAAAAAGCGGCATAACAAACAACCCCCTTTCCGGAGAAACCCCACGGGGGAAAAGAAAATCAAAACCGTCCGGAATGAAGAAAGGCAGGGGAGAAGGGGAAAGAGAATGCAAAAGTACTgga               >  2:426257/1‑138 (MQ=255)
                                                                                                                               ggCGGGTAAAAAAACGTGGCATAACAAACAACTCCTTGTCCGGAGAACCCCAACGGGGAAAACGAAATTCAATACCGTCGGTAATGACGAAAGGCAGGTGAGAAGGGCAACGAGACTGCAAATGTACTGAATTTCATGTa      <  2:608688/140‑1 (MQ=255)
                                                                                                                               ggCGGGTAAAAAAACGTGGCATAACAAACAACTCCTTGTCCGGAGAACCCCAACGGGGAAAACGAAATTCAATACCGTCGGTAATGACGAAAGGCAGGTGAGAAGGGCAACGAGACTGCAAATGTACTGAATTTCATGTa      <  2:872107/140‑1 (MQ=255)
                                                                                                                               ggCGGGTAAAAAAACGTGGCATAACAAACAACTCCTTGTCCGGAGAACCCCAACGGGGAAAACGAAATTCAATACCGTCGGTAATGACGAAAGGCAGGTGAGAAGGGCAACGAGACTGCAAATGTACTGAATTTCATGTa      <  1:192671/140‑1 (MQ=255)
                                                                                                                                 cGGGTAAAAAAACGTGGGATAACAAACAACTCCTTGTCCGGAGAACCCCAACGGGGAAAACGAAATTCAATGCCGTCGGTAATGACGAAAGGCAGGTGAGAAGGGCAACGAGACTGCAAATGTACTGAATTTCATGTAAt    <  2:911237/140‑1 (MQ=255)
                                                                                                                                 cGGGTAAAAAAACGTGGCATAACAAACAACTCCTTGTCCGGAGAACCCCAACGGGGAAAACGAAATTCAATACCGTCGGTAATGACGAAAGGCAGGTGAGAAGGGCAACGAGACTGCAAATGTACTGAATTTCATGTAAt    <  2:38383/140‑1 (MQ=255)
                                                                                                                                    ggTAAAACCACGTGGCATAATAAACACTTCCTTGTCCGGACAACCCCATCGGGGAAAACGAAATTCAATACCGTCGGTAATGACGAAAGGCAGGTGAGAAGGGCAACGAGACTGCAAATGTACTGAATTTCATGTAATCa  <  1:636649/140‑1 (MQ=255)
                                                                                                                                      |                                                                                                                                          
GCCAGACTGACATCGGCATTAAAATCACGTACCACATTAATGATCCCAGGCTGGATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAACTCCTTGTCCGGAGAACCCCAACGGGGAAAACGAAATTCAATACCGTCGGTAATGACGAAAGGCAGGTGAGAAGGGCAACGAGACTGCAAATGTACTGAATTTCATGTAATCA  >  NC_000913/4568371‑4568641

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