Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,413,7350AT43.3% 130.2 / 61.2 90Y2* (TAT→TAAydaCDUF1187 family protein, Rac prophage; putative double‑strand break reduction protein
Reads supporting (aligned to +/- strand):  ref base A (2/49);  new base T (39/0);  total (41/49)
Fisher's exact test for biased strand distribution p-value = 1.75e-23
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTTACTTTCTGCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCATAACCGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGAAATGGCTTTCACCGTTACTGATGG  >  NC_000913/1413592‑1413885
                                                                                                                                                |                                                                                                                                                            
tttACTTTCTGCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATg                                                                                                                                                        <  1:1049261/150‑1 (MQ=255)
  tACTTTCTGCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTAAATGaa                                                                                                                                                      >  1:446328/1‑150 (MQ=255)
   aCTTTCTGCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAac                                                                                                                                                     <  1:552098/150‑1 (MQ=255)
       tCTGCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCt                                                                                                                                                <  1:614687/150‑1 (MQ=255)
        cTGCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTTATTTTTTATATTaaaaacccc                                                                                                                                               >  2:51705/1‑144 (MQ=255)
        cTGCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGCAAATTTTTATATGAACACCTcc                                                                                                                                              >  1:985852/1‑151 (MQ=255)
           ctctctGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTATATGAAAAcccc                                                                                                                                               >  2:593767/1‑145 (MQ=255)
            tctctGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTTCATGAACACCCCCatt                                                                                                                                           >  2:1334838/1‑150 (MQ=255)
              tctGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTAAATTAAAACcccccctttt                                                                                                                                        >  1:539891/1‑145 (MQ=255)
              tctGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAAATTTTTAAAAGAAAACCCCCcttatt                                                                                                                                       >  1:586099/1‑151 (MQ=255)
               ctGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTTAATGAACACCCCCcttatt                                                                                                                                       >  1:1355200/1‑150 (MQ=255)
               ctGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCatta                                                                                                                                          <  2:1364635/148‑1 (MQ=255)
                 ggAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTTTATTAAAACCCCCcttattatt                                                                                                                                   >  2:1347397/1‑151 (MQ=255)
                 ggAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTAAATTAACAACTCCcttattttt                                                                                                                                   >  1:672599/1‑151 (MQ=255)
                 ggAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCattattat                                                                                                                                    <  2:402192/150‑1 (MQ=255)
                 ggAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCattatt                                                                                                                                       <  2:832436/148‑1 (MQ=255)
                      cGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTTAATTTTTTCATTAATACcccccctttttt                                                                                                                                      >  2:561736/1‑137 (MQ=255)
                      cGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTTCATGAACACCCCCCCTTTTAttt                                                                                                                                  >  1:962519/1‑147 (MQ=255)
                       gCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGt                                                                                                                             <  1:470238/151‑1 (MQ=255)
                        ccTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTAAATGAAAAACCCCCTTTTTTTTTTCAGTg                                                                                                                           >  1:1002538/1‑151 (MQ=255)
                        ccTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTg                                                                                                                           <  2:468047/151‑1 (MQ=255)
                        ccTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTg                                                                                                                           <  1:45300/151‑1 (MQ=255)
                         cTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTACATGAAAAACTCCCTTATTATTTTCAGTgg                                                                                                                          >  1:1273937/1‑151 (MQ=255)
                         cTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTAAATGAAAACCCCCCTTTTTTTTTTCTGTgg                                                                                                                          >  1:1227904/1‑151 (MQ=255)
                            cttcttTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGtt                                                                                                                        <  2:148665/150‑1 (MQ=255)
                              tcttTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTTGATTACAACCTCCCTTATTTTTTTTCGTGGTTc                                                                                                                       >  2:510793/1‑149 (MQ=255)
                              tcttTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTc                                                                                                                       <  1:1209295/149‑1 (MQ=255)
                                 tttttACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTTAATGAAAACCCCCATTTTTTTTTTCCGTGGGTTGtttt                                                                                                                  >  1:1199113/1‑150 (MQ=255)
                                 tttttACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTACAATAATACCCCCCTTTTTTTTTTCAGTGGttgttttt                                                                                                                  >  1:109674/1‑145 (MQ=255)
                                 tttttACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGttt                                                                                                                   <  1:403066/150‑1 (MQ=255)
                                 tttttACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGttt                                                                                                                   <  2:925349/150‑1 (MQ=255)
                                    ttACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTAtt                                                                                                                <  1:668763/150‑1 (MQ=255)
                                     tACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTTAATTAAAACCTCCATTATTTTTTTCCTTGGGTCGttttttt                                                                                                               >  2:961810/1‑149 (MQ=255)
                                       ccTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCAt                                                                                                            <  1:387471/151‑1 (MQ=255)
                                       ccTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCAt                                                                                                            <  2:1132496/151‑1 (MQ=255)
                                       ccTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCAt                                                                                                            <  2:331448/151‑1 (MQ=255)
                                       ccTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAAATTTTTTCATTAACACCCCCCTTATTTTTTTCCGTGGTTCGttttttt                                                                                                               >  1:1027670/1‑147 (MQ=255)
                                         tGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGCAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTAAATGAAAAACCCCCTTTTTTTTTTCAGTGGGTCGTTTAtt                                                                                                                >  1:3107/1‑145 (MQ=255)
                                          gagaGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCtt                                                                                                         <  2:153804/151‑1 (MQ=255)
                                           agagCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTATATTAAAACCCCCCTTATTTTTTTTTTTGGTTTGTTTTTTTCATCttt                                                                                                        >  2:515865/1‑151 (MQ=255)
                                           agagCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCttt                                                                                                        <  1:788746/151‑1 (MQ=255)
                                            gagCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTc                                                                                                       <  2:351995/151‑1 (MQ=255)
                                               cATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTACAATTAAACCCCCCATATTATTTTCCATTGGTCGTTTATTTCATCtttt                                                                                                       >  1:144301/1‑147 (MQ=255)
                                                aTTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGt                                                                                                   <  2:687153/151‑1 (MQ=255)
                                                   ttttCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGc                                                                                                 <  1:1107194/150‑1 (MQ=255)
                                                     ttCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGctt                                                                                               <  1:1282206/150‑1 (MQ=255)
                                                       cGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTTAAAGAAAACCCCCCATTTTATTTTCAATGGGTTGTTTTTTTCCTCtttttggtggttttt                                                                                            >  2:1009254/1‑139 (MQ=255)
                                                       cGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGcttctt                                                                                            <  2:971498/151‑1 (MQ=255)
                                                         cATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCtttt                                                                                          <  1:7870/151‑1 (MQ=255)
                                                            tCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAAATTTTTATATTAAAACCCCCCCTTTTTTTTTCCGTGGGTTGGTTTTTTCCTTTTTTTGGGGCttttttttt                                                                                        >  2:448535/1‑149 (MQ=255)
                                                               gATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCa                                                                                       <  2:550306/148‑1 (MQ=255)
                                                               gATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTACGTTTATTCCATCTTTCGAGTGCTTCTTTTTCa                                                                                       <  1:1205258/148‑1 (MQ=255)
                                                                aTTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTTATTTTTTATATTAACACCACCCTTATTTTTTTCCACGGGTTGTTTTTTTTCTCTTTTGAGTGGttttttttttctt                                                                                    >  1:1227879/1‑145 (MQ=255)
                                                                   tCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTcca                                                                                <  1:852272/150‑1 (MQ=255)
                                                                        aGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTTCATCTTTCGAGTGCTTCTTTTTCACTTCCACCATa                                                                          <  1:1083835/150‑1 (MQ=255)
                                                                         gTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCATAAc                                                                        <  2:159304/151‑1 (MQ=255)
                                                                         gTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCATAAc                                                                        <  1:1222927/151‑1 (MQ=255)
                                                                          tttAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTTAAATAAAACCCCCCTTTTTTTTTTCAGTGGGTTGTTTTTTTTATCTTTCTTGTGCTTTTTTTTTTCTTCCACCATaa                                                                         >  2:704903/1‑149 (MQ=255)
                                                                           ttAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTTAATGAAGACCCCCCTTATTATTTTCAGGGGGTCGTTTTTTTCCTTTTTTTAGTGGTTCTTTTTTTCTTTCTCCATaaac                                                                       >  1:383524/1‑148 (MQ=255)
                                                                           ttAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCATAAcc                                                                       <  1:254795/150‑1 (MQ=255)
                                                                              gATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCATAACCGGtt                                                                   <  2:820435/151‑1 (MQ=255)
                                                                               aTTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTTATTTTTTAAATGTAAACCCCCACTTTTATTTTCATGTGTTTGTTTTTTTCATTTTTTTGGTGGTTTTTTTTTACTTtcacccctaaccg                                                                      >  2:1399494/1‑140 (MQ=255)
                                                                                ttttGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAAATTTTTACATTAAAACCTCCCTTATTTTTTTCAGTGGGTTGTTTTTTTCATCTTTTTAGGGGTTTTTTTTTACTTCCACCCTAAcc                                                                       >  1:957737/1‑145 (MQ=255)
                                                                                 tttGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCATAACCGGTTCg                                                                 <  2:962498/150‑1 (MQ=255)
                                                                                 tttGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCATAACCGGTTCg                                                                 <  1:1375698/150‑1 (MQ=255)
                                                                                   tGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTTATTTTTTTAAATAAAACCCCCATTATTTTTTTCAGTGGGTTGTTTATTTCATCTTTCTTGTGCTTTTTTTTTTCTTCTACCATAAACCGGTTgggg                                                               >  2:711805/1‑151 (MQ=255)
                                                                                   tGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTAAATTAACACCCCCCTTTTTTTTTTCCGTGGGTTGTTTTTTTCATCTTTTGCGGGGTTCTTTTTTACTTTCCCCCTAAACCGGTCgggg                                                              >  1:130154/1‑150 (MQ=255)
                                                                                     aaTATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTTCATGAAAAACCCCCTTTTTTTTTTCTGTGGGTTGGTTTTTTCCTCTTTTTAGTGGTTCTTTTTTTCTTTCCCCCTAACCGGGTCTGGGtt                                                            >  1:1393570/1‑151 (MQ=255)
                                                                                      atatCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCATAACCGGTTCGGGATTc                                                           <  1:340135/151‑1 (MQ=255)
                                                                                      atatCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGATATTTTTTATATTTTAACCCCCCATTTTATTTTTAATGGGTCGTTTTTTTTCTCTTTTTAGTGCTTTTTTTTTTCTTTCTCCATAACCgggttg                                                                 >  1:676208/1‑141 (MQ=255)
                                                                                        atCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTAAATGAAAACCCCCCTTTTTTTTTTCAGTGGTTTGTTTTTTTCCTCTTTTTAGTGCTTCTTTTTTTCTTTCACCCTAACCCGGTCGGGGtt                                                            >  1:145121/1‑148 (MQ=255)
                                                                                        atCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCATAACCGGTTCGGGATTCg                                                          <  2:980561/150‑1 (MQ=255)
                                                                                             gTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTAAATGAACAACTCCCTTATTATTTTCAGTGGGTTGGTTTTTTCCTCTTTTTAGGGCTTTTTTTTTTCTTTCCCCCTAAACCGGTCGGGGGTCGccttc                                                     >  1:1013399/1‑146 (MQ=255)
                                                                                                  gTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCATAACCGGTTCGGGATTCGCATCCGTTaca                                               <  2:472904/151‑1 (MQ=255)
                                                                                                     agtagGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAAATTTTTTAAAAGACACCCCCCTTTTTATTTTCTTGGGGTCGTTTTTTTTATTTTTTTGGTGGTTTTTTTTTACTTTCCCCCTAAACCGGTCCGGGTTTCCGTCCGTTTCACt                                             >  1:639285/1‑150 (MQ=255)
                                                                                                      gtagGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCATAACCGGTTCGGGATTCGCATCCGTTACACTTc                                           <  2:53509/151‑1 (MQ=255)
                                                                                                        agGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTTAATTAAAACCCCCCTTATTTTTTTCATTGGGTTGGTTTTTTTATATTTTTAGGGCTTTTTTTTTTCTTTCACCACAACCCGGTCGGGGTTTGCATTCGTTTCACtt                                            >  2:460328/1‑148 (MQ=255)
                                                                                                             tGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTACATTAACACCCCCCTTATTTTTTTCCGTGGGTTGGTTTTTTCCTCTTTTCAGTGCTTCTTTTTTTCTTTCACCCTAACCCGGTCGGGGTTTGCCTCCCTTTCACTTccc                                         >  1:171395/1‑144 (MQ=255)
                                                                                                               ccACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCATAACCGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTaa                                  <  2:807536/151‑1 (MQ=255)
                                                                                                                cACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCATAACCGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAat                                 <  1:749929/151‑1 (MQ=255)
                                                                                                                         tttCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCATAACCGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGaa                         <  2:1385259/150‑1 (MQ=255)
                                                                                                                           tCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCATAACCGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGAAAt                       <  2:682185/150‑1 (MQ=255)
                                                                                                                                 gTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCATAACCGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGAAATGGCTTTc                <  1:961831/151‑1 (MQ=255)
                                                                                                                                 gTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCATAACCGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGAAATGGCTTTc                <  1:1140021/151‑1 (MQ=255)
                                                                                                                                   ggCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCATAACCGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGAAATGGCTTTCa               <  1:651238/150‑1 (MQ=255)
                                                                                                                                   ggCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCATAACCGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGAAATGGCTTTCAc              <  1:593780/151‑1 (MQ=255)
                                                                                                                                    gCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCATAACCGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGAAATGGCTTTCAcc             >  2:629801/1‑151 (MQ=255)
                                                                                                                                      ggTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCATAACCGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGAAATGGCTTTCACCGt           <  1:892192/151‑1 (MQ=255)
                                                                                                                                          aaTTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCATAACCGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGAAATGGCTTTCACCGTTa         <  1:745119/150‑1 (MQ=255)
                                                                                                                                            ttttATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCATAACCGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGAAATGGCTTTCACCGTTACTg      >  1:58496/1‑151 (MQ=255)
                                                                                                                                               tataCATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCATAACCGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGAAATGGCTTTCACCGTTACTGATg   >  2:1224069/1‑151 (MQ=255)
                                                                                                                                                ataCATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCATAACCGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGAAATGGCTTTCACCGTTACTGATgg  <  1:1108796/151‑1 (MQ=255)
                                                                                                                                                ataCATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCATAACCGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGAAATGGCTTTCACCGTTACTGATgg  <  2:1393544/151‑1 (MQ=255)
                                                                                                                                                ataCATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCATAACCGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGAAATGGCTTTCACCGTTACTGATg   <  1:510856/150‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                            
TTTACTTTCTGCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCATAACCGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGAAATGGCTTTCACCGTTACTGATGG  >  NC_000913/1413592‑1413885

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