Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,413,7350AT46.9% 63.7 / 50.4 64Y2* (TAT→TAAydaCDUF1187 family protein, Rac prophage; putative double‑strand break reduction protein
Reads supporting (aligned to +/- strand):  ref base A (8/26);  new base T (30/0);  total (38/26)
Fisher's exact test for biased strand distribution p-value = 3.83e-11
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTTTTACTTTCTGCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCATAACCGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGAAATGGCTTTCACCGTTACTG  >  NC_000913/1413590‑1413881
                                                                                                                                                  |                                                                                                                                                      
gTTTTACTTTCTGCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACAt                                                                                                                                                   <  2:546216/151‑1 (MQ=255)
gTTTTACTTTCTGCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACAt                                                                                                                                                   <  1:978200/151‑1 (MQ=255)
   ttACTTTCTGCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAAATTTTTACa                                                                                                                                                    >  2:838230/1‑147 (MQ=255)
    tACTTTCTGCTCTCTGGAAACGCCTGCTTTTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGaa                                                                                                                                                <  1:1184194/150‑1 (MQ=255)
    tACTTTCTGCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGaa                                                                                                                                                <  2:849190/150‑1 (MQ=255)
      cTTTCTGCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAaca                                                                                                                                              <  1:1003808/150‑1 (MQ=255)
            gCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCa                                                                                                                                        <  1:615615/150‑1 (MQ=255)
                tctGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTTCATGAACACCTcccatttt                                                                                                                                   >  2:295108/1‑145 (MQ=255)
                 ctGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTATATTAAAaccacccctttt                                                                                                                                   >  1:797824/1‑144 (MQ=255)
                           cTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGGAAATTTTTAAATTAAAACCCCCCTTATTTTTTTCCGTgg                                                                                                                      >  2:1032269/1‑151 (MQ=255)
                              cttcttTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAAATTTTTTAATGAAAACCCCCCTTTTTTTTTTCAGTgggt                                                                                                                    >  2:238879/1‑148 (MQ=255)
                                    ttttACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTTTAATAACACCTCCCTTTTTTTTTTTATTTGTTCGttttt                                                                                                             >  1:562289/1‑149 (MQ=255)
                                     tttACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTAt                                                                                                             <  2:109606/150‑1 (MQ=255)
                                      ttACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTTAATTTTTTAATTTAAAACTCCCTTTTTTTTTTCTGTGGTTTGttttttt                                                                                                           >  1:411641/1‑150 (MQ=255)
                                      ttACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACAAGAAAACCCCCCTTATTATTTTCAGTGGGTCGtttttt                                                                                                            >  1:135293/1‑150 (MQ=255)
                                       tACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTTATTTTTTAAAAGAACACCTCCACTTTTATTTTCAGAGGGTTGttttttt                                                                                                           >  1:368267/1‑149 (MQ=255)
                                       tACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTc                                                                                                           <  2:772727/150‑1 (MQ=255)
                                          cTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTTCATTAAAACCTCCCTTTTTTTTTTCCGTGGTTCGtttttt                                                                                                            >  1:995860/1‑146 (MQ=255)
                                           tGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTTAATTTTTAAATTAAAGCACCCCTTATTCTTTTCACTGGGTGGtttttttcttat                                                                                                      >  2:170026/1‑145 (MQ=255)
                                            gagaGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTAAATTAACACCCCCCTTTTTTTTTTCCGTGGTTTGTTTTTTTCCTCtt                                                                                                     >  2:974978/1‑151 (MQ=255)
                                            gagaGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCt                                                                                                      <  2:575882/150‑1 (MQ=255)
                                              gagCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTTTAATAAAACCCCCATTATTATTTTCAATGGTTTTTTTTTTTCATCttt                                                                                                    >  2:854869/1‑150 (MQ=255)
                                              gagCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTc                                                                                                   <  2:1161390/151‑1 (MQ=255)
                                                gCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTAAATGAAAACCCCCCTTATTTTTTTCTTTGGGTCGTTTTTTTCATCttt                                                                                                    >  2:667299/1‑148 (MQ=255)
                                                 cATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTAAAAGAAAACCCCCCTTATTATTTTCCGTGGGTTGTTTTTTTCATCttt                                                                                                    >  2:1039854/1‑147 (MQ=255)
                                                  aTTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAAATTTTTTAAAGAAAACCCCCCCTTTTTTTTTCAGTGGTTTGGTTTTTTCATCttt                                                                                                    >  2:262391/1‑146 (MQ=255)
                                                    tttttCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTg                                                                                             <  1:770550/150‑1 (MQ=255)
                                                         cGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGcttctt                                                                                       <  2:534053/151‑1 (MQ=255)
                                                         cGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGcttctt                                                                                       <  2:1112592/151‑1 (MQ=255)
                                                          gCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTATATTAACACCTCCCTTATTTTTTCTAGTGGGTCGTTTTTTTCTTCTTTCGATTGCtttttt                                                                                      >  1:573734/1‑151 (MQ=255)
                                                          gCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCttt                                                                                      <  1:639784/151‑1 (MQ=255)
                                                            aTTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATTAAAACATCTCTTATTTTTTTTAGTGGTTGGTTTTTTTCCTCTTTTCTGTGCtttttttt                                                                                    >  1:700006/1‑151 (MQ=255)
                                                                 gATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACtt                                                                               <  1:664632/151‑1 (MQ=255)
                                                                      cGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTccac                                                                           <  1:1200995/150‑1 (MQ=255)
                                                                            tttAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATATAAGAACAACTCCCTTATTTTTTTCAGTGGTTTGTTTTTTTCCTCTTTCTGGTGCTTCTTTTTTACTTTCACCCTaa                                                                     >  2:830791/1‑149 (MQ=255)
                                                                                gATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTATATTAACACCCCCCCTTTTTTTTTCCGTGGTTTGTTTTTTTCATCTTTTTTGTGCTTCTTTTTTTCTTTCACCCCAACCCGtt                                                               >  1:1217270/1‑151 (MQ=255)
                                                                                gATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTATATGAAAACCTTCCTTTTTTTTTTTTGTGGGTTGTTTTTTTTCATTTTCGTGTGCTTTTTTTTTTCTTCTCCCATAAcccggt                                                               >  1:971459/1‑147 (MQ=255)
                                                                                gATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCATAACCGGtt                                                               <  2:387995/151‑1 (MQ=255)
                                                                                 aTTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCATAACCGGTTc                                                              <  1:51136/151‑1 (MQ=255)
                                                                                     tGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTATAAGAAAACCCCCCTTATTTTTTTCCATGGTGTGTTTTTTTTTTCTTTTTAGGGCTTTTTTTTTTCTTTCACCCCAAcccggttggg                                                           >  1:774591/1‑142 (MQ=255)
                                                                                      gAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAAATTTTTAAAAAAACAACCCCATTTTTTTTTTCAGGGGGTTGTTTTTTTCCTCTTTTGCGTGCTTCTTTTTTACTTTCCCCATAAACCGGTTggg                                                            >  2:359296/1‑149 (MQ=255)
                                                                                       aaTATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTTATTTTTTAAATTAAGACCCCCCTTATTTTTTTTCAGGGGTTGTTTTTTTATTTTTTTTAGCTCTTTTTTTTTTCTTTCACCATAACCgacttggg                                                           >  1:288722/1‑141 (MQ=255)
                                                                                         tatCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTAAATTAAAACCTCCCTTATTTTTTTCAGTGGTTTGGTTTTTTCATCTTTTTAGTGGTTCTTTTTTACTTTCACCATAAACCGTTTggg                                                           >  2:960819/1‑146 (MQ=255)
                                                                                          atCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTTTAAGTAAAACTCCATTTTTTTTTTCAATGGGTTGTTTTTTTCCTCTTTTGGGTGGTTTCTTTTTTACTTCACCCTAAACCGTTCgggg                                                          >  1:1209147/1‑145 (MQ=255)
                                                                                           tCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTACATTAAAACCCCCCTTTTTTTTTTCCGAGGGTCGTTTATTTCATCTTTTTTGTGCCTCTTTTTTTCTTTCCCCCTAAACAGGTCGGGGTTctcc                                                    >  1:966966/1‑148 (MQ=255)
                                                                                              tGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTACATTAACAACCCCCTTATTTTTTTCAGTGGTTTGTTTTTTTCCTCTTTTGAGTGCTTCTTTTTTACTTCCCCCCTAAACGGTTCGGGGTTTGCATcc                                                 >  2:1181344/1‑151 (MQ=255)
                                                                                               gTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGGAATTTTTTTCAAGAACACCCCCCTTTTTATTTTCAGGGGGTTGTTTTTTCCCTCTTTCGTGGGCTTCTTTTTTCCTTTCCCCATAACCGGCTCGGGGTTTGCATcc                                                 >  2:143109/1‑150 (MQ=255)
                                                                                                 ccAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCATAACCGGTTCGGGATTCGCATCCGtt                                              <  2:140363/151‑1 (MQ=255)
                                                                                                 ccAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCATAACCGGTTCGGGATTCGCATCCGt                                               <  2:1133648/150‑1 (MQ=255)
                                                                                                  cAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATATTTTAAAAGAAAACCTCCCCTATTTATTTTCGTTGGTTTTTCATTCAATCTTTCTAGTGCTTCTTTTTTTCTTTCTCCTTAAACGGACTGGGAtttcgattccttt                                             >  1:990021/1‑140 (MQ=255)
                                                                                                    gTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTTAATTAACACCCCCCTTTTTATTTTCAGTGGGTTGTTTTTTTCCTCTTTTTAGTGGTTTTTTTTTTCTTCCCCCCTAACCGGTTCgggggttgcctcc                                                 >  2:1003576/1‑135 (MQ=255)
                                                                                                      tagtagGAGTGCCACCTTCCTTTTCAATAGTGGCGGGAATTTTTTTAAAGAACACCCCCCTTTTTTTTTTTAGTGGGTCGTTTTTTTCATCTTTTTCATGGTTCTTTTTTTCTTTCCCCATTACCGGTTTGGGGtttttctctctt                                              >  1:61496/1‑136 (MQ=255)
                                                                                                                gCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCATAACCGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGt                                <  2:763840/150‑1 (MQ=255)
                                                                                                                gCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCATAACCGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTa                               <  2:417381/151‑1 (MQ=255)
                                                                                                                        ccttTTCAATAGTGGCGGTAAATTTTTAAATAAAAAACACCCTTTTTATTTTCAGTGGTTTGTTTTTTTCTTCTTTTTGCTAGTTTTTTTTTTATTTCACCACCACCGGTTCTGGGTTCGTACCCGTTTCACTTTTCTCCGTTatat                           >  1:792047/1‑147 (MQ=255)
                                                                                                                           tttCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCATAACCGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGa                      <  2:446567/149‑1 (MQ=255)
                                                                                                                              cAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCATAACCGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGAAATg                  <  1:860755/150‑1 (MQ=255)
                                                                                                                                  aGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCATAACCGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGAAATGGCttt             <  2:107111/151‑1 (MQ=255)
                                                                                                                                    tGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCATAACCGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGAAATGGCTTTc            >  2:1129775/1‑150 (MQ=255)
                                                                                                                                    tGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCATAACCGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGAAATGGCTTTc            >  1:338007/1‑150 (MQ=255)
                                                                                                                                        ggTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCATAACCGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGAAATGGCTTTCACCGt       >  2:258797/1‑151 (MQ=255)
                                                                                                                                             aTTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCATAACCGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGAAATGGCTTTCACCGTTACt   <  2:70703/151‑1 (MQ=255)
                                                                                                                                              ttttATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCATAACCGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGAAATGGCTTTCACCGTTACTg  >  2:1100270/1‑151 (MQ=255)
                                                                                                                                               tttATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCATAACCGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGAAATGGCTTTCACCGTTACTg  >  2:247546/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                      
GTTTTACTTTCTGCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCATAACCGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGAAATGGCTTTCACCGTTACTG  >  NC_000913/1413590‑1413881

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