Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,508,6620TA41.1% 170.9 / 44.6 112E384D (GAA→GATyhfWphosphopentomutase‑related metalloenzyme superfamily protein
Reads supporting (aligned to +/- strand):  ref base T (20/46);  new base A (44/2);  total (64/48)
Fisher's exact test for biased strand distribution p-value = 3.58e-13
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 114 total aligned reads displayed.

CACGTTGTTTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTGATCCGCCATCAC  >  NC_000913/3508513‑3508811
                                                                                                                                                       |                                                                                                                                                     
cACGTTGTTTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAAtt                                                                                                                                                      <  1:751332/151‑1 (MQ=255)
  cGTTGTTTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTca                                                                                                                                                    <  1:113429/151‑1 (MQ=255)
     tgttTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGGGGGAAAAATaaca                                                                                                                                                  >  4:428129/1‑150 (MQ=255)
     tgttTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGGAAAAATcaca                                                                                                                                                  >  3:117251/1‑150 (MQ=255)
      gttTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTcacaca                                                                                                                                                <  1:939303/151‑1 (MQ=255)
      gttTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTcaca                                                                                                                                                  >  2:1246974/1‑149 (MQ=255)
           tAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGGAAAAAttaaa                                                                                                                                                   >  4:488096/1‑141 (MQ=255)
            aaaGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGAGGGGCGGGGGAAAAATTACACACgg                                                                                                                                             >  1:969703/1‑148 (MQ=255)
             aaGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGc                                                                                                                                         <  2:381503/151‑1 (MQ=255)
              aGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGGAAAAATCacacac                                                                                                                                               >  2:650969/1‑144 (MQ=255)
                   gagagaTCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGAGCGGAAAAATTcacaccagctaagccccc                                                                                                                                    >  1:995451/1‑137 (MQ=255)
                    agagaTCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGGAAAAATTACAAACGGGAGCCCCcac                                                                                                                                   >  2:1120835/1‑150 (MQ=255)
                       gaTCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCa                                                                                                                               <  1:1302669/151‑1 (MQ=255)
                       gaTCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGGGGCGGGCGGGAAAAATAACAACCGGGAGgcccccc                                                                                                                                   >  1:348417/1‑145 (MQ=255)
                       gaTCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGAGGCGGGCGGGAAAAATAAAACACgggcggccccccaac                                                                                                                                >  1:827955/1‑137 (MQ=255)
                         tCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGCGGGGCGCGGGAAAAAACACACACGGTAGccccccc                                                                                                                                   >  1:1103952/1‑143 (MQ=255)
                          cAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGGGGCGCGCGGGAAAAATAAAAAACGGGAGCCCCCACATCa                                                                                                                               >  1:377438/1‑148 (MQ=255)
                           aaTATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGaaa                                                                                                                           <  4:159001/151‑1 (MQ=255)
                           aaTATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGCGTGGGGCGGGAAAAATAAAACCCGGGAcccccccccctct                                                                                                                              >  1:36445/1‑141 (MQ=255)
                            atatGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAg                                                                                                                          <  1:1313776/151‑1 (MQ=255)
                             tatGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGGAAAAATTACAAACCGGAGCCCCCCCATCAGaa                                                                                                                            >  1:825564/1‑148 (MQ=255)
                             tatGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGCGGGGCGCGGGAAAAAATACACACCGGAGACCCCCCATccaaaa                                                                                                                           >  1:1228760/1‑144 (MQ=255)
                              atGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGGAAAATTCACACACCGGACCCCCCACCTccaaaaac                                                                                                                         >  1:1294478/1‑143 (MQ=255)
                                   gATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCgtggtg                                                                                                                   <  3:173558/151‑1 (MQ=255)
                                     tGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGCGCGGGAAAAATCACACACGGTAGCCCCCCCCTCAGAAAGCgggg                                                                                                                     >  2:1651423/1‑147 (MQ=255)
                                      gCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGGAAAAATTACACACGGTAGCCCCCCCATCAGAAAACGTGGTgcg                                                                                                                 >  3:563510/1‑150 (MQ=255)
                                      gCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGGAAAAATCACACACGGGAGCCCCCCCATCAGAAAACGGGGGgcg                                                                                                                 >  4:453602/1‑150 (MQ=255)
                                         cATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGGAAAAATTACAAACCGGAGACCCCCCATCAGAAAACgtgggggg                                                                                                                 >  2:152361/1‑143 (MQ=255)
                                          aTGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGGGGGGGAAAAATTAAAAACCGGAGCCCCCACATCAGaaaagggggtggg                                                                                                                 >  2:614220/1‑136 (MQ=255)
                                            aggCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGCGCGGGAAAAATTACACACCGGAGCCCCCCCATCCGAAAGAGGGGGGGGCACAcc                                                                                                           >  1:918273/2‑150 (MQ=255)
                                             ggCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGa                                                                                                         <  3:92143/151‑1 (MQ=255)
                                                gAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGCGCGGGGAAAATCACACACGGGAGCCCCCCCATCAGAAAAAGGGGTGCGCACACCGa                                                                                                         >  1:1494621/1‑148 (MQ=255)
                                                gAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCTCCCACATCAGAAAGCGTGGTGCGCACACCGAGCt                                                                                                      <  2:1646047/151‑1 (MQ=255)
                                                 aGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGCGCGGGAAAAATTACACACGGGAACCCCCACACCCGAAAGCGGGGGGCGCACACCCAGATg                                                                                                     >  1:937094/1‑151 (MQ=255)
                                                 aGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTg                                                                                                     <  2:446910/151‑1 (MQ=255)
                                                   ggTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCg                                                                                                   <  3:413185/151‑1 (MQ=255)
                                                     tgtAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGCGCGGGAAAAATTACACACGGTAGCCCCCCCATCAGAAAGCGGGGGGGGCCAACCCAGCcgcg                                                                                                   >  3:205723/1‑149 (MQ=255)
                                                     tgtAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAAATCACACCCGGGAGCCCCCCCATCAGAAAACGTGGGGCGCACACCGAGCTg                                                                                                     >  2:1630068/1‑147 (MQ=255)
                                                     tgtAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGATGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGt                                                                                                  >  2:458466/1‑150 (MQ=255)
                                                      gtAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAg                                                                                                <  2:852452/151‑1 (MQ=255)
                                                        aaCCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGGGGGAAAAATCACAAACGGTAGCCCCCCCACCAGAAAGCGGGGGGGGCACACCGAGCTGCGGAg                                                                                                >  1:259629/1‑149 (MQ=255)
                                                        aaCCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAga                                                                                               <  1:1503282/150‑2 (MQ=255)
                                                         aCCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGa                                                                                             <  2:543481/151‑1 (MQ=255)
                                                         aCCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGa                                                                                             <  2:386357/151‑1 (MQ=255)
                                                         aCCCGTCGGATCAAAACTCATAGGGTGTCTCATGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGAGGCGGCGGGCGGGAAAAAAAAAACACGGGAGGCACCACCACAAAAAAAGTGGGGGGGACAcccgaccgagtag                                                                                                >  1:1299212/1‑137 (MQ=255)
                                                          cccGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGGAAAAATCACACACGGGAGCCCCCCCATCAGAAAACGGGGGGCGCCCACCCAGCTGCGTAGCggt                                                                                            >  4:504863/1‑149 (MQ=255)
                                                           ccGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGat                                                                                            <  3:301668/150‑1 (MQ=255)
                                                            cGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGCGCGGGAAAAATCAAACACGGTAACCCCCCCCTCCAAAAACGGGGGGCGCCCACCGAGCCGCGTCGCGatc                                                                                           >  1:1673544/1‑150 (MQ=255)
                                                            cGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGatca                                                                                          <  1:134819/151‑1 (MQ=255)
                                                             gTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGAGGCGCGCGGGAAAAATTACACACGGTAGCCCCCACAACAGAAAACGGGGTGCGCACACCCAGCTGCGGAGCGgtca                                                                                          >  2:76164/1‑150 (MQ=255)
                                                               cGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGGAAAAATAACACACGGGAACCCCCACATCAGAAAAAGGGGGGGGCACACCGAGCTGCGGAGCGATCATcc                                                                                       >  1:189808/1‑151 (MQ=255)
                                                               cGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATcc                                                                                       <  3:60337/151‑1 (MQ=255)
                                                               aggATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTTGCGCGCGGGAAAATTCACACACGGGAGCCCCCACATCAGAAAACGTGGGGCGCACACCGAGCTGCGTAGCGATCATcc                                                                                       >  2:354824/2‑151 (MQ=255)
                                                                  aTCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCtt                                                                                    <  2:603013/151‑1 (MQ=255)
                                                                     aaaaCTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGc                                                                                  <  2:703300/150‑1 (MQ=255)
                                                                       aaCTCATAGGGTGTCTCCTGCAAAGCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGa                                                                               <  2:960035/151‑1 (MQ=255)
                                                                            aTAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATaaa                                                                           <  2:83080/150‑1 (MQ=255)
                                                                            aTAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAAc                                                                          <  1:1019157/151‑1 (MQ=255)
                                                                              aGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCAAACAAGGGTGGCCCCCCACCAAGAAAAGTGGGGCGCCCACCGAGCCGCGGAACGATCATCCCCTTCTGGTaaa                                                                           >  2:1301096/1‑148 (MQ=255)
                                                                                ggTGTCTCCTGCAAACCGGAGGTAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCACGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCa                                                                        <  2:1138002/149‑1 (MQ=255)
                                                                                  tgtCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGCGCGGGAAAAATTAAAAAAGGTAGCCCCCACACCAGAAAACGTGGTGGGCACACCCCGCTGCGGAGCGATCCTCCCTTGCTGATAAACCaac                                                                      >  4:305765/1‑147 (MQ=255)
                                                                                   gtCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCAc                                                                    <  1:624160/150‑1 (MQ=255)
                                                                                   gtCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACt                                                                   <  1:1271464/151‑1 (MQ=255)
                                                                                     ctcCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCCCGGAAAAATTCAAACACGGTAGCCCCCCCATCAGAAAAAGGGGGGCGGACACCCAGCTGCGCAGAGAACATCCCCTTCTGGTAAAACAGCACttg                                                                 >  1:1210563/1‑149 (MQ=255)
                                                                                       ccTGAAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAAATCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTCGCTGATAAACAAGCACTGGca                                                               <  1:665597/151‑1 (MQ=255)
                                                                                        cTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGGGCGGGAAAAATCAAAAACAGGATGCCCCACAACCAAAAGAGTGGTGCGCACACCGCGCTGCGTAGCGATCATCCCCTTCTGATAAAACAACACTGGcac                                                              >  1:1625120/1‑151 (MQ=255)
                                                                                              aCCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTc                                                        <  4:189614/151‑1 (MQ=255)
                                                                                                 ggagggagGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTcgcg                                                     <  2:783287/151‑1 (MQ=255)
                                                                                                  gagggaggAAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGCGCGGGAAAAATCAAACACGGTAGCCCCCCCATCAGAAAACGGGGGGCGCCCACCCAGCCGCGTAGCGATCATCCCCTGCTGATAAAACAACACCGGGACCACCTCCCgg                                                    >  2:1336959/1‑151 (MQ=255)
                                                                                                   agggaggAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGGAAAATTCACACACCGTAGCCCCCCCATCAGAAAACGGGGGGCGCACACCCAGCTGCGTAGCGATCATCCCTTTCTGATAAAACAGCACTGGCACCACCTCCCggg                                                   >  2:1515346/1‑151 (MQ=255)
                                                                                                   agggaggAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAAATCACACACCGTAGCCCCCACATCAGAAAGCGTGGGGCGCACAACGAGCTGCGTAACGATCATCCCTTGCTGGTAAACCAGCACTGGCACCACCTCGCggg                                                   >  2:1435085/1‑151 (MQ=255)
                                                                                                   agggaggAAAGAAAAGAGCGACCATTTTGTGGCGGCGGGGCGCGGGAAAAATCACAAACGGGAGCCCCCCCCACAGAAAAGGTGGGGGGCGCAACCAGCCGCGAGGTGATCATCCCCTGCTGTTAAAACAGCCCTGGACCCACTTAGCggg                                                   >  1:1140759/1‑151 (MQ=255)
                                                                                                    gggaggAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGGAAAAATCACACACGGGAGCCCCCACCTCAGAAAGCGTGGTGCGCACACCGAGCCGCGGAGCGATCATCCCCTGCTGATAAAACAGCACTGGCACCACCTCGCgggg                                                  >  2:351149/1‑150 (MQ=255)
                                                                                                     ggaggaAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGCGCGGGAAAAATTACACACGGGAACCCCCCCATCAGAAAACGGGGGGCGCACAACGAGGCGCGGAGCGATCATCCCCTGCTGATAAAACAACACCGGCACCCCCTCGCggg                                                   >  2:338173/1‑149 (MQ=255)
                                                                                                       aggaAAGAAAAGAGCGACCATTTTGTGGCGGTGGTGCGGGGGAAAAATTAAAAAAGGGAAGCCCCCCCTCCAAAAACGGGGGGGGGGCAACGAGCCGCGGAAGGATCATCCCTTTCTTATAAAACAGAACTGGcaccaccttccggg                                                   >  4:360992/1‑139 (MQ=255)
                                                                                                        ggaAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGCGCGGGAAAAATTACACACGGTAACCCCCCCAACAGAAAAAGGGGGGCGCACACCGAGCTGCGTAGCGATCATCCCCTGCTGATAAAACAGCACTGGCACCACTTCGCGGGGAtgg                                              >  3:367641/1‑151 (MQ=255)
                                                                                                          aaagaaaagaGCGACCATTTTGTGGCGGTGGGGCGCGGGAAAAATTACACACGGGAGCCCCCCCATCAGAAAACGTGGGGGGGACACCCAGCCGCGGAGCGATCATCCCCTGCTGATAAAACAGCACCGGCCCCCCCTCGCgggggtgg                                              >  2:1639839/1‑144 (MQ=255)
                                                                                                                aagaGCGTCCATTTTGTGGCGGTGACGCGCGGAAAAATTCACAGACGGTAGCCCCCACAACAGATAGCGTGGTGCGTACACCGAGCTGCGTAGCGGTCATCCCTTGCTGATAAACCAGCCCTGGCACCACTTCGCGGGTATGGTGGCtgtg                                      <  1:1324245/151‑1 (MQ=255)
                                                                                                                   agCGACCATTTTGTGGCGGTGGCGCGCGGGAAAATTCACACACGGTAGCCCCCACATCAGAAAACGGGGGGGGCACACCCAGCTGCGTAGCGATCATCCCCTTCTGATAAAACAGCACCGGCACCACCTCGCGGGGATGGGGGGTGTGAcc                                   >  2:919027/1‑151 (MQ=255)
                                                                                                                        ccATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATgg                              <  2:21683/151‑1 (MQ=255)
                                                                                                                             ttGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCgg                          <  1:809487/150‑1 (MQ=255)
                                                                                                                                 ggcggtggcgCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGt                     <  2:801909/151‑1 (MQ=255)
                                                                                                                                    ggtggcgCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGTTCGGATCGTTg                   >  1:1228002/1‑150 (MQ=255)
                                                                                                                                    ggtggcgCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGc                  >  1:505075/1‑151 (MQ=255)
                                                                                                                                      tggcgCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGGATGGTGGCTGTGACCAATGGTCGGATCGTTGcc                 >  2:141287/1‑150 (MQ=255)
                                                                                                                                       ggcgCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGt               <  3:567000/151‑1 (MQ=255)
                                                                                                                                       ggcgCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGt               <  2:1294628/151‑1 (MQ=255)
                                                                                                                                        gcgcgcGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTg              <  1:1090072/151‑1 (MQ=255)
                                                                                                                                        gcgcgcGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTg              <  2:805677/151‑1 (MQ=255)
                                                                                                                                        gcgcgcGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTg              >  2:63736/1‑151 (MQ=255)
                                                                                                                                           cgcgGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTGATc           <  2:1624646/151‑1 (MQ=255)
                                                                                                                                                ggAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTGATCCg         >  4:467718/1‑150 (MQ=255)
                                                                                                                                                 gAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTGATCCGcc       >  2:589534/1‑151 (MQ=255)
                                                                                                                                                 gAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTGATCCGcc       <  1:194930/151‑1 (MQ=255)
                                                                                                                                                  aaaaaTTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAACCCAGCACTGGCACCTCTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTGATCCGCCa      <  1:346373/151‑1 (MQ=255)
                                                                                                                                                  aaaaaTTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTGATCCGCCa      <  2:636308/151‑1 (MQ=255)
                                                                                                                                                   aaaaTTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTGATCCGCCAt     <  1:838134/151‑1 (MQ=255)
                                                                                                                                                    aaaTTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTGATCCGCCATc    >  2:173088/1‑151 (MQ=255)
                                                                                                                                                     aaTTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTGATCCGCCATc    >  1:4337/1‑150 (MQ=255)
                                                                                                                                                      aTTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTGATCCGCCATCAc  <  2:810068/151‑1 (MQ=255)
                                                                                                                                                       |                                                                                                                                                     
CACGTTGTTTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTGATCCGCCATCAC  >  NC_000913/3508513‑3508811

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