Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 286,206 G→C 5.9% G271A (GGC→GCC)  yagG → CP4‑6 prophage; putative sugar transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913286,2060GC 5.9% 232.9 / 5.3 68G271A (GGC→GCC) yagGCP4‑6 prophage; putative sugar transporter
Reads supporting (aligned to +/- strand):  ref base G (28/36);  new base C (2/2);  total (30/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.14e-01

TGCTGCTGGGCAACAGCCAGTGGCGCATCATGTGCGCGTTCAAGATGATGGCGACCTGCTCCAACGTGGTGCGCGGCGGGGCGACGCTCTACTTCGTGAAATACGTGATGGATCACCCGGAGTTGGCGACCCAGTTTTTACTTTACGGCAGCCTCGCCACCATGTTCGGCTCGCTTTGCTCCTCACGCCTGCTGGGCCGCTTCGACCGCGTCACCGCCTTCAAGTGGATCATCGTCGCCTACTCGCTGATCAGCCTGCTGATTTTCGTCACCCCGGCGGAGCACATCGCGC  >  NC_000913/286059‑286349
                                                                                                                                                   |                                                                                                                                               
tgctgctgGGCTACAGCCAGTGGCGCATCATGTGCGCGTTCAAGATGATGGCGCCCTGCTCCAACGGGGTGCGCGGCGGGGCGACGCTCTACTTCGGGAAATACGTGATGGATCACCCGGAGTTGGCGACCCAGTTTTTACTTTACGGc                                                                                                                                                <  4:228497/149‑1 (MQ=255)
   tgctgGGCAACAGCCAGTGGCGCATCATGTGCGCGTTCAAGATGATGGCGACCTGCTCCAACGTGGTGCGCGGCGGGGCGACGCTCTACTTCGTGAAATACGTGATGGATCACCCGGAGTTGGCGACCCAGTTTTTACTTTACGGCAGc                                                                                                                                             >  3:289634/1‑149 (MQ=255)
    gctgGGCAACAGACAGTGGCGAATCATGTGCGCGTTCAAGATGATGGAGACCTGCTACAACGTGGTGCTCGGCGTGGCGACGATCTAATTCGTGAAATACGTGATGGATCACCCGGAGTTGGCAAACCAGTGTTTACTTTACGGCAgac                                                                                                                                            >  8:49831/1‑147 (MQ=255)
        ggCAACAGCCAGTGGCGCATCATGTGCGCGTTCAAGATGATGGCGACCTGCTCCAACGTGGTGCGCGGCGGGGCGACGCTCTACTTCGTGAAATACGTGATGGATCACCCGGAGTTGGCGACCCAGTTTTTACTTTACGGCAGCCTc                                                                                                                                          <  6:45316/147‑1 (MQ=255)
             cAGCCAGTGGCGCATCATGTGCGCGTTCAAGATGATGGCGACCTGCTCCAACGTGGTGCGCGGCGGGGCGACGCTCTACTTCGTGAAATACGTGATGGATCACCCGGAGTTGGCGACCCAGTTTTTACTTTACGGCAGCCTCGccacca                                                                                                                                   >  6:314928/1‑149 (MQ=255)
                ccAGTGGCGCATCATGTGCGCGTTCAAGATGATGGCGACCTGCTCCAACGTGGTGCGCGGCGGGGCGACGCTCTACTTCGTGAAATACGTGATGGATCACCCGGAGTTGGCGACCCAGTTTTTACTTTACGGCAGCCTCGCCACCATGt                                                                                                                                <  7:359404/149‑1 (MQ=255)
                ccAGTGGCGCATCATGTGCGCGTTCAAGATGATGGCGACCTGCTCCAACGTGGTGCGCGGCGGGGCGACGCTCTACTTCGTGAAATACGTGATGGATCACCCGGAGTTGGCGACCCAGTTTTTACTTTACGGCAGCCTCGCCACCATGt                                                                                                                                <  5:315090/149‑1 (MQ=255)
                ccAGTGGCGCATCATGTGCGCGTTCAAGATGATGGCGACCTGCTCCAACGTGGTGCGCGGCGGGGCGACGCTCTACTTCGTGAAATACGTGATGGATCACCCGGAGTTGGCGACCCAGTTTTTACTTTACGGCAGCCTCGCCACCATGt                                                                                                                                >  3:310205/1‑149 (MQ=255)
                  aGTGGCGCATCATGTGCGCGTTCAAGATGATGGCGACCTGCTCCAACGTGGTGCGCGGCGGGGCGACGCTCTACTTCGTGAAATACGTGATGGATCACCCGGAGTTGGCGACCCAGTTTTTACTTTACGGCAGCCTCGCCACCATGTTc                                                                                                                              >  8:16155/1‑149 (MQ=255)
                    tGGCGCATCATGTGCGCGTTCAAGATGATGGCGACCTGCTCCAACGTGGTGCTCGGCGGGGCGACGCTCTACTTCGTGAAATACGTGATGGATCACCCGGAGTTGGCGACCCAGTTTTTACTTTACGGCAGCCTCGCCACCATGTTCgg                                                                                                                            <  1:9256/149‑1 (MQ=255)
                    tGGCGCATCATGTGCGCGTTCAAGATGATGGCGACCTGCTCCAACGTGGTGCGCGGCGGGGCGACGCTCTACTTCGTGAAATACGTGATGGATCACCCGGAGTTGGCGACCCAGTTTTTACTTTACgg                                                                                                                                                 >  8:97672/1‑128 (MQ=255)
                    tGGCGCATCATGTGCGCGTTCAAGATGATGGCGACCTGCTCCAACGTGGTGCGCGGCGGGGCGACGCTCTACTTCGTGAAATACGTGATGGATCACCCGGAGTTGGCGACCCAGTTTTTACTTTACgg                                                                                                                                                 <  7:97672/128‑1 (MQ=255)
                     ggCGCATCATGTGCGCGTTCAAGATGATGGCGACCTGCTCCAACGTGGTGCGCGGCGGGGCGACGCTCTACTTCGTGAAATACGTGATGGATCACCCGGAGTTGGCGACCCAGTTTTTACTTTACGGCAGCCTCGCCACCATGTTCGGc                                                                                                                           <  2:82734/149‑1 (MQ=255)
                          atcatGTGCGCGTTCAAGATGATGGCGACCTGCTCCAACGTGGTGCGCGGCGGGGCGACGCTCTACTTCGTGAAATACGTGATGGATCACCCGGAGTTGGCGACCCAGTTTTTACTTTACGGCAGCCTCGCCACCATGTTCGGCTCg                                                                                                                        <  1:50300/147‑1 (MQ=255)
                          atcatGTGCGCGTTCAAGATGATGGCGACCTGCTCCAACGTGGTGCGCGGCGGGGCGACGCTCTACTTCGTGAAATACGTGATGGATCACCCGGAGTTGGCGACCCAGTTTTTACTTTACGCCAGc                                                                                                                                             <  1:145423/126‑1 (MQ=255)
                          atcatGTGCGCGTTCAAGATGATGGCGACCTGCTCCAACGTGGTGCGCGGCGGGGCGACGCTCTACTTCGTGAAATACGTGATGGATCACCCGGAGTTGGCGACCCAGTTTTTACTTTACGCCAGc                                                                                                                                             >  2:140331/1‑126 (MQ=255)
                          atcatGTGCGCGTTCAAGATGATGGCGACCTGCTCCAACGTGGTGCGCGGCGGGGCGACGCTCTACTTCGTGAAATACGTGATGGATCACCCGGAGTTGGCGACCCAGTTTTTACTTTACGCCAGc                                                                                                                                             <  1:140331/126‑1 (MQ=255)
                          atcatGTGCGCGTTCAAGATGATGGCGACCTGCTCCAACGTGGTGCGCGGCGGGGCGACGCTCTACTTCGTGAAATACGTGATGGATCACCCGGAGTTGGCGACCCAGTTTTTACTTTACGCCAGc                                                                                                                                             >  2:145423/1‑126 (MQ=255)
                          atcatGTGCGCGTTCAAGATGATGGCGACCTGCTCCAACGTGGAGCGCGGCGGGGCGACGCTCTACTTCGTGAAATACGTGATGGATCACCCGGAGTTGGCGACCCAGTTTTTACTTTACGGCAGCCTCGCCACCATGTTCGGCTCg                                                                                                                        >  2:50300/1‑147 (MQ=255)
                             atGTGCGCGTTCAAGATGATGGCGACCTGCTCCAACGTGGTGCGCGGCGGGGCGACGCTCTACTTCGTGAAATACGTGATGGATCACCCGGAGTTGGCGACCCAGTTTTTACTTTACGGCAGCCTCGCCACCATGTTCGGCTCGCTTTg                                                                                                                   <  6:63688/149‑1 (MQ=255)
                                tgCGCGTTCAAGATGATGGCGACCTGCTCCAACGTGGTGCGCGGCGGGGCGACGCTCTACTTCGTGAAATACGTGATGGATCACCCGGAGTTGGCGACCCAGTTTTTACTTTACGGCAGCCTCGCCACCATGTTCGGCTCGCTTTGctc                                                                                                                <  4:266030/149‑1 (MQ=255)
                                tgCGCGTTCAAGATGATGGCGACCTGCTCCAACGTGGTGCGCGGCGGGGCGACGCTCTACTTCGTGAAATACGTGATGGATCACCCGGAGTTGGCGACCCAGTTTTTACTTTACGGCAGCCTCGCCACCATGTTCGGCTCGCTTTGct                                                                                                                 <  3:177175/148‑1 (MQ=255)
                                      ttCAAGATGATGGCGACCTGCTCCAACGTGGTGCGCGGCGGGGCGACGCTCTACTTCGTGAAATACGTGATGGATCACCCGGAGTTGGCGACCCAGTTTTTACTTTACGGCAGCCTCGCCACCATGTTCGGCTCGCTTTGCTCCTCACg                                                                                                          >  2:181838/1‑149 (MQ=255)
                                        cAAGATGATGGCGACCTGCTCCAACGTGGTGCGCGGCGGGGCGACGCTCTACTTCGTGAAATACGTGATGGATCACCCGGAGTTGGCGACCCAGTTTTTACTTTACGGCAGCCTCGCCACCATGTTCGGCTCGCTTTGCTCCTCACGc                                                                                                         <  7:16155/148‑1 (MQ=255)
                                              gatgGCGACCTGCTCCAACGTGGTGCGCGGCGGGGCGACGCTCTACTTCGTGAAATACGTGATGGATCACCCGGAGTTGGCGACCCAGTTTTTACTTTACGGCAGCCTCGCCACCATGTTCGGCTCGCTTTGCTCCTCACGCCTGCTgg                                                                                                  >  8:240216/1‑149 (MQ=255)
                                               atgGCGACCTGCTCCAACGTGGTGCGCGGCGGGGCGACGCTCTACTTCGTGAAATACGTGATGGATCACCCGGAGTTGGCGACCCAGTTTTTACTTTACGGCAGCCTCGCCACCATGTTCGGCTCGCTTTGCTCCTCACGCCTGCTgg                                                                                                  <  1:181838/148‑1 (MQ=255)
                                                tgGCGACCTGCTCCAACGTGGTGCGCGGCGGGGCGACGCTCTACTTCGTGAAATACGTGATGGATCACCCGGAGTTGGCGACCCAGTTTTTACTTTACGGCAGCCTCGCCACCATGTTCGGCTCGCTTTGCTCCTCACGCCTGCTGGGc                                                                                                <  3:344497/149‑1 (MQ=255)
                                                   cGACCTGCTCCAACGTGGTGCGCGGCGGGGCGACGCTCTACTTCGTGAAATACGTGATGGATCACCCGGAGTTGGCGACCCAGTTTTTACTTTACGGCAGCCTCGCCACCATGTTCGGCTCGCTTTGCTCCTCACGCCTGCTGGGCCGc                                                                                             >  2:214716/1‑149 (MQ=255)
                                                      ccTGCTCCAACGTGGTGCGCGGCGGGGCGACGCTCTACTTCGTGAAATACGTGATGGATCACCCGGAGTTGGCGACCCAGTTTTTACTTTACGGCAGCCTCGCCACCATGTTCGGCTCGCTTTGCTCCTCACGCCTGCTGGGCCGCTTc                                                                                          <  5:314928/149‑1 (MQ=255)
                                                        tGCTCCAACGTGGTGCGCGGCGGGGCGACGCTCTACTTCGTGAAATACGTGATGGATCACCCGGAGTTGGCGACCCAGTTTTTACTTTACGGCAGCCTCGCCACCATGTTCGGCTCGCTTTGCTCCTCACGCCTGCTGGGCCGCTTCGa                                                                                        >  6:337140/1‑149 (MQ=255)
                                                             cAACGTGGTGCGCGTCGGGGCGACGCTCGACTTCGTGAAATACGTGATGGATCACCCGGAGTTGGCGACCCAGTTTTTACTTTACGGCAGCCTCGCCACCATGTTCGGCTCGCTTTGCTCCTCACGCCTGCTGGACCGCTTCGACcgc                                                                                    <  4:28918/148‑1 (MQ=255)
                                                             cAACGTGGTGCGCGGCGGGGCGACGCTCTACTTCGTGAAATACGTGATGGATCCCCCGGAGTTGGCGCCCCAGTTTTTACTTTACGGCAGCCTCGCCACCATGTTCGGCTCGCTTTGCTCCTCACGCCTGCTGGGCCGCTTCGACcgcg                                                                                   >  7:176898/1‑149 (MQ=255)
                                                             cAACGTGGTGCGCGGCGGGGCGACGCTCTACTTCGTAAAATACGTGATGGATCACCCGGAGTTGGCGACCCAGTTTTTACTTTACGGCAGCCTCGCCACCATGTTCGGCTCGCTTTGCTCCTCACGCCTGCTGGGCCGCTTCGACcgcg                                                                                   >  7:207272/1‑149 (MQ=255)
                                                               aCGTGGTGCGCGGCGGGGCGACGCTCTACTTCGTGAAATACGTGATGGATCACCCGGAGTTGGCGACCCAGTTTTTACTTTACGGCAGCCTCGCCACCATGTTCGGCTCGCTTTGCTCCTCACGCCTGCTGGGCCGCTTCGACCGCGTc                                                                                 <  3:70591/149‑1 (MQ=255)
                                                                       cgcgGCGGGGCGACGCTCTACTTCGTGAAATACGTGATGGATCACCCGGAGTTGGCGACCCAGTTTTTACTTTACGGCAGCCTCGCCACCATGTTCGGCTCGCTTTGCTCCTCACGCCTGCTGGGCCGCTTCGACCGCGTCACCGCCtt                                                                         <  6:324542/149‑1 (MQ=255)
                                                                         cggcggGGCGACGCTCTACTTCGTGAAATACGTGATGGATCACCCGGAGTTGGCGACCCAGTTTTTACTTTACGGCAGCCTCGCCACCATGTTCGGCTCGCTTTGCTCCTCACGCCTGCTGGGCCGCTTCGACCGCGTCACCGCCTTc                                                                        >  8:275377/1‑148 (MQ=255)
                                                                               ggCGACGCTCTACTTCGTGAAATACGTGATGGATCACCCGGAGTTGGCGACCCAGGTTTTACTTTACGGCAGCCTCGCCACCATGTTCGGCTCGCTTTGCTCCTCACGCCTGCTGGGCCGCTTCGACCGCGTCACCGCCTTCAAGTGGa                                                                 <  8:33671/149‑1 (MQ=255)
                                                                                 cGACGCTCTACTTCGTGAAATACGTGATGGATCACCCGGAGTTGGCGACCCAGTTTTTACTTTACGGCAGCCTCGCCACCATGTTCGGCTCGCTTTGCTCCTCACGCCTGCTGGGCCGCTTCGACCGCGTCACCGCCTTCAAGTGGatc                                                               <  3:290189/149‑1 (MQ=255)
                                                                                 cGACGCTCTACTTCGTGAAATACGTGATGGATCACCCGGAGTTGGCGACCCAGTTTTTACTTTACGGCAGCCTCGCCACCATGTTCGGCTCGCTTTGCTCCTCACGCCTGCTGGGCCGCTTCGACCGCGTCACCGCCTTCAAGTGGatc                                                               <  2:30218/149‑1 (MQ=255)
                                                                                  gACGCTCTACTTCGTGAAATACGTGATGGATCACCCGGAGTTGGCGACCCAGTTTTTACTTTACGGCAGCCTCGCCACCATGTTCGGCTCGCTTTGCTCCTCACGCCTGCTGGGCCGCTTCGACCGCGTCACCGCCTTCAAGTGGatca                                                              <  3:135551/149‑1 (MQ=255)
                                                                                   aCGCTCTACTTCGTGAAATACGTGATGGATCACCCGGAGTTGGCGACCCAGTTTTTACTTTACGGCAGCCTCGCCACCATGTTCGGCTCGCTTTGCTCCTCACGCCTGCTGGGCCGCTTCGACCGCGTCACCGCCTTCAAGTGGatcat                                                             <  5:113607/149‑1 (MQ=255)
                                                                                     gCTCTACTTCGTGAAATACGTGATGGATCACCCGGAGTTGGCGACCCAGTTTTTACTTTACGGCAGCCTCGCCACCATGTTCGGCTCGCTTTGCTCCTCACGCCTGCTGGGCCGCTTCGACCGCGTCACCGCCTTCAAGTGGATCAtcg                                                           >  1:69283/1‑149 (MQ=255)
                                                                                           cTTCGTGAAATACGTGATGGATCACCCGGAGTTGGCGACCCAGTTTTTACTTTACGGCAGCCTCGCCACCATGTTCGGCTCGCTTTGCTCCTCTCGCCTGCTGGGCCGCTTCGACCGCGTCACCGCCTTCAAGTGGATCATCGTCGcc                                                      <  4:310205/148‑1 (MQ=255)
                                                                                           cTTCGTGAAATACGTGATGGATCACCCGGAGTTGGCGACCCAGTTTTTACTTTACGGCAGCCTCGCCACCATGTTCGGCTCGCTTTGCTCCTCACGCCTGCTGGGCCGCTTCGACCGCGTCACCGCCTTCAAGTGGATCATCGTCGc                                                       <  4:289634/147‑1 (MQ=255)
                                                                                           cTTCGTGAAATACGTGATGGATCACCCGGAGTTGGCGACCCAGTTTTTACTTTACGGCAGCCTCGCCACCATGTTCGGCTCGCTTTGCTCCTCACGCCTGCTGGGCCGCTTCGACCGCGTCACCGCCTTCAAGTGGATCATCGTCGc                                                       <  2:69283/147‑1 (MQ=255)
                                                                                            ttCGTGAAATACGTGATGGATCACCCGGAGTTGGCGACCCAGTTTTTACTTTACGGCAGCCTCGCCACCATGTTCGGCTCGCTTTGCTCCTCACGCCTGCTGGGCCGCTTCGACCGCGTCACCGCCTTCAAGTGGatca                                                              >  1:69483/1‑139 (MQ=255)
                                                                                            ttCGTGAAATACGTGATGGATCACCCGGAGTTGGCGACCCAGTTTTTACTTTACGGCAGCCTCGCCACCATGTTCGGCTCGCTTTGCTCCTCACGCCTGCTGGGCCGCTTCGACCGCGTCACCGCCTTCAAGTGGatca                                                              <  2:69483/139‑1 (MQ=255)
                                                                                            ttCGTGAAATACGTGATGGATCACCCGGAGTTGGCGACCCAGTTTTTACTTTACGGCAGCCTCGCCACCATGTTCGGCTCGCTTTGCTCCTCACGCCTGCTGGGCCGCTTCGACCGCGTCACCGCCTTCAAGTGGATCATCGTCGCCTa                                                    >  7:373424/1‑149 (MQ=255)
                                                                                             tCGTGAAATACGTGATGGATCACCCGGAGTTGGCGACCCAGTTTTTACTTTACGGCAGCCTCGCCACCATGTTCGGCTCGCTTTGCTCCTCACGCCTGCTGGGCCGCTTCGACCGCGTCAc                                                                               <  3:87291/121‑1 (MQ=255)
                                                                                             tCGTGAAATACGTGATGGATCACCCGGAGTTGGCGACCCAGTTTTTACTTTACGGCAGCCTCGCCACCATGTTCGGCTCGCTTTGCTCCTCACGCCTGCTGGGCCGCTTCGACCGCGTCAc                                                                               >  4:87291/1‑121 (MQ=255)
                                                                                                tGAAATACGTGATGGATCACCCGGAGTTGGCGACCCAGTTTTTACTTTACGGCAGCCTCGCCACCATGTTCGGCTCGCTTTGCTCCTCACGCCTGCTGGGCCGCTTCGACCGCGTCACCGCCTTCAAGTGGATCATCGTCGCCTACTCg                                                <  1:214716/149‑1 (MQ=255)
                                                                                                 gAAATACGTGATGGATCACCCGGAGTTGGCGACCCAGTTTTTACTTTACGGCAGCCTCGCCACCATGTTCGGCTCGCTTTGCTCCTCACGCCTGCTGGGCCGCTTCGACCGCGTCACCGCCTTCAAGTGGATCATCGTCGCCTACTCGc                                               <  3:128605/149‑1 (MQ=255)
                                                                                                       cGTGATGGATCACCCGGAGTTGGCGACCCAGTTTTTACTTTACGGCAGCCTCGCCACCATGTTCGGCTCGCTTTGCTCCTCACGCCTGCTGGGCCGCTTCGACCGCGTCACCGCCTTCAAGTGGATCATCGTCGCCTACTCGCTGATc                                          <  5:337140/148‑1 (MQ=255)
                                                                                                        gTGATGGATCCCCCGGAGTTGGCGCCCCAGGTTTTACTTGACGGCCGCCTCGCCACCATGTTCGGCTCGCTTGTCTCCACAAGCCGGCTGGGCCGCTGCGACCGCGGCCCCGCCGTCAAGTGGATCATCGTCGCCTCCGCGCTGCTCAg                                        >  5:53485/1‑149 (MQ=255)
                                                                                                         tGATGGATCACCCGGAGTTGGCGACCCAGTTTTTACTTTACGGCAGCCTCGCCACCATGTTCGGCTCGCTTTGCTCCTCACGCCTGCTGGGCCGCTTCGACCGCGTCACCGCCTTCAAGAGGATCATCGTCGCCTACTCGCTGATCAGc                                       >  4:75940/1‑149 (MQ=255)
                                                                                                          gATGGATCACCCGGAGTTGGCGACCCAGTTTTTACTTTACGGCAGCCTCGCCACCATGTTCGGCTCGCTTTGCTCCTCACGCCTGCTGGGCCGCTTCGACCGCGTCACCGCCTTCAAGAGGATCATCGTCGCCTACTCGCTGATCAGcc                                      >  5:169924/1‑149 (MQ=255)
                                                                                                             ggATCACCCGGAGTTGGCGACCCAGTTTTTACTTTACGGCAGCCTCGCCACCATGTTCGGCTCGCTTTGCTCCTCACGCCTGCTGGGCCGCTTCGACCGCGTCACCGCCTTCAAGTGGATCATCGTCGCCTACTCGCTGATCAGCctgc                                   >  8:24397/1‑149 (MQ=255)
                                                                                                              gATCACCCGGAGTTGGCGACCCAGTTTTTACTTTACGGCAGCCTCGCCACCATGTTCGGCTCGCTTTGCTCCTCACGCCTGCTGGGCCGCTTCGACCGCGTCACCGCCTTCAAGTGGATCATCGTCGCCTACTCGCTGATCAGCctgc                                   <  3:60785/148‑1 (MQ=255)
                                                                                                                tCACCCGGAGTTGGCGACCCAGTTTTTACTTTACGGCAGCCTCGCCACCATGTTCGGCTCGCTTTGCTCCTCACGCCTGCTGGGCCGCTTCGACCGCGTCACCGCCTTCAAGTGGATCATCGTCGCCTACTCGCTGATCAGCCTGCTGa                                >  4:88381/1‑149 (MQ=255)
                                                                                                                     cGGAGTTGGCGACCCAGTTTTTACTTTACGGCAGCCTCGCCACCATGTTCGGCTCGCTTTGCTCCTCACGCCTGCTGGGCCGCTTCGACCGCGTCACCGCCTTCAAGTGGATCATCGTCGCCTACTCGCTGATCAGCCTGCTGATTTTc                           <  8:207272/149‑1 (MQ=255)
                                                                                                                                aCCCAGTTTTTACTTTACGGCAGCCTCGCCACCATGTTCGGCTCGCTTTGCTCCTCACGCCTGCTGGGCCGCTTCGACCGCGTCACCGCCTTCAAGTGGATCATCGTCGCCTACTCGCTGATCAGCCTGCTGATTTTCGTCACCCcggc                >  7:141507/1‑149 (MQ=255)
                                                                                                                                aCCCAGTTTTTACTTTACGGCAGCCTCGCCACCATGTTCGGCTCGCTTTGCTCCTCACGCCTGCTGGGCCGCTTCGACCGCGTCACCGCCTTCAAGTGGATCATCGTCGCCTACTCGCTGATCAGCCTGCTGATTTTCGTCACCCcggc                >  8:223974/1‑149 (MQ=255)
                                                                                                                                  ccAGTTTTTACTTTACGGCAGCCTCGCCACCATGTTCGGCTCGCTTTGCTCCTCACGCCTGCTGGGCCGCTTCGACCGCGTCACCGCCTTCAAGTGGATCATCGTCGCCTACTCGCTGATCAGCCTGCTGATTTTCGTCACCCcggcgg              <  7:275377/149‑1 (MQ=255)
                                                                                                                                      tttttACTTTACGGCAGCCTCGCCACCATGTTCGGCTCGCTTTGCTCCTCACGCCTGCTGGGCCGCTTCGACCGCGTCACCGCCTTCAAGTGGATCATCGTCGCCTACTCGCTGATCAGcc                                      >  5:5924/1‑121 (MQ=255)
                                                                                                                                      tttttACTTTACGGCAGCCTCGCCACCATGTTCGGCTCGCTTTGCTCCTCACGCCTGCTGGGCCGCTTCGACCGCGTCACCGCCTTCAAGTGGATCATCGTCGCCTACTCGCTGATCAGcc                                      <  6:5924/121‑1 (MQ=255)
                                                                                                                                           actttacGGCAGCCTCGCCACCATGTTCGGCTCGCTTTGCTCCTCACGCCTGCTGGGCCGCTTCGACCGCGTCACCGCCTTCAAGTGGATCATCGTCGCCTACTCGCTGATCAGCCTGCTGATTTTCGTCACCCCGGCGGAGCACATc      <  8:176898/148‑1 (MQ=255)
                                                                                                                                            ctttacGGCAGCCTCGCCACCATGTTCGGCTCGCTTTGCTCCTCACGCCTGCTGGGCCGCTTCGACCGCGTCACCGCCTTCAAGTGGATCATCGTCGCCTACTCGCTGATCAGCCTGCTGATTTTCGTCACCCCGGCGGAGCACATcgc    >  5:282401/1‑149 (MQ=255)
                                                                                                                                             tttacGGCAGCCTCGCCACCATGTTCGGCTCGCTTTGCTCCTCACGCCTGCTGGGCCGCTTCGACCGCGTCACCGCCTTCAAGTGGATCATCGTCGCCTACTCGCTGATCAGGCTGCTGATTTTCGTCACCCCGGCGGAGCACATcgcg   >  1:308220/1‑149 (MQ=255)
                                                                                                                                             tttacGGCAGCCTCGCCACCATGTTCGGCTCGCTTTGCTCCTCACGCCTGCTGGGCCGCTTCGACCGCGTCACCGCCTTCAAGTGGATCATCGTCGCCTACTCGCTGATCAGCCTGCTGATTTTCGTCACCCCGGCGGAGCACATcgcg   >  4:249279/1‑149 (MQ=255)
                                                                                                                                               tacGGCAGCCTCGCCACCATGTTCGGCTCGCTTTGCTCCTCACGCCTGCTGGGCCGCTTCGACCGCGTCACCGCCTTCAAGTGGATCATCGTCGCCTACTCGCTGATCAGCCTGCTGATTTTCGTCACCCCGGCGGAGCACATCgcgc  <  7:24397/148‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                               
TGCTGCTGGGCAACAGCCAGTGGCGCATCATGTGCGCGTTCAAGATGATGGCGACCTGCTCCAACGTGGTGCGCGGCGGGGCGACGCTCTACTTCGTGAAATACGTGATGGATCACCCGGAGTTGGCGACCCAGTTTTTACTTTACGGCAGCCTCGCCACCATGTTCGGCTCGCTTTGCTCCTCACGCCTGCTGGGCCGCTTCGACCGCGTCACCGCCTTCAAGTGGATCATCGTCGCCTACTCGCTGATCAGCCTGCTGATTTTCGTCACCCCGGCGGAGCACATCGCGC  >  NC_000913/286059‑286349

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