Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 284,450 A→G E417G (GAG→GGG)  yagF → CP4‑6 prophage; dehydratase family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913284,4500AG98.5% 224.0 / ‑5.0 65E417G (GAG→GGG) yagFCP4‑6 prophage; dehydratase family protein
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base G (34/30);  total (34/31)
Fisher's exact test for biased strand distribution p-value = 4.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

ATGCTCCACCTGCGCGACCTCGGCCTGCTGCATCTGGACGCCATGACCGTGACCGGCCAGACGGTGGGCGAGAACCTTGAATGGTGGCAGGCGTCCGAGCGCCGGGCGCGCTTCCGCCAGTGCCTGCGCGAGCAGGACGGCGTAGAGCCGGATGACGTGATCCTGCCGCCGGAGAAGGCAAAAGCGAAAGGGCTGACCTCGACGGTCTGCTTCCCGACGGGCAACATCGCTCCGGAAGGTTCGGTGATCAAGGCCACGGCGATCGACCCGTCGGTGGTGGGCGAAGATGGC  >  NC_000913/284305‑284595
                                                                                                                                                 |                                                                                                                                                 
aTGCTCCACCTGCGCGACCTCGGCCTGCTGCATCTGGACGCCATGACCGTGACCGGCCAGACGGTGGGCGAGAACCTTGAATGGTGGCAGGCGTCCGAGCGCCGGGCGCGCTTCCGCCAGTGCCTGCGCGAGCAGGACGGCGTAGGGCCgg                                                                                                                                              >  1:372359/1‑151 (MQ=255)
 tGCTCCACCTGCGCGACCTCGGCCTGCTGCATCTGGACGCCATGACCGTGACCGGCCAGACGGTGGGCGAGAACCTTGAATGGTGGCAGGCGTCCGAGCGCCGGGCGCGCTTCCGCCAGTGCCTGCGCGAGCAGGACGGCGTAGGGCCgg                                                                                                                                              >  1:189812/1‑150 (MQ=255)
  gCTCCACCTGCGCGACCTCGGCCTGCTGCATCTGGACGCCATGACCGTGACCGGCCAGACGGTGGGCGAGAACCTTGAATGGTGGCAGGCGTCCGAGCGCCGGGCGCGCTTCCGCCAGTGCCTGCGCGAGCAGGACGGCGTAGGGCCGGa                                                                                                                                             <  2:51951/150‑1 (MQ=255)
  gCTCCACCTGCGCGACCTCGGCCTGCTGCATCTGGACGCCATGACCGTGACCGGCCAGACGGTGGGCGAGAACCTTGAATGGTGGCAGGCGTCCGAGCGCCGGGCGCGCTTCCGCCAGTGCCTGCGCGAGCAGGACGGCGTAGAGCCGGAt                                                                                                                                            <  1:342384/151‑1 (MQ=255)
      cACCTGCGCGACCTCGGCCTGCTGCATCTGGACGCCATGACCGTGACCGGCCAGACGGTGGGCGAGAACCTTGAATGGTGGCAGGCGTCCGAGCGCCGGGCGCGCTTCCGCCAGTGCCTGCGCGAGCAGGACGGCGTAGGGCCGGATGACg                                                                                                                                        >  1:326621/1‑151 (MQ=255)
        ccTGCGCGACCTCGGCCTGCTGCATCTGGACGCCATGACCGTGACCGGCCAGACGGTGGGCGAGAACCTTGAATGGTGGCAGGCGTCCGAGCGCCGGGCGCGCTTCCGCCAGTGCCTGCGCGAGCAGGACGGCGTAGGGCCGGATGACGTg                                                                                                                                      >  1:239103/1‑151 (MQ=255)
          tGCGCGACCTCGGCCTGCTGCATCTGGACGCCATGACCGTGACCGGCCAGACGGTGGGCGAGAACCTTGAATGGTGGCAGGCGTCCGAGCGCCGGGCGCGCTTCCGCCAGTGCCTGCGCGAGCAGGACGGCGTAGGGCCGGATGACGTGAt                                                                                                                                    >  1:8647/1‑151 (MQ=255)
           gcgcGACCTCGGCCTGCTGCATCTGGACGCCATGACCGTGACCGGCCAGACGGTGGGCGAGAACCTTGAATGGTGGCAGGCGTCCGAGCGCCGGGCGCGCTTCCGCCAGTGCCTGCGCGAGCAGGACGGCGTAGGGCCGGATGACGTGa                                                                                                                                     <  2:472365/149‑1 (MQ=255)
               gACCTCGGCCTGCTGCATCTGGACGCCATGACCGTGACCGGCCAGACGGTGGGCGAGAACCTTGAATGGTGGCAGGCGTCCGAGCGCCGGGCGCGCTTCCGCCAGTGCCTGCGCGAGCAGGACGGCGTAGGGCCGGATGACGTGATCCTgc                                                                                                                               >  1:584581/1‑151 (MQ=255)
                 ccgcGGCCTGCTGCATCTGGACGCCATGACCGTGACCGGCCAGACGGTGGGCGAGAACCTTGAATGGTGGCAGGCGTCCGAGCGCCGGGCGCGCTTCCGCCAGTGCCTGCGCGAGCAGGACGGCGTAGGGCCGGATGACGTGATCCTgccg                                                                                                                             >  1:46177/4‑151 (MQ=255)
                 ccTCGGCCTGCTGCATCTGGACGCCATGACCGTGACCGGCCAGACGGTGGGCGAGAACCTTGAATGGTGGCAGGCGTCCGAGCGCCGGGCGCGCTTCCGCCAGTGCCTGCGCGAGCAGGACGGCGTAGGGCCGGATGACGTGATCCTgccg                                                                                                                             >  1:626456/1‑151 (MQ=255)
                     ggCCTGCTGCATCTGGACGCCATGACCGTGACCGGCCAGACGGTGGGCGAGAACCTTGAATGGTGGCAGGCGTCCGAGCGCCGGGCGCGCTTCCGCCAGTGCCTGCGCGAGCAGGACGGCGTAGGGCCGGATGACGTGATCCTGCCGCCgg                                                                                                                         >  1:319507/1‑151 (MQ=255)
                        ctgctgCATCTGGACGCCATGACCGTGACCGGCCAGACGGTGGGCGAGAACCTTGAATGGTGGCAGGCGTCCGAGCGCCGGGCGCGCTTCCGCCAGTGCCTGCGCGAGCAGGACGGCGTAGGGCCGGATGACGTGATCCTGCCGCCgg                                                                                                                         >  1:55260/1‑148 (MQ=255)
                          gctgcATCTGGACGCCATGACCGTGACCGGCCAGACGGTGGGCGAGAACCTTGAATGGTGGCAGGCGTCCGAGCGCCGGGCGCGCTTCCGCCAGTGCCTGCGCGAGCAGGACGGCGTAGGGCCGGATGACGTGATCCTGCCGCCGGAGaa                                                                                                                     <  2:126135/150‑1 (MQ=255)
                          gctgcATCTGGACGCCATGACCGTGACCGGCCAGACGGTGGGCGAGAACCTTGAATGGTGGCAGGCGTCCGAGCGCCGGGCGCGCTTCCGCCAGTGCCTGCGCGAGCAGGACGGCGTAGGGCCGGATGACGTGATCCTGCCGCCGGAGaa                                                                                                                     <  1:3820/150‑1 (MQ=255)
                           ctgcATCTGGACGCCATGACCGTGACCGGCCAGACGGTGGGCGAGAACCTTGAATGGTGGCAGGCGTCCGAGCGCCGGGCGCGCTTCCGCCAGTGCCTGCGCGAGCAGGACGGCGTAGGGCCGGATGACGTGATCCTGCCGCCGGAGAAgg                                                                                                                   >  1:200850/1‑151 (MQ=255)
                             gcATCTGGACGCCATGACCGTGACCGGCCAGACGGTGGGCGAGAACCTTGAATGGTGGCAGGCGTCCGAGCGCCGGGCGCGCTTCCGCCAGTGCCTGCGCGAGCAGGACGGCGTAGGGCCGGATGACGTGATCCTGCCGCCGGAGAAgg                                                                                                                   <  2:18194/149‑1 (MQ=255)
                             gcATCTGGACGCCATGACCGTGACCGGCCAGACGGTGGGCGAGAACCTTGAATGGTGGCAGGCGTCCGAGCGCCGGGCGCGCTTCCGCCAGTGCCTGCGCGAGCAGGACGGCGTAGGGCCGGATGACGTGATCCTGCCGCCGGAGAAGGc                                                                                                                  <  1:243734/150‑1 (MQ=255)
                                  tGGACGCCATGACCGTGACCGGCCAGACGGTGGGCGAGAACCTTGAATGGTGGCAGGCGTCCGAGCGCCGGGCGCGCTTCC‑CCAGTGCCTGCGCGAGCAGGACGGCGTAGGGCCGGATGACGTGATCCTGCCGCCGGAGAAGGCAAAAGCg                                                                                                           >  1:326870/1‑151 (MQ=255)
                                      cGCCATGACCGTGACCGGCCAGACGGTGTGCGAGAACCTTGAATGGTGGCAGGCGTCCGAGCGCCGGGCGCGCTTCCGCCAGTGCCTGCGCGAGCAGGACGGCGTAGGGCCGGATGACGTGATCCTGCCGCCGGAGAAGGCAAAAGCGaa                                                                                                         >  1:404851/1‑150 (MQ=255)
                                         cATGACCGTGACCGGCTAGACGGTGGGCGAGAACCTTGAATGGTGGCAGGCGTCCGAGCGCCGGGCGCGCTTCCGCCAGTGCCTGCGCGAGCAGGACGGCGTAGGGCCGGATGACGTGATCCTGCCGCCGGAGAAGGCAAAAGCGAAAggg                                                                                                     <  1:577028/151‑1 (MQ=255)
                                          aTGACCGTGACCGGCCAGACGGTGGGCGAGAACCTTGAATGGTGGCAGGCGTCCGAGCGCCGGGCGCGCTTCCGCCAGTGCCTGCGCGAGCAGGACGGCGTAGGGCCGGATGACGTGATCCTGCCGCCGGAGAAGGCAAAAGCGAAAGGGc                                                                                                    >  2:641978/1‑151 (MQ=255)
                                                gtgaccgGCCAGACGGTGGGCGAGAACCTTTAATGGTGGCAGGCGTCCGAGCGCCGGGCGCGCTTCCGCCAGTGCCTGCGCGAGCAGGACGGCGTAGGGCCGGATGACGTGATCCTGCCGCCGGAGAAGGCAAAAGCGAAAGGGCTGACCt                                                                                              <  2:36059/151‑1 (MQ=255)
                                                gtgaccgGCCAGACGGTGGGCGAGAACCTTGAATGGTGGCAGGCGTCCGAGCGCCGGGCGCGCTTCCGCCAGTGCCTGCGCGAGCAGGACGGCGTAGGGCCGGATGACGTGATCCTGCCGCCGGAGAAGGCAAAAGCGAAAGGGCTGACCt                                                                                              >  1:648994/1‑151 (MQ=255)
                                                  gaccgGTCAGACGGTGGGCGAGAACCTTGAATGGTGGCAGGCGTCCGAGCGCCGGGCGCGCTTCCGCCAGTGCCTGCGCGAGCAGGACGGCGTAGGGCCGGATGACGTGATCCTGCCGCCGGAGAAGGCAAAAGCGAAAGGGCTGACCTc                                                                                             <  2:59842/150‑1 (MQ=255)
                                                       gCCAGACGGTGGGCGAGAACCTTGAATGGTGGCAGGCGTCCGAGCGCCGGGCGCGCTTCCGCCAGTGCCTGCGCGAGCAGGACGGCGTAGGGCCGGATGACGTGATCCTGCCGCCGGAGAAGGCAAAAGCGAAAGGGCTGACCTCGACGGt                                                                                       >  2:369676/1‑151 (MQ=255)
                                                          aGACGGTGGGCGAGAACCTTGAATGGTGGCAGGCGTCCGAGCGCCGGGCGCGCTTCCGCCAGTGCCTGCGCGAGCAGGACGGCGTAGGGCCGGATGACGTGATCCTGCCGCCGGAGAAGGCAAAAGCGAAAGGGCTGACCTCGACGGTCTg                                                                                    >  2:281825/1‑151 (MQ=255)
                                                          aGACGGTGGGCGAGAACCTTGAATGGTGGCAGGCGTCCGAGCGCCGGGCGCGCTTCCGCCAGTGCCTGCGCGAGCAGGACGGCGTAGGGCCGGATGACGTGATCCTGCCGCCGGAGAAGGCAAAAGCGAAAGGGCTGACCTCGACGGTCTg                                                                                    >  1:494829/1‑151 (MQ=255)
                                                              ggTGGGCGAGAACCTTGAATGGTGGCAGGCGTCCGAGCGCCGGGCGCGCTTCCGCCAGTGCCTGCGCGAGCAGGACGGCGTAGGGCCGGATGACGTGATCCTGCCGCCGGAGAAGGCAAAAGCGAAAGGGCTGACCTCGACGGTCTGCTTc                                                                                <  2:263066/151‑1 (MQ=255)
                                                              ggTGGGCGAGAACCTTGAATGGTGGCAGGCGTCCGAGCGCCGGGCGCGCTTCCGCCAGTGCCTGCGCGAGCAGGACGGCGTAGGGCCGGATGACGTGATCCTGCCGCCGGAGAAGGCAAAAGCGAAAGGGCTGACCTCGACGGTCTGCTTc                                                                                <  1:513272/151‑1 (MQ=255)
                                                               gTGGGCGAGAACCTTGAATGGTGGCAGGCGTCCGAGCGCCGGGCGCGCTTCCGCCAGTGCCTGCGCGAGCAGGACGGCGTAGGGCCGGATGACGTGATCCTGCCGCCGGAGAAGGCAAAAGCGAAAGGGCTGACCTCGACGGTCTGCTTcc                                                                               <  2:36478/151‑1 (MQ=255)
                                                               gTGGGCGAGAACCTTGAATGGTGGCAGGCGTCCGAGCGCCGGGCGCGCTTCCGCCAGTGCCTGCGCGAGCAGGACGGCGTAGGGCCGGATGACGTGATCCTGCCGCCGGAGAAGGCAAAAGCGAAAGGGCTGACCTCGACGGTCTGCTTcc                                                                               >  1:16805/1‑151 (MQ=255)
                                                                tGGGCGAGAACCTTGAATGGTGGCAGGCGTCCGAGCGCCGGGCGCGCTTCCGCCAGTGCCTGCGCGAGCAGGACGGCGTAGGGCCGGATGACGTGATCCTGCCGCCGGAGAAGGCAAAAGCGAAAGGGCTGACCTCGACGGTCTGCTTccc                                                                              >  1:136202/1‑151 (MQ=255)
                                                                    cGAGAACCTTGAATGGTGGCAGGCGTCCGAGCGCCGGGCGCGCTTCCGCCAGTGCCTGCGCGAGCAGGACGGCGTAGGGCCGGATGACGTGATCCTGCCGCCGGAGAAGGCAAAAGCGAAAGGGCTGACCTCGACGGTCTGCTTCCCGACg                                                                          >  1:133742/1‑151 (MQ=255)
                                                                    cGAGAACCTTGAATGGTGGCAGGCGTCCGAGCGCCGGGCGCGCTTCCGCCAGTGCCTGCGCGAGCAGGACGGCGTAGGGCCGGATGACGTGATCCTGCCGCCGGAGAAGGCAAAAGCGAAAGGGCTGACCTCGACGGTCTGCTTCCCGACg                                                                          >  2:218659/1‑151 (MQ=255)
                                                                     gagaACCTTGAATGGTGGCAGGCGTCCGAGCGCCGGGCGCGCTTCCGCCAGTGCCTGCGCGAGCAGGACGGCGTAGGGCCGGATGACGTGATCCTGCCGCCGGAGAAGGCAAAAGCGAAAGGGCTGACCTCGACGGTCTGCTTCCCGACgg                                                                         >  1:29444/1‑151 (MQ=255)
                                                                         aCCTTGAATGGTGGCAGGCGTCCGAGCGCCGGGCGCGCTTCCGCCAGTGCCTGCGCGAGCAGGACGGCGTAGGGCCGGATGACGTGATCCTGCCGCCGGAGAAGGCAAAAGCGAAAGGGCTGACCTCGACGGTCTGCTTCCCGACGGGCa                                                                      >  2:183118/1‑150 (MQ=255)
                                                                             tGAATGGTGGCAGGCGTCCGAGCGCCGGGCGCGCTTCCGCCAGTGCCTGCGCGAGCAGGACGGCGTAGGGCCGGATGACGTGATCCTGCCGCCGGAGAAGGCAAAAGCGAAAGGGCTGACCTCGACGGTCTGCTTCCCGACGGGCAACa                                                                   <  1:453018/149‑1 (MQ=255)
                                                                              gAATGGTGGCAGGCGTCCGAGCGCCGGGCGCGCTTCCGCCAGTGCCTGCGCGAGCAGGACGGCGTAGGGCCGGATGACGTGATCCTGCCGCCGGAGAAGGCAAAAGCGAAAGGGCTGACCTCGACGGTCTGCTTCCCGACGGGCAACATCg                                                                <  1:640966/151‑1 (MQ=255)
                                                                               aaTGGTGGCAGGCGTCCGAGCGCCGGGCGCGCTTCCGCCAGTGCCTGCGCGAGCAGGACGGCGTAGGGCCGGATGACGTGATCCTGCCGCCGGAGAAGGCAAAAGCGAAAGGGCTGACCTCGACGGTCTGCTTCCCGACGGGCAACATCg                                                                >  1:453338/1‑150 (MQ=255)
                                                                                aTGGTGGCAGGCGTCCGAGCGCCGGGCGCGCTTCCGCCAGTGCCTGCGCGAGCAGGACGGCGTAGGGCCGGATGACGTGATCCTGCCGCCGGAGAAGGCAAAAGCGAAAGGGCTGACCTCGACGGTCTGCTTCCCGACGGGCAACATCGCt                                                              <  2:578152/151‑1 (MQ=255)
                                                                                  ggtggCAGGCGTCCGAGCGCCGGGCGCGCTTCCGCCAGTGCCTGCGCGAGCAGGACGGCGTAGGGCCGGATGACGTGATCCTGCCGCCGGAGAAGGCAAAAGCGAAAGGGCTGACCTCGACGGTCTGCTTCCCGACGGGCAACATCGCTcc                                                            <  1:367854/151‑1 (MQ=255)
                                                                                  gggggcaggcaTCCGAGGGCGGGGCGCGGTGCCGCCAGGGCCTGCGCGAGCAGGACGGCGTAGGGCCGGATGACGTGATCCTGCCGCCGGAGAAGGCAAAAGCGAAAGGGCAGACGTCGGCGGTCTGGTTCCCGACGGGCAAGATCGCTcc                                                            <  2:129898/148‑1 (MQ=255)
                                                                                       cAGGCGTCCGAGCGCCGGGCGCGCTTCCGCCAGTGCCTGCGCGAGCAGGACGGCGTAGGGCCGGATGACGTGATCCTGCCGCCGGAGAAGGCAAAAGCGAAAGGGCTGACCTCGACGGTCTGCTTCCCGACGGGCAACATCGCTCCGGa                                                         >  1:454867/1‑149 (MQ=255)
                                                                                       cAGGCGTCCGAGCGCCGGGCGCGCTTCCGCCAGTGCCTGCGCGAGCAGGACGGCGTAGGGCCGGATGACGTGATCCTGCCGCCGGAGAAGGCAAAAGCGAAAGGGCTGACCTCGACGGTCTGCTTCCCGACGGGCAACATCGCTCCGGa                                                         >  1:251750/1‑149 (MQ=255)
                                                                                          gCGTCCGAGCGCCGGGCGCGCTTCCGCCAGTGCCTGCGCGAGCAGGACGGCGTAGGGCCGGATGACGTGATCCTGCCGCCGGAGAAGGCAAAAGCGAAAGGGCTGACCTCGACGGTCTGCTTCCCGACGGGCAACATCGCTCCGGAAGGtt                                                    <  1:260933/151‑1 (MQ=255)
                                                                                               cGAGCGCCGGGCGCGCTTCCGCCAGTGCCTGCGCGAGCAGGACGGCGTAGGGCCGGATGACGTGATCCTGCCGCCGGAGAAGGCAAAAGCGAAAGGGCTGACCTCGACGGTCTGCTTCCCGACGGGCAACATCGCTCCGGAAGGTTCGGTg                                               >  2:215788/1‑151 (MQ=255)
                                                                                                       gggCGCGCTTCCGCCAGTGCCTGCGCGAGCAGGACGGCGTAGGGCCGGATGACGTGATCCTGCCGCCGGAGAAGGCAAAAGCGAAAGGGCTGACCTCGACGGTCTGCTTCCCGACGGGCAACATCGCTCCGGAAGGTTCGGTGATCAAgg                                        <  2:454855/150‑1 (MQ=255)
                                                                                                         gcgcgcTTCCGCCAGTGCCTGCGCGAGCAGGACGGCGTAGGGCCGGATGACGTGATCCTGCCGCCGGAGAAGGCAAAAGCGAAAGGGCTGACCTCGACGGTCTGCTTCCCGACGGGCAACATCGCTCCGGAAGGTTCGGTGATCAAGGcc                                      >  2:550735/1‑150 (MQ=255)
                                                                                                               ttCCGCCAGTGCCTGCGCGAGCAGGACGGCGTAGGGCCGGATGACGTGATCCTGCCGCCGGAGAAGGCAAAAGCGAAAGGGCTGACCTCGACGGTCTGCTTCCCGACGGGCAACATCGCTCCGGAAGGTTCGGTGATCAAGGCCACGGCg                                >  2:408844/1‑150 (MQ=255)
                                                                                                                tCCGCCAGTGCCTGCGCGAGCAGGACGGCGTAGGGCCGGATGACGTGATCCTGCCGCCGGAGAAGGCAAAAGCGAAAGGGCTGACCTCGACGGTCTGCTTCCCGACGGGCAACATCGCTCCGGAAGGTTCGGTGATCAAGGCCACGGCGa                               <  1:527008/150‑1 (MQ=255)
                                                                                                                 ccGCCAGTGCCTGCGCGAGCAGGACGGCGTAGGGCCGGATGACGTGATCCTGCCGCCGGAGAAGGCAAAAGCGAAAGGGCTGACCTCGACGGTCTGCTTCCCGACGGGCAACATCGCTCCGGAAGGTTCGGTGATCAAGGCCACGGc                                 <  2:569087/147‑1 (MQ=255)
                                                                                                                 ccGCCAGTGCCTGCGCGAGCAGGACGGCGTAGGGCCGGATGACGTGATCCTGCCGCCGGAGAAGGCAAAAGCGAAAGGGCTGACCTCGACGGTCTGCTTCCCGACGGGCAACATCGCTCCGGAAGGTTCGGTGATCAAGGCCACGGCGATc                             <  2:523231/151‑1 (MQ=255)
                                                                                                                 ccGCCAGTGCCGGCGCGAGCAGGACGGCGTAGGGCCGGATGACGTGATCCTGCCGCCGGAGAAGGCAAAAGCGAAAGGGCTGACCTCGACGGTCTGCTTCCCGACGGGCAACATCGCTCCGGAAGGTTCGGTGATCAAGGCCACGGCGAt                              <  2:650182/150‑1 (MQ=255)
                                                                                                                  cGCCAGTGCCTGCGCGAGCAGGACGGCGTAGGGCCGGATGACGTGATCCTGCCGCCGGAGAAGGCAAAAGCGAAAGGGCTGACCTCGACGGTCTGCTTCCCGACGGGCAACATCGCTCCGGAAGGTTCGGTGATCAAGGCCACGGCGATCg                            <  2:71332/151‑1 (MQ=255)
                                                                                                                      aGTGCCTGCGCGAGCAGGACGGCGTAGGGCCGGATGACGTGATCCTGCCGCCGGAGAAGGCAAAAGCGAAAGGGCTGACCTCGACGGTCTGCTTCCCGACGGGCAACATCGCTCCGGAAGGTTCGGTGATCAAGGCCACGGCGATCGAccc                        >  1:383048/1‑151 (MQ=255)
                                                                                                                         gCCTGCGCGAGCAGGACGGCGTAGGGCCGGATGACGTGATCCTGCCGCCGGAGAAGGCAAAAGCGAAAGGGCTGACCTCGACGGTCTGCTTCCCGACGGGCAACATCGCTCCGGAAGGTTCGGTGATCAAGGCCACGGCGATCGACCCGt                      <  1:85982/150‑1 (MQ=255)
                                                                                                                                cgAGCAGGACGGCGTAGGGCCGGATGACGTGATCCTGCCGCCGGAGAAGGCAAAAGCGAAAGGGCTGACCTCGACGGTCTGCTTCCCGACGGGCAACATCGCTCCGGAAGGTTCGGTGATCAAGGCCACGGCGATCGACCCGTCGgtggtg              >  1:145501/1‑151 (MQ=255)
                                                                                                                                  aGCAGGACGGCGTAGGGCCGGATGACGTGATCCTGCCGCCGGAGAAGGCAAAAGCGAAAGGGCTGACCTCGACGGTCTGCTTCCCGACGGGCAACATCGCTCCGGAAGGTTCGGTGATCAAGGCCACGGCGATCGACCCGTCGGTGGTggg            >  1:385698/1‑151 (MQ=255)
                                                                                                                                   gCAGGACGGCGTAGGGCCGGATGACGTGATCCTGCCGCCGGAGAAGGCAAAAGCGAAAGGGCTGACCTCGACGGTCTGCTTCCCGACGGGCAACATCGCTCCGGAAGGTTCGGTGATCAAGGCCACGGCGATCGACCCGTCGGTGGTggg            <  2:616681/150‑1 (MQ=255)
                                                                                                                                     aGGACGGCGTAGGGCCGGATGACGTGATCCTGCCGCCGGAGAAGGCAAAAGCGAAAGGGCTGACCTCGACGGTCTGCTTCCCGACGGGCAACATCGCTCCGGAAGGTTCGGTGATCAAGGCCACGGCGATCGACCCGTCGGTGGTGGGCGa         <  1:429533/151‑1 (MQ=255)
                                                                                                                                      ggACGGCGTAGGGCCGGATGACGTGATCCTGCCGCCGGAGAAGGCAAAAGCGAAAGGGCTGACCTCGACGGTCTGCTTCCCGACGGGCAACATCGCTCCGGAAGGTTCGGTGATCAAGGCCACGGCGATCGACCCGTCGGTGGTGGGCGaa        <  2:267864/151‑1 (MQ=255)
                                                                                                                                        aCGGCGTAGGGCCGGATGACGTGATCCTGCCGCCGGAGAAGGCAAAAGCGAAAGGGCTGACCTCGACGGTCTGCTTCCCGACGGGCAACATCGCTCCGGAAGGTTCGGTGATCAAGGCCACGGCGATCGACCCGTCGGTGGTGGGCGAAGa      <  2:153966/151‑1 (MQ=255)
                                                                                                                                          ggCGTAGGGCCGGATGACGTGATCCTGCCGCCGGAGAAGGCAAAAGCGAAAGGGCTGACCTCGACGGTCTGCTTCCCGACGGGCAACATCGCTCCGGAAGGTTCGGTGATCAAGGCCACGGCGATCGACCCGTCGGTGGTGGGCGAAGATg    >  1:295559/1‑151 (MQ=255)
                                                                                                                                            cGTAGGGCCGGATGACGTGATCCTGCCGCCGGAGAAGGCAAAAGCGAAAGGGCTGACCTCGACGGTCTGCTTCCCGACGGGCAACATCGCTCCGGAAGGTTCGGTGATCAAGGCCACGGCGATCGACCCGTCGGTGGTGGGCGAAGATGGc  <  1:482548/151‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                 
ATGCTCCACCTGCGCGACCTCGGCCTGCTGCATCTGGACGCCATGACCGTGACCGGCCAGACGGTGGGCGAGAACCTTGAATGGTGGCAGGCGTCCGAGCGCCGGGCGCGCTTCCGCCAGTGCCTGCGCGAGCAGGACGGCGTAGAGCCGGATGACGTGATCCTGCCGCCGGAGAAGGCAAAAGCGAAAGGGCTGACCTCGACGGTCTGCTTCCCGACGGGCAACATCGCTCCGGAAGGTTCGGTGATCAAGGCCACGGCGATCGACCCGTCGGTGGTGGGCGAAGATGGC  >  NC_000913/284305‑284595

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