Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,413,7350AT45.6% 102.3 / 51.2 79Y2* (TAT→TAAydaCDUF1187 family protein, Rac prophage; putative double‑strand break reduction protein
Reads supporting (aligned to +/- strand):  ref base A (4/39);  new base T (36/0);  total (40/39)
Fisher's exact test for biased strand distribution p-value = 2.53e-18
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGTTTTACTTTCTGCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCATAACCGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGAAATGGCTTTCACCGTTACTGATGG  >  NC_000913/1413589‑1413885
                                                                                                                                                     |                                                                                                                                                           
aGTTTTACTTTCTGCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGTGGTAATTTTATACa                                                                                                                                                        <  2:1010138/151‑1 (MQ=255)
    ttACTTTCTGCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTatatta                                                                                                                                                     >  2:201120/1‑148 (MQ=255)
      aCTTTCTGCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAac                                                                                                                                                   <  2:788984/150‑1 (MQ=255)
        tttCTGCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAAATTTTTATATTAaca                                                                                                                                                  >  1:781814/1‑149 (MQ=255)
         ttCTGCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAacac                                                                                                                                                 <  2:870517/149‑1 (MQ=255)
         ttCTGCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAacac                                                                                                                                                 <  2:739716/149‑1 (MQ=255)
         ttCTGCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACAcc                                                                                                                                                <  2:1255193/150‑1 (MQ=255)
           cTGCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTTAATTTTTTAATGaaa                                                                                                                                                   >  1:962464/1‑144 (MQ=255)
            tGCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTAAATGAAAACcccc                                                                                                                                             >  2:758263/1‑150 (MQ=255)
             gCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTCAATTTTTTCATGAaaac                                                                                                                                                 >  2:876186/1‑145 (MQ=255)
               tctctGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTTAATTAACACCCCCctt                                                                                                                                          >  2:248528/1‑150 (MQ=255)
                ctctGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACGCCTCCatt                                                                                                                                          <  1:1477826/149‑1 (MQ=255)
                 tctGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTCATTTTTTATATTaaaaacccccctatt                                                                                                                                       >  2:1372603/1‑138 (MQ=255)
                    ggAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCattattatt                                                                                                                                    <  2:373482/151‑1 (MQ=255)
                    ggAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCattattatt                                                                                                                                    <  2:1448890/151‑1 (MQ=255)
                    ggAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCattatt                                                                                                                                       <  1:538254/148‑1 (MQ=255)
                     gAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTATATTAAAACCCCCCTtttttttt                                                                                                                                   >  1:94546/1‑151 (MQ=255)
                      aaaCGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTGAATTTTTTAAATAAAAACCCCCtttttttttt                                                                                                                                   >  1:126818/1‑151 (MQ=255)
                          gCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTAAATATAAAACCCCACTATTATTTTCac                                                                                                                               >  2:2912/1‑149 (MQ=255)
                           ccTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTTAATTTTTTAATGAACAACTCCCTttttttttt                                                                                                                                  >  2:783051/1‑145 (MQ=255)
                               cttcttTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGtt                                                                                                                          <  2:288629/150‑1 (MQ=255)
                                ttcttTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTATATTAAAACCTCCCTTATTTTTTTTTGTGGtt                                                                                                                          >  2:480420/1‑149 (MQ=255)
                                    tttttACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGt                                                                                                                       <  2:864518/148‑1 (MQ=255)
                                    tttttACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGGAATTTTTTAAATATAAACCCCCATTTTTTTTTTCAATgggttggtt                                                                                                                     >  2:1164520/1‑143 (MQ=255)
                                     ttttACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTCAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTTATTTTTTACATGAAAAACCCCCTTATTATTTTTTGTGGTTCGttttt                                                                                                                   >  1:153143/1‑149 (MQ=255)
                                          ccTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTTCATTAACACCCCCCTTATTATTTTCACTGGTTTGttttttt                                                                                                                >  2:1134123/1‑147 (MQ=255)
                                             gagaGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTAAAAGAACACCTCCCCTTTTTTTTTCCGTGGttttttttttttctc                                                                                                             >  2:37629/1‑136 (MQ=255)
                                              agagCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCttt                                                                                                         <  1:835333/151‑1 (MQ=255)
                                              agagCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCttt                                                                                                         <  1:1379626/151‑1 (MQ=255)
                                                agCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGAAATTTTTTTTATTGAGACCCCCCTTATTATTTCCTGCGGATTGttttttttttcttt                                                                                                         >  1:625674/1‑141 (MQ=255)
                                                 gCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGa                                                                                                      <  2:450295/151‑1 (MQ=255)
                                                 gCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGa                                                                                                      <  1:1097789/151‑1 (MQ=255)
                                                  cATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTTATTTTTTATAATAAAAACCCCCTTTTTTTTTTCCTTGGGTTGTTTTTTTCATCtttt                                                                                                        >  1:1529662/1‑147 (MQ=255)
                                                    ttttttCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTTCAATAAAAACCCCCTTATTTTTTTCAGCGGGTCTTTTTTTTCTTCtttttt                                                                                                      >  1:687993/1‑145 (MQ=255)
                                                        ttCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTTTATAAAAACCTCCCTTATTATTTTCACTGGGTTGTTTTTTTCCTTTTTTTAGTGctt                                                                                                >  2:1228472/1‑150 (MQ=255)
                                                        ttCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTACATTAACACCCCCCTTTTTTTTTTCAGTGGTTTGTTTTTTTCATTTTTTTTGTGctt                                                                                                >  2:91238/1‑150 (MQ=255)
                                                            cATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCtttt                                                                                           <  2:120028/151‑1 (MQ=255)
                                                             aTTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCttttt                                                                                          <  2:1301358/151‑1 (MQ=255)
                                                                  gATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCa                                                                                        <  2:662478/148‑1 (MQ=255)
                                                                  gATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACt                                                                                     <  1:15486/150‑1 (MQ=255)
                                                                    tttCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTCATTTTTTTTATAAAAAACACCCTTTTTTTTTTCACTGGGTTGTTTTTTTCATCTTTTTCGTTCtttttttt                                                                                          >  1:1221614/1‑144 (MQ=255)
                                                                      tCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTTATTTTTTATATTAAAAACCCCCCTTTTTTTTTCAGTGGGTTGTTTTTTTTTTCTTTTTAGTGGTTCttttt                                                                                          >  1:325804/1‑142 (MQ=255)
                                                                      tCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTcca                                                                                <  2:1100960/150‑1 (MQ=255)
                                                                       cGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTccac                                                                              <  2:1122478/150‑1 (MQ=255)
                                                                        gTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTCAATTTTTATAATAACACCCCCACTTTTATTTTCAGGGGGTTGTTTTTTTCCTCCTTTGGGTGGTTTCTTTTTACTTtcac                                                                              >  1:1318917/1‑149 (MQ=255)
                                                                         ttAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTAAAATAAAAACCCCCTTTTTTTTTTCAGTGGGTTGTTTTTTTCATCTTTTTAGTGGttttttttttctt                                                                                     >  2:115223/1‑139 (MQ=255)
                                                                          tAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTATAATAAAACCCCCCATATTTTTTTCATTGGTTTGGTTTTTTTCTCTTTTGAGTGGTTCTTTTTTTCTTccccc                                                                             >  2:660738/1‑148 (MQ=255)
                                                                          tAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTACAAGAACACCCCCCTTATTATTTTCAGTGGTTCGTTTTTTTCCTCTTTTCAGTGCTTTTTTTTTTCTTccccc                                                                             >  2:1299109/1‑148 (MQ=255)
                                                                          tAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCAt                                                                           <  1:667374/150‑1 (MQ=255)
                                                                           aGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCATa                                                                          <  2:475124/150‑1 (MQ=255)
                                                                                aGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCATAACCgg                                                                     <  2:1547508/150‑1 (MQ=255)
                                                                                 gATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCATAACCGGtt                                                                   <  1:816351/151‑1 (MQ=255)
                                                                                 gATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGGAATTTTTTTAAATAAAAACCCCACTATTTTTTTCCGTGGGTCGGTTATTTCCTTTTTTCGGGGCTTCTTTTTTTCTTTCACCACAAACgggt                                                                   >  2:1331843/1‑149 (MQ=255)
                                                                                    tttGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTAAATGAACACCCCCCTTATTATTTTCCGCGGGTTGTTTTTTTTCTCTTTTTAGTGCTTCTTTTTTTCTTCCACCCTAAACgggt                                                                   >  2:468136/1‑146 (MQ=255)
                                                                                          tatCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTTATTTTTTAAATGAACACCTCCCTTATTATTTTCATTGGGTTGTTTATTTCATCTTTTTTATGCTTCTTTTTTACTTCCACCCTAACCGGtgtggggtt                                                            >  1:163972/1‑141 (MQ=255)
                                                                                            tCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAAATTTTTTCATAAACAACTCCCTTATTTTTTTCAATGGTTCGGTTATTTCTTCTTTTTAGTGGTTTTTTTTTACTTACCCCCCAACCCGTCCGGGGTTTGcc                                                        >  1:1093219/1‑150 (MQ=255)
                                                                                              ttGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAAATTTTTTAATTAAAACCCCCCCTTTTTTTTTCCCTGGGTTGTTTTTTTCATTTTTTGAGTGGTTCTTTTTTTCTTTCCCCATAACCCGGTTGGGGTTctcctc                                                      >  1:698133/1‑146 (MQ=255)
                                                                                              ttGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAAATTTTTAAATAAAAAACCCCCCTATTTTTTTCAGTGGGTCGTTTTTTTCATCTTTTTAGTGGTTCTTTTTTTCTTTCCCCCTAAACCGGTTGGGGtttgcat                                                        >  1:627061/1‑146 (MQ=255)
                                                                                                gTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTATATTAAAAACTCCCCTATTATTTTCAATGGTTTGTTTTTTTCATCTTTTTAGGGCTTTTTTTTTTCTTTCACCCAAAACCGTTTGGGGGttgcat                                                        >  2:1359372/1‑144 (MQ=255)
                                                                                                gTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCATAACCGGTTCGGGATTCGCATCCg                                                    <  1:1420461/151‑1 (MQ=255)
                                                                                                gTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCATAACCGGTTCGGGATTCGCATCCg                                                    <  1:444348/151‑1 (MQ=255)
                                                                                                      ttAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTTAATTAAAACCCCCCCTTTTTTTTTTATTGGGTCGTTTTTTTTTTCTTTTTTGTGCTTTTTTTTTTTTTTCCCCATAAACCGTTTgggggtttcatccctt                                                  >  1:117196/1‑134 (MQ=255)
                                                                                                      ttAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAAATTTTTAAATAAAAACCCCCCTTTTTTTTTTCTATGGGTTGTTTTTTTCCTCTTTTTTGTGGTTTTTTTTTTCTTCCCCCCAAAACCGGTCGGGGTTCGCATccccttc                                                >  2:429532/1‑144 (MQ=255)
                                                                                                       tagtagGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCATAACCGGTTCGGGATTCGCATCCGTTacac                                              <  2:786036/150‑1 (MQ=255)
                                                                                                        agtagGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCATAACCGGTTCGGGATTCGCATCCGTTACACt                                             <  1:1027918/150‑1 (MQ=255)
                                                                                                         gtagGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCATAACCGGTTCGGGATTCGCATCCGTTACACTTc                                           <  2:42165/151‑1 (MQ=255)
                                                                                                                     ccttccttTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCATAACCGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAatat                               <  1:843620/151‑1 (MQ=255)
                                                                                                                         ccttTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCATAACCGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCt                            <  2:396864/150‑1 (MQ=255)
                                                                                                                            tttCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCATAACCGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGaa                         <  2:1321662/150‑1 (MQ=255)
                                                                                                                              tCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCATAACCGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGAAAt                       <  2:962567/150‑1 (MQ=255)
                                                                                                                              tCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCATAACCGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGAAAt                       <  2:1292223/150‑1 (MQ=255)
                                                                                                                               cAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCATAACCGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGAAATg                      <  1:351132/150‑1 (MQ=255)
                                                                                                                                     tGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCATAACCGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGAAATGGCTTTc                >  2:44446/1‑150 (MQ=255)
                                                                                                                                      ggCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCATAACCGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGAAATGGCTTTCAc              >  2:694659/1‑151 (MQ=255)
                                                                                                                                         ggTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCATAACCGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGAAATGGCTTTCACCGt           <  1:339441/151‑1 (MQ=255)
                                                                                                                                         ggTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCATAACCGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGAAATGGCTTTCACCGt           <  1:248498/151‑1 (MQ=255)
                                                                                                                                          gTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCATAACCGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGAAATGGCTTTCACCGtt          >  2:818717/1‑151 (MQ=255)
                                                                                                                                              aaTTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCATAACCGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGAAATGGCTTTCACCGTTa         <  2:94557/150‑1 (MQ=255)
                                                                                                                                                 ttttATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCATAACCGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGAAATGGCTATCACCGTTACTg      >  1:1183527/1‑151 (MQ=255)
                                                                                                                                                    tataCATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCATAACCGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGAAATGGCTTTCACCGTTACTGAt    <  2:1375810/150‑1 (MQ=255)
                                                                                                                                                     ataCATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCATAACCGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGAAATGGCTTTCACCGTTACTGATgg  <  2:471396/151‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                           
AGTTTTACTTTCTGCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCATAACCGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGAAATGGCTTTCACCGTTACTGATGG  >  NC_000913/1413589‑1413885

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: