Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,080,8560TG40.3% 125.6 / 10.0 73intergenic (+43/‑500)putP/efeUproline:sodium symporter/ferrous iron permease (pseudogene)
Reads supporting (aligned to +/- strand):  ref base T (9/34);  new base G (29/0);  total (39/34)
Fisher's exact test for biased strand distribution p-value = 1.87e-12
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTAGTTTGCTGGGTAAAGCGCCGTCAGCGGCGATGCAAAAACGCTTTGCCGAGGCCGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTGTTTGGCTTCAGCAGCGGGTTGCGCTCCCTTAATGTGCCTCGCCATATAAATTGAATGGTGCAGGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCTGTCCCGGCCTTCGGGGAACGCTTCAGCGATTTTGACGCCACCAACACCCGAGCTGTTATTATGT  >  NC_000913/1080710‑1081004
                                                                                                                                                   |                                                                                                                                                     
ttAGTTTGCTGGGTAAAGCGCCGTCAGCGGCGATGCAAAAACGCTTTGCCGAGGCCGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTGTTTGGCTTCAGCAGCGGGTTgc                                                                                                                                                    <  2:569898/150‑1 (MQ=255)
  aGTTTGCTGGGTAAAGCGCCGTCAGCGGCGATGCAAAAACGCTTTGCCGAGGCCGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGCTTTGTTTTGGTTTCGCAGCGGGTTgcgc                                                                                                                                                  >  2:9051/1‑150 (MQ=255)
     ttGCTGGGTAAAGCGCCGTCAGCGGCGATGCAAAAACGCTTTGCCGAGGCCGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTGTTTTGCTTTAACAGCGGGGTGCGCTccc                                                                                                                                              >  1:891273/1‑151 (MQ=255)
        cTGGGTAAAGCGCCGTCAGCGGCGATGCAAAAACGCTTTGCCGAGGCCGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTTTTTTGGTTCAGCCGCGGGGTGCGCccccct                                                                                                                                            >  2:415948/1‑148 (MQ=255)
          gggTAAGGCGCCGTCAGCGGTGATGCAAAAACTCTTTGCCGAGGCCGATGCTCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTGTTTGGCTTCAGCAGCGGGTTGCGCTCCCTTaa                                                                                                                                          <  2:28828/150‑1 (MQ=255)
                aGCGCCGTCAGCGGCGATGCAAAAACGCTTTGCCGAGGCCGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTGTTTGGCTTCAGCAGCGGGTTGCGCTCCCTTAATGTGCCt                                                                                                                                   <  1:2616/151‑1 (MQ=255)
                 gcgcCGTCAGCGGCGATGCAAAAACGCTTTGCCGAGGCCGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTGTTTGGCTTCAGCAGCGGGTTGCGCTCCCTTAATGTTCCTc                                                                                                                                  <  2:997384/151‑1 (MQ=255)
                     cGTCAGCGGCGATGCAAAAACGCTTTGCCGAGGCCGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTGTTTGGCTTCAGCAGCGGGTTGCGCTCCCTTAATGTGCCTCGCCa                                                                                                                              <  1:1301695/151‑1 (MQ=255)
                      gTCAGCGGCGATGCAAAAACGCTTTGCCGAGGCCGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTGTTTGGCTTCAGCAGCGGGGTGCGCTCCCTTAATGTGCCTCGcc                                                                                                                               >  1:709996/1‑149 (MQ=255)
                       tCAGCGGCGATGCAAAAACGCTTTGCCGAGGCCGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTTTTTTGGTTCAGCCGCGGGGTGCGCTCCCTTTATGTGCCTCGcc                                                                                                                               >  1:853667/1‑148 (MQ=255)
                        cAGCGGCGATGCAAAAACGCTTTGCCGAGGCCGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTGTTTGGCTTCAGCAGCGGGTTGCGCTCCCTTAATGTGCCTCGCCata                                                                                                                            <  1:308366/150‑1 (MQ=255)
                         aGCGGCGATGCAAAAACGCTTTGCCGAGGCCGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTTTTTGGCTTCAGCCGCGGGGTGCGCCCCCTTAATGTGCCTCCCCatat                                                                                                                           >  2:808537/1‑150 (MQ=255)
                          gcggcgATGCAAAAACGCTTTGCCGAGGCCGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTGTTTGGCTTCAGCAGCGGGTTGCGCTCCCTTAATGTGCCTCGCCATATaa                                                                                                                         <  1:515602/151‑1 (MQ=255)
                           cggcgATGCAAAAACGCTTTGCCGAGGCCGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTTCTTCGCTTTAGTCAGCGGGCGCGCCCCCTTTAttcgccccccc                                                                                                                               >  2:1008424/1‑134 (MQ=255)
                              cgATGCAAAAACGCTTTGCCGAGGCCGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGCTTTTTTTTGGTTTAGCAGCGGGTTGGGCTCCCCTTATGTGGCCTGGCATATaaaat                                                                                                                      >  2:515876/1‑148 (MQ=255)
                               gATGCAAAAACGCTTTGCCGAGGCCGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTGTTTGGCTTCAGCAGCGGGTTGCGCTCCCTTAATGTGCCTCGCCATATAAATTGa                                                                                                                    <  2:160865/151‑1 (MQ=255)
                               gATGCAAAAACGCTTTGCCGAGGCCGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTGTTTGGCTTCAGCAGCGGGTTGCGCTCCCTTAATGTGCCTCGCCATATAAATTGa                                                                                                                    <  1:847808/151‑1 (MQ=255)
                                   cAAAAACGCTTTGCCGAGGCCGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTGTTTTGCCTCAGCAGCGGGGTTCGCTCCCCTAATGTGTCTCGCCCTATAAATTGAAttg                                                                                                                >  1:430989/1‑149 (MQ=255)
                                     aaaaCGCTTTGCCGAGGCCGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTGTTTGGCTTCAGCAGCGGGTTGCGCTCCCTTAATGTGCCTCGCCATATAAATTGAATGGTg                                                                                                              <  2:660122/151‑1 (MQ=255)
                                      aaaCGCTTTGCCGAGGCCGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTGTTTGGCTTCAGCAGCGGGTTGCGCTCCCTTAATGTGCCTCGCCATATAAATTGAATGGTGc                                                                                                             <  1:1211633/151‑1 (MQ=255)
                                       aaCGCTTTGCCGAGGCCGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTGTTTGGCTTCAGCAGCGGGTTGCGCTCCCTTAATGTGCCTCGCCATATAAATTGAATGGTGCa                                                                                                            <  2:550084/151‑1 (MQ=255)
                                        aCGCTTTGCCGAGGCCGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTTTTTGGCCTCAGCAGCGGGGTTCGCTCCCCTTATGTGCCTCGCCATATAAAATGAATAGTg                                                                                                              >  1:61227/1‑148 (MQ=255)
                                           cTTTGCCGAGGCCGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTTTTTTGCTTTAGCAGCGGGGTGCGCTCCCCTTATTTGCCCCCCCCTATAAAATTAATGGGGCAggg                                                                                                         >  1:406248/1‑150 (MQ=255)
                                             ttGCCGAGGCCGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTGTTTGGCTTCAGCAGCGGGTTGCGCTCCCTTAATGTGCCTCGCCATATAAATTGAATGGTGCAGGGAg                                                                                                       <  1:232642/150‑1 (MQ=255)
                                               gCCGAGGCCGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTGTTTGGCTTCAGCAGCGGGTTGCGCTCCCTTAATGTGCCTCGCCATATAAATTGAATGGTGCAGGGAgcgc                                                                                                    <  2:1303280/151‑1 (MQ=255)
                                               gCCGAGGCCGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTGTTTGGCTTCAGCAGCGGGTTGCGCTCCCTTAATGTGCCTCGCCATATAAATTGAATGGTGCAGGGAgcgc                                                                                                    <  2:253093/151‑1 (MQ=255)
                                                ccGAGGCCGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTGTTTGGCTTCAGCAGCGGGGTGCGCTCCCCTTATTTGCCTCCCCCTATAAATTGAATGGTGCAGGGAggg                                                                                                     >  2:81862/1‑147 (MQ=255)
                                                  gAGGCCGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTTTTTTGCTTCAGCAGCGGGGTGCGCCCCCTTAATGTGCCTCTCCATATAAATTTAATGGGGGAGGGGgcgc                                                                                                    >  2:513798/1‑148 (MQ=255)
                                                  gAGGCCGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTTTTTGGCTTCAGCAGCGGGGTGCGCTCCCTTAATGTGCCTCCCCCTATAAATTGAATGGTGCAGGGAgcgc                                                                                                    >  1:704369/1‑148 (MQ=255)
                                                  gAGGCCGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTGTTTGGCTTCAGCAGCGGGTTGCGCTCCCTTAATGTGCCTCGCCATATAAATTGAATGGTGCAGGGAGCGCGCa                                                                                                 <  2:938219/151‑1 (MQ=255)
                                                    ggCCGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTGTTTGGCTTCAGCAGCGGGTTGCGCTCCCTTAATGTGCCTCGCCATATAAATTGAATGGTGCAGGGAGCGCGCAgg                                                                                               <  2:1180757/151‑1 (MQ=255)
                                                       cGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTGTTTGGGTTCAGCCACGGGGTGCGCTCCCTTAATGTGCCTCGCCATATAAATTGAAAGGGGCAGGGAGCGCGCAggggg                                                                                            >  1:1196555/1‑151 (MQ=255)
                                                       cGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTGTTTGGCTTCAGCAGCGGGTTGCGCTCCCTTAATGTGCCTCGCCATATAAATTGAATGGTGCAGGGAGCGCGCAggggg                                                                                            <  2:741987/151‑1 (MQ=255)
                                                       cGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTGTTTGGCTTCAGCAGCGGGTTGCGCTCCCTTAATGTGCCTCGCCATATAAATTGAATGGTGCAGGGAGCGCGCAggggg                                                                                            <  1:702721/151‑1 (MQ=255)
                                                        gATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTGTTTGGCTTCAGCAGCGGGTTGCGCTCCCTTAATGTGCCTCGCCATATAAATTGAATGGTGCAGGGAGCGCGCAGGGggc                                                                                           <  2:5355/151‑1 (MQ=255)
                                                           gcgcACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTGTTTGGCTTCAGCAGCGGGGTGCGCTCCCTTAATGTGCCTCGCCCTATAAAATGAATGGTGCAGGGAGCGCGCAGGGggcgg                                                                                         >  1:616504/1‑150 (MQ=255)
                                                                   tCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTTTTTTGCTTCAGCAACGGGGTGCGCTCCCTTAATGTGGCTCGCCCTATAAATTGAATGGGGCAGGGGGGGCGCAGGGGGCGGCCCATcgcc                                                                                >  1:698354/1‑151 (MQ=255)
                                                                    cATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTTTTTTGCTTTCGCAGCCGGGTGCGCTCCCTTAATGTGCCCCGCCATATAAATTGGATGGTGCAGGGGGCGCGCCGGGGGCGGCCAATcgcc                                                                                >  2:310146/1‑150 (MQ=255)
                                                                     aTTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTTTTTGGCTTCAGCAGCGGGGTGCGCTCCCCTAATGTGCCTCGCCATATAAATTTAAAGGGGGAGGGAGCGCGCAGGGGGGGGCCAATcgccgc                                                                              >  2:25864/1‑151 (MQ=255)
                                                                      ttCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTGTTTTGCTTCAGCAGCGGGGTGCGCTCCCCTTATTTGCCCCGCCCTATAAATTGAATGGGGCAGGGGGCGGGCAGGGGGGGGGCAAAcggcccc                                                                             >  2:1225294/1‑146 (MQ=255)
                                                                       tCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTTTTTTGCTTCAGCAGCGGGGTGCGCTCCCTTAATTTGCCTCGCCCTATAAATTGGATGGGGCCGGGGGCGCGCAGGGGGGGGGCCAAccccgcc                                                                             >  2:88005/1‑150 (MQ=255)
                                                                        cGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTGTTTTGCTTCCGCAACGGGGTGCGCCCCCTTAATGTGGCTCCCCCTATAAAATTAATTGGGCAGGGAGCGCGCAGGGGGGGGCCAAAcgccgcccc                                                                           >  2:419537/1‑149 (MQ=255)
                                                                           cTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTGTTTGGCCTCCGCCGCGGGGTGCGCCCCCTTTAAGTGCCTCGCCCTATAAAATGAAAGGGGGAGGGGGCGGGCAGGGGGCGGCCCATTCCCGCCGcccc                                                                        >  2:494104/1‑151 (MQ=255)
                                                                            tCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTTTTTTGCTTCTGCCGCGGGTTTCGCTCCCCTAAATTGGCTCGCCCTATTAAATGTATGGGGGAGGGGGCGCGCAGGGGGCGGCCAATCGCcccccccccc                                                                       >  2:1282345/1‑151 (MQ=255)
                                                                            tCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTTTTTGGCTTCAGCAGCGGGGTGCGCCCCCTTAATGTGCCTCGCCATATAAATTGAATGGTGCAGGGAGGGCGCAGGGGGCGGCCAATCgcccccgccccc                                                                       >  2:606867/1‑151 (MQ=255)
                                                                            tCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTGTTTGGCTTCAGCAGCGGGGTGCGCCCCCCTAATGTGGCTCGCCATATAAATTTAATGGGGCAGGGAGCGCGCCGGGGGCGGCCCATCGCcccccccccc                                                                       >  2:1024281/1‑151 (MQ=255)
                                                                              cgccgTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTGTTTGGCTTCAGCAGCGGGTTGCGCTCCCTTAATGTGCCTCGCCATATAAATTGAATGGTGCAGGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCctg                                                                     <  2:683460/151‑1 (MQ=255)
                                                                                       ggTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTTTTTGGCTTCCGCAGCGGGTTGGGGTCCCCTAATGTGCCTCCGCATATAAATTGAATGGGGCAGGGAGCCCGCAGGGGGCGGCCCATCGCCGCCCCCCCccgcgtttccgg                                                            >  2:731630/1‑143 (MQ=255)
                                                                                          tGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTGTTTGGCTTTAGCAGCGGGGTGGGCTCCCTTAATGTGCCCCGCCCTATAAATTGAATGGGGGAGGGGGGGCGCCGGGGGGGGCCAATCCCCCCCGCCCCCCGCTGTCCCgggc                                                          >  2:652094/1‑148 (MQ=255)
                                                                                            cAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTGTTTGGCTTCAGCAGCGGGTTGCGCTCCCTTAATGTGCCTCGCCATATAAATTGAATGGTGCAGGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCTGTCCCGGCCTTc                                                      <  1:198230/151‑1 (MQ=255)
                                                                                                  aGCTAAGGGACTTAGCCTGCGGCGGTTTTGTTTGGCTTCAGCAGCGGGGTGCGCCCCCTTAATGTGCCTCGCCATATAAATTGAATGGGGCAGGGAGCGCGCAGGGGGCGGCCAAACGCCGCCCCCCCCCGCTGTCCCGGCCTTCgggg                                                  >  1:1167634/1‑149 (MQ=255)
                                                                                                    cTAAGGGACTTAGCCTGCGGCGGGTTTTTTTTGCTTCACCAGCGGGGTGCCCGCCCTTAATGTGCCTCGCCCTATAAATTTTATGGTGCAGGGGGGCCGCAGGGGGCGGGCCAACGCCGCCGCCCCCCGCTGTCCCGGCCCTCgggg                                                  >  1:8971/1‑147 (MQ=255)
                                                                                                     tAAGGGACTTAGCCTGCGGCGGTTTTGTTTGGCTTCAGCAGCGGGTTGCGCTCCCTTAATGTGCCTCGCCATATAAATTGAATGGTGCAGGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCTGTCCCGGCCTTCGGGGAAc                                               <  2:1160164/149‑1 (MQ=255)
                                                                                                      aaGGGACTTAGCCTGCGGCGGTTTTGTTTGGCTTCAGCAGCGGGTTGCGCTCCCTTAATGTGCCTCGCCATATAAATTGAATGGTGCAGGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCTGTCCCGGCCTTCGGGGAACGc                                             <  1:636249/150‑1 (MQ=255)
                                                                                                       aGGGACTTAGCCTGCGGCGGTTTTTTTTGGGTTCAACAACGGGGTGCGCTCCCCTAATGTGCCTCCCCCTATAAATTTGATGGTGCAGGGAGCGCGCAGGGGGCGGGCAAACGCCGCCCCCCCCCGCTGTCCCCGCCTTCgggggaagcct                                           >  2:388442/1‑144 (MQ=255)
                                                                                                          gACTTAGCCTGCGGCGGTTTTGTTTGGCTTCAGCAGCGGGGTGGGCTCCCTTAATGTGCCTCGCCATATAAATTGAATGGGGGAGGGAGCGCGCCGGGGGGGGGCCAACGCCCCCGCCCCCTCCTGTCCCCGCCCTCGGGGGACGc                                             >  2:1056775/1‑146 (MQ=255)
                                                                                                             ttAGCCTGCGGCGGTTTTTTTTTGGTTCAACAACGGGGTGCGGTCCCTTTAAATGGCTCGCCCTTTTAATTTAATGGTGCATGGAGCGCGCAGGGGGCGGCCAAACGCCCCCCCCCCCCGCTGTCCCCGCCTTTCGGGGAACCCtacggg                                       >  2:900083/1‑145 (MQ=255)
                                                                                                                    gcggcgGTTTTGTTTGGCTTCAGCAGCGGGTTGCGCTCCCTTAATGTGCCTCGCCATATAAATTGAATGGTGCAGGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCTGTCCCGGGCTTCGGGGAACGCTTCAGCGATTTTGAc                              <  1:1267461/151‑1 (MQ=255)
                                                                                                                    gcggcgGTTTTGTTTGGCTTCAGCAGCGGGTTGCGCTCCCTTAATGTGCCTCGCCATATAAATTGAATGGTGCAGGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCTGTCCCGGCCTTCGGGGAACGCTTCAGCGATTTTGAc                              <  1:824499/151‑1 (MQ=255)
                                                                                                                     cggcggTTTTTTTTGGGTTCAGCAACGGGTTGCGCCCCCCTAATGTGCCTCGCCCTATAAATTGAATGGTGCAGGGGGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCCGCTGTCCCGGCCCTCGGGGGACGCTTCAGCGGTTTTGAc                              >  2:1524/1‑150 (MQ=255)
                                                                                                                     cggcggTTTTGTTTTGCTTTAGCAGCGGGGTGCCCCCCCTTAATGGGCCTCGCCCTATAAATTGAATGGTGCAGGGAGCGCGCAGGGGGCGGCCAAACGCCGCCGCCCCCCGCTGTCCCGGCCTTCGGGGGACGCCTCAGCGAattt                                 >  2:1290320/1‑147 (MQ=37)
                                                                                                                      ggcggTTTTGTTTGGCTTCAGCAGCGGGTTGCGCTCCCTTAATGTGCCTCGCCATATAAATTGAATGGTGCAGGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCTGTCCCGGCCTTCGGGGAACGCTTCAGCGATTTTGACGc                            <  1:1302970/151‑1 (MQ=255)
                                                                                                                         ggTTTTGTTTGGCTTCAGCAGCGGGTTGCGCTCCCTTAATGTGCCTCGCCATATAAATTGAATGGTGCAGGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCTGTCCCGGCCTTCGGGGAACGCTTCAGCGATTTTGACGccac                         >  2:1147289/1‑151 (MQ=255)
                                                                                                                          gTTTTGTTTGGCTTCAGCAGCGGGTTGCGCTCCCTTAATGTGCCTCGCCATATAAATTGAATGGTGCAGGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCTGTCCCGGCCTTCGGGGAACGCTTCAGCGATTTTGACGccacc                        >  2:159101/1‑151 (MQ=255)
                                                                                                                              tgtttgGCTTCAGCAGCGGGTTGCGCTCCCTTAATGTGCCTCGCCATATAAATTGAATGGTGCAGGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCTGTCCCGGCCTTCGGGGAACGCTTCAGCGATTTTGACGCCACCAac                     >  2:335402/1‑150 (MQ=255)
                                                                                                                                  tgGCTTCAGCAGCGGGTTGCGCTCCCTTAATGTGCCTCGCCATATAAATTGAATGGTGCAGGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCTGTCCCGGCCTTCGGGGAACGCTTCAGCGATTTTGACGCCACCAACAccc                 <  2:1035490/150‑1 (MQ=255)
                                                                                                                                   ggCTTCAGCAGCGGGTTGCGCTCCCTTAATGTGCCTCGCCATATAAATTGAATGGTGCAGGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCTGTCCCGGCCTTCGGGGAACGCTTCAGCGATTTTGACGCCACCAACACCCg                <  2:810174/150‑1 (MQ=255)
                                                                                                                                       tCAGCAGCGGGTTGCGCTCCCTTAATGTGCCTCGCCATATAAATTGAATGGTGCAGGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCTGTCCCGGCCTTCGGGGAACGCTTCAGCGATTTTGACGCCACCAACACCCGAGc             <  2:662100/149‑1 (MQ=255)
                                                                                                                                       tCAGCAGCGGGTTGCGCTCCCTTAATGTGCCTCGCCATATAAATTGAATGGTGCAGGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCTGTCCCGGCCTTCGGGGAACGCTTCAGCGATTTTGACGCCACCAACACCCGAGc             <  2:215496/149‑1 (MQ=255)
                                                                                                                                         agcagcGGGTTGCGCTCCCTTAATGTGCCTCGCCATATAAATTGAATGGTGCAGGGAGCGCGCAGGGGGCGACCAATCGCCGCCGCCCCCTGCTGTCCCGGCCTTCGGGGAACGCTTCAGCGATTTTGACGCCACCAACACCCGAGCTGtt         <  2:215268/151‑1 (MQ=255)
                                                                                                                                          gcagcGGGTTGCGCTCCCTTAATGTGCCTCGCCATATAAATTGAATGGTGCAGGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCTGTCCCGGCCTTCGGGGAACGCTTCAGCGATTTTGACGCCACCAACACCCGAGCTGtta        <  2:913408/151‑1 (MQ=255)
                                                                                                                                            cagcGGGTTGCGCTCCCTTAATGTGCCTCGCCATATAAATTGAATGGTGCAGGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCTGTCCCGGCCTTCGGGGAACGCTTCAGCGATTTTGACGCCACCAACACCCGAGCTGttat       >  2:73351/1‑151 (MQ=255)
                                                                                                                                                 ggTTGCGCTCCCTTAATGTGCCTCGCCATATAAATTGAATGGTGCAGGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCTGTCCCGGCCTTCGGGGAACGCTTCAGCGATTTTGACGCCACCAACACCCGAGCTGTTATTATGt  <  2:881811/151‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
TTAGTTTGCTGGGTAAAGCGCCGTCAGCGGCGATGCAAAAACGCTTTGCCGAGGCCGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTGTTTGGCTTCAGCAGCGGGTTGCGCTCCCTTAATGTGCCTCGCCATATAAATTGAATGGTGCAGGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCTGTCCCGGCCTTCGGGGAACGCTTCAGCGATTTTGACGCCACCAACACCCGAGCTGTTATTATGT  >  NC_000913/1080710‑1081004

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