Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913589,9340TG47.1% 98.7 / 28.4 71T341P (ACC→CCC) nfrAbacteriophage N4 receptor, outer membrane subunit
Reads supporting (aligned to +/- strand):  ref base T (1/36);  new base G (33/0);  total (35/36)
Fisher's exact test for biased strand distribution p-value = 3.39e-19
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CATCAGTTGCCAGGTTAGTTGATCCAGGCGGGTAAGATTTGCCGGTTCTTGCTGATACAGCAATCGTGCCAGACGCAGAGCTTCAGCCTTGTTACGGGTCGCCACGCTGACAGCATAACGCTCCTCAAGCATTTCATTGGCGGGGAGGGTGGCGAGCAG‑TTTTTGCGCTGCGTCGTA‑CTGACCTTCTTTTAACAGCACCGGTAGCGTCGCGCCAACAACATACTGGC‑GGTTG‑TCGGCAAACTGTACCGTATAATTC‑GCCAACGCCTGAACGGGGTTAGCGCTGTATTT  >  NC_000913/589785‑590072
                                                                                                                                                     |                                                                                                                                               
cATCAGTTGCCAGGTTAGTTGATCCAGGCGGGTAAGATTTGCCGGTTCTTGCTGATACAGCAATCGTGCCAGACGCAGAGCTTCAGCCTTGTTACGGGTCGCCACGCTGACAGCATAACGCTCCTCAAGCATTTCATTGGCGGGGAGGGTg                                                                                                                                                <  2:1328695/151‑1 (MQ=255)
 aTCAGTTGCCAGGTTAGTTGATCCAGGCGGGTAAGAATTTCCGGGTCTTGCTGATACCGCAATCGTGCCACACGCCGAGCTTCAGCCTTGTTACGGGTCGCCCCGCTGACAGCCTAACGCTCCTCAAGCATTTCATTGGCGgggggggggg                                                                                                                                               >  1:1043051/1‑151 (MQ=255)
       tGCCAGGTTAGTTGATCCAGGCGGGTAAGATTTTCCGGGTCTTGCTGATAAAGCAATCGTGCCCGACGCAGAGCTTCAGCCTTTTTACGGGTCGCCACGCTGACAGCATAACGCTCCTCAAGCATTTCATTGGCGggggggggggggg                                                                                                                                            >  1:2855/1‑145 (MQ=255)
         ccAGGTTAGTTGATCCAGGCGGGTAAGATTTGCCGGTTCTTGCTGATACAGCAATCGTGCCAGACGCAGAGCTTCAGCCTTGTTACGGGTCGCCACGCTGACAGCATAACGCTCCTCAAGCATTTCATTGGCGGGGAGGGTGGCGAGCAG‑t                                                                                                                                      <  1:1515702/151‑1 (MQ=255)
         ccAGGTTAGTTGATCCAGGCGGGTAAGAATTGCCGGTTCTTGCCGATACAGCAATCGTGCCAGACGGAGAGCTTCAGCCTTGGTACGGGGTGCCACGCCGACAGCATAAAGCTCCTCAAGAATTTTCATGGCGgggggggggggggggag                                                                                                                                        >  1:962791/1‑143 (MQ=255)
          cAGGTTAGTTGATCCAGGCGGGTAAGATTTTCCGGTTCTTGCTGATACAGCAATCGTGCCAGACGCAGAGCTTCAGCCTTGTTACGGGTCGCCACGCTGACAGCATAACGCTCCTCAAGCATTTCATTGGCGgggggggggggggggggg                                                                                                                                       >  2:604004/1‑142 (MQ=255)
          cAGGTTAGTTGATCCAGGCGGGTAAGATTTGCCGGTTCTTGCTGATACAGCAATCGTGCCAGACGCAGAGCTTCAGCCTTGTTACGGGTCGCCACGCTGACAGCATAACGCTCCTCAAGCATTTCATTGGCGGGGAGGGTGGCGAGCAG‑tt                                                                                                                                     <  1:609980/151‑1 (MQ=255)
              ttAGTTGATCCAGGCGGGTTAAATTTGCCGGGTCTTGCTGATACAGCAATCGTGCCCGACGCAGAGCTTCAGCCTTGTTACGGGGCGCCACGCTGACAGCATAACGCTCCTCAAGCATTTCATTGGCGgggggggggggggggggggtttt                                                                                                                                  >  1:462931/1‑138 (MQ=255)
                    gATCCAGGCGGGTAAGATTTGCCGGTTCTTGCTGATACAGCAATCGTGCCAGACGCAGAGCTTCAGCCTTGTTACGGGTCGCCACGCTGACAGCATAACGCTCCTCAAGCATTTCATTGGCGGGGAGGGTGGCGAGCAG‑TTTTTGCGCTGc                                                                                                                           <  2:335829/151‑1 (MQ=255)
                      tCCAGGCGGGTAAGATTTGCCGGTTCTTGCTGATACAGCAATCGTGCCAGACGCAGAGCTTCAGCCTTGTTACGGGTCGCCACGCTGACAGCATAACGCTCCTCAAGCATTTCATTGGCGggggggggggggggggggtttttt                                                                                                                                 >  1:405533/1‑130 (MQ=255)
                        cAGGCGGGTAAGATTTGCCGGTTCTTGCTGATACAGCAATCGTGCCAGACGCAGAGCTTCAGCCTTGTTACGGGTCGCCACGCTGACAGCATAACGCTCCTCAAGCATTTCATTGGCGGGGGGGGGGGGGGGGAG‑gtttttt                                                                                                                                >  2:1128119/1‑140 (MQ=255)
                          ggCGGGTAAGATTTGCCGGTTCTTGCTGATACAGCAATCGTGCCAGACGCAGAGCTTCAGCCTTGTTACGGGTCGCCACGCTGACAGCATAACGCTCCTCAAGCATTTCATTGGCGGGGAGGGTGGCGAGCAG‑TTTTTGCGCTGcgtcgt                                                                                                                      <  1:396542/150‑1 (MQ=255)
                          ggCGGGTAAGATTTGCCGGTTCTTGCTGATACAGCAATCGTGCCAGACGCAGAGCTTCAGCCTTGTTACGGGTCGCCACGCTGACAGCATAACGCTCCTCAAGCATTTCATTGGCGGGGAGGGTGGCGAGCAG‑TTTTTGCGCTGCGTCGTa                                                                                                                     <  1:745776/151‑1 (MQ=255)
                              ggTAAGATTTGCCGGTTCTTGCTGATACAGCAATCGTGCCAGACGCAGAGCTTCAGCCTTGTTACGGGTCGCCACGCTGACAGCATAACGCTCCTCAAGCATTTCATTGGCGGGGAGGGTGGCGAGCAG‑TTTTTGCGCTGCGTCGTA‑CTGa                                                                                                                <  2:1435199/151‑1 (MQ=255)
                               gTAAGATTTGCCGGTTCTTGCTGATACAGCAATCGTGCCAGACGCAGAGCTTCAGCCTTGTTACGGGTCGCCACGCTGACAGCATAACGCTCCTCAAGCATTTCATTGGCGGGGGGGGGGGGGGGGAG‑GTTTTTCGCTGCGTCGTA‑Ct                                                                                                                  >  2:834710/1‑148 (MQ=255)
                               gTAAGATTTGCCGGTTCTTGCTGATACAGCAATCGTGCCAGACGCAGAGCTTCAGCCTTGTTACGGGTCGCCACGCTGACAGCATAACGCTCCTCAAGCATTTCATTGGCGGGGGGGGGGGGGGGGAG‑GTTTTTCGCCGCGTCGTA‑CTg                                                                                                                 >  1:1232998/1‑149 (MQ=255)
                                tAAGATTTGCCGGTTCTTGCTGATACAGCAATCGTGCCAGACGCAGAGCTTCAGCCTTGTTACGGGTCGCCACGCTGACAGCATAACGCTCCTCAAGCATTTCATTGGCGGGGAGGGGGGCGGGCCG‑GTTTTTCGGTGCGTCGTa                                                                                                                     >  2:1002763/1‑145 (MQ=255)
                                  aGATTTGCCGGTTCTTGCTGATACAGCAATCGTGCCAGACGCAGAGCTTCAGCCTTGTTACGGGTCGCCACGCTGACAGCATAACGCTCCTCAAGCATTTCATTGGCGGGGGGGGGGGCGGGGAGATTTTTGC‑CTGCGTCGTA‑cctaacct                                                                                                            >  2:1102555/1‑144 (MQ=255)
                                  aGATTTGCCGGTTCTTGCTGATACAGCAATCGTGCCAGACGCAGAGCTTCAGCCTTGTTACGGGTCGCCACGCTGACAGCATAACGCTCCTCAAGCATTTCATTGGCGGGGAGGGTGGCGAGCAG‑TTTTTGCGCTGCGTCGTA‑CTGACctt                                                                                                            <  1:814981/151‑1 (MQ=255)
                                   gATTTGCCGGTTCTTGCTGATACAGCAATCGTGCCAGACGCAGAGCTTCAGCCTTGTTACGGGTCGCCACGCTGACAGCATAACGCTCCTCAAGCATTTCATTGGCGGGGGGGGGGGGGAGCAG‑TTTTTTCGCTGCGTCGTA‑CTGACctt                                                                                                            >  2:553489/1‑150 (MQ=255)
                                    aTTTGCCGGTTCTTGCTGATACAGCAATCGTGCCAGACGCAGAGCTTCAGCCTTGTTACGGGTCGCCACGCTGACAGCATAACGCTCCTCAAGCATTTCATTGGCGGGGGGGGGGGGGGGGAG‑GTTTTTCGCCGCGGCGGA‑CTTAActt                                                                                                            >  1:1268970/1‑149 (MQ=255)
                                    aTTTGCCGGTTCTTGCTGATACAGCAATCGTGCCAGACGCAGAGCTTCAGCCTTGTTACGGGTCGCCACGCTGACAGCATAACGCTCCTCAAGCATTTCATTGGCGGGGGGGGGGGGGGGCAG‑TTTTTTTGGTGCGGCGTA‑Cttgccct                                                                                                            >  1:1018030/1‑143 (MQ=255)
                                    aTTTGCCGGTTCTTGCTGATACAGCAATCGTGCCAGACGCAGAGCTTCAGCCTTGTTACGGGTCGCCACGCTGACAGCATAACGCTCCTCAAGCATTTCATTGGCGGGGAGGGTGGCGAGCAG‑TTTTTGCGCTGCGTCGTA‑CTGACcttc                                                                                                           <  1:569076/150‑1 (MQ=255)
                                        gCCGGTTCTTGCTGATACAGCAATCGTGCCAGACGCAGAGCTTCAGCCTTGTTACGGGTCGCCACGCTGACAGCATAACGCTCCTCAAGCATTTCATTGGCGGGGAGGGTGGCGAGCAG‑TTTTTGCGCTGCGTCGTA‑CTGACCTTCTTTTa                                                                                                      <  1:1348494/151‑1 (MQ=255)
                                        gCCGGTTCTTGCTGATACAGCAATCGTGCCAGACGCAGAGCTTCAGCCTTGTTACGGGTCGCCACGCTGACAGCATAACGCTCCTCAAGCATTTCATTGGCGGGGAGGGGGGCGAGCCG‑GTTTTTCGCTGCGTCGGA‑CTGACCTTCttttt                                                                                                      >  2:280229/1‑150 (MQ=255)
                                           ggTTCTTGCTGATACAGCAATCGTGCCAGACGCAGAGCTTCAGCCTTGTTACGGGTCGCCACGCTGACAGCATAACGCTCCTCAAGCATTTCATTGGCGGGGGGGGGGGGGAGGAG‑GTTTTGCGCTGCGTCGGA‑CTGACCCTCttttta                                                                                                     >  1:1235938/1‑147 (MQ=255)
                                           ggTTCTTGCTGATACAGCAATCGTGCCAGACGCAGAGCTTCAGCCTTGTTACGGGTCGCCACGCTGACAGCATAACGCTCCTCAAGCATTTCATTGGCGGGGAGGGTGGCGAGCAG‑TTTTTGCGCTGCGTCGTA‑CTGACCTTCTTTTAACa                                                                                                   <  1:117467/151‑1 (MQ=255)
                                             ttCTTGCTGATTCAGCAATCGTGCCAGACGCAGACCTTCAGCCTTGTTACGGGTCGACACGCTGACAACATCACGCTCCTCAAGCATTTCATTGGCGGGGAGGGTGGCGCGCAG‑TTTTTGCGCTGCGTCGTA‑CTGCCCTTCTTTTAACa                                                                                                   <  1:506480/149‑1 (MQ=255)
                                              tCTTGCTGATACAGCAATCGTGCCAGACGCAGAGCTTCAGCCTTGTTACGGGTCGCCACGCTGACAGCATAACGCTCCTCAAGCATTTCATTGGCGGGGGGGGGGGGGAGGAG‑GTTTTTCGCTGCGGCGTA‑CTGACCTTCtttttaaaac                                                                                                 >  2:1482743/1‑144 (MQ=255)
                                               cTTGCTGATACAGCAATCGTGCCAGACGCAGAGCTTCAGCCTTGTTACGGGTCGCCACGCTGACAGCATAACGCTCCTCAAGCATTTCATTGGCGGGGGGGGGGGGGGGGAG‑TTTTTTCGCTTCGTCGGT‑CTGACCTTTTTTTTACAacac                                                                                               >  1:367561/1‑151 (MQ=255)
                                               cTTGCTGATACAGCAATCGTGCCAGACGCAGAGCTTCAGCCTTGTTACGGGTCGCCACGCTGACAGCATAACGCTCCTCAAGCATTTCATTGGCGGGGGGGGGGGGGGGGAG‑GTTTTTTGCTGCGGCGTA‑CTGGCCTTCtttttaaaaca                                                                                                >  2:450236/1‑143 (MQ=255)
                                                ttGCTGATACAGCAATCGTGCCAGACGCAGAGCTTCAGCCTTGTTACGGGTCGCCACGCTGACAGCATAACGCTCCTCAAGCATTTCATTGGCGGGGAGGGGGGGGAGCAG‑GTTTTGCGCTGCGCCGGA‑CTTACCTTTTTTTAACAGCAcc                                                                                              >  2:661272/1‑151 (MQ=255)
                                                 tGCTGATACAGCAATCGTGCCAGACGCAGAGCTTCAGCCTTGTTACGGGTCGCCACGCTGACAGCATAACGCTCCTCAAGCATTTCATTGGCGGGGAGGGTGGCGAGCAG‑TTTTTGCGCTGCGTCGTA‑CTGACCTTCTTTTAACAGCAc                                                                                               <  1:621479/149‑1 (MQ=255)
                                                  gCTGATACAGCAATCGTGCCAGACGCAGAGCTTCAGCCTTGTTACGGGTCGCCACGCTGACAGCATAACGCTCCTCAAGCATTTCATTGGCGGGGGGGGGGGGGAGGAG‑ATTTTTCGCTGCGTCGTG‑CTGACCTTCTTTTTAAAGCAACg                                                                                             >  2:1134875/1‑150 (MQ=255)
                                                  gCTGATACAGCAATCGTGCCAGACGCAGAGCTTCAGCCTTGTTACGGGTCGCCACGCTGACAGCATAACGCTCCTCAAGCATTTCATTGGCGGGGAGGGTGGCGAGCAG‑TTTTTGCGCTGCGTCGTA‑CTGACCTTCTTTTAACAGCACCgg                                                                                            <  2:137135/151‑1 (MQ=255)
                                                    tGATACAGCAATCGTGCCAGACGCAGAGCTTCAGCCTTGTTACGGGTCGCCACGCTGACAGCATAACGCTCCTCAAGCATTTCATTGGCGGGGGGGGGGGGGAGCAG‑TTTTTGCGCTGCGGCGTA‑CTGAACTTCTTTTTACAGCACCgg                                                                                            >  1:626707/1‑149 (MQ=255)
                                                     gATACAGCAATCGTGCCAGACGCAGAGCTTCAGCCTTGTTACGGGTCGCCACGCTGACAGCATAACGCTCCTCAAGCATTTCATTGGCGGGGAGGGTGGCGAGCAG‑TTTTTGCGCTGCGTCGTA‑CTGACCTTCTTTTAACAGCACCGGTa                                                                                          <  1:415122/150‑1 (MQ=255)
                                                     gATACAGCAATCGTGCCAGACGCAGAGCTTCAGCCTTGTTACGGGTCGCCACGCTGACAGCATAACGCTCCTCAAGCATTTCATTGGCGGGGAGGGTGGCGAGCAG‑TTTTTGCGCTGCGTCGTA‑CTGACCTTCTTTTAACAGCACCGGTa                                                                                          <  2:1029272/150‑1 (MQ=255)
                                                     gATACAGCAATCGTGCCAGACGCAGAGCTTCAGCCTTGTTACGGGTCGCCACGCTGACAGCATAACGCTCCTCAAGCATTTCATTGGCGGGGAGGGTGGCGAGCAG‑TTTTTGCGCTGCGTCGTA‑CTGACCTTCTTTTAACAGCACCGGTAg                                                                                         <  1:1316813/151‑1 (MQ=255)
                                                     gATACAGCAATCGTGCCAGACGCAGAGCTTCAGCCTTGTTACGGGTCGCCACGCTGACAGCATAACGCTCCTCAAGCATTTCATTGGCGGGGAGGGTGGCGAGCAG‑TTTTTGCGCTGCGTCGTA‑CTGACCTTCTTTTAACAGCACCGGTAg                                                                                         <  2:1105368/151‑1 (MQ=255)
                                                        aCAGCAATCGTGCCAGACGCAGAGCTTCAGCCTTGTTACGGGTCGCCACGCTGACAGCATAACGCTCCTCAAGCATTTCATTGGCGGGGAGGGTGGCGAGCAG‑TTTTTGCGCTGCGTCGTA‑CTGACCTTCTTTTAACAGCACCGGTAGCGt                                                                                      <  2:784223/151‑1 (MQ=255)
                                                         cAGCAATCGTGCCAGACGCAGAGCTTCAGCCTTGTTACGGGTCGCCACGCTGACAGCATAACGCTCCTCAAGCATTTCATTGGCGGGGGGGGGGGCGGGCAG‑GTTTTTCGCCGCGTCGGG‑CTGACCTTTTTTTTACAGCACCgggtgcg                                                                                       >  1:9214/1‑144 (MQ=255)
                                                          aGCAATCGTGCCAGACGCAGAGCTTCAGCCTTGTTACGGGTCGCCACGCTGACAGCATAACGCTCCTCAAGCATTTCATTGGCGGGGGGGGGGGGGAGCCG‑TTTTTTTGCCGCGTCGTA‑CTGTCCTTCTTTTTAAAGCACCGGGAGCg                                                                                       >  2:1545751/1‑148 (MQ=255)
                                                           gCAATCGTGCCAGACGCAGAGCTTCAGCCTTGTTACGGGTCGCCACGCTGACAGCATAACGCTCCTCAAGCATTTCATTGGCGGGGAGGGTGGCGAGCAG‑TTTTTGCGCTGCGTCGTA‑CTGACCTTCTTTTAACAGCACCGGTAGCGTcgc                                                                                   <  2:1029470/151‑1 (MQ=255)
                                                           gCAATCGTGCCAGACGCAGAGCTTCAGCCTTGTTACGGGTCGCCACGCTGACAGCATAACGCTCCTCAAGCATTTCATTGGCGGGGAGGGTGGCGAGCAG‑TTTTTGCGCTGCGTCGTA‑CTGACCTTCTTTTAACAGCACCGGTAGCGTcgc                                                                                   <  2:26938/151‑1 (MQ=255)
                                                            cAATCGTGCCAGACGCAGAGCTTCAGCCTTGTTACGGGTCGCCACGCTGACAGCATAACGCTCCTCAAGCATTTCATTGGCGGGGAGGGGGGGGGGCAG‑TTTTTTCGCTGGGTCGGA‑ATGGACTTCTTTTTAAAACACCCGTAGGGTcgc                                                                                   >  2:797734/1‑150 (MQ=255)
                                                              aTCGTGCCAGACGCAGAGCTTCAGCCTTGTTACGGGTCGCCACGCTGACAGCATAACGCTCCTCAAGCATTTCATTGGCGGGGGGGGGGGCGAGCAG‑GTTTTTCGCTGCGTCGTA‑CTGACCTTCTTTTTACAACACCCGTAGCGTCGCGcc                                                                                >  2:1545987/1‑151 (MQ=255)
                                                                 gTGCCAGACGCAGAGCTTCAGCCTTGTTACGGGTCGCCACGCTGACAGCATAACGCTCCTCAAGCATTTCATTGGCGGGGGGGGGGGGGGGAAG‑GTTTTTCGCTGCCGCGGT‑CTGGCCCTCTTTTTACAACACCGGGAGCGGCGCGccc                                                                               >  1:536987/1‑148 (MQ=255)
                                                                 gTGCCAGACGCAGAGCTTCAGCCTTGTTACGGGTCGCCACGCTGACAGCATAACGCTCCTCAAGCATTTCATTGGCGGGGGGGGGGGCGGGCAG‑GTTTTTCGCTGCGTCGTA‑CTGCCCCTTTTTTTAAAGCACCGGGAGCGTCGCGcc                                                                                >  2:1158069/1‑148 (MQ=255)
                                                                         cGCAGAGCTTCAGCCTTGTTACGGGTCGCCACGCTGACAGCATAACGCTCCTCAAGCATTTCATTGGCGGGGAGGGTGGCGAGCAG‑TTTTTGCGCTGCGTCGTA‑CTGACCTTCTTTTAACAGCACCGGTAGCGTCGCGCCAACAACATACt                                                                     <  1:222230/151‑1 (MQ=255)
                                                                          gCAGAGCTTCAGCCTTGTTACGGGTCGCCACGCTGACAGCATAACGCTCCTCAAGCATTTCATTGGCGGGGGGGGGGGGGGGGAG‑TTTTTTCGCTGCGGCGTG‑CTGACCCTCTTTTTACAGCACCCGGAGCGGCGCGCccacaa                                                                           >  1:962878/1‑144 (MQ=255)
                                                                           cAGAGCTTCAGCCTTGTTACGGGTCGCCACGCTGACAGCATAACGCTCCTCAAGCATTTCATTGGCGGGGAGGGTGGCGAGCAG‑TTTTTGCGCTGCGTCGTA‑CTGACCTTCTTTTAACAGCACCGGTAGCGTCGCGCCAACAACATACTgg                                                                   <  2:1005078/151‑1 (MQ=255)
                                                                             gagCTTCAGCCTTGTTACGGGTCGCCACGCTGACAGCATAACGCTCCTCAAGCATTTCATTGGCGGGGGGGGGGGCGGGCAG‑TTTTTTCGCTGCGTCGTA‑CTGACCTTCTTTTTACAGCACCGGGAGCGGCGCGCCCACAACATACTGGC‑g                                                                >  2:508275/1‑151 (MQ=255)
                                                                             gagCTTCAGCCTTGTTACGGGTCGCCACGCTGACAGCATAACGCTCCTCAAGCATTTCATTGGCGGGGAGGGTGGCGAGCAG‑TTTTTGCGCTGCGTCGTA‑CTGACCTTCTTTTAACAGCACCGGTAGCGTCGCGCCAACAACATACTGGC‑g                                                                <  2:487548/151‑1 (MQ=255)
                                                                               gCTTCAGCCTTGTTACGGGTCGCCACGCTGACAGCATAACGCTCCTCAAGCATTTCATTGGCGGGGGGGGGGGCGGGGAG‑GTTTTTTGCTGGGGCGCA‑GTTACCCTCTTTTTACAACACCGGGAGCGTGGCGCCCACAACAAACTGGC‑gg                                                               >  2:1506606/1‑150 (MQ=255)
                                                                                cTTCAGCCTTGTTACGGGTCGCCACGCTGACAGCATAACGCTCCTCAAGCATTTCATTGGCGGGGGGGGGGGGGGGCAG‑GGTTTTTGCTGCGTCGTA‑CTGTCCTTCTTTTTTAAGCAACCGTGGCGTCGCGCCAACAAAACACTTGC‑ggg                                                              >  1:1059727/1‑149 (MQ=255)
                                                                                cTTCAGCCTTGTTACGGGTCGCCACGCTGACAGCATAACGCTCCTCAAGCATTTCATTGGCGGGGGGGGGGGGGAGCAG‑GTTTTTTGCCGCGTCGTG‑CCGTCCCTCTTTTTTAAACACCGGGAGGGTCGCGCCCAAAAAAAACTGGC‑ggg                                                              >  1:1387663/1‑149 (MQ=255)
                                                                                 ttCAGCCTTGTTACGGGTCGCCACGCTGACAGCATAACGCTCCTCAAGCATTTCATTGGCGGGGAGGGGGGCGAGAAG‑TTTTTTCGCTGCGGCGTA‑CTGACCTTCTTTTAACAACACCGGTAGCGGCGCGCCAAAAACATACTGGC‑gg                                                               >  2:224647/1‑148 (MQ=255)
                                                                                     gCCTTGTTACGGGTCGCCACGCTGACAGCATAACGCTCCTCAAGCATTTCATTGGCGGGGGGGGGGGGGGGCAG‑GTTTTTTGCTGCGGCGTA‑CTGACCTTTTTTTTACAGCACCCGGAGCGGCGCGCCCAAAACAAACTGGC‑ggg                                                              >  1:799154/1‑144 (MQ=255)
                                                                                     gCCTTGTTACGGGTCGCCACGCTGACAGCATAACGCTCCTCAAGCATTTCATTGGCGGGGGGGGGGGCGGGCAG‑GTTTTGCGCCGCGTCGGA‑CTGACCTTCTTTTTACAACACCGGTAGCGTCGCGCCAACAACATACTGGC‑ggg                                                              >  1:726678/1‑144 (MQ=255)
                                                                                     gCCTTGTTACGGGTCGCCACGCTGACAGCATAACGCTCCTCAAGCATTTCATTGGCGGGGAGGGCGGTGAGGAG‑GTTTTGAGCCGCGGCGGT‑CCGACCTTCTTTTATCAGACACGGTGGCGGTCCGCcaacacactaatgggg                                                                 >  1:13929/1‑132 (MQ=255)
                                                                                      ccTTTTTACGGGTCGCCACGCTGACAGCATAACGCTCCTCAAGCATTTCATTGGCGGGGGGGGGGGGGAGGCG‑TTTTTTCGCTGCCTCCTA‑CTGACCTTCTTTTAACAGCACCGGTAGCGTCGCGCCCCCAACATACTGGC‑gg                                                               >  2:489567/1‑143 (MQ=255)
                                                                                       cTTGTTACGGGTCGCCACGCTGACAGCATAACGCTCCTCAAGCATTTCATTGGCGGGGAGGGTGGCGAGCAG‑TTGTTGCGCTGCGTCGTA‑CTGACCTTCTTTTAACAGCACCGGTAGCGTCGCGCCAACAACATACTGGC‑GGTTG‑TCGGc                                                      <  1:1102495/150‑1 (MQ=255)
                                                                                               gggTCGCCACGCTGACAGCATAACGCTCCTCAAGCATTTCATTGGCGGGGGGGGGGGCGGGCAG‑GTTTTTTGCTGCGTCGTA‑CTGACCTTCTTTTTAAAGCACCGGTAGCGTCGCGCCAACAACATAATGGG‑GGGTggtgggggaa                                                   >  2:1044864/1‑137 (MQ=255)
                                                                                                   cGCCACGCTGACAGCATAACGCTCCTCAAGCATTTCATTGGCGGGGAGGGTGGCGAGCAG‑TTTTTGCGCTGCGTCGTA‑CTGACCTTCTTTTAACAGCACCGGTAGCGTCGCGCCAACAACATACTGGC‑GGTTG‑TCGGCAAACTGTACCGta                                         <  2:289377/151‑1 (MQ=255)
                                                                                                    gCCACGCTGACAGCATAACGCTCCTCAAGCATTTCATTGGCGGGGAGGGTGGCGAGCAG‑TTTTTGCGCTGCGTCGTA‑CTGACCTTCTTTTAACAGCACCGGTAGCGTCGCGCCAACAACATACTGGC‑GGTTG‑TCGGCAAACTGTACCGtat                                        <  2:841413/151‑1 (MQ=255)
                                                                                                           tGACAGCATAACGCTCCTCAAGCATTTCATTGGCGGGGGGGGGGGGGGGCAG‑CTTTTTCGCTTGGTCGTATCTACCCTTTTTTTAAAACAACCGG‑GGCGGCGCGCCAACCACATACTGGCGGGGTG‑TTGGCaaaaattacgccgt                                        >  1:1434188/1‑134 (MQ=255)
                                                                                                            gACAGCATAGCGCTCCTCAAGCATTTCATTGGCGGGGAGGGTGGCGAGCAG‑TTTTTGCGCTGCGTCGTA‑CTGACCTTCTTTTAACAGCACCGGTAGCGTCGCGCCAACAACATACTGGC‑GGTTG‑TCGGCAAACTGTACCGTATAATTC‑Gc                                <  1:816270/150‑1 (MQ=255)
                                                                                                                gCATAACGCTCCTCAAGCATTTCATTGGCGGGGAGGGTGGCGAGCAG‑TTTTTGCGCTGCGTCGTA‑CTGACCTTCTTTTAACAGCACCGGTAGCGTCGCGCCAACAACATACTGGC‑GGTTG‑TCGGCAAACTGTACCGTATAATTC‑GCCAACg                           <  2:1196488/151‑1 (MQ=255)
                                                                                                                 cATAACGCTCCTCAAGCATTTCATTGGCGGGGAGGGTGGCGAGCAG‑TTTTTGCGCTGCGTCGTA‑CTGACCTTCTTTTAACAGCACCGGTAGCGTCGCGCCAACAACATACTGGC‑GGTTG‑TCGGCAAACTGTACCGTATAATTC‑GCCAACGc                          <  2:1413335/151‑1 (MQ=255)
                                                                                                                     aCGCTCCTCAAGCATTTCATTGGCGGGGAGGGTGGCGAGCAG‑TTTTTGCGCTGCGTCGTA‑CTGACCTTCTTTTAACAGCACCGGTAGCGTCGCGCCAACAACATACTGGC‑GGTTG‑TCGGCAAACTGTACCGTATAATTC‑GCCAACGCCTGa                      <  1:406749/151‑1 (MQ=255)
                                                                                                                          cctcAAGCATTTCATTGGCGGGGAGGGTGGCGAGCAG‑TTTTTGCGCTGCGTCGTA‑CTGACCTTCTTTTAACAGCACCGGTAGCGTCGCGCCAACAACATACTGGC‑GGTTG‑TCGGCAAACTGTACCGTATAATTC‑GCCAACGCCTGAACgg                  <  1:1506516/150‑1 (MQ=255)
                                                                                                                               aGCATTTCATTGGCGGGGAGGGTGGCGAGCAG‑TTTTTGCGCTGCGTCGTA‑CTGACCTTCTTTTAACAGCACCGGTAGCGTCGCGCCAACAACATACTGGC‑GGTTG‑TCGGCAAACTGTACCGTATAATTC‑GCCAACGCCTGAACGGGGTTAg            <  2:818853/151‑1 (MQ=255)
                                                                                                                                     tCATTGGCGGGGAGGGTGGCGAGCAG‑TTTTTGCGCTGCGTCGTA‑CTGACCTTCTTTTAACAGCACCGGTAGCGTCGCGCCAACAACATACTGGC‑GGTTG‑TCGGCAAACTGTACCGTATAATTC‑GCCAACGCCTGAACGGGGTTAGCGCTg       <  1:563973/150‑1 (MQ=255)
                                                                                                                                      cATTGGCGGGGAGGGTGGCGAGCAG‑TTTTTGCGCTGCGTCGTA‑CTGACCTTCTTTTAACAGCACCGGTAGCGTCGCGCCAACAACATACTGGC‑GGTTG‑TCGGCAAACTGTACCGTATAATTC‑GCCAACGCCTGAACGGGGTTAGCGCTGt      <  1:91236/150‑1 (MQ=255)
                                                                                                                                      cATTGGCGGGGAGGGTGGCGAGCAG‑TTTTTGCGCTGCGTCGTA‑CTGACCTTCTTTTAACAGCACCGGTAGCGTCGCGCCAACAACATACTGGC‑GGTTG‑TCGGCAAACTGTACCGTATAATTC‑GCCAACGCCTGAACGGGGTTAGCGCTGt      <  2:1228607/150‑1 (MQ=255)
                                                                                                                                           gCGGGGAGGGTGGCGAGCAG‑TTTTTGCGCTGCGTCGTA‑CTGACCTTCTTTTAACAGCACCGGTAGCGTCGCGCCAACAACATACTGGC‑GGTTGTTCGGCAAA‑TCTACCGTATAATTCACCCAACACCGGAAGgggggtgcggcggctttt  >  1:219174/1‑136 (MQ=255)
                                                                                                                                                     |                                                                                                                                               
CATCAGTTGCCAGGTTAGTTGATCCAGGCGGGTAAGATTTGCCGGTTCTTGCTGATACAGCAATCGTGCCAGACGCAGAGCTTCAGCCTTGTTACGGGTCGCCACGCTGACAGCATAACGCTCCTCAAGCATTTCATTGGCGGGGAGGGTGGCGAGCAG‑TTTTTGCGCTGCGTCGTA‑CTGACCTTCTTTTAACAGCACCGGTAGCGTCGCGCCAACAACATACTGGC‑GGTTG‑TCGGCAAACTGTACCGTATAATTC‑GCCAACGCCTGAACGGGGTTAGCGCTGTATTT  >  NC_000913/589785‑590072

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