Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,413,7350AT40.9% 173.4 / 64.5 110Y2* (TAT→TAAydaCDUF1187 family protein, Rac prophage; putative double‑strand break reduction protein
Reads supporting (aligned to +/- strand):  ref base A (9/56);  new base T (45/0);  total (54/56)
Fisher's exact test for biased strand distribution p-value = 3.30e-22
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.87e-01

Only 100 of 113 total aligned reads displayed.

TTACTTTCTGCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCATAACCGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGAAATGGCTTTCACCGTTACTGAT  >  NC_000913/1413593‑1413883
                                                                                                                                               |                                                                                                                                                       
ttACTTTCTGCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTCATTTTATTCa                                                                                                                                                     >  2:644279/1‑147 (MQ=255)
 tACTTTCTGCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGaa                                                                                                                                                 <  1:1759294/150‑1 (MQ=255)
 tACTTTCTGCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGaa                                                                                                                                                 <  1:285442/150‑1 (MQ=255)
 tACTTTCTGCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGaa                                                                                                                                                 <  1:434198/150‑1 (MQ=255)
    tttCTGCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTATATTTTTTAATGAac                                                                                                                                                >  2:201347/1‑148 (MQ=255)
    tttCTGCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTATATTAACAcc                                                                                                                                            >  2:19653/1‑151 (MQ=255)
     ttCTGCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTACATGAACAcc                                                                                                                                            >  2:1413888/1‑150 (MQ=255)
     ttCTGCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAacac                                                                                                                                             <  1:1886361/149‑1 (MQ=255)
      tCTGCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTTAATGAaca                                                                                                                                              >  1:1808602/1‑147 (MQ=255)
         gCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCat                                                                                                                                       <  1:1547777/151‑1 (MQ=255)
           tctctGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTTTATGAacaacccc                                                                                                                                         >  1:1536055/1‑142 (MQ=255)
           tctctGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTTAATTAAAACCTCCctt                                                                                                                                      >  2:572324/1‑150 (MQ=255)
              ctGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCattat                                                                                                                                    <  1:1245222/149‑1 (MQ=255)
                ggAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCattattatt                                                                                                                                <  2:1194415/151‑1 (MQ=255)
                ggAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCattatt                                                                                                                                   <  1:39794/148‑1 (MQ=255)
                 gAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTAttt                                                                                                                               <  1:1430519/151‑1 (MQ=255)
                   aaCGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTATATGAACACCCCCCTTAtttttt                                                                                                                               >  2:1831534/1‑149 (MQ=255)
                    aCGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTTAATGAACACCCCCCTTATTAttt                                                                                                                               >  2:1040364/1‑148 (MQ=255)
                    aCGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCa                                                                                                                            <  2:1794616/151‑1 (MQ=255)
                      gCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGt                                                                                                                          <  2:254687/151‑1 (MQ=255)
                       ccTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTAAATGAAAACCCCCCTTTTTTTTTTCAGTg                                                                                                                         >  2:1718918/1‑151 (MQ=255)
                        cTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTTTATGAACACCCCCCTTTTTTTTTTCAATgg                                                                                                                        >  1:834361/1‑151 (MQ=255)
                          gCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTgg                                                                                                                        <  1:679166/149‑1 (MQ=255)
                             tcttTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTAAATGAACAACTCCCATAtttttttccgggggt                                                                                                                      >  2:410055/1‑139 (MQ=255)
                               ttttttACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTTCATTAAAAACCCCCTTTTTTTTTTCAGTgg                                                                                                                        >  1:1892957/1‑144 (MQ=255)
                               ttttttACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTATATGTACAACCCCCTTTTTATTTTTATTGGtttcgtt                                                                                                                 >  1:1103517/1‑146 (MQ=255)
                               ttttttACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAAAACCTCCATTATTTTTTTCAGTGgtttgttt                                                                                                                 >  2:1233483/1‑151 (MQ=255)
                                tttttACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTATATTAAAACCCCCCTTTTTTTTTTCTGTGGGTTGttt                                                                                                                 >  1:1256408/1‑150 (MQ=255)
                                tttttACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGttt                                                                                                                 <  1:1507427/150‑1 (MQ=255)
                                tttttACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGt                                                                                                                   <  1:1785497/148‑1 (MQ=255)
                                  tttACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTAt                                                                                                              <  1:970537/150‑1 (MQ=255)
                                   ttACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTACATTAAAAACCCCCTTTTTTTTTTCCGTGGTTCCtttttt                                                                                                             >  2:1489869/1‑150 (MQ=255)
                                      ccTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTTAATAAACACCCCCCTTTTTTTTTTCAGTGGTTTGtttttt                                                                                                             >  2:605825/1‑147 (MQ=255)
                                        tGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTAAAAGAACACCCCCCTTATTTTTTTCCGTGGTTCGTTTAtt                                                                                                             >  2:325976/1‑145 (MQ=255)
                                         gagaGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAGTTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCt                                                                                                       <  2:372414/150‑1 (MQ=255)
                                             gCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTTATTTTTTAAATGAAACCCCCCCTTATTTTTTTCAGTGGGTCTTTTTTTTCATCttt                                                                                                     >  2:1473191/1‑148 (MQ=255)
                                             gCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTAAATTTAAACCCCCATTTTTTTTTTCTGTGGGTTGtttttttcctttttttt                                                                                                  >  1:1526073/1‑140 (MQ=255)
                                               aTTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCTAGt                                                                                                <  1:1716018/151‑1 (MQ=255)
                                               aTTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGt                                                                                                <  2:1527817/151‑1 (MQ=255)
                                                ttttttCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTTAAAGAAAAACCCCATTATTTTTTTCAGTGGGTTGTTTTTTTCCTCTTTTTAGTg                                                                                               >  1:1555681/1‑151 (MQ=255)
                                                  ttttCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTTTATTAAAACCCCCCTTTTTTTTTTCAGTGGGTTGTTTTTTTTATCttt                                                                                                     >  1:1101927/1‑143 (MQ=255)
                                                  ttttCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTTCATTAAAACCCCCCTTTTTTTTTTCTCTGGTTTGTTTTTTTCATCTTTTGAGt                                                                                                >  1:574301/1‑148 (MQ=255)
                                                   tttCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTAAATGAACACCTCCATTATTTTTTTCAGTGGTTCGGTTTTTTTATCTTTTTAGTGctt                                                                                            >  1:525056/1‑151 (MQ=255)
                                                   tttCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGct                                                                                             <  2:1268397/150‑1 (MQ=255)
                                                   tttCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGGAAATTTTTTAATAAAAACCTCCATTATTTTTTTCAGTTGTTCGTTTTTTTCATCTTTcttatggtt                                                                                            >  2:1596252/1‑143 (MQ=255)
                                                      cGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGcttctt                                                                                         <  2:281194/151‑1 (MQ=255)
                                                       gCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATTCATGCACACCTCCATTATTTTTTTCAGTGGGTTGTTTTTTTCATCTTTTGTGTGCTTCttt                                                                                        >  2:1872099/1‑151 (MQ=255)
                                                        cATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCtttt                                                                                       <  2:689926/151‑1 (MQ=255)
                                                         aTTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTAGTTCGTTTATTCCATCTTTCGAGTGCTTCttttt                                                                                      <  1:822348/151‑1 (MQ=255)
                                                           tCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTc                                                                                     <  2:1790210/150‑1 (MQ=255)
                                                              gATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCAc                                                                                   <  2:1365316/149‑1 (MQ=255)
                                                               aTTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTAAATGAACACCCCCCTTTTTATTTTCCGTGGGTCGTTTTTTTCATCTTTCCAGTGGTTTTTTTTTACtt                                                                                 >  2:1475833/1‑150 (MQ=255)
                                                                 ttCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTAAATGAACACCCCCCTTATTTTTTTCAGTGGGTCGGTTTTTTCATCTTTTTAGGGCTTTTTTTTTACtt                                                                                 >  2:1791936/1‑148 (MQ=255)
                                                                   cGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTTTATTAACACCCCCCTTCTTTTTTTCAGTGGGTTGTTTTTTTCTTCTTTTTAGTGCTTCtttttt                                                                                     >  1:1309615/1‑141 (MQ=255)
                                                                   cGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTTAATTAAAACCCCCCTTATTTTTTTCAGTGGTTTGTTTTTTTTATCTTTTTAGTGGTTCTTTTTTTCTTccacc                                                                            >  1:296921/1‑151 (MQ=255)
                                                                    gTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTccacca                                                                           <  1:1132971/151‑1 (MQ=255)
                                                                       aGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCATaa                                                                       <  2:56044/151‑1 (MQ=255)
                                                                       aGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCATaa                                                                       <  1:652521/151‑1 (MQ=255)
                                                                       aGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCATa                                                                        <  1:657637/150‑1 (MQ=255)
                                                                         tttAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTACATAAAAACCACCATTTTTTTTTTCACAGGTTCGTTTTTTTCATCTTTTTAGTGCttttttttttctt                                                                                 >  2:1388924/1‑135 (MQ=255)
                                                                          ttAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTAAAATAAAACCCCCCTTATTTTTTTCAATGGGTTGTTTTTTTCATCTTTTTAGTGGTTTTTTTTTTCTTTCACCATaaac                                                                     >  1:351889/1‑148 (MQ=255)
                                                                          ttAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCATAAcc                                                                     <  1:1904946/150‑1 (MQ=255)
                                                                           tAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTATATGAAAACCCCCCCTATTTTTTTCAGTGGGTCGTTTTTTTCATCTTTTTAGTGCTTCTTTTTTTCTTtcaccccaaaccg                                                                    >  1:673321/1‑143 (MQ=255)
                                                                              aTTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCATAACCGGTTc                                                                <  2:1799003/151‑1 (MQ=255)
                                                                              aTTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCATAACCGGTTc                                                                <  2:164108/151‑1 (MQ=255)
                                                                                tttGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTTCAAGAACACCTCCCTTTTTTTTTTCCGTGGTTTGTTTTTTTCCTCTTTTGAGGGGTTCTTTTTTTCTTCCACCCTAAACCGGt                                                                  >  2:481169/1‑148 (MQ=255)
                                                                                  tGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTAAATGAAAACCTCCCTTATTTTTTTCACTGGGTCGTTTTTTTCCTCTTTCTCGTGCTTTTTTTTTACTTTCCCCATAACCCGTTCtgg                                                             >  1:427031/1‑150 (MQ=255)
                                                                                  tGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCATAACCGGTTCggg                                                             <  2:493443/150‑1 (MQ=255)
                                                                                   gAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCATAACCGGTTCGGGAt                                                           <  2:629886/151‑1 (MQ=255)
                                                                                    aaTATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCATAACCGGTTCGGGAtt                                                          <  2:1017961/151‑1 (MQ=255)
                                                                                    aaTATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAAATTTTTTCATGAAAAACTCCCTTATTTTTTTCAGTGGGTCGTTTTTTTCCTCTTTTGGGGTGTTTTTTTTTACTTTCACCCTAAACGGTTCgggggt                                                          >  2:196967/1‑148 (MQ=255)
                                                                                     atatCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTTAATAAAAACACCCCTTTTTTTTTTCAGAGGGTCGTTTTTTTCCTCTTTTTTGTGCTTCTTTTTTTTTTTCACCCTAACCGGGTCGGGCtt                                                          >  1:1136077/1‑150 (MQ=255)
                                                                                     atatCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTACATGAACAACCCCCTTATTTTTTTCAGTGGGTCGTTTTTTTCATTTTTTTGGTGGTTTTTTTTTCCTTCCCCCCTAAACGGTTCggg                                                             >  2:679179/1‑147 (MQ=255)
                                                                                      tatCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTACAACAAAAACCCCCTTTTTTTTTTCCGTGGGTTGCTTTTTTCCTCTTTCGCGAGGCTTTTTTTTACTTTCCCCCTAAACCGTTCgggg                                                            >  1:375801/1‑146 (MQ=255)
                                                                                           tGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTTAATGAAAACCCCCCTTTTTTTTTTCCGTGGGTTGTTTTTTTCATTTTTTGTGGGCTTCTTTTTTACTTCCCCCCTAACCCGGTCGGGGTTCTCATcc                                                   >  2:1521616/1‑151 (MQ=255)
                                                                                           tGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCATAACCGGTTCGGGATTCGCATc                                                    <  1:476602/150‑1 (MQ=255)
                                                                                               cAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTATATGTACAACCCCCCTTTTTTTTTCTGGGGGTTGTTTTTTTCCTTTTTCCCTTGCTTTTTTTTTTTTTCCCCCATTACCCGTGCGGGGtttggttcccgt                                                >  1:853693/1‑140 (MQ=255)
                                                                                               cAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACAATAACACCCCCATTATTATTTTCAGTGGTTCGTTTATTTCATCTTTCGTGTGCTTTTTTTTTACTTCCCCCACAAACGGGTCGGGGTTCGCATcc                                                   >  2:1586467/1‑147 (MQ=255)
                                                                                                    agtagGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTAAATGAACACCTCCCCTTTTATTTTCCATGGGTTGTTTTTTTTCTCTTTTGTGTTGTTCTTTTTTTTTTTCCCCCTAAACGGTTCGGGGtttggctccctt                                                >  1:41123/1‑135 (MQ=255)
                                                                                                       agGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTACATTAACAACACCACTATTTTTTTCAGGGGGTTGTTTATTTCAATCTTTTAGGGGTTTTTTTTTACTTTCTCCACAACCAGGTCGGGGTTCGCAACCGTTaccccttt                                        >  2:852627/1‑144 (MQ=255)
                                                                                                         gAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCATAACCGGTTCGGGATTCGCATCCGTTACACTTCGCTc                                     <  2:418515/151‑1 (MQ=255)
                                                                                                          aGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTAAATGAAAACCCCCCTTATTTTTTTCCGTGGGTTGTTTTTTTCCTCTTTTTAGTGCTTCTTTTTTTCTTTCCCCCTAAACCGTTTGGGGTTCGCCTCCGTTTCACTtctctc                                     >  1:446870/1‑150 (MQ=255)
                                                                                                           gTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCATAACCGGTTCGGGATTCGCATCCGTTACACTTCGCTCgg                                   <  2:462361/151‑1 (MQ=255)
                                                                                                                    tccttTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCATAACCGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATc                            <  1:56290/149‑1 (MQ=255)
                                                                                                                     ccttTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCATAACCGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTg                         <  1:572223/151‑1 (MQ=255)
                                                                                                                         ttCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCATAACCGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGaaa                      <  2:335376/150‑1 (MQ=255)
                                                                                                                          tCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCATAACCGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGAAAt                     <  1:1027674/150‑1 (MQ=255)
                                                                                                                            aaTAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCATAACCGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGAAATGGc                  <  2:1753606/151‑1 (MQ=255)
                                                                                                                             aTAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCATAACCGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGAAATGGCt                 <  2:286503/151‑1 (MQ=255)
                                                                                                                               aGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGTGTGCTTCTTTTTCACTTCCACCATAACCGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGAAATGGCttt               <  2:480787/151‑1 (MQ=255)
                                                                                                                               aGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCATAACCGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGAAATGGCtt                <  2:875611/150‑1 (MQ=255)
                                                                                                                                  ggCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCATAACCGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGAAATGGCTTTCa             <  1:1181801/150‑1 (MQ=255)
                                                                                                                                  ggCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCATAACCGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGAAATGGCTTTCAc            >  2:1107705/1‑151 (MQ=255)
                                                                                                                                       tAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCATAACCGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGAAATGGCTTTCACCGtt        >  1:255976/1‑150 (MQ=255)
                                                                                                                                        aaTTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCATAACCGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGAAATGGCTTTCACCGTTa       <  1:1763435/150‑1 (MQ=255)
                                                                                                                                          aTTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCATAACCGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGAAATGGCTTTCACCGTTACt     <  2:1895689/151‑1 (MQ=255)
                                                                                                                                            tttATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCATAACCGGTTCGGGATTCGCCTCCGTTACACTTCGCTCGGTAATATCCTGAAATGGCTTTCACCGTTACt     >  2:283980/1‑149 (MQ=255)
                                                                                                                                            tttATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCATAACCGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGAAATGGCTTTCACCGTTACTg    >  1:1376791/1‑150 (MQ=255)
                                                                                                                                              tataCATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCATAACCGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGAAATGGCTTTCACCGTTACTGAt  <  2:350963/150‑1 (MQ=255)
                                                                                                                                              tataCATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCATAACCGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGAAATGGCTTTCACCGTTACTGAt  <  2:1187743/150‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                                       
TTACTTTCTGCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCATAACCGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGAAATGGCTTTCACCGTTACTGAT  >  NC_000913/1413593‑1413883

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