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

CATCAGTTGCCAGGTTAGTTGATCCAGGCGGGTAAGATTTGCCGGTTCTTGCTGATACAGCAATCGTGCCAGACGCAGAGCTTCAGCCTTGTTACGGGTCGCCACGCTGACAGCATAACGCTCCTCAAGCATTTCATTGGCGGGGAGGGTGGCGAGCAGTTTTTGCGCTGCGTCGTACTGACCTTCTTTTAACAGCACCGGTAGCGTCGCGCCAACAACATACTGGCGGTTGTCGGCAAACTGTACCGTATAATTCGCCAACGCCTGAACGGGGTTAGCGCTGTATTT  >  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)
         ccAGGTTAGTTGATCCAGGCGGGTAAGATTTGCCGGTTCTTGCTGATACAGCAATCGTGCCAGACGCAGAGCTTCAGCCTTGTTACGGGTCGCCACGCTGACAGCATAACGCTCCTCAAGCATTTCATTGGCGGGGAGGGTGGCGAGCAGt                                                                                                                                      <  1:1515702/151‑1 (MQ=255)
         ccAGGTTAGTTGATCCAGGCGGGTAAGAATTGCCGGTTCTTGCCGATACAGCAATCGTGCCAGACGGAGAGCTTCAGCCTTGGTACGGGGTGCCACGCCGACAGCATAAAGCTCCTCAAGAATTTTCATGGCGgggggggggggggggag                                                                                                                                        >  1:962791/1‑143 (MQ=255)
          cAGGTTAGTTGATCCAGGCGGGTAAGATTTTCCGGTTCTTGCTGATACAGCAATCGTGCCAGACGCAGAGCTTCAGCCTTGTTACGGGTCGCCACGCTGACAGCATAACGCTCCTCAAGCATTTCATTGGCGgggggggggggggggggg                                                                                                                                       >  2:604004/1‑142 (MQ=255)
          cAGGTTAGTTGATCCAGGCGGGTAAGATTTGCCGGTTCTTGCTGATACAGCAATCGTGCCAGACGCAGAGCTTCAGCCTTGTTACGGGTCGCCACGCTGACAGCATAACGCTCCTCAAGCATTTCATTGGCGGGGAGGGTGGCGAGCAGtt                                                                                                                                     <  1:609980/151‑1 (MQ=255)
              ttAGTTGATCCAGGCGGGTTAAATTTGCCGGGTCTTGCTGATACAGCAATCGTGCCCGACGCAGAGCTTCAGCCTTGTTACGGGGCGCCACGCTGACAGCATAACGCTCCTCAAGCATTTCATTGGCGgggggggggggggggggggtttt                                                                                                                                  >  1:462931/1‑138 (MQ=255)
                    gATCCAGGCGGGTAAGATTTGCCGGTTCTTGCTGATACAGCAATCGTGCCAGACGCAGAGCTTCAGCCTTGTTACGGGTCGCCACGCTGACAGCATAACGCTCCTCAAGCATTTCATTGGCGGGGAGGGTGGCGAGCAGTTTTTGCGCTGc                                                                                                                           <  2:335829/151‑1 (MQ=255)
                      tCCAGGCGGGTAAGATTTGCCGGTTCTTGCTGATACAGCAATCGTGCCAGACGCAGAGCTTCAGCCTTGTTACGGGTCGCCACGCTGACAGCATAACGCTCCTCAAGCATTTCATTGGCGggggggggggggggggggtttttt                                                                                                                                 >  1:405533/1‑130 (MQ=255)
                        cAGGCGGGTAAGATTTGCCGGTTCTTGCTGATACAGCAATCGTGCCAGACGCAGAGCTTCAGCCTTGTTACGGGTCGCCACGCTGACAGCATAACGCTCCTCAAGCATTTCATTGGCGGGGGGGGGGGGGGGGAGgtttttt                                                                                                                                >  2:1128119/1‑140 (MQ=255)
                          ggCGGGTAAGATTTGCCGGTTCTTGCTGATACAGCAATCGTGCCAGACGCAGAGCTTCAGCCTTGTTACGGGTCGCCACGCTGACAGCATAACGCTCCTCAAGCATTTCATTGGCGGGGAGGGTGGCGAGCAGTTTTTGCGCTGcgtcgt                                                                                                                      <  1:396542/150‑1 (MQ=255)
                          ggCGGGTAAGATTTGCCGGTTCTTGCTGATACAGCAATCGTGCCAGACGCAGAGCTTCAGCCTTGTTACGGGTCGCCACGCTGACAGCATAACGCTCCTCAAGCATTTCATTGGCGGGGAGGGTGGCGAGCAGTTTTTGCGCTGCGTCGTa                                                                                                                     <  1:745776/151‑1 (MQ=255)
                              ggTAAGATTTGCCGGTTCTTGCTGATACAGCAATCGTGCCAGACGCAGAGCTTCAGCCTTGTTACGGGTCGCCACGCTGACAGCATAACGCTCCTCAAGCATTTCATTGGCGGGGAGGGTGGCGAGCAGTTTTTGCGCTGCGTCGTACTGa                                                                                                                <  2:1435199/151‑1 (MQ=255)
                               gTAAGATTTGCCGGTTCTTGCTGATACAGCAATCGTGCCAGACGCAGAGCTTCAGCCTTGTTACGGGTCGCCACGCTGACAGCATAACGCTCCTCAAGCATTTCATTGGCGGGGGGGGGGGGGGGGAGGTTTTTCGCTGCGTCGTACt                                                                                                                  >  2:834710/1‑148 (MQ=255)
                               gTAAGATTTGCCGGTTCTTGCTGATACAGCAATCGTGCCAGACGCAGAGCTTCAGCCTTGTTACGGGTCGCCACGCTGACAGCATAACGCTCCTCAAGCATTTCATTGGCGGGGGGGGGGGGGGGGAGGTTTTTCGCCGCGTCGTACTg                                                                                                                 >  1:1232998/1‑149 (MQ=255)
                                tAAGATTTGCCGGTTCTTGCTGATACAGCAATCGTGCCAGACGCAGAGCTTCAGCCTTGTTACGGGTCGCCACGCTGACAGCATAACGCTCCTCAAGCATTTCATTGGCGGGGAGGGGGGCGGGCCGGTTTTTCGGTGCGTCGTa                                                                                                                     >  2:1002763/1‑145 (MQ=255)
                                  aGATTTGCCGGTTCTTGCTGATACAGCAATCGTGCCAGACGCAGAGCTTCAGCCTTGTTACGGGTCGCCACGCTGACAGCATAACGCTCCTCAAGCATTTCATTGGCGGGGGGGGGGGCGGGGAGATTTTTGCCTGCGTCGTAcctaacct                                                                                                            >  2:1102555/1‑144 (MQ=255)
                                  aGATTTGCCGGTTCTTGCTGATACAGCAATCGTGCCAGACGCAGAGCTTCAGCCTTGTTACGGGTCGCCACGCTGACAGCATAACGCTCCTCAAGCATTTCATTGGCGGGGAGGGTGGCGAGCAGTTTTTGCGCTGCGTCGTACTGACctt                                                                                                            <  1:814981/151‑1 (MQ=255)
                                   gATTTGCCGGTTCTTGCTGATACAGCAATCGTGCCAGACGCAGAGCTTCAGCCTTGTTACGGGTCGCCACGCTGACAGCATAACGCTCCTCAAGCATTTCATTGGCGGGGGGGGGGGGGAGCAGTTTTTTCGCTGCGTCGTACTGACctt                                                                                                            >  2:553489/1‑150 (MQ=255)
                                    aTTTGCCGGTTCTTGCTGATACAGCAATCGTGCCAGACGCAGAGCTTCAGCCTTGTTACGGGTCGCCACGCTGACAGCATAACGCTCCTCAAGCATTTCATTGGCGGGGGGGGGGGGGGGGAGGTTTTTCGCCGCGGCGGACTTAActt                                                                                                            >  1:1268970/1‑149 (MQ=255)
                                    aTTTGCCGGTTCTTGCTGATACAGCAATCGTGCCAGACGCAGAGCTTCAGCCTTGTTACGGGTCGCCACGCTGACAGCATAACGCTCCTCAAGCATTTCATTGGCGGGGGGGGGGGGGGGCAGTTTTTTTGGTGCGGCGTACttgccct                                                                                                            >  1:1018030/1‑143 (MQ=255)
                                    aTTTGCCGGTTCTTGCTGATACAGCAATCGTGCCAGACGCAGAGCTTCAGCCTTGTTACGGGTCGCCACGCTGACAGCATAACGCTCCTCAAGCATTTCATTGGCGGGGAGGGTGGCGAGCAGTTTTTGCGCTGCGTCGTACTGACcttc                                                                                                           <  1:569076/150‑1 (MQ=255)
                                        gCCGGTTCTTGCTGATACAGCAATCGTGCCAGACGCAGAGCTTCAGCCTTGTTACGGGTCGCCACGCTGACAGCATAACGCTCCTCAAGCATTTCATTGGCGGGGAGGGTGGCGAGCAGTTTTTGCGCTGCGTCGTACTGACCTTCTTTTa                                                                                                      <  1:1348494/151‑1 (MQ=255)
                                        gCCGGTTCTTGCTGATACAGCAATCGTGCCAGACGCAGAGCTTCAGCCTTGTTACGGGTCGCCACGCTGACAGCATAACGCTCCTCAAGCATTTCATTGGCGGGGAGGGGGGCGAGCCGGTTTTTCGCTGCGTCGGACTGACCTTCttttt                                                                                                      >  2:280229/1‑150 (MQ=255)
                                           ggTTCTTGCTGATACAGCAATCGTGCCAGACGCAGAGCTTCAGCCTTGTTACGGGTCGCCACGCTGACAGCATAACGCTCCTCAAGCATTTCATTGGCGGGGGGGGGGGGGAGGAGGTTTTGCGCTGCGTCGGACTGACCCTCttttta                                                                                                     >  1:1235938/1‑147 (MQ=255)
                                           ggTTCTTGCTGATACAGCAATCGTGCCAGACGCAGAGCTTCAGCCTTGTTACGGGTCGCCACGCTGACAGCATAACGCTCCTCAAGCATTTCATTGGCGGGGAGGGTGGCGAGCAGTTTTTGCGCTGCGTCGTACTGACCTTCTTTTAACa                                                                                                   <  1:117467/151‑1 (MQ=255)
                                             ttCTTGCTGATTCAGCAATCGTGCCAGACGCAGACCTTCAGCCTTGTTACGGGTCGACACGCTGACAACATCACGCTCCTCAAGCATTTCATTGGCGGGGAGGGTGGCGCGCAGTTTTTGCGCTGCGTCGTACTGCCCTTCTTTTAACa                                                                                                   <  1:506480/149‑1 (MQ=255)
                                              tCTTGCTGATACAGCAATCGTGCCAGACGCAGAGCTTCAGCCTTGTTACGGGTCGCCACGCTGACAGCATAACGCTCCTCAAGCATTTCATTGGCGGGGGGGGGGGGGAGGAGGTTTTTCGCTGCGGCGTACTGACCTTCtttttaaaac                                                                                                 >  2:1482743/1‑144 (MQ=255)
                                               cTTGCTGATACAGCAATCGTGCCAGACGCAGAGCTTCAGCCTTGTTACGGGTCGCCACGCTGACAGCATAACGCTCCTCAAGCATTTCATTGGCGGGGGGGGGGGGGGGGAGTTTTTTCGCTTCGTCGGTCTGACCTTTTTTTTACAacac                                                                                               >  1:367561/1‑151 (MQ=255)
                                               cTTGCTGATACAGCAATCGTGCCAGACGCAGAGCTTCAGCCTTGTTACGGGTCGCCACGCTGACAGCATAACGCTCCTCAAGCATTTCATTGGCGGGGGGGGGGGGGGGGAGGTTTTTTGCTGCGGCGTACTGGCCTTCtttttaaaaca                                                                                                >  2:450236/1‑143 (MQ=255)
                                                ttGCTGATACAGCAATCGTGCCAGACGCAGAGCTTCAGCCTTGTTACGGGTCGCCACGCTGACAGCATAACGCTCCTCAAGCATTTCATTGGCGGGGAGGGGGGGGAGCAGGTTTTGCGCTGCGCCGGACTTACCTTTTTTTAACAGCAcc                                                                                              >  2:661272/1‑151 (MQ=255)
                                                 tGCTGATACAGCAATCGTGCCAGACGCAGAGCTTCAGCCTTGTTACGGGTCGCCACGCTGACAGCATAACGCTCCTCAAGCATTTCATTGGCGGGGAGGGTGGCGAGCAGTTTTTGCGCTGCGTCGTACTGACCTTCTTTTAACAGCAc                                                                                               <  1:621479/149‑1 (MQ=255)
                                                  gCTGATACAGCAATCGTGCCAGACGCAGAGCTTCAGCCTTGTTACGGGTCGCCACGCTGACAGCATAACGCTCCTCAAGCATTTCATTGGCGGGGGGGGGGGGGAGGAGATTTTTCGCTGCGTCGTGCTGACCTTCTTTTTAAAGCAACg                                                                                             >  2:1134875/1‑150 (MQ=255)
                                                  gCTGATACAGCAATCGTGCCAGACGCAGAGCTTCAGCCTTGTTACGGGTCGCCACGCTGACAGCATAACGCTCCTCAAGCATTTCATTGGCGGGGAGGGTGGCGAGCAGTTTTTGCGCTGCGTCGTACTGACCTTCTTTTAACAGCACCgg                                                                                            <  2:137135/151‑1 (MQ=255)
                                                    tGATACAGCAATCGTGCCAGACGCAGAGCTTCAGCCTTGTTACGGGTCGCCACGCTGACAGCATAACGCTCCTCAAGCATTTCATTGGCGGGGGGGGGGGGGAGCAGTTTTTGCGCTGCGGCGTACTGAACTTCTTTTTACAGCACCgg                                                                                            >  1:626707/1‑149 (MQ=255)
                                                     gATACAGCAATCGTGCCAGACGCAGAGCTTCAGCCTTGTTACGGGTCGCCACGCTGACAGCATAACGCTCCTCAAGCATTTCATTGGCGGGGAGGGTGGCGAGCAGTTTTTGCGCTGCGTCGTACTGACCTTCTTTTAACAGCACCGGTa                                                                                          <  1:415122/150‑1 (MQ=255)
                                                     gATACAGCAATCGTGCCAGACGCAGAGCTTCAGCCTTGTTACGGGTCGCCACGCTGACAGCATAACGCTCCTCAAGCATTTCATTGGCGGGGAGGGTGGCGAGCAGTTTTTGCGCTGCGTCGTACTGACCTTCTTTTAACAGCACCGGTa                                                                                          <  2:1029272/150‑1 (MQ=255)
                                                     gATACAGCAATCGTGCCAGACGCAGAGCTTCAGCCTTGTTACGGGTCGCCACGCTGACAGCATAACGCTCCTCAAGCATTTCATTGGCGGGGAGGGTGGCGAGCAGTTTTTGCGCTGCGTCGTACTGACCTTCTTTTAACAGCACCGGTAg                                                                                         <  1:1316813/151‑1 (MQ=255)
                                                     gATACAGCAATCGTGCCAGACGCAGAGCTTCAGCCTTGTTACGGGTCGCCACGCTGACAGCATAACGCTCCTCAAGCATTTCATTGGCGGGGAGGGTGGCGAGCAGTTTTTGCGCTGCGTCGTACTGACCTTCTTTTAACAGCACCGGTAg                                                                                         <  2:1105368/151‑1 (MQ=255)
                                                        aCAGCAATCGTGCCAGACGCAGAGCTTCAGCCTTGTTACGGGTCGCCACGCTGACAGCATAACGCTCCTCAAGCATTTCATTGGCGGGGAGGGTGGCGAGCAGTTTTTGCGCTGCGTCGTACTGACCTTCTTTTAACAGCACCGGTAGCGt                                                                                      <  2:784223/151‑1 (MQ=255)
                                                         cAGCAATCGTGCCAGACGCAGAGCTTCAGCCTTGTTACGGGTCGCCACGCTGACAGCATAACGCTCCTCAAGCATTTCATTGGCGGGGGGGGGGGCGGGCAGGTTTTTCGCCGCGTCGGGCTGACCTTTTTTTTACAGCACCgggtgcg                                                                                       >  1:9214/1‑144 (MQ=255)
                                                          aGCAATCGTGCCAGACGCAGAGCTTCAGCCTTGTTACGGGTCGCCACGCTGACAGCATAACGCTCCTCAAGCATTTCATTGGCGGGGGGGGGGGGGAGCCGTTTTTTTGCCGCGTCGTACTGTCCTTCTTTTTAAAGCACCGGGAGCg                                                                                       >  2:1545751/1‑148 (MQ=255)
                                                           gCAATCGTGCCAGACGCAGAGCTTCAGCCTTGTTACGGGTCGCCACGCTGACAGCATAACGCTCCTCAAGCATTTCATTGGCGGGGAGGGTGGCGAGCAGTTTTTGCGCTGCGTCGTACTGACCTTCTTTTAACAGCACCGGTAGCGTcgc                                                                                   <  2:1029470/151‑1 (MQ=255)
                                                           gCAATCGTGCCAGACGCAGAGCTTCAGCCTTGTTACGGGTCGCCACGCTGACAGCATAACGCTCCTCAAGCATTTCATTGGCGGGGAGGGTGGCGAGCAGTTTTTGCGCTGCGTCGTACTGACCTTCTTTTAACAGCACCGGTAGCGTcgc                                                                                   <  2:26938/151‑1 (MQ=255)
                                                            cAATCGTGCCAGACGCAGAGCTTCAGCCTTGTTACGGGTCGCCACGCTGACAGCATAACGCTCCTCAAGCATTTCATTGGCGGGGAGGGGGGGGGGCAGTTTTTTCGCTGGGTCGGAATGGACTTCTTTTTAAAACACCCGTAGGGTcgc                                                                                   >  2:797734/1‑150 (MQ=255)
                                                              aTCGTGCCAGACGCAGAGCTTCAGCCTTGTTACGGGTCGCCACGCTGACAGCATAACGCTCCTCAAGCATTTCATTGGCGGGGGGGGGGGCGAGCAGGTTTTTCGCTGCGTCGTACTGACCTTCTTTTTACAACACCCGTAGCGTCGCGcc                                                                                >  2:1545987/1‑151 (MQ=255)
                                                                 gTGCCAGACGCAGAGCTTCAGCCTTGTTACGGGTCGCCACGCTGACAGCATAACGCTCCTCAAGCATTTCATTGGCGGGGGGGGGGGGGGGAAGGTTTTTCGCTGCCGCGGTCTGGCCCTCTTTTTACAACACCGGGAGCGGCGCGccc                                                                               >  1:536987/1‑148 (MQ=255)
                                                                 gTGCCAGACGCAGAGCTTCAGCCTTGTTACGGGTCGCCACGCTGACAGCATAACGCTCCTCAAGCATTTCATTGGCGGGGGGGGGGGCGGGCAGGTTTTTCGCTGCGTCGTACTGCCCCTTTTTTTAAAGCACCGGGAGCGTCGCGcc                                                                                >  2:1158069/1‑148 (MQ=255)
                                                                         cGCAGAGCTTCAGCCTTGTTACGGGTCGCCACGCTGACAGCATAACGCTCCTCAAGCATTTCATTGGCGGGGAGGGTGGCGAGCAGTTTTTGCGCTGCGTCGTACTGACCTTCTTTTAACAGCACCGGTAGCGTCGCGCCAACAACATACt                                                                     <  1:222230/151‑1 (MQ=255)
                                                                          gCAGAGCTTCAGCCTTGTTACGGGTCGCCACGCTGACAGCATAACGCTCCTCAAGCATTTCATTGGCGGGGGGGGGGGGGGGGAGTTTTTTCGCTGCGGCGTGCTGACCCTCTTTTTACAGCACCCGGAGCGGCGCGCccacaa                                                                           >  1:962878/1‑144 (MQ=255)
                                                                           cAGAGCTTCAGCCTTGTTACGGGTCGCCACGCTGACAGCATAACGCTCCTCAAGCATTTCATTGGCGGGGAGGGTGGCGAGCAGTTTTTGCGCTGCGTCGTACTGACCTTCTTTTAACAGCACCGGTAGCGTCGCGCCAACAACATACTgg                                                                   <  2:1005078/151‑1 (MQ=255)
                                                                             gagCTTCAGCCTTGTTACGGGTCGCCACGCTGACAGCATAACGCTCCTCAAGCATTTCATTGGCGGGGGGGGGGGCGGGCAGTTTTTTCGCTGCGTCGTACTGACCTTCTTTTTACAGCACCGGGAGCGGCGCGCCCACAACATACTGGCg                                                                >  2:508275/1‑151 (MQ=255)
                                                                             gagCTTCAGCCTTGTTACGGGTCGCCACGCTGACAGCATAACGCTCCTCAAGCATTTCATTGGCGGGGAGGGTGGCGAGCAGTTTTTGCGCTGCGTCGTACTGACCTTCTTTTAACAGCACCGGTAGCGTCGCGCCAACAACATACTGGCg                                                                <  2:487548/151‑1 (MQ=255)
                                                                               gCTTCAGCCTTGTTACGGGTCGCCACGCTGACAGCATAACGCTCCTCAAGCATTTCATTGGCGGGGGGGGGGGCGGGGAGGTTTTTTGCTGGGGCGCAGTTACCCTCTTTTTACAACACCGGGAGCGTGGCGCCCACAACAAACTGGCgg                                                               >  2:1506606/1‑150 (MQ=255)
                                                                                cTTCAGCCTTGTTACGGGTCGCCACGCTGACAGCATAACGCTCCTCAAGCATTTCATTGGCGGGGGGGGGGGGGGGCAGGGTTTTTGCTGCGTCGTACTGTCCTTCTTTTTTAAGCAACCGTGGCGTCGCGCCAACAAAACACTTGCggg                                                              >  1:1059727/1‑149 (MQ=255)
                                                                                cTTCAGCCTTGTTACGGGTCGCCACGCTGACAGCATAACGCTCCTCAAGCATTTCATTGGCGGGGGGGGGGGGGAGCAGGTTTTTTGCCGCGTCGTGCCGTCCCTCTTTTTTAAACACCGGGAGGGTCGCGCCCAAAAAAAACTGGCggg                                                              >  1:1387663/1‑149 (MQ=255)
                                                                                 ttCAGCCTTGTTACGGGTCGCCACGCTGACAGCATAACGCTCCTCAAGCATTTCATTGGCGGGGAGGGGGGCGAGAAGTTTTTTCGCTGCGGCGTACTGACCTTCTTTTAACAACACCGGTAGCGGCGCGCCAAAAACATACTGGCgg                                                               >  2:224647/1‑148 (MQ=255)
                                                                                     gCCTTGTTACGGGTCGCCACGCTGACAGCATAACGCTCCTCAAGCATTTCATTGGCGGGGGGGGGGGGGGGCAGGTTTTTTGCTGCGGCGTACTGACCTTTTTTTTACAGCACCCGGAGCGGCGCGCCCAAAACAAACTGGCggg                                                              >  1:799154/1‑144 (MQ=255)
                                                                                     gCCTTGTTACGGGTCGCCACGCTGACAGCATAACGCTCCTCAAGCATTTCATTGGCGGGGGGGGGGGCGGGCAGGTTTTGCGCCGCGTCGGACTGACCTTCTTTTTACAACACCGGTAGCGTCGCGCCAACAACATACTGGCggg                                                              >  1:726678/1‑144 (MQ=255)
                                                                                     gCCTTGTTACGGGTCGCCACGCTGACAGCATAACGCTCCTCAAGCATTTCATTGGCGGGGAGGGCGGTGAGGAGGTTTTGAGCCGCGGCGGTCCGACCTTCTTTTATCAGACACGGTGGCGGTCCGCcaacacactaatgggg                                                                 >  1:13929/1‑132 (MQ=255)
                                                                                      ccTTTTTACGGGTCGCCACGCTGACAGCATAACGCTCCTCAAGCATTTCATTGGCGGGGGGGGGGGGGAGGCGTTTTTTCGCTGCCTCCTACTGACCTTCTTTTAACAGCACCGGTAGCGTCGCGCCCCCAACATACTGGCgg                                                               >  2:489567/1‑143 (MQ=255)
                                                                                       cTTGTTACGGGTCGCCACGCTGACAGCATAACGCTCCTCAAGCATTTCATTGGCGGGGAGGGTGGCGAGCAGTTGTTGCGCTGCGTCGTACTGACCTTCTTTTAACAGCACCGGTAGCGTCGCGCCAACAACATACTGGCGGTTGTCGGc                                                      <  1:1102495/150‑1 (MQ=255)
                                                                                               gggTCGCCACGCTGACAGCATAACGCTCCTCAAGCATTTCATTGGCGGGGGGGGGGGCGGGCAGGTTTTTTGCTGCGTCGTACTGACCTTCTTTTTAAAGCACCGGTAGCGTCGCGCCAACAACATAATGGGGGGTggtgggggaa                                                   >  2:1044864/1‑137 (MQ=255)
                                                                                                   cGCCACGCTGACAGCATAACGCTCCTCAAGCATTTCATTGGCGGGGAGGGTGGCGAGCAGTTTTTGCGCTGCGTCGTACTGACCTTCTTTTAACAGCACCGGTAGCGTCGCGCCAACAACATACTGGCGGTTGTCGGCAAACTGTACCGta                                         <  2:289377/151‑1 (MQ=255)
                                                                                                    gCCACGCTGACAGCATAACGCTCCTCAAGCATTTCATTGGCGGGGAGGGTGGCGAGCAGTTTTTGCGCTGCGTCGTACTGACCTTCTTTTAACAGCACCGGTAGCGTCGCGCCAACAACATACTGGCGGTTGTCGGCAAACTGTACCGtat                                        <  2:841413/151‑1 (MQ=255)
                                                                                                           tGACAGCATAACGCTCCTCAAGCATTTCATTGGCGGGGGGGGGGGGGGGCAGCTTTTTCGCTTGGTCGTATCTACCCTTTTTTTAAAACAACCGGGGCGGCGCGCCAACCACATACTGGCGGGGTGTTGGCaaaaattacgccgt                                        >  1:1434188/1‑134 (MQ=255)
                                                                                                            gACAGCATAGCGCTCCTCAAGCATTTCATTGGCGGGGAGGGTGGCGAGCAGTTTTTGCGCTGCGTCGTACTGACCTTCTTTTAACAGCACCGGTAGCGTCGCGCCAACAACATACTGGCGGTTGTCGGCAAACTGTACCGTATAATTCGc                                <  1:816270/150‑1 (MQ=255)
                                                                                                                gCATAACGCTCCTCAAGCATTTCATTGGCGGGGAGGGTGGCGAGCAGTTTTTGCGCTGCGTCGTACTGACCTTCTTTTAACAGCACCGGTAGCGTCGCGCCAACAACATACTGGCGGTTGTCGGCAAACTGTACCGTATAATTCGCCAACg                           <  2:1196488/151‑1 (MQ=255)
                                                                                                                 cATAACGCTCCTCAAGCATTTCATTGGCGGGGAGGGTGGCGAGCAGTTTTTGCGCTGCGTCGTACTGACCTTCTTTTAACAGCACCGGTAGCGTCGCGCCAACAACATACTGGCGGTTGTCGGCAAACTGTACCGTATAATTCGCCAACGc                          <  2:1413335/151‑1 (MQ=255)
                                                                                                                     aCGCTCCTCAAGCATTTCATTGGCGGGGAGGGTGGCGAGCAGTTTTTGCGCTGCGTCGTACTGACCTTCTTTTAACAGCACCGGTAGCGTCGCGCCAACAACATACTGGCGGTTGTCGGCAAACTGTACCGTATAATTCGCCAACGCCTGa                      <  1:406749/151‑1 (MQ=255)
                                                                                                                          cctcAAGCATTTCATTGGCGGGGAGGGTGGCGAGCAGTTTTTGCGCTGCGTCGTACTGACCTTCTTTTAACAGCACCGGTAGCGTCGCGCCAACAACATACTGGCGGTTGTCGGCAAACTGTACCGTATAATTCGCCAACGCCTGAACgg                  <  1:1506516/150‑1 (MQ=255)
                                                                                                                               aGCATTTCATTGGCGGGGAGGGTGGCGAGCAGTTTTTGCGCTGCGTCGTACTGACCTTCTTTTAACAGCACCGGTAGCGTCGCGCCAACAACATACTGGCGGTTGTCGGCAAACTGTACCGTATAATTCGCCAACGCCTGAACGGGGTTAg            <  2:818853/151‑1 (MQ=255)
                                                                                                                                     tCATTGGCGGGGAGGGTGGCGAGCAGTTTTTGCGCTGCGTCGTACTGACCTTCTTTTAACAGCACCGGTAGCGTCGCGCCAACAACATACTGGCGGTTGTCGGCAAACTGTACCGTATAATTCGCCAACGCCTGAACGGGGTTAGCGCTg       <  1:563973/150‑1 (MQ=255)
                                                                                                                                      cATTGGCGGGGAGGGTGGCGAGCAGTTTTTGCGCTGCGTCGTACTGACCTTCTTTTAACAGCACCGGTAGCGTCGCGCCAACAACATACTGGCGGTTGTCGGCAAACTGTACCGTATAATTCGCCAACGCCTGAACGGGGTTAGCGCTGt      <  1:91236/150‑1 (MQ=255)
                                                                                                                                      cATTGGCGGGGAGGGTGGCGAGCAGTTTTTGCGCTGCGTCGTACTGACCTTCTTTTAACAGCACCGGTAGCGTCGCGCCAACAACATACTGGCGGTTGTCGGCAAACTGTACCGTATAATTCGCCAACGCCTGAACGGGGTTAGCGCTGt      <  2:1228607/150‑1 (MQ=255)
                                                                                                                                           gCGGGGAGGGTGGCGAGCAGTTTTTGCGCTGCGTCGTACTGACCTTCTTTTAACAGCACCGGTAGCGTCGCGCCAACAACATACTGGCGGTTGTTCGGCAAATCTACCGTATAATTCACCCAACACCGGAAGgggggtgcggcggctttt  >  1:219174/1‑136 (MQ=255)
                                                                                                                                                     |                                                                                                                                               
CATCAGTTGCCAGGTTAGTTGATCCAGGCGGGTAAGATTTGCCGGTTCTTGCTGATACAGCAATCGTGCCAGACGCAGAGCTTCAGCCTTGTTACGGGTCGCCACGCTGACAGCATAACGCTCCTCAAGCATTTCATTGGCGGGGAGGGTGGCGAGCAGTTTTTGCGCTGCGTCGTACTGACCTTCTTTTAACAGCACCGGTAGCGTCGCGCCAACAACATACTGGCGGTTGTCGGCAAACTGTACCGTATAATTCGCCAACGCCTGAACGGGGTTAGCGCTGTATTT  >  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: