Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 265,090 A→G A39A (GCT→GCC yafX ← CP4‑6 prophage; uncharacterized protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913265,0900AG98.9% 321.5 / NA 90A39A (GCT→GCCyafXCP4‑6 prophage; uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base G (50/39);  minor base T (0/1);  total (50/40)
Fisher's exact test for biased strand distribution p-value = 4.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCATGTCTCATCGTTGTCAGGCTCGGGGGACATAAATGCGCCGCCATTGCTGAGGGTGTAGAACGACCAGATACCACCGCTGTAGCCCTCACAGAAGCGGTCCATCCAGGCGAAGATATGCGGCTCCAGGAGTAGCCACTGCGGGATAGCGCCAAAGTACTGTGGCCAGAAATCGATACGCTGTTCATCGGGGACCGGCGTGACGGTCAGTTCAAATTCGGGTTGGTTAGCGGGTGCGAGGTCGTGCTGCGTCTGTGTTGTCATGGGTATGTCTCCGTCAATAAAAACGCCAGCG  >  NC_000913/264943‑265237
                                                                                                                                                   |                                                                                                                                                   
ccATGTCTCATCGTTGTCAGGCTCGGGGGACATAAATGCGCCGCCATTGCTGAGGGTGTAGAACGACCAGATACCACCGCTGTAGCCCTCACAGAAGCGGTCCATCCAGGCGAAGATATGCGGCTCCAGGAGTAGCCACTGCGGGATGgcg                                                                                                                                                  <  2:609722/151‑1 (MQ=255)
   tGTCTCATCGTTGTCAGGCTCGGGGGACATAAATGCGCCGCCATTGCTGAGGGTGTAGAACGACCAGATACCACCGCTGTAGCCCTCACAGAAGCGGTCCATCCAGGCGAAGATATGCGGCTCCAGGAGTAGCCACTGCGGGTTGgcgc                                                                                                                                                 >  1:1132827/1‑149 (MQ=255)
   tGTCTCATCGTTGTCAGGCTCGGGGGACATAAATGCGCCGCCATTGCTGAGGGTGTAGAACGACCAGATACCACCGCTGTAGCCCTCACAGAAGCGGTCCATCCAGGCGAAGATATGCGGCTCCAGGAGTAGCCACTGCGGGATGGCGcc                                                                                                                                                >  2:497746/1‑150 (MQ=255)
    gTCTCATCGTTGTCAGGCTCGGGGGACATAAATGCGCCGCCATTGCTGAGGGTGTAGAACGACCAGATACCACCGCTGTAGCCCTCACAGAAGCGGTCCATCCAGGCGAAGATATGCGGCTCCAGGAGTAGCCACTGCGGGATGGCGCCa                                                                                                                                               >  2:985918/1‑150 (MQ=255)
    gTCTCATCGTTGTCAGGCTCGGGGGACATAAATGCGCCGCCATTGCTGAGGGTGTAGAACGACCAGATACCACCGCTGTAGCCCTCACAGAAGCGGTCCATCCAGGCGAAGATATGCGGCTCCAGGAGTAGCCACTGCGGGATGGCGCCa                                                                                                                                               <  2:1421051/150‑1 (MQ=255)
        cATCGTTGTCAGGCTCGGGGGACATAAATGCGCCGCCATTGCTGAGGGTGTAGAACGACCAGATACCACCGCTGTAGCCCTCACAGAAGCGGTCCATCCAGGCGAAGATATGCGGCTCCAGGAGTAGCCACTGCGGGATGGCGCCAAAg                                                                                                                                            <  1:115840/149‑1 (MQ=255)
         aTCGTTGTCAGGCTCGGGGGACATAAATGCGCCGCCATTGCTGAGGGTGTAGAACGACCAGATACCACCGCTGTAGCCCTCACAGAAGCGGTCCATCCAGGCGAAGATATGCGGCTCCAGGAGTAGCCACTGCGGGATGGCGCCaaa                                                                                                                                             <  2:1337820/147‑1 (MQ=255)
         aTCGTTGTCAGGCTCGGGGGACATAAATGCGCCGCCATTGCTGAGGGTGTAGAACGACCAGATACCACCGCTGTAGCCCTCACAGAAGCGGTCCATCCAGGCGAAGATATGCGGCTCCAGGAGTAGCCACTGCGGGATGGCGCCAAAGTAc                                                                                                                                         >  2:1349026/1‑151 (MQ=255)
              tGTCAGGCTCGGGGGACATAAATGCGCCGCCATTGCTGAGGGTGTAGAACGACCAGATACCACCGCTGTAGCCCTCACAGAAGCGGTCCATCCAGGCGAAGATATGCGGCTCCAGGAGTAGCCACTGCGGGATGGCGCCAAAGTACtgtg                                                                                                                                     <  1:84129/150‑1 (MQ=255)
                 cAGGCTCGGGGGACATAAATGCGCCGCCATTGCTGAGGGTGTAGAACGACCAGATACCACCGCTGTAGCCCTCACAGAAGCGGTCCATCCAGGCGAAGATATGCGGCTCCAGGAGTAGCCACTGCGGGATGGCGCCAAAGTACTGTGGcc                                                                                                                                  >  1:942378/1‑150 (MQ=255)
                 cAGGCTCGGGGGACATAAATGCGCCGCCATTGCTGAGGGTGTAGAACGACCAGATACCACCGCTGTAGCCCTCACAGAAGCGGTCCATCCAGGCGAAGATATGCGGCTCCAGGAGTAGCCACTGCGGGATGGCGCCAAAGTACTGTGGcc                                                                                                                                  >  1:570403/1‑150 (MQ=255)
                  aGGCTCGGGGGACATAAATGCGCCGCCATTGCTGAGGGTGTAGAACGACCAGATACCACCGCTGTAGCCCTCACAGAAGCGGTCCATCCAGGCGAAGATATGCGGCTCCAGGAGTAGCCACTGCGGGATGGCGCCAAAGTACTGTGGcc                                                                                                                                  >  2:355454/1‑149 (MQ=255)
                  aGGCTCGGGGGACATAAATGCGCCGCCATTGCTGAGGGTGTAGAACGACCAGATACCACCGCTGTAGCCCTCACAGAAGCGGTCCATCCAGGCGAAGATATGCGGCTCCAGGAGTAGCCACTGCGGGATGGCGCCAAAGTACTGTGGCCAg                                                                                                                                <  1:1034356/151‑1 (MQ=255)
                   ggCTCGGGGGACATAAATGCGCCGCCATTGCTGAGGGTGTAGAACGACCAGATACCACCGCTGTAGCCCTCACAGAAGCGGTCCATCCAGGCGAAGATATGCGGCTCCAGGAGTAGCCACTGCGGGATGGCGCCAAAGTACTGTGGCCAGa                                                                                                                               <  2:683130/151‑1 (MQ=255)
                    gCTCGGGGGACATAAATGCGCCGCCATTGCTGAGGGTGTAGAACGACCAGATACCACCGCTGTAGCCCTCACAGAAGCGGTCCATCCAGGCGAAGATATGCGGCTCCAGGAGTAGCCACTGCGGGATGGCGCCAAAGTACTGTGGCCAGa                                                                                                                               >  2:78649/1‑150 (MQ=255)
                      tCGGGGGACATAAATGCGCCGCCATTGCTGAGGGTGTAGAACGACCAGATACCACCGCTGTAGCCCTCACAGAAGCGGTCCATCCAGGCGAAGATATGCGGCTCCAGGAGTAGCCACTGCGGGATGGCGCCAAAGTACTGTGGCCAGAAAt                                                                                                                            >  2:413854/1‑151 (MQ=255)
                        gggggACATAAATGCGCCGCCATTGCTGAGGGTGTAGAACGACCAGATACCACCGCTGTAGCCCTCACAGAAGCGGTCCATCCAGGCGAAGATATGCGGCTCCAGGAGTAGCCACTGCGGGATGGCGCCAAAGTACTGTGGCCAGaaa                                                                                                                             <  1:1101721/148‑1 (MQ=255)
                         ggggACATAAATGCGCCGCCATTGCTGAGGGTGTAGAACGACCAGATACCACCGCTGTAGCCCTCACAGAAGCGGTCCATCCAGGCGAAGATATGCGGCTCCAGGAGTAGCCACTGCGGGATGGCGCCAAAGTACTGTGGCCAGAAATCg                                                                                                                          >  2:981219/1‑150 (MQ=255)
                            gACATAAATGCGCCGCCATTGCTGAGGGTGTAGAACGACCAGATACCACCGCTGTAGCCCTCACAGAAGCGGTCCATCCAGGCGAAGATATGCGGCTCCAGGAGTAGCCACTGCGGGATGGCGCCAAAGTACTGTGGCCAGAAATCGATAc                                                                                                                      >  2:1142577/1‑151 (MQ=255)
                                tAAATGCGCCGCCATTGCTGAGGGTGTAGAACGACCAGATACCACCGCTGTAGCCCTCACAGAAGCGGTCCATCCAGGCGAAGATATGCGGCTCCAGGAGTAGCCACTGCGGGATGGCGCCAAAGTACTGTGGCCAGAAATCGATACGCTg                                                                                                                  >  1:1339159/1‑151 (MQ=255)
                                tAAATGCGCCGCCATTGCTGAGGGTGTAGAACGACCAGATACCACCGCTGTAGCCCTCACAGAAGCGGTCCATCCAGGCGAAGATATGCGGCTCCAGGAGTAGCCACTGCGGGATGGCGCCAAAGTACTGTGGCCAGAAATCGATACGCTg                                                                                                                  >  1:763058/1‑151 (MQ=255)
                                 aaaTGCGCCGCCATTGCTGAGGGTGTAGAACGACCAGATACCACCGCTGTAGCCCTCACAGAAGCGGTCCATCCAGGCGAAGATATGCGGCTCCAGGAGTAGCCACTGCGGGATGGCGCCAAAGTACTGTGGCCAGAAATCGATACGCTGt                                                                                                                 >  1:616521/1‑151 (MQ=255)
                                  aaTGCGCCGCCATTGCTGAGGGTGTAGAACGACCAGATACCACCGCTGTAGCCCTCACAGAAGCGGTCCATCCAGGCGAAGATATGCGGCTCCAGGAGTAGCCACTGCGGGATGGCGCCAAAGTACTGTGGCCAGAAATCGATACGCTGtt                                                                                                                >  1:727450/1‑151 (MQ=255)
                                  aaTGCGCCGCCATTGCTGAGGGTGTAGAACGACCAGATACCACCGCTGTAGCCCTCACAGAAGCGGTCCATCCAGGCGAAGATATGCGGCTCCAGGAGTAGCCACTGCGGGATGGCGCCAAAGTACTGTGGCCAGAAATCGATACGCTGtt                                                                                                                <  1:1211107/151‑1 (MQ=255)
                                  aaTGCGCCGCCATTGCTGAGGGTGTAGAACGACCAGATACCACCGCTGTAGCCCTCACAGAAGCGGTCCATCCAGGCGAAGATATGCGGCTCCAGGAGTAGCCACTGCGGGATGGCGCCAAAGTACTGTGGCCAGAAATCGATACGCTGtt                                                                                                                <  2:973122/151‑1 (MQ=255)
                                     gcgcCGCCATTGCTGAGGGTGTAGAACGACCAGATACCACCGCTGTAGCCCTCACAGAAGCGGTCCATCCAGGCGAAGATATGCGGCTCCAGGAGTAGCCACTGCGGGATGGCGCCAAAGTACTGTGGCCAGAAATCGATACGCTGTTCa                                                                                                              <  2:927656/150‑1 (MQ=255)
                                     gcgcCGCCATTGCTGAGGGTGTAGAACGACCAGATACCACCGCTGTAGCCCTCACAGAAGCGGTCCATCCAGGCGAAGATATGCGGCTCCAGGAGTAGCCACTGCGGGATGGCGCCAAAGTACTGTGGCCAGAAATCGATACGCTGTTCAt                                                                                                             >  2:90489/1‑151 (MQ=255)
                                         cgccATTGCTGAGGGTGTAGAACGACCAGATACCACCGCTGTAGCCCTCACAGAAGCGGTCCATCCAGGCGAAGATATGCGGCTCCAGGAGTAGCCACTGCGGGATGGCGCCAAAGTACTGTGGCCAGAAATCGATACGCTGTTCATCggg                                                                                                         <  1:1344959/151‑1 (MQ=255)
                                           ccATTGCTGAGGGTGTAGAACGACCAGATACCACCGCTGTAGCCCTCACAGAAGCGGTCCATCCAGGCGAAGATATGCGGCTCCAGGAGTAGCCACTGCGGGATGGCGCCAAAGTACTGTGGCCAGAAATCGATACGCTGTTCa                                                                                                              >  1:1342817/1‑144 (MQ=255)
                                           ccATTGCTGAGGGTGTAGAACGACCAGATACCACCGCTGTAGCCCTCACAGAAGCGGTCCATCCAGGCGAAGATATGCGGCTCCAGGAGTAGCCACTGCGGGATGGCGCCAAAGTACTGTGGCCAGAAATCGATACGCTGTTCATCgggg                                                                                                        >  1:444824/1‑150 (MQ=255)
                                              ttGCTGAGGGTGTAGAACGACCAGATACCACCGCTGTAGCCCTCACAGAAGCGGTCCATCCAGGCGAAGATATGCGGCTCCAGGAGTAGCCACTGCGGGATGGCGCCAAAGTACTGTGGCCAGAAATCGATACGCTGTTCATCGGGGACCg                                                                                                    >  1:395271/1‑151 (MQ=255)
                                               tGCTGAGGGTGTAGAACGACCAGATACCACCGCTGTAGCCCTCACAGAAGCGGTCCATCCAGGCGAAGATATGCGGCTCCAGGAGTAGCCACTGCGGGATGGCGCCAAAGTACTGTGGCCAGAAATCGATACGCTGTTCATCGGGGACCgg                                                                                                   >  2:1084997/1‑151 (MQ=255)
                                               tGCTGAGGGTGTAGAACGACCAGATACCACCGCTGTAGCCCTCACAGAAGCGGTCCATCCAGGCGAAGATATGCGGCTCCAGGAGTAGCCACTGCGGGATGGCGCCAAAGTACTGTGGCCAGAAATCGATACGCTGTTCATCGGGGACCg                                                                                                    <  1:815316/150‑1 (MQ=255)
                                                  tGAGGGTGTAGAACGACCAGATACCACCGCTGTAGCCCTCACAGAAGCGGTCCATCCAGGCGAAGATATGCGGCTCCAGGAGTAGCCACTGCGGGATGGCGCCAAAGTACTGTGGCCAGAAATCGATACGCTGTTCATCGGGGACCGGCg                                                                                                 >  2:516695/1‑150 (MQ=255)
                                                  tGAGGGTGTAGAACGACCAGATACCACCGCTGTAGCCCTCACAGAAGCGGTCCATCCAGGCGAAGATATGCGGCTCCAGGAGTAGCCACTGCGGGATGGCGCCAAAGTACTGTGGCCAGAAATCGATACGCTGTTCATCGGGGACCGGCGt                                                                                                >  1:1349512/1‑151 (MQ=255)
                                                   gAGGGTGTAGAACGACCAGATACCACCGCTGTAGCCCTCACAGAAGCGGTCCATCCAGGCGAAGATATGCGGCTCCAGGAGTAGCCACTGCGGGATGGCGCCAAAGTACTGTGGCCAGAAATCGATACGCTGTTCATCGGGGACCGGCGTg                                                                                               >  2:349494/1‑151 (MQ=255)
                                                   gAGGGTGTAGAACGACCAGATACCACCGCTGTAGCCCTCACAGAAGCGGTCCATCCAGGCGAAGATATGCGGCTCCAGGAGTAGCCACTGCGGGATGGCGCCAAAGTACTGTGGCCAGAAATCGATACGCTGTTCATCGGGGACCGGCGTg                                                                                               <  1:453692/151‑1 (MQ=255)
                                                    aGGGTGTAGAACGACCAGATACCACCGCTGTAGCCCTCACAGAAGCGGTCCATCCAGGCGAAGATATGCGGCTCCAGGAGTAGCCACTGCGGGATGGCGCCAAAGTACTGTGGCCAGAAATCGATACGCTGTTCATCGGGGACCGGCGTg                                                                                               >  1:857373/1‑150 (MQ=255)
                                                     gggTGTAGAACGACCAGATACCACCGCTGTAGCCCTCACAGAAGCGGTCCATCCAGGCGAAGATATGCGGCTCCAGGAGTAGCCACTGCGGGATGGCGCCAAAGTACTGTGGCCAGAAATCGATACGCTGTTCATCGGGGACCGGCGTGAc                                                                                             <  1:834031/151‑1 (MQ=255)
                                                       gtgtAGAACGACCAGATACCACCGCTGTAGCCCTCTCAGAAGCGGTCCATCCAGGCGAAGATATGCGGCTCCATGAGTAGCCACTGCGGGATGGCGCCAAAGTACTGTGTCCAGAAATCGATACGATGTTCATCGGGGACCGGCGTGACg                                                                                            <  2:1161176/150‑1 (MQ=255)
                                                       gtgtAGAACGACCAGATACCACCGCTGTAGCCCTCACAGAAGCGGTCCATCCAGGCGAAGATATGCGGCTCCAGGAGTAGCCACTGCGGGATGGCGCCAAAGTACTGTGGCCAGAAATCGATACGCTGTTCATCGGGGACCGGCGTGACgg                                                                                           <  2:112195/151‑1 (MQ=255)
                                                        tgtAGAACGACCAGATACCACCGCTGTAGCCCTCACAGAAGCGGTCCATCCAGGCGAAGATATGCGGCTCCAGGAGTAGCCACTGCGGGATGGCGCCAAAGTACTGTGGCCAGAAATCGATACGCTGTTCATCGGGGACCGGCGTGACGGt                                                                                          >  1:671650/1‑151 (MQ=255)
                                                         gtAGAACGACCAGATACCACCGCTGTAGCCCTCACAGAAGCGGTCCATCCAGGCGAAGATATGCGGCTCCAGGAGTAGCCACTGCGGGATGGCGCCAAAGTACTGTGGCCAGAAATCGATACGCTGTTCATCGGGGACCGGCGTGACGGTc                                                                                         <  2:416390/151‑1 (MQ=255)
                                                            gAACGACCAGATACCACCGCTGTAGCCCTCACAGAAGCGGTCCATCCAGGCGAAGATATGCGGCTCCAGGAGTAGCCACTGCGGGATGGCGCCAAAGTACTGTGGCCAGAAATCGATACGCTGTTCATCGGGGACCGGCGTGACGGTCAGt                                                                                      >  1:929064/1‑151 (MQ=255)
                                                                gACCAGATACCACCGCTGTAGCCCTCACAGAAGCGGTCCATCCAGGCGAAGATATGCGGCTCCAGGAGTAGCCACTGCGGGATGGCGCCAAAGTACTGTGGCCAGAAATCGATACGCTGTTCATCGGGGACCGGCGTGACGGTCAGTTc                                                                                    >  1:171297/1‑149 (MQ=255)
                                                                gACCAGATACCACCGCTGTAGCCCTCACAGAAGCGGTCCATCCAGGCGAAGATATGCGGCTCCAGGAGTAGCCACTGCGGGATGGCGCCAAAGTACTGTGGCCAGAAATCGATACGCTGTTCATCGGGGACCGGCGTGACGGTCAGTTCa                                                                                   >  1:1407265/1‑150 (MQ=255)
                                                                  ccAGATACCACCGCTGTAGCCCTCACAGAAGCGGTCCATCCAGGCGAAGATATGCGGCTCCAGGAGTAGCCACTGCGGGATGGCGCCAAAGTACTGTGGCCAGAAATCGATACGCTGTTCATCGGGGACCGGCGTGACGGTCAGTTCAAAt                                                                                >  2:699284/1‑151 (MQ=255)
                                                                   cAGATACCACCGCTGTAGCCCTCACAGAAGCGGTCCATCCAGGCGAAGATATGCGGCTCCAGGAGTAGCCACTGCGGGATGGCGCCAAAGTACTGTGGCCAGAAATCGATACGCTGTTCATCGGGGACCGGCGTGACGGTCAGTTCAAAtt                                                                               <  1:317101/151‑1 (MQ=255)
                                                                    aGATACCACCGCTGTAGCCCTCACAGAAGCGGTCCATCCAGGCGAAGATATGCGGCTCCAGGAGTAGCCACTGCGGGATGGCGCCAAAGTACTGTGGCCAGAAATCGATACGCTGTTCATCGGGGACCGGCGTGACGGTCAGTTCAAAtt                                                                               >  2:839079/1‑150 (MQ=255)
                                                                    aGATACCACCGCTGTAGCCCTCACAGAAGCGGTCCATCCAGGCGAAGATATGCGGCTCCAGGAGTAGCCACTGCGGGATGGCGCCAAAGTACTGTGGCCAGAAATCGATACGCTGTTCATCGGGGACCGGCGTGACGGTCAGTTCAAATTc                                                                              >  2:459371/1‑151 (MQ=255)
                                                                      aTACCACCGCTGTAGCCCTCACAGAAGCGGTCCATCCAGGCGAAGATATGCGGCTCCAGGAGTAGCCACTGCGGGATGGCGCCAAAGTACTGTGGCCAGAAATCGATACGCTGTTCATCGGGGACCGGCGTGACGGTCAGTTCAAATTCgg                                                                            <  1:1316949/151‑1 (MQ=255)
                                                                         ccaccGCTGTAGCCCTCACAGAAGCGGTCCATCCAGGCGAAGATATGCGGCTCCAGGAGTAGCCACTGCGGGATGGCGCCAAAGTACTGTGGCCAGAAATCGATACGCTGTTCATCGGGGACCGGCGTGACGGTCAGTTCAAATTCGggtt                                                                         >  1:1014756/1‑151 (MQ=255)
                                                                             cGCTGTAGCCCTCACAGAAGCGGTCCATCCAGGCGAAGATATGCGGCTCCAGGAGTAGCCACTGCGGGATGGCGCCAAAGTACTGTGGCCAGAAATCGATACGCTGTTCATCGGGGACCGGCGTGACGGTCAGTTCAAATTCGggttggtt                                                                     <  2:823985/151‑1 (MQ=255)
                                                                              gCTGTAGCCCTCACAGAAGCGGTCCATCCAGGCGAAGATATGCGGCTCCAGGAGTAGCCACTGCGGGATGGCGCCAAAGTACTGTGGCCAGAAATCGATACGCTGTTCATCGGGGACCGGCGTGACGGTCAGTTCAAATTCGggttggtt                                                                     >  2:482381/1‑150 (MQ=255)
                                                                                tGTAGCCCTCACAGAAGCGGTCCATCCAGGCGAAGATATGCGGCTCCAGGAGTAGCCACTGCGGGATGGCGCCAAAGTACTGTGGCCAGAAATCGATACGCTGTTCATCGGGGACCGGCGTGACGGTCAGTTCAAATTCGGGTTGGTTAGc                                                                  >  1:916240/1‑151 (MQ=255)
                                                                                 gTAGCCCTCACAGAAGCGGTCCATCCAGGCGAAGATATGCGGCTCCAGGAGTAGCCACTGCGGGATGGCGCCAAAGTACTGTGGCCAGAAATCGATACGCTGTTCATCGGGGACCGGCGTGACGGTCAGTTCAAATTCGGGTTGGTTAGc                                                                  <  2:171275/150‑1 (MQ=255)
                                                                                  tAGCCCTCACAGAAGCGGTCCATCCAGGCGAAGATATGCGGCTCCAGGAGTAGCCACTGCGGGATGGCGCCAAAGTACTGTGGCCAGAAATCGATACGCTGTTCATCGGGGACCGGCGTGACGGTCAGTTCAAATTCGGGTTGGTTAGCg                                                                 <  1:314124/150‑1 (MQ=255)
                                                                                   aGCCCTCACAGAAGCGGTCCATCCAGGCGAAGATATGCGGCTCCAGGAGTAGCCACTGCGGGATGGCGCCAAAGTACTGTGGCCAGAAATCGATACGCTGTTCATCGGGGACCGGCGTGACGGTCAGTTCAAATTCGGGTTGGTTAGCggg                                                               >  2:191476/1‑151 (MQ=255)
                                                                                       cTCACAGAAGCGGTCCATCCAGGCGAAGATATGCGGCTCCAGGAGTAGCCACTGCGGGATGGCGCCAAAGTACTGTGGCCAGAAATCGATACGCTGTTCATCGGGGACCGGCGTGACGGTCAGTTCAAATTCGGGTTGGTTAGCGGGTGCg                                                           >  1:406384/1‑151 (MQ=255)
                                                                                       cTCACAGAAGCGGTCCATCCAGGCGAAGATATGCGGCTCCAGGAGTAGCCACTGCGGGATGGCGCCAAAGTACTGTGGCCAGAAATCGATACGCTGTTCATCGGGGACCGGCGTGACGGTCAGTTCAAATTCGGGTTGGTTAGCGGGTGCg                                                           >  2:730078/1‑151 (MQ=255)
                                                                                         cacaGAAGCGGTCCATCCAGGCGAAGATATGCGGCTCCAGGAGTAGCCACTGCGGGATGGCGCCAAAGTACTGTGGCCAGAAATCGATACGCTGTTCATCGGGGACCGGCGTGACGGTCAGTTCAAATTCGGGTTGGTTAGCGGGTGCGAg                                                         <  2:588808/151‑1 (MQ=255)
                                                                                         cacaGAAGCGGTCCATCCAGGCGAAGATATGCGGCTCCAGGAGTAGCCACTGCGGGATGGCGCCAAAGTACTGTGGCCAGAAATCGATACGCTGTTCATCGGGGACCGGCGTGACGGTCAGTTCAAATTCGGGTTGGTTAGCGGGTGCGAg                                                         <  1:1199395/151‑1 (MQ=255)
                                                                                          acaGAAGCGGTCCATCCAGGCGAAGATATGCGGCTCCAGGAGTAGCCACTGCGGGATGGCGCCAAAGTACTGTGGCCAGAAATCGATACGCTGTTCATCGGGGACCGGCGTGACGGTCAGTTCAAATTCGGGTTGGTTAGCGGGTGCGAgg                                                        <  1:1295907/151‑1 (MQ=255)
                                                                                               aGCGGTCCATCCAGGCGAAGATATGCGGCTCCAGGAGTAGCCACTGCGGGATGGCGCCAAAGTACTGTGGCCAGAAATCGATACGCTGTTCATCGGGGACCGGCGTGACGGTCAGTTCAAATTCGGGTTGGTTAGCGGGTGCGAGGTCGtg                                                   >  2:92874/1‑151 (MQ=255)
                                                                                                gCGGTCCATCCAGGCGAAGATATGCGGCTCCAGGAGTAGCCACTGCGGGATGGCGCCAAAGTACTGTGGCCAGAAATCGATACGCTGTTCATCGGGGACCGGCGTGACGGTCAGTTCAAATTCGGGTTGGTTAGCGGGTGCGAGGTCGtgc                                                  >  1:703255/1‑151 (MQ=255)
                                                                                                      catccaGGCGAAGATATGCGGCTCCAGGAGTAGCCACTGCGGGATGGCGCCAAAGTACTGTGGCCAGAAATCGATACGCTGTTCATCGGGGACCGGCGTGACGGTCAGTTCAAATTCGGGTTGGTTAGCGGGTGCGAGGTCGTGCTGCGTc                                            <  1:986012/151‑1 (MQ=255)
                                                                                                      catccaGGCGAAGATATGCGGCTCCAGGAGTAGCCACTGCGGGATGGCGCCAAAGTACTGTGGCCAGAAATCGATACGCTGTTCATCGGGGACCGGCGTGACGGTCAGTTCAAATTCGGGTTGGTTAGCGGGTGCGAGGTCGTGCTGCGTc                                            <  1:520849/151‑1 (MQ=255)
                                                                                                      catccaGGCGAAGATATGCGGCTCCAGGAGTAGCCACTGCGGGATGGCGCCAAAGTACTGTGGCCAGAAATCGATACGCTGTTCATCGGGGACCGGCGTGACGGTCAGTTCAAATTCGGGTTGGTTAGCGGGTGCGAGGTCGTGCTGCGTc                                            >  1:1062042/1‑151 (MQ=255)
                                                                                                       atccaGGCGAAGATATGCGGCTCCAGGAGTAGCCACTGCGGGATGGCGCCAAAGTACTGTGGCCAGAAATCGATACGCTGTTCAGCGGGGACCGGCGTTACGGTCAGTTCAAATTCGGGTTGGTTAGCGGGTGCGAGGGCGTGCTGCGTCt                                           >  2:671114/1‑151 (MQ=255)
                                                                                                        tccaGGCGAAGATATGCGGCTCCAGGAGTAGCCACTGCGGGATGGCGCCAAAGTACTGTGGCCAGAAATCGATACGCTGTTCATCGGGGACCGGCGTGACGGTCAGTTCAAATTCGGGTTGGTTAGCGGGTGCGAGGTCGTGCTGCGTCtg                                          >  1:146207/1‑151 (MQ=255)
                                                                                                             gCGAAGATATGCGGCTCCAGGAGTAGCCACTGCGGGATGGCGCCAAAGTACTGTGGCCAGAAATCGATACGCTGTTCATCGGGGACCGGCGTGACGGTCAGTTCAAATTCGGGTTGGTTAGCGGGTGCGAGGTCGTGCTGCGTCTGtgttg                                     <  2:174475/151‑1 (MQ=255)
                                                                                                               gAAGATATGCGGCTCCAGGAGTAGCCACTGCGGGATGGCGCCAAAGTACTGTGGCCAGAAATCGATACGCTGTTCATCGGGGACCGGCGTGACGGTCAGTTCAAATTCGGGTTGGTTAGCGGGTGCGAGGTCGTGCTGCGTCTGTGTTGTc                                   <  1:355502/151‑1 (MQ=255)
                                                                                                                   atatGCGGCTCCAGGAGTAGCCACTGCGGGATTGCGCCAAAGTACTGTGGCCAGAAATCGATACGCTGTTCATCGGGGACCGGCGTGACGGTCAGTTCAAATTCGGGTTGGTTAGCGGGTGCGAGGTCGTGCTGCGTCTGTGTTGTCATgg                               <  1:1189826/151‑1 (MQ=255)
                                                                                                                     atGCGGCTCCAGGAGTAGCCACTGCGGGATGGCGCCAAAGTACTGTGGCCAGAAATCGATACGCTGTTCATCGGGGACCGGCGTGACGGTCAGTTCAAATTCGGGTTGGTTAGCGGGTGCGAGGTCGTGCTGCGTCTGTGTTGTCATGGGt                             <  2:200711/151‑1 (MQ=255)
                                                                                                                     atGCGGCTCCAGGAGTAGCCACTGCGGGATGGCGCCAAAGTACTGTGGCCAGAAATCGATACGCTGTTCATCGGGGACCGGCGTGACGGTCAGTTCAAATTCGGGTTGGTTAGCGGGTGCGAGGTCGTGCTGCGTCTGTGTTGTCATGGGt                             >  1:573732/1‑151 (MQ=255)
                                                                                                                      tGCGGCTCCAGGAGTAGCCACTGCGGGATGGCGCCAAAGTACTGTGGCCAGAAATCGATACGCTGTTCATCGGGGACCGGCGTGACGGTCAGTTCAAATTCGGGTTGGTTAGCGGGTGCGAGGTCGTGCTGCGTCTGTGTTGTCATGGGt                             >  1:1137171/1‑150 (MQ=255)
                                                                                                                       gCGGCTCCAGGAGTAGCCACTGCGGGATGGCGCCAAAGTACTGTGGCCAGAAATCGATACGCTGTTCATCGGGGACCGGCGTGACGGTCAGTTCAAATTCGGGTTGGTTAGCGGGTGCGAGGTCGTGCTGCGTCTGTGTTGTCATGGGTa                            <  1:1250568/150‑1 (MQ=255)
                                                                                                                             ccAGGAGTAGCCACTGCGGGATGGCGCCAAAGTACTGTGGCCAGAAATCGATACGCTGTTCATCGGGGACCGGCGTGACGGTCAGTTCAAATTCGGGTTGGTTAGCGGGTGCGAGGTCGTGCTGCGTCTGTGTTGTCATGGGTATGtct                       <  2:1455485/149‑1 (MQ=255)
                                                                                                                             ccAGGAGTAGCCACTGCGGGATGGCGCCAAAGTACTGTGGCCAGAAATCGATACGCTGTTCATCGGGGACCGGCGTGACGGTCAGTTCAAATTCGGGTTGGTTAGCGGGTGCGAGGTCGTGCTGCGTCTGTGTTGTCATGGGTATGTCTcc                     >  1:424366/1‑151 (MQ=255)
                                                                                                                             ccAGGAGTAGCCACTGCGGGATGGCGCCAAAGTACTGTGGCCAGAAATCGATACGCTGTTCATCGGGGACCGGCGTGACGGTCAGTTCAAATTCGGGTTGGTTAGCGGGTGCGAGGTCGTGCTGCGTCTGTGTTGTCATGGGTATGTCTcc                     <  1:1134918/151‑1 (MQ=255)
                                                                                                                              cAGGAGTAGCCACTGCGGGATGGCGCCAAAGTACTGTGGCCAGAAATCGATACGCTGTTCATCGGGGACCGGCGTGACGGTCAGTTCAAATTCGGGTTGGTTAGCGGGTGCGAGGTCGTGCTGCGTCTGTGTTGTCATGGGTATGTCTcc                     >  2:890138/1‑150 (MQ=255)
                                                                                                                                   gTAGCCACTGCGGGATGGCGCCAAAGTACTGTGGCCAGAAATCGATACGCTGTTCATCGGGGACCGGCGTGACGGTCAGTTCAAATTCGGGTTGGTTAGCGGGTGCGAGGTCGTGCTGCGTCTGTGTTGTCATGGGTATGTCTCCGTCAAt               <  2:1370379/151‑1 (MQ=255)
                                                                                                                                    tAGCCACTGCGGGATGGCGCCAAAGTACTGTGGCCAGAAATCGATACGCTGTTCATCGGGGACCGGCGTGACGGTCAGTTCAAATTCGGGTTGGTTAGCGGGTGCGAGGTCGTGCTGCGTCTGTGTTGTCATGGGTATGTCTCCGTCAAt               >  1:377740/1‑150 (MQ=255)
                                                                                                                                     aGCCACTGCGGGATGGCGCCAAAGTACTGTGGCCAGAAATCGATACGCTGTTCATCGGGGACCGGCGTGACGGTCAGTTCAAATTCGGGTTGGTTAGCGGGTGCGAGGTCGTGCTGCGTCTGTGTTGTCATGGGTATGTCTCCGTCAATa              >  1:318090/1‑150 (MQ=255)
                                                                                                                                       ccACTGCGGGATGGCGCCAAAGTACTGTGGCCAGAAATCGATACGCTGTTCATCGGGGACCGGCGTGACGGTCAGTTCAAATTCGGGTTGGTTAGCGGGTGCGAGGTCGTGCTGCGTCTGTGTTGTCATGGGTATGTCTCCGTCAATaaa            >  1:1330151/1‑150 (MQ=255)
                                                                                                                                         aCTGCGGGATGGCGCCAAAGTACTGTGGCCAGAAATCGATACGCTGTTCATCGGGGACCGGCGTGACGGTCAGTTCAAATTCGGGTTGGTTAGCGGGTGCGAGGTCGTGCTGCGTCTGTGTTGTCATGGGTATGTCTCCGTCAATAAAAAc         >  2:1343247/1‑151 (MQ=255)
                                                                                                                                           tGCGGGATGGCGCCAAAGTACTGTGGCCAGAAATCGATACGCTGTTCATCGGGGACCGGCGTGACGGTCAGTTCAAATTCGGGTTGGTTAGCGGGTGCGAGGTCGTGCTGCGTCTGTGTTGTCATGGGTATGTCTCCGTCAATaaaaa          <  1:1397541/148‑1 (MQ=255)
                                                                                                                                            gCGGGATGGCGCCAAAGTACTGTGGCCAGAAATCGATACGCTGTTCATCGGGGACCGGCGTGACGGTCAGTTCAAATTCGGGTTGGTTAGCGGGTGCGAGGTCGTGCTGCGTCTGTGTTGTCATGGGTATGTCTCCGTCAATAAAAACGcc      <  1:497804/151‑1 (MQ=255)
                                                                                                                                              gggATGGCGCCAAAGTACTGTGGCCAGAAATCGATACGCTGTTCATCGGGGACCGGCGTGACGGTCAGTTCAAATTCGGGTTGGTTAGCTGGTGCGAGGTCGTGCTGCGTCTGTGTTGTCATGGGTATGTCTCCGTCAATAAAAACGCCAg    <  2:918717/151‑1 (MQ=255)
                                                                                                                                                gATGGCGCCAAAGTACTGTGGCCAGAAATCGATACGCTGTTCATCGGGGACCGGCGTGACGGTCAGTTCAAATTCGGGTTGGTTAGCGGGTGCGAGGTCGTGCTGCGTCTGTGTTGTCATGGGTATGTCTCCGTCAATAAAAACGCCAgcg  <  2:857285/151‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
CCATGTCTCATCGTTGTCAGGCTCGGGGGACATAAATGCGCCGCCATTGCTGAGGGTGTAGAACGACCAGATACCACCGCTGTAGCCCTCACAGAAGCGGTCCATCCAGGCGAAGATATGCGGCTCCAGGAGTAGCCACTGCGGGATAGCGCCAAAGTACTGTGGCCAGAAATCGATACGCTGTTCATCGGGGACCGGCGTGACGGTCAGTTCAAATTCGGGTTGGTTAGCGGGTGCGAGGTCGTGCTGCGTCTGTGTTGTCATGGGTATGTCTCCGTCAATAAAAACGCCAGCG  >  NC_000913/264943‑265237

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