Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 276,186 G→A 5.9% A288T (GCT→ACT)  mmuP → CP4‑6 prophage; putative S‑methylmethionine transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913276,1860GA 5.9% 222.7 / 3.0 68A288T (GCT→ACT) mmuPCP4‑6 prophage; putative S‑methylmethionine transporter
Reads supporting (aligned to +/- strand):  ref base G (32/32);  new base A (2/2);  total (34/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.89e-01

TACTACCATCGCGCGACTGATTATTTTCTTTATCGGCACCGTGTTTGTGCTGGCAGCGCTGATCCCGATGCAGCAGGTGGGCGTGGAGAAAAGCCCGTTTGTGCTGGTATTTGAGAAAGTAGGGATCCCGTACGCCGCTGATATTTTTAACTTCGTGATCCTGACGGCTATTCTTTCTGCAGCGAACTCCGGGTTATATGCCTCCGGGCGCATGCTGTGGTCGTTGTCGAATGAACGTACGCTACCGGCCTGTTTTGCGCGAGTAACGAAAAACGGC  >  NC_000913/276050‑276326
                                                                                                                                        |                                                                                                                                            
tactacCATCGCGCGACTGATTATTTTCTTTATCGGCACCGTGTGTGTGCTGGCAGCGCTGATCCCGATGCAGCAGGTGGGCGTGGAGAAAAGCCCGTTTGTGCTGGTATTTGAGAAAGTAGGGATCCCGTACGCCGCTGATAtttt                                                                                                                                    >  3:142751/1‑147 (MQ=255)
    acCATCGCGCGACTGATTATTTTCTTTATCGGCACCGTGTTTGTGCTGGCAGCGCTGATCCCGATGCAGCAGGTGGGCGTGGAGAAAAGCCCGTTTGTGCTGGTATTTGAGAAAGTAGGGATCCCGTACGCCGCTGATATTTTTAACt                                                                                                                               >  2:268134/1‑148 (MQ=255)
    acCATCGAGCGACTGCTTATTTTCTTTATCGGCACAGTGTTTGTGCTGGCAGCGCTGATCCCGATGCAGCAGTTGGGCGTGGAGAAAAGCGCGTTTGTGCTGGTATTTGAGAAAGTAGGGATCCCGTAAGCCGCTGATATTTTTAACt                                                                                                                               >  6:239183/1‑148 (MQ=255)
     ccATCGCGCGACTGATTATTTTCTTTATCGGCACCGTGTTTGTGCTGGCAGCGCTGATCCCGATGCAGCAGGTGGGCGTGGAGAAAAGCCCGTTTGTGCTGGTATTTGAGAAAGTAGGGATCCCGTACGCCGCTGATATTTTTAACtt                                                                                                                              >  3:109415/1‑148 (MQ=255)
           cgcgACTGATTATTTTCTTTATCGGCACCGTGTTTGTGCTGGCAGCGCTGATCCCGATGCAGCAGGTGGGCGTGGAGAAAAGCCCGTTTGTGCTGGTATTTGAGAAAGTAGGGATCCCGTACGCCGCTGATATTTTTAACTTCGTGAt                                                                                                                        <  1:268134/148‑1 (MQ=255)
             cgACTGATTTTTTTCTTTATCGGCACCGTGTTTGTGCTGGCAGCGCTGATCCCTATGCAGCAGGTGGGCGTGGAGAAAAGCCCGTTTGTGCTGGTATTTGAGAAAGTAGGGATCCCGTACGCCGCTGATATTTTTAACTTCGTGATcc                                                                                                                      <  1:296768/148‑1 (MQ=255)
             cgACTGATTATTTTCTTTATCGGCACCGTGTTTGTGCTGGCAGCGCTGATCCCGATGCAGCAGGTGGGCGTGGAGAAAAGCCCGTTTGTGCTGGTATTTGAGAAAGTAGGGATCCCGTACGCCGCTGATATTTTTAACTTCGTGATcc                                                                                                                      >  3:263303/1‑148 (MQ=255)
                     tattTTCTTTATCGGCACCGTGTTTGTGCTGGCAGCGCTGATCCCGATGCAGCAGGTGGGCGTGGAGAAAAGCCCGTTTGTGCTGGTATTTGAGAAAGTAGGGATCCCGTACGCCGCTGATATTTTTAACTTCGTGATCCTGACGGc                                                                                                               >  6:137489/1‑147 (MQ=255)
                        tttCTTTATCGGCACCGTGTTTGTGCTGGCAGCGCTGATCCCGATGCAGCAGGTGGGCGTGGAGAAAAGCCCGTTTGTGCTGGTATTTGAGAAAGTAGGGATCCCGTACGCCGCTGATATTTTTAACTTCGTGATCCTGACGGCTAtt                                                                                                           >  8:125747/1‑148 (MQ=255)
                          tCTTTATCGGCATCGTGTTTGTGCTGGCAGCGCTGATCCCGATGCAGCAGGTGGGCGTGGAGAAAAGCCCGTTTGTGCTGGTATTTGAGAAAGTAGGGATCCCGTACGCCGCTGATATTTTTAACTTCGTGATCCTg                                                                                                                    <  1:154789/137‑1 (MQ=255)
                          tCTTTATCGGCAACGTGTTTGTGCTGGCAGCGCTGATCCCGATGCAGCAGGTGGGCGTGGAGAAAAGCCCGTTTGTGCTGGTATTTGAGAAAGTAGGGATCCCGTACGCCGCTGATATTTTTAACTTCGTGATCCTg                                                                                                                    >  2:154789/1‑137 (MQ=255)
                            tttATCGGCACCGTGTTTGTGCTGGCAGCGCTGATCCCGATGCAGCAGGTGGGCGTGGAGAGAAGCCCGTTTGTGCTGGTATTTGAGAAAGTAGGGATCCCGTACGCCGCTGATATTTTTAACTTCGTGATCCTGACGGCTAt                                                                                                            >  8:85275/1‑143 (MQ=255)
                            tttATCGGCACCGTGTTTGTGCTGGCAGCGCTGATCCCGATGCAGCAGGTGGGCGTGGAGAGAAGCCCGTTTGTGCTGGTATTTGAGAAAGTAGGGATCCCGTACGCCGCTGATATTTTTAACTTCGTGATCCTGACGGCTAt                                                                                                            <  7:85275/143‑1 (MQ=255)
                            tttATCGGCACCGTGTTTGTGCTGGCAGCGCTGATCCCGATGCAGCAGGTGGGCGTGGAGAAAAGCCCGTTTGTGCTGGTATTTGAGAAAGTAGGGATCCCGTACGCCGCTGATATTTTTAACTTCGTGATCCTGACGGCTAttctt                                                                                                        >  8:985/1‑147 (MQ=255)
                              tATCGGCACCGTGTTTGTGCTGGCAGCGCTGATCCCGATGCAGCAGGTGGGCGTGGAGAAAAGCCCGTTTGTGCTGGTATTTGAGAAAGTAGGGATCCCGTACGCCGCTGATATTTTTAACTTCGTGATCCTGACGGCTAttctttc                                                                                                      <  7:985/147‑1 (MQ=255)
                                             tgtgCTGGCAGCGCTGATCCCGATGCAGCAGGTGGGCGTGGAGAAAAGCCCGTTTGTGCTGGTATTTGAGAAAGTAGGGATCCCGTACGCCGCTGATATTTTTAACTTCGTGATCCTGACGGCTATTCTTTCTGCAGCGAACTCCgg                                                                                       <  5:137489/147‑1 (MQ=255)
                                              gtgCTGGCAGCGCTGATCCCGATGCAGCAGGTGGGCGTGGAGAAAAGCCCGTTTGTGCTGGTATTTGAGAAAGTAGGGATCCCGTACGCCGCTGATATTTTTAACTTCGTGATCCTGACGGCTATTCTTTCTGCAGCGAACTCCggg                                                                                      <  8:195909/147‑1 (MQ=255)
                                                 cTGGCAGCGCTGATCCCGATGCAGCAGGTGGGCGTGGAGAAAAGCCCGTTTGTGCTGGTATTTGAGAAAGTAGGGATCCCGTACGCCGCTGATATTTTTAACTTCGTGATCCTGACGGCTATTCTTTCTGCAGCGAACTCCGGGTtat                                                                                  <  8:215367/148‑1 (MQ=255)
                                                 cTGGCAGCGCTGATCCCGATGCAGCAGGTGGGCGTGGAGAAAAGCCCGTTTGTGCTGGTATTTGAGAAAGTAGGGATCCCGTACGCCGCTGATATTTTTAACTTCGTGATCCTGACGGCTATTCTTTCTGCAGCGAACTCCGGGTtat                                                                                  <  7:202741/148‑1 (MQ=255)
                                                    gCAGCGCTGATCCCGATGCAGCAGGTGGGCGGGGAGAATAGCCCGTTTGTGCTGGTATTTGAGAAAGTAGGGATCCCGTACGCCGCTGATATTTTTAACTTCGTGATCCTGACGGCTATTCTTTCTGCAGCGAACTCCGGGTtata                                                                                 <  8:166248/146‑1 (MQ=255)
                                                        cgcTGATCCCGATGCAGCAGGTGGGCGTGGAGAAAAGCCCGTTTGTGCTGGTATTTGAGAAAGTAGGGATCCCGTACGCCGCTGATATTTTTAACTTCGTGATCCTGACGGCTATTCTTTCTGCAGCGAACTCCGGGTTATATGCCTc                                                                           >  6:251038/1‑148 (MQ=255)
                                                        cgcTGATCCCGATGCAGCAGGTGGGCGTGGAGAAAAGCCCGTTTGTGCTGGTATTTGAGAAAGTAGGGATCCCGTACGCCGCTGATATTTTTAACTTCGTGATCCTGACGGCTATTCTTTCTGCAGCGAACTCCGGGTTATATGCCTc                                                                           <  5:251038/148‑1 (MQ=255)
                                                            gATCCCGATGCAGCAGGTGGGCGTGGAGAAAAGCCCGTTTGTGCTGGTATTTGAGAAAGTAGGGATCCCGTACGCCGCTGATATTTTTAACTTCGTGATCCTGACGGCTATTCTTTCTGCAGCGAACTCCGGGTTATATGCCTCCggg                                                                       <  7:266649/148‑1 (MQ=255)
                                                               accGATGCAGCAGGTGGCCGTGGAGAAAAGCCCGTTTGAGCGTGTATTTGAGAAAGTAGGGATCCCGTAAGCCGCTGATATTTTTATCTTCGTGATCCTGACGGCTAATCTTTCTGCAGCTAACTCCGGGTTATATGCCTCCGGgcgc                                                                    <  5:239183/147‑1 (MQ=255)
                                                                    tGCAGCAGGTGGGCGTGGAGAAAAGCCCGTTTGTGCTGGTATTTGAGAAAGTAGGGATCCCGTACGCCGCTGATATTTTTAACTTCGTGATCCTGACGGCTATTCTTTCTGCAGCGAACTCCGGGTtata                                                                                 <  7:112594/130‑1 (MQ=255)
                                                                    tGCAGCAGGTGGGCGTGGAGAAAAGCCCGTTTGTGCTGGTATTTGAGAAAGTAGGGATCCCGTACGCCGCTGATATTTTTAACTTCGTGATCCTGACGGCTATTCTTTCTGCAGCGAACTCCGGGTtata                                                                                 >  8:112594/1‑130 (MQ=255)
                                                                      cagcagCTGGGCGTGGAGAAAAGCCCGTTTGTGCTGGTATTTGAGAAAGTAGGGATCCCGTACGCCGCTGATATTTTTAACTTCGTGATCCTGACGGCTATTCTTTCTGCAGCGAACTCCGGGTTATATGCCTCCGGGCGCATGCtgt                                                             <  7:125747/148‑1 (MQ=255)
                                                                        gcagGTGGGCGTGGAGAAAAGCCCGTTTGTGCTGGTATTTGAGAAAGTAGGGATCCCGTACGCCGCTGATATTTTTAACTTCGTGATCCTGACGGCTATTCTTTCTGCAGCGAACTCCGGGTTATATGCCTCCGGGCGCATGCTGTgg                                                           >  5:82748/1‑148 (MQ=255)
                                                                         cagGTGGGCGTGGAGAAAAGCCCGTTTGTGCTGGTATTTGAGAAAGTAGGGATCCCGTACGCCGCTGATATTTTTAACTTCGTGATCCTGACGGCTATTCTTTCTGCAGCGAACTCCGGGTTATATGCCTCCGGGCGCATGCTGTGGt                                                          >  5:240714/1‑148 (MQ=255)
                                                                              gggCGTGGAGAAAAGCCCGTTTGTGCTGGTATTTGAGAAAGTAGGGATCCCGTACGCCGCTGATATTTTTAACTTCGTGATCCTGACGGCTATTCTTTCTGCAGCGAACTCCGGGTTATATGCCTCCGGGCGCATGCTGTGGTCGTTg                                                     >  6:101117/1‑148 (MQ=255)
                                                                               ggCGTGGATAAACGCCCGGTTGGGCTGGGATTTGAGAAAGTAGGGATCCCGTCCGCCGCTGATATTTTTAACTTCGTGATCCTGACGGCTATTCTTTCTGCAGCGAACTCCGGGTTATATGCCTCCGGGCGCATGCTGTGGTCGTTGt                                                    <  1:252080/148‑1 (MQ=255)
                                                                                   tGGAGAAAAGCCCGTTTGTGCTGGTATTTGAGAAAGTAGGGATCCCGTACGCCGCTGATATTTTTAACTTCGTGATCCTGACGGCTATTCTTTCTGCAGCGAACTCCGGGTTATATGCCTCCGGGCGCATGCTGTGGTCGTTGt                                                    <  1:274483/144‑1 (MQ=255)
                                                                                   tGGAGAAAAGCCCGTTTGTGCTGGTATTTGAGAAAGTAGGGATCCCGTACGCCGCTGATATTTTTAACTTCGTGATCCTGACGGCTATTCTTTCTGCAGCGAACTCCGGGTTATATGCCTCCGGGCGCATGCTGTGGTCGTTGt                                                    >  2:274483/1‑144 (MQ=255)
                                                                                   tGGAGAAAAGCCCGTTTGTGCTGGTATTTGAGAAAGTAGGGATCCCGTACGCCGCTGATATTTTTAACTTCGTGATCCTGACGGCTATTCTTTCTGCAGCGAACTCCGGGTTATATGCCTCCGGGCGCATGCTGTGGTCGTTGTCGaa                                                <  4:263303/148‑1 (MQ=255)
                                                                                    ggAGAAAAGCCCGTTTGTGCTGGTATTTGAGAAAGTAGGGATCCCGTACGCCGCTGATATTTTTAACTTCGTGATCCTGACGGCTATTCTTTCTGCAGCGAACTCCGGGTTATATGCCTCCGGGCGCATGCTGTGGTCGTTGTCGa                                                 <  7:3593/146‑1 (MQ=255)
                                                                                    ggAGAAAAGCCCGTTTGTGCTGGTATTTGAGAAAGTAGGGATCCCGTACGCCGCTGATATTTTTAACTTCGTGATCCTGACGGCTATTCTTTCTGCAGCGAACTCCGGGTTATATGCCTCCGGGCGCATGCTGTGGTCGTTGTCGa                                                 >  8:3593/1‑146 (MQ=255)
                                                                                      agaAAAGCCCGTTTGTGCTGGTATTTGAGAAAGTAGTGATCCCGTACGCCGCTGATATTTTTAACTTCGTGATCCTGACGGCTATTCTTTCTGCAGCGAACTCCGGGTTATATGCCTCCGGGCGCATGCTGTGGTCGTTGTCGAATGa                                             <  5:101117/148‑1 (MQ=255)
                                                                                      agaAAAGCCCGTTTGTGCTGGTATTTGAGAAAGTAGGGATCCCGTACGCCGCTGATATTTTTAACTTCGTGATCCTGACGGCTATTCTTTCTGCAGCGAACTCCGGGTTATATGCCTCCGGGCGCATGCTGTGGTCGTTGTCGa                                                 >  1:170117/1‑144 (MQ=255)
                                                                                      agaAAAGCCCGTTTGTGCTGGTATTTGAGAAAGTAGGGATCCCGTACGCCGCTGATATTTTTAACTTCGTGATCCTGACGGCTATTCTTTCTGCAGCGAACTCCGGGTTATATGCCTCCGGGCGCATGCTGTGGTCGTTGTCGa                                                 <  2:170117/144‑1 (MQ=255)
                                                                                        taaaGCGCGTTTGTGCTGGTATTTGAGAATTTAGTGATCCCGTACGCCGCTGATATTTTTAACTTCGTGATCCTGACGGCTATTCTTTCTGCAGCGAACTCCGGGTTATATGCCTCCGGGCGCATGCtgtg                                                            <  1:50457/130‑1 (MQ=255)
                                                                                        aaaaGCCCGTTTGTGCTGGTATTTGAGAAAGTAGGGATCCCGTACGCCGCTGATATTTTTAACTTCGTGATCCTGACGGCTATTCTTTCTGCAGCGAACTCCGGGTTATATGCCTCCGGGCGCATGCTGTGGTCGTTGTCGAATGAAc                                           >  2:203249/1‑148 (MQ=255)
                                                                                         aaaGCCCGTTTGTGCTGGTATTTGAGAAAGTAGGGATCCCGTACGCCGCTGATATTTTTAACTTCGTGATCCTGACGGCTATTCTTTCTGCAGCGAACTCCGGGTTATATGCCTCCGGGCGCATGCTGTGGTCGTTGTCGAATGAACg                                          >  4:301598/1‑148 (MQ=255)
                                                                                         aaaGCCCGTTTGTGCTGGTATTTGAGAAAGTAGGGATCCCGTACGCCACTGATATTTTTAACTTCGTGATCCTGACGGCTATTCTTTCTGCAGCGAACTCCGGGTTATATGCCTCCGGgcgc                                                                    <  8:65282/122‑1 (MQ=255)
                                                                                         aaaGCCCGTTTGTGCTGGTATTTGAGAAAGTAGGGATCCCGTACGCCACTGATATTTTTAACTTCGTGATCCTGACGGCTATTCTTTCTGCAGCGAACTCCGGGTTATATGCCTCCGGgcgc                                                                    <  8:69788/122‑1 (MQ=255)
                                                                                         aaaGCCCGTTTGTGCTGGTATTTGAGAAAGTAGGGATCCCGTACGCCACTGATATTTTTAACTTCGTGATCCTGACGGCTATTCTTTCTGCAGCGAACTCCGGGTTATATGCCTCCGGgcgc                                                                    >  7:65282/1‑122 (MQ=255)
                                                                                         aaaGCCCGTTTGTGCTGGTATTTGAGAAAGTAGGGATCCCGTACGCCACTGATATTTTTAACTTCGTGATCCTGACGGCTATTCTTTCTGCAGCGAACTCCGGGTTATATGCCTCCGGgcgc                                                                    >  7:69788/1‑122 (MQ=255)
                                                                                            gCCCGTTTGTGCTGGTATTTGAGAAAGTAGGGATCCCGTACGCCGCTGATATTTTTAACTTCGTGATCCTGACGGCTATTCTTTCTGCAGCGAACTCCGGGTTATATGCCTCCGGGCGCATGCTGTGGTCGTTGTc                                                   <  3:176740/136‑1 (MQ=255)
                                                                                            gCCCGTTTGTGCTGGTATTTGAGAAAGTAGGGATCCCGTACGCCGCTGATATTTTTAACTTCGTGATCCTGACGGCTATTCTTTCTGCAGCGAACTCCGGGTTATATGCCTCCGGGCGCATGCTGTGGTCGTTGTc                                                   <  7:275165/136‑1 (MQ=255)
                                                                                            gCCCGTTTGTGCTGGTATTTGAGAAAGTAGGGATCCCGTACGCCGCTGATATTTTTAACTTCGTGATCCTGACGGCTATTCTTTCTGCAGCGAACTCCGGGTTATATGCCTCCGGGCGCATGCTGTGGTCGTTGTc                                                   >  8:275165/1‑136 (MQ=255)
                                                                                            gCCCGTTTGTGCTGGTATTTGAGAAAGTAGGGATCCCGTACGCCGCTGATATTTTTAACTTCGTGATCCTGACGGCTATTCTTTCTGCAGCGAACTCCGGGTTATATGCCTCCGGGCGCATGCTGTGGTCGTTGTc                                                   >  4:176740/1‑136 (MQ=255)
                                                                                              ccGTTTGTGCTGGTATTTGAGAAAGTAGGGCTCCCGTACGCCGCTGATATTTTTAACTTAGTGATCCTGACGGCTATTCTTTCTGCAGCGAACTCCGGGTTATATGCCTCCGGGCGCATGCTGTGGTCGTTGTCGAATGAACGTACGc                                     >  1:86954/1‑148 (MQ=255)
                                                                                              ccGTTTGTGCTGGTATTTGAGAAAGTAGGGATCCCGTACGCCGCTGATATTTTTAACTTCGTGATCCTGACg                                                                                                                 <  8:284687/72‑1 (MQ=255)
                                                                                              ccGTTTGTGCTGGTATTTGAGAAAGTAGGGATCCCGTACGCCGCTGATATTTTTAACTTCGTGATCCTGACg                                                                                                                 >  7:284687/1‑72 (MQ=255)
                                                                                                         ggTATTTGAGAAAGTAGGGATCCCGTACGCCGCTGATATTTTTAACTTCGTGATCCTGACGGCTATTCTTTCTGCAGCGAACTCCGGGTtata                                                                                 >  4:51991/1‑93 (MQ=255)
                                                                                                         ggTATTTGAGAAAGTAGGGATCCCGTACGCCGCTGATATTTTTAACTTCGTGATCCTGACGGCTATTCTTTCTGCAGCGAACTCCGGGTtata                                                                                 <  3:51991/93‑1 (MQ=255)
                                                                                                         ggTATTTGAGAAAGTAGGGATCCCGTACGCCGCTGATATTTTTAACTTCGTGATCCTGACGGCTATTCTTTCTGCAGCGAACTCCGGGTTATATGCCTCCGGGCGCATGCTGTGGTCGTTATCGAATGAACGTACGCTACCGGCCTGt                          >  7:272868/1‑148 (MQ=255)
                                                                                                          gTATTTGAGAAAGTAGGGATCCCGTACGCCGCTGATATTTTTAACTTCGTGATCCTGACGGCTATTCTTTCTGCAGCGAACTCCGGGTTATATGCCTCCGGGCGCATGCTGTGGTCGTTGTCGAATGAACGTACGCTACCGGCCTGtt                         >  8:325451/1‑148 (MQ=255)
                                                                                                                 agaAAGTAGGGATCCCGTACGCCGCTGATATTTTTAACTTCGTGATCCTGACGGCTATTCTTTCTGCAGCGAACTCCGGGTTATATGCCTCCGGGCGCATGCTGTGGTCGTTGTCGAATGAACGTACGCTACCGGCCTGTTTTGcgcg                  <  3:301598/148‑1 (MQ=255)
                                                                                                                  gaAAGTAGGGATCCCGTACGCCGCTGATATTTTTAACTTCGTGATCCTGACGGCTATTCTTTCTGCAGCGAACTCCGGGTTATATGCCTCCGGGCGCATGCTGTGGTCGTTGTCGAATGAACGTACGCTACCGGCCTGTTTTgcgc                   >  1:132854/1‑146 (MQ=255)
                                                                                                                  gaAAGTAGGGATCCCGTACGCCGCTGATATTTTTAACTTCGTGATCCTGACGGCTATTCTTTCTGCAGCGAACTCCGGGTTATATGCCTCCGGGCGCATGCTGTGGTCGTTGTCGAATGAACGTACGCTACCGGCCTGTTTTgcgc                   >  2:177411/1‑146 (MQ=255)
                                                                                                                  gaAAGTAGGGATCCCGTACGCCGCTGATATTTTTAACTTCGTGATACTGACGGCTATTCTTTCTGCAGCGAACTCAGGGTTATATGCCTCCGGGCGCATGCTGTGGTCGTTGTCGAATGAACGTACGCTACCGGCCTGTTTTgcgc                   >  7:115614/1‑146 (MQ=255)
                                                                                                                       tAGGGATCCCGTACGCCGCTGATATTTTTAACTTCGTGATCCTGACGGCTATTCTTTCTGCAGCGAACTCCGGGTTATATGCCTCCGGGCGCATGCTGTGGTCGTTGTCGAATGAACGTACGCTACCGGCCTGTTTTGcgcg                  >  6:206695/1‑142 (MQ=255)
                                                                                                                       tAGGGATCCCGTACGCCGCTGATATTTTTAACTTCGTGATCCTGACGGCTATTCTTTCTGCAGCGAACTCCGGGTTATATGCCTCCGGGCGCATGCTGTGGTCGTTGTCGAATGAACGTACGCTACCGGCCTGTTTTGcgcg                  <  5:206695/142‑1 (MQ=255)
                                                                                                                          ggATCCCGTACGCCGCTGATATTTTTAACTTCGTGATCCTGACGGCTATTCTTTCTGCAGCGAACTCCGGGTTATATGCCTCCGGGCGCATGCTGTGGTCGTTGTCGAATGAACGTACGCTACCGGCCTGTTTTGCGCGAGTAACGaa         <  6:126580/148‑1 (MQ=255)
                                                                                                                            aTCCCGTACGCCGCTGATATTTTTAACTTCGTGATCCTGACGGCTATTCTTTCTGCAGCGAACTCCGGGTTATATGCCTCCGGGCGCATGCTGTGGTCGTTGTCGAATGAACGTACGCTACCgg                               <  6:252209/124‑1 (MQ=255)
                                                                                                                            aTCCCGTACGCCGCTGATATTTTTAACTTCGTGATCCTGACGGCTATTCTTTCTGCAGCGAACTCCGGGTTATATGCCTCCGGGCGCATGCTGTGGTCGTTGTCGAATGAACGTACGCTACCgg                               >  5:252209/1‑124 (MQ=255)
                                                                                                                                cGTACGCCGCTGATATTTTTAACTTCGTGATCCTGACGGCTATTCTTTCTGCAGCGAACTCCGGGTTATATGCCTCCGGGCGCATGCTGTGGTCGTTGTCGAATGAACGTACGCTACCGGCCTGTTTTGCGCGAGTAACGAAAAACgg   <  4:109415/148‑1 (MQ=255)
                                                                                                                                 gTACCCCGCAGATATTTTTAACTTCGTGATCCTGACGGCTAATCTATCTGCAGCGAACTCCGGGATATATGCCTCCGGGCTCATGCTGTGGTCGTTGTCGAATTAACGTACGCTACCTGCATGTTTTGCGCGAGTAACGAAAAACGGc  <  6:230975/148‑1 (MQ=255)
                                                                                                                                        |                                                                                                                                            
TACTACCATCGCGCGACTGATTATTTTCTTTATCGGCACCGTGTTTGTGCTGGCAGCGCTGATCCCGATGCAGCAGGTGGGCGTGGAGAAAAGCCCGTTTGTGCTGGTATTTGAGAAAGTAGGGATCCCGTACGCCGCTGATATTTTTAACTTCGTGATCCTGACGGCTATTCTTTCTGCAGCGAACTCCGGGTTATATGCCTCCGGGCGCATGCTGTGGTCGTTGTCGAATGAACGTACGCTACCGGCCTGTTTTGCGCGAGTAACGAAAAACGGC  >  NC_000913/276050‑276326

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: ‑