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

TTAGTTTGCTGGGTAAAGCGCCGTCAGCGGCGATGCAAAAACGCTTTGCCGAGGCCGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTGTTTGGCTTCAG‑CAGCGGGTTGCGCTCCCTTAATGTGCCTCGCCATATAAATTGAATGGTGCAGGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCTGTCCCGGCC‑TTCGGGGAACGCTTCAGCGATTTTGACGCCACCAACACCCGAGCTGTTATTATGT  >  NC_000913/1080710‑1081004
                                                                                                                                                   |                                                                                                                                                     
ttAGTTTGCTGGGTAAAGCGCCGTCAGCGGCGATGCAAAAACGCTTTGCCGAGGCCGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTGTTTGGCTTCAG‑CAGCGGGTTgc                                                                                                                                                    <  2:569898/150‑1 (MQ=255)
  aGTTTGCTGGGTAAAGCGCCGTCAGCGGCGATGCAAAAACGCTTTGCCGAGGCCGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGCTTTGTTTTGGTTTCG‑CAGCGGGTTgcgc                                                                                                                                                  >  2:9051/1‑150 (MQ=255)
     ttGCTGGGTAAAGCGCCGTCAGCGGCGATGCAAAAACGCTTTGCCGAGGCCGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTGTTTTGCTTTAA‑CAGCGGGGTGCGCTccc                                                                                                                                              >  1:891273/1‑151 (MQ=255)
        cTGGGTAAAGCGCCGTCAGCGGCGATGCAAAAACGCTTTGCCGAGGCCGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTTTTTTGGTTCAG‑CCGCGGGGTGCGCccccct                                                                                                                                            >  2:415948/1‑148 (MQ=255)
          gggTAAGGCGCCGTCAGCGGTGATGCAAAAACTCTTTGCCGAGGCCGATGCTCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTGTTTGGCTTCAG‑CAGCGGGTTGCGCTCCCTTaa                                                                                                                                          <  2:28828/150‑1 (MQ=255)
                aGCGCCGTCAGCGGCGATGCAAAAACGCTTTGCCGAGGCCGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTGTTTGGCTTCAG‑CAGCGGGTTGCGCTCCCTTAATGTGCCt                                                                                                                                   <  1:2616/151‑1 (MQ=255)
                 gcgcCGTCAGCGGCGATGCAAAAACGCTTTGCCGAGGCCGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTGTTTGGCTTCAG‑CAGCGGGTTGCGCTCCCTTAATGTTCCTc                                                                                                                                  <  2:997384/151‑1 (MQ=255)
                     cGTCAGCGGCGATGCAAAAACGCTTTGCCGAGGCCGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTGTTTGGCTTCAG‑CAGCGGGTTGCGCTCCCTTAATGTGCCTCGCCa                                                                                                                              <  1:1301695/151‑1 (MQ=255)
                      gTCAGCGGCGATGCAAAAACGCTTTGCCGAGGCCGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTGTTTGGCTTCAG‑CAGCGGGGTGCGCTCCCTTAATGTGCCTCGcc                                                                                                                               >  1:709996/1‑149 (MQ=255)
                       tCAGCGGCGATGCAAAAACGCTTTGCCGAGGCCGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTTTTTTGGTTCAG‑CCGCGGGGTGCGCTCCCTTTATGTGCCTCGcc                                                                                                                               >  1:853667/1‑148 (MQ=255)
                        cAGCGGCGATGCAAAAACGCTTTGCCGAGGCCGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTGTTTGGCTTCAG‑CAGCGGGTTGCGCTCCCTTAATGTGCCTCGCCata                                                                                                                            <  1:308366/150‑1 (MQ=255)
                         aGCGGCGATGCAAAAACGCTTTGCCGAGGCCGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTTTTTGGCTTCAG‑CCGCGGGGTGCGCCCCCTTAATGTGCCTCCCCatat                                                                                                                           >  2:808537/1‑150 (MQ=255)
                          gcggcgATGCAAAAACGCTTTGCCGAGGCCGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTGTTTGGCTTCAG‑CAGCGGGTTGCGCTCCCTTAATGTGCCTCGCCATATaa                                                                                                                         <  1:515602/151‑1 (MQ=255)
                           cggcgATGCAAAAACGCTTTGCCGAGGCCGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTTCTTCGCTTTAGTCAGCGGG‑CGCGCCCCCTTTAttcgccccccc                                                                                                                               >  2:1008424/1‑134 (MQ=255)
                              cgATGCAAAAACGCTTTGCCGAGGCCGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGCTTTTTTTTGGTTTAG‑CAGCGGGTTGGGCTCCCCTTATGTGGCCTGGCATATaaaat                                                                                                                      >  2:515876/1‑148 (MQ=255)
                               gATGCAAAAACGCTTTGCCGAGGCCGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTGTTTGGCTTCAG‑CAGCGGGTTGCGCTCCCTTAATGTGCCTCGCCATATAAATTGa                                                                                                                    <  2:160865/151‑1 (MQ=255)
                               gATGCAAAAACGCTTTGCCGAGGCCGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTGTTTGGCTTCAG‑CAGCGGGTTGCGCTCCCTTAATGTGCCTCGCCATATAAATTGa                                                                                                                    <  1:847808/151‑1 (MQ=255)
                                   cAAAAACGCTTTGCCGAGGCCGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTGTTTTGCCTCAG‑CAGCGGGGTTCGCTCCCCTAATGTGTCTCGCCCTATAAATTGAAttg                                                                                                                >  1:430989/1‑149 (MQ=255)
                                     aaaaCGCTTTGCCGAGGCCGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTGTTTGGCTTCAG‑CAGCGGGTTGCGCTCCCTTAATGTGCCTCGCCATATAAATTGAATGGTg                                                                                                              <  2:660122/151‑1 (MQ=255)
                                      aaaCGCTTTGCCGAGGCCGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTGTTTGGCTTCAG‑CAGCGGGTTGCGCTCCCTTAATGTGCCTCGCCATATAAATTGAATGGTGc                                                                                                             <  1:1211633/151‑1 (MQ=255)
                                       aaCGCTTTGCCGAGGCCGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTGTTTGGCTTCAG‑CAGCGGGTTGCGCTCCCTTAATGTGCCTCGCCATATAAATTGAATGGTGCa                                                                                                            <  2:550084/151‑1 (MQ=255)
                                        aCGCTTTGCCGAGGCCGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTTTTTGGCCTCAG‑CAGCGGGGTTCGCTCCCCTTATGTGCCTCGCCATATAAAATGAATAGTg                                                                                                              >  1:61227/1‑148 (MQ=255)
                                           cTTTGCCGAGGCCGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTTTTTTGCTTTAG‑CAGCGGGGTGCGCTCCCCTTATTTGCCCCCCCCTATAAAATTAATGGGGCAggg                                                                                                         >  1:406248/1‑150 (MQ=255)
                                             ttGCCGAGGCCGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTGTTTGGCTTCAG‑CAGCGGGTTGCGCTCCCTTAATGTGCCTCGCCATATAAATTGAATGGTGCAGGGAg                                                                                                       <  1:232642/150‑1 (MQ=255)
                                               gCCGAGGCCGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTGTTTGGCTTCAG‑CAGCGGGTTGCGCTCCCTTAATGTGCCTCGCCATATAAATTGAATGGTGCAGGGAgcgc                                                                                                    <  2:1303280/151‑1 (MQ=255)
                                               gCCGAGGCCGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTGTTTGGCTTCAG‑CAGCGGGTTGCGCTCCCTTAATGTGCCTCGCCATATAAATTGAATGGTGCAGGGAgcgc                                                                                                    <  2:253093/151‑1 (MQ=255)
                                                ccGAGGCCGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTGTTTGGCTTCAG‑CAGCGGGGTGCGCTCCCCTTATTTGCCTCCCCCTATAAATTGAATGGTGCAGGGAggg                                                                                                     >  2:81862/1‑147 (MQ=255)
                                                  gAGGCCGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTTTTTTGCTTCAG‑CAGCGGGGTGCGCCCCCTTAATGTGCCTCTCCATATAAATTTAATGGGGGAGGGGgcgc                                                                                                    >  2:513798/1‑148 (MQ=255)
                                                  gAGGCCGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTTTTTGGCTTCAG‑CAGCGGGGTGCGCTCCCTTAATGTGCCTCCCCCTATAAATTGAATGGTGCAGGGAgcgc                                                                                                    >  1:704369/1‑148 (MQ=255)
                                                  gAGGCCGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTGTTTGGCTTCAG‑CAGCGGGTTGCGCTCCCTTAATGTGCCTCGCCATATAAATTGAATGGTGCAGGGAGCGCGCa                                                                                                 <  2:938219/151‑1 (MQ=255)
                                                    ggCCGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTGTTTGGCTTCAG‑CAGCGGGTTGCGCTCCCTTAATGTGCCTCGCCATATAAATTGAATGGTGCAGGGAGCGCGCAgg                                                                                               <  2:1180757/151‑1 (MQ=255)
                                                       cGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTGTTTGGGTTCAG‑CCACGGGGTGCGCTCCCTTAATGTGCCTCGCCATATAAATTGAAAGGGGCAGGGAGCGCGCAggggg                                                                                            >  1:1196555/1‑151 (MQ=255)
                                                       cGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTGTTTGGCTTCAG‑CAGCGGGTTGCGCTCCCTTAATGTGCCTCGCCATATAAATTGAATGGTGCAGGGAGCGCGCAggggg                                                                                            <  2:741987/151‑1 (MQ=255)
                                                       cGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTGTTTGGCTTCAG‑CAGCGGGTTGCGCTCCCTTAATGTGCCTCGCCATATAAATTGAATGGTGCAGGGAGCGCGCAggggg                                                                                            <  1:702721/151‑1 (MQ=255)
                                                        gATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTGTTTGGCTTCAG‑CAGCGGGTTGCGCTCCCTTAATGTGCCTCGCCATATAAATTGAATGGTGCAGGGAGCGCGCAGGGggc                                                                                           <  2:5355/151‑1 (MQ=255)
                                                           gcgcACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTGTTTGGCTTCAG‑CAGCGGGGTGCGCTCCCTTAATGTGCCTCGCCCTATAAAATGAATGGTGCAGGGAGCGCGCAGGGggcgg                                                                                         >  1:616504/1‑150 (MQ=255)
                                                                   tCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTTTTTTGCTTCAG‑CAACGGGGTGCGCTCCCTTAATGTGGCTCGCCCTATAAATTGAATGGGGCAGGGGGGGCGCAGGGGGCGGCCCATcgcc                                                                                >  1:698354/1‑151 (MQ=255)
                                                                    cATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTTTTTTGCTTTCG‑CAGCCGGGTGCGCTCCCTTAATGTGCCCCGCCATATAAATTGGATGGTGCAGGGGGCGCGCCGGGGGCGGCCAATcgcc                                                                                >  2:310146/1‑150 (MQ=255)
                                                                     aTTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTTTTTGGCTTCAG‑CAGCGGGGTGCGCTCCCCTAATGTGCCTCGCCATATAAATTTAAAGGGGGAGGGAGCGCGCAGGGGGGGGCCAATcgccgc                                                                              >  2:25864/1‑151 (MQ=255)
                                                                      ttCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTGTTTTGCTTCAG‑CAGCGGGGTGCGCTCCCCTTATTTGCCCCGCCCTATAAATTGAATGGGGCAGGGGGCGGGCAGGGGGGGGGCAAAcggcccc                                                                             >  2:1225294/1‑146 (MQ=255)
                                                                       tCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTTTTTTGCTTCAG‑CAGCGGGGTGCGCTCCCTTAATTTGCCTCGCCCTATAAATTGGATGGGGCCGGGGGCGCGCAGGGGGGGGGCCAAccccgcc                                                                             >  2:88005/1‑150 (MQ=255)
                                                                        cGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTGTTTTGCTTCCG‑CAACGGGGTGCGCCCCCTTAATGTGGCTCCCCCTATAAAATTAATTGGGCAGGGAGCGCGCAGGGGGGGGCCAAAcgccgcccc                                                                           >  2:419537/1‑149 (MQ=255)
                                                                           cTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTGTTTGGCCTCCG‑CCGCGGGGTGCGCCCCCTTTAAGTGCCTCGCCCTATAAAATGAAAGGGGGAGGGGGCGGGCAGGGGGCGGCCCATTCCCGCCGcccc                                                                        >  2:494104/1‑151 (MQ=255)
                                                                            tCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTTTTTTGCTTCTG‑CCGCGGGTTTCGCTCCCCTAAATTGGCTCGCCCTATTAAATGTATGGGGGAGGGGGCGCGCAGGGGGCGGCCAATCGCcccccccccc                                                                       >  2:1282345/1‑151 (MQ=255)
                                                                            tCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTTTTTGGCTTCAG‑CAGCGGGGTGCGCCCCCTTAATGTGCCTCGCCATATAAATTGAATGGTGCAGGGAGGGCGCAGGGGGCGGCCAATCgcccccgccccc                                                                       >  2:606867/1‑151 (MQ=255)
                                                                            tCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTGTTTGGCTTCAG‑CAGCGGGGTGCGCCCCCCTAATGTGGCTCGCCATATAAATTTAATGGGGCAGGGAGCGCGCCGGGGGCGGCCCATCGCcccccccccc                                                                       >  2:1024281/1‑151 (MQ=255)
                                                                              cgccgTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTGTTTGGCTTCAG‑CAGCGGGTTGCGCTCCCTTAATGTGCCTCGCCATATAAATTGAATGGTGCAGGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCctg                                                                     <  2:683460/151‑1 (MQ=255)
                                                                                       ggTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTTTTTGGCTTCCG‑CAGCGGGTTGGGGTCCCCTAATGTGCCTCCGCATATAAATTGAATGGGGCAGGGAGCCCGCAGGGGGCGGCCCATCGCCGCCCCCCCccgcgtttccgg                                                            >  2:731630/1‑143 (MQ=255)
                                                                                          tGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTGTTTGGCTTTAG‑CAGCGGGGTGGGCTCCCTTAATGTGCCCCGCCCTATAAATTGAATGGGGGAGGGGGGGCGCCGGGGGGGGCCAATCCCCCCCGCCCCCCGCTGTCCCgggc                                                          >  2:652094/1‑148 (MQ=255)
                                                                                            cAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTGTTTGGCTTCAG‑CAGCGGGTTGCGCTCCCTTAATGTGCCTCGCCATATAAATTGAATGGTGCAGGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCTGTCCCGGCC‑TTc                                                      <  1:198230/151‑1 (MQ=255)
                                                                                                  aGCTAAGGGACTTAGCCTGCGGCGGTTTTGTTTGGCTTCAG‑CAGCGGGGTGCGCCCCCTTAATGTGCCTCGCCATATAAATTGAATGGGGCAGGGAGCGCGCAGGGGGCGGCCAAACGCCGCCCCCCCCCGCTGTCCCGGCC‑TTCgggg                                                  >  1:1167634/1‑149 (MQ=255)
                                                                                                    cTAAGGGACTTAGCCTGCGGCGGGTTTTTTTTGCTTCAC‑CAGCGGGGTGCCCGCCCTTAATGTGCCTCGCCCTATAAATTTTATGGTGCAGGGGGGCCGCAGGGGGCGGGCCAACGCCGCCGCCCCCCGCTGTCCCGGCC‑CTCgggg                                                  >  1:8971/1‑147 (MQ=255)
                                                                                                     tAAGGGACTTAGCCTGCGGCGGTTTTGTTTGGCTTCAG‑CAGCGGGTTGCGCTCCCTTAATGTGCCTCGCCATATAAATTGAATGGTGCAGGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCTGTCCCGGCC‑TTCGGGGAAc                                               <  2:1160164/149‑1 (MQ=255)
                                                                                                      aaGGGACTTAGCCTGCGGCGGTTTTGTTTGGCTTCAG‑CAGCGGGTTGCGCTCCCTTAATGTGCCTCGCCATATAAATTGAATGGTGCAGGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCTGTCCCGGCC‑TTCGGGGAACGc                                             <  1:636249/150‑1 (MQ=255)
                                                                                                       aGGGACTTAGCCTGCGGCGGTTTTTTTTGGGTTCAA‑CAACGGGGTGCGCTCCCCTAATGTGCCTCCCCCTATAAATTTGATGGTGCAGGGAGCGCGCAGGGGGCGGGCAAACGCCGCCCCCCCCCGCTGTCCCCGCC‑TTCgggggaagcct                                           >  2:388442/1‑144 (MQ=255)
                                                                                                          gACTTAGCCTGCGGCGGTTTTGTTTGGCTTCAG‑CAGCGGGGTGGGCTCCCTTAATGTGCCTCGCCATATAAATTGAATGGGGGAGGGAGCGCGCCGGGGGGGGGCCAACGCCCCCGCCCCCTCCTGTCCCCGCC‑CTCGGGGGACGc                                             >  2:1056775/1‑146 (MQ=255)
                                                                                                             ttAGCCTGCGGCGGTTTTTTTTTGGTTCAA‑CAACGGGGTGCGGTCCCTTTAAATGGCTCGCCCTTTTAATTTAATGGTGCATGGAGCGCGCAGGGGGCGGCCAAACGCCCCCCCCCCCCGCTGTCCCCGCCTTTCGGGGAACCCtacggg                                       >  2:900083/1‑145 (MQ=255)
                                                                                                                    gcggcgGTTTTGTTTGGCTTCAG‑CAGCGGGTTGCGCTCCCTTAATGTGCCTCGCCATATAAATTGAATGGTGCAGGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCTGTCCCGGGC‑TTCGGGGAACGCTTCAGCGATTTTGAc                              <  1:1267461/151‑1 (MQ=255)
                                                                                                                    gcggcgGTTTTGTTTGGCTTCAG‑CAGCGGGTTGCGCTCCCTTAATGTGCCTCGCCATATAAATTGAATGGTGCAGGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCTGTCCCGGCC‑TTCGGGGAACGCTTCAGCGATTTTGAc                              <  1:824499/151‑1 (MQ=255)
                                                                                                                     cggcggTTTTTTTTGGGTTCAG‑CAACGGGTTGCGCCCCCCTAATGTGCCTCGCCCTATAAATTGAATGGTGCAGGGGGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCCGCTGTCCCGGCC‑CTCGGGGGACGCTTCAGCGGTTTTGAc                              >  2:1524/1‑150 (MQ=255)
                                                                                                                     cggcggTTTTGTTTTGCTTTAG‑CAGCGGGGTGCCCCCCCTTAATGGGCCTCGCCCTATAAATTGAATGGTGCAGGGAGCGCGCAGGGGGCGGCCAAACGCCGCCGCCCCCCGCTGTCCCGGCC‑TTCGGGGGACGCCTCAGCGAattt                                 >  2:1290320/1‑147 (MQ=37)
                                                                                                                      ggcggTTTTGTTTGGCTTCAG‑CAGCGGGTTGCGCTCCCTTAATGTGCCTCGCCATATAAATTGAATGGTGCAGGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCTGTCCCGGCC‑TTCGGGGAACGCTTCAGCGATTTTGACGc                            <  1:1302970/151‑1 (MQ=255)
                                                                                                                         ggTTTTGTTTGGCTTCAG‑CAGCGGGTTGCGCTCCCTTAATGTGCCTCGCCATATAAATTGAATGGTGCAGGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCTGTCCCGGCC‑TTCGGGGAACGCTTCAGCGATTTTGACGccac                         >  2:1147289/1‑151 (MQ=255)
                                                                                                                          gTTTTGTTTGGCTTCAG‑CAGCGGGTTGCGCTCCCTTAATGTGCCTCGCCATATAAATTGAATGGTGCAGGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCTGTCCCGGCC‑TTCGGGGAACGCTTCAGCGATTTTGACGccacc                        >  2:159101/1‑151 (MQ=255)
                                                                                                                              tgtttgGCTTCAG‑CAGCGGGTTGCGCTCCCTTAATGTGCCTCGCCATATAAATTGAATGGTGCAGGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCTGTCCCGGCC‑TTCGGGGAACGCTTCAGCGATTTTGACGCCACCAac                     >  2:335402/1‑150 (MQ=255)
                                                                                                                                  tgGCTTCAG‑CAGCGGGTTGCGCTCCCTTAATGTGCCTCGCCATATAAATTGAATGGTGCAGGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCTGTCCCGGCC‑TTCGGGGAACGCTTCAGCGATTTTGACGCCACCAACAccc                 <  2:1035490/150‑1 (MQ=255)
                                                                                                                                   ggCTTCAG‑CAGCGGGTTGCGCTCCCTTAATGTGCCTCGCCATATAAATTGAATGGTGCAGGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCTGTCCCGGCC‑TTCGGGGAACGCTTCAGCGATTTTGACGCCACCAACACCCg                <  2:810174/150‑1 (MQ=255)
                                                                                                                                       tCAG‑CAGCGGGTTGCGCTCCCTTAATGTGCCTCGCCATATAAATTGAATGGTGCAGGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCTGTCCCGGCC‑TTCGGGGAACGCTTCAGCGATTTTGACGCCACCAACACCCGAGc             <  2:662100/149‑1 (MQ=255)
                                                                                                                                       tCAG‑CAGCGGGTTGCGCTCCCTTAATGTGCCTCGCCATATAAATTGAATGGTGCAGGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCTGTCCCGGCC‑TTCGGGGAACGCTTCAGCGATTTTGACGCCACCAACACCCGAGc             <  2:215496/149‑1 (MQ=255)
                                                                                                                                         ag‑cagcGGGTTGCGCTCCCTTAATGTGCCTCGCCATATAAATTGAATGGTGCAGGGAGCGCGCAGGGGGCGACCAATCGCCGCCGCCCCCTGCTGTCCCGGCC‑TTCGGGGAACGCTTCAGCGATTTTGACGCCACCAACACCCGAGCTGtt         <  2:215268/151‑1 (MQ=255)
                                                                                                                                          g‑cagcGGGTTGCGCTCCCTTAATGTGCCTCGCCATATAAATTGAATGGTGCAGGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCTGTCCCGGCC‑TTCGGGGAACGCTTCAGCGATTTTGACGCCACCAACACCCGAGCTGtta        <  2:913408/151‑1 (MQ=255)
                                                                                                                                            cagcGGGTTGCGCTCCCTTAATGTGCCTCGCCATATAAATTGAATGGTGCAGGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCTGTCCCGGCC‑TTCGGGGAACGCTTCAGCGATTTTGACGCCACCAACACCCGAGCTGttat       >  2:73351/1‑151 (MQ=255)
                                                                                                                                                 ggTTGCGCTCCCTTAATGTGCCTCGCCATATAAATTGAATGGTGCAGGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCTGTCCCGGCC‑TTCGGGGAACGCTTCAGCGATTTTGACGCCACCAACACCCGAGCTGTTATTATGt  <  2:881811/151‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
TTAGTTTGCTGGGTAAAGCGCCGTCAGCGGCGATGCAAAAACGCTTTGCCGAGGCCGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTGTTTGGCTTCAG‑CAGCGGGTTGCGCTCCCTTAATGTGCCTCGCCATATAAATTGAATGGTGCAGGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCTGTCCCGGCC‑TTCGGGGAACGCTTCAGCGATTTTGACGCCACCAACACCCGAGCTGTTATTATGT  >  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: ‑