Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,770,319 C→G 7.4% P206R (CCG→CGG)  yfjR → CP4‑57 prophage; putative DNA‑binding transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,770,3190CG 7.4% 160.2 / 3.9 65P206R (CCG→CGG) yfjRCP4‑57 prophage; putative DNA‑binding transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base C (30/30);  new base G (3/2);  total (33/32)
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.77e-01

CTCATCCACCAGCAGACCTGGTATCTGGTGGCTGAACACGAAGGGCATATCGCCGTATTCACACTTGATGAAATCCATCTGATTCAGCCTTTGCAGGAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTATCCAGCATTCACTGCTTACGTTTGTTCCGGCCGACAGCCCACCGGAATAGCGCAGGCGTTGTTATCAACCCGGCAACATGGAGGAGCCCAGTGCCCGTTATTGCCATT  >  NC_000913/2770177‑2770462
                                                                                                                                              |                                                                                                                                               
cTCATCCACCAGCAGACCTGGTATCTGGTGGCTGCCCACGAAGGGCATATCGCCGTATTCACACTTGATCAAATCCATCTGATTCAGCCTTTACAGGAGACTTGTCGCCGCAATGACAGGCTCTGTCGCCTGGTTGAAGACCCGGTCtt                                                                                                                                           >  1:115214/1‑149 (MQ=255)
cTCATCCACCAGCAGACCTGGTATCTGGTGGCTGAACACGAAGGGCATATCGCCGTATTCACACTTGATGAATTCCATCTGATTCAGCCTTTGCAGGAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCCgtgctt                                                                                                                                           >  2:84294/1‑144 (MQ=255)
 tCATCCACCAGCAGACCTGGTATCTGGTGGCTGAACACGACGGGCATATCGCCGTATTCACACTTGATGAAATCCATCTGATTCAGCCTTTGCAGGATACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAATACCCGGTCttc                                                                                                                                          >  2:434045/1‑149 (MQ=255)
   aTCCACCAGCAGACCTGGTATCTGGTGGCTGAACACGAAGGGCATATCGCCGTATTCACACTTGATGAAATCCATCTGATTCAGCCTTTGCAGGAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCttcat                                                                                                                                        <  1:950062/149‑1 (MQ=255)
   aTCCACCAGCAGACCTGGTATCTGGTGGCTGAACACGAAGGGCATATCGCCGTATTCACACTTGATGAAATCCATCTGATTCAGCCTTTGCAGGAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCttcat                                                                                                                                        >  2:369234/1‑149 (MQ=255)
   aTCCACCAGCAGACCTGGTATCTGGTGGCTGAACACGAAGGGCATATCGCCGTATTCACACTTGATGAAATCCATCTGATTCAGCCTTTGCAGGAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCttcat                                                                                                                                        >  1:672317/1‑149 (MQ=255)
   aTCCACCAGCAGACCTGGTATCTGGTGGCTGAACACGAAGGGCATATCGCCGTATTCACACTTGATGAAATCCATCTGATTCACCCTTTGCAGGAGACTTTTCGCAGCAATGACAGTCTGTGTCGCATCGTTGAAGACACGGTCttcat                                                                                                                                        >  1:937414/1‑149 (MQ=255)
     ccaccaGCAGACCTGGTATCTGGTGGCTGAACACGAAGGGCATATCGCCGTATTCACACTTGATGAAATCCATCTGATTCAGCCTTTGCAGGAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCttcattc                                                                                                                                      <  1:716583/149‑1 (MQ=255)
       accaGCAGACCTGGTATCTGGTGGCTGAACACGAAGTGCATATCGCCGTATTCACACTTGATGAAATCCATCTGATTCAGCCTTTGCAGGAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCTTCATTCAg                                                                                                                                    <  2:1004892/149‑1 (MQ=255)
       accaGCAGACCTGGTATCTGGTGGCTGAACACGAAGGGCATATCGCCGTATTCACACTTGATGAAATCCATCTGATTCAGCCTTTGCAGGAGCCTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCttcattcac                                                                                                                                    >  2:852845/1‑148 (MQ=255)
       accaGCAGACCTGGTATCTGGTGGCTGAACACGAAGGGCATATCGCCGTATTCACACTTGATGAAATCCAACTGATTCAGCCTTTGCAGGAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCTTCATTCAg                                                                                                                                    >  1:30264/1‑149 (MQ=255)
        ccaGCAGACCTGGTATCTGGTGGCTGAACACGAAGGGCATATCGCCGTATTCACACTTGATGAAATCCATCTGATTCAGCCTTTGCAGGAGCCTTTTCGCCGCAATGCCAGTCTGTGTCGCCTGGTTGAAGACCCGGTCTTCATTCAgg                                                                                                                                   >  1:139727/1‑149 (MQ=255)
         cagcagACCTGGTATCTGGTGGCTGAACACGAAGGGCATATCGCCGTATTCACACTTGATGAAATCCATCTGATTCAGCCTTTGCAGGAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGAATGAAGACCGGGTCTTCAACCAGGc                                                                                                                                  <  2:13948/149‑1 (MQ=255)
            cagACCTGGTATCTGGTGGCTGAACACGAAGGGCATATCGCCGTATTCACACTTGATGAAATCCATCTGATTCAGCCTTTGCAGGAGACTGTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCTTCATTCAGGCCtt                                                                                                                               >  2:813520/1‑149 (MQ=255)
                          ggtggCTGAACACGAAGGGCATATCGCCGTATTCACACTTGATGAAATCCTTCTGATTCAGCCTTTGCAGGAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCt                                                                                                                 <  2:737466/149‑1 (MQ=255)
                          ggtggCTGAACACGAAGGGCATATCGCCGTATTCACACTTGATGAAATCCATCTGATTCAGCCTTTGCAGGAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCt                                                                                                                 <  1:119442/149‑1 (MQ=255)
                            tggCTGAACACGAAGGGCATATCGCCGTATTCACACTTGATGAAATCCATCTGATTCAGCCTTTGCAGGAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCttt                                                                                                               <  1:569342/149‑1 (MQ=255)
                              gCTGAACACGAAGGGCATATCGCCGTATTCACACTTGATGAAATCCATCTGATTCAGCCTTTGCAGGAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGTCCCGGTCTTCATTCAGGCCTTACCCCCTTTTCGCTTTAt                                                                                                             >  1:667556/1‑149 (MQ=255)
                                tGAACACGAAGGGCATATCGCCGTATTCACACTTGATGAAATCCATCTGATTCAGCCTTTGCAGGAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTATcc                                                                                                           >  2:100598/1‑149 (MQ=255)
                                  aaCACGAAGGGCATATCGCCGTATTCACACTTGATGAAATCCATCTGATTCAGCCTTTGCAGGAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTATCCAg                                                                                                         <  2:1001299/149‑1 (MQ=255)
                                      cGAAGGGCATATCGCCGTATTCACACTTGATGAAATCCATCTGCTTCAGCCTTTGCAGGAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTATCCAGCAtt                                                                                                     <  2:220038/149‑1 (MQ=255)
                                        aaGGGCATATCGCCGTATTCACACTTGATGAAATCCATCTGATTCAGCCTTTGCAGGAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCTTCATTCAGGCCTTACACCATTTTCGCTTTATCCAGCATTCa                                                                                                   >  2:322619/1‑149 (MQ=255)
                                         aGGGCATATCGCCGTATTCACACTTGATGAAATCCATCTGATTCAGCCTTTGCAGGAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTATCCAGCATTCAc                                                                                                  <  2:443528/149‑1 (MQ=255)
                                           ggCATATCGCCGTATTCACACTTGATGAAATCCATCTGATTCAGCCTTTGCAGGAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGCCCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTATCCAGCATTCACTg                                                                                                <  2:616150/149‑1 (MQ=255)
                                                atCGCCGTATTCACACTTGATGAAATCCATCTGATTCAGCCTTTGCAGGAGACTTTTCGCCGCAATGACAGTCTTTGTCGCCTGGATGAAGACCCGGTCTACATTCAGGCCTTACCCCATTTTCGCTTTATCCAGCATTCACTGCATAc                                                                                           >  1:967902/1‑149 (MQ=255)
                                                atCGCCGTATTCACACTTGATGAAATCCATCTGATTCAGCCTTTGCAGGAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCTTCATTCAGGCGTTACCCCATTTTCGCTTTATCCAGCATTCACTGCTTAc                                                                                           >  1:398308/1‑149 (MQ=255)
                                                atCGCCGTATTCACACTTGATGAAATCCATCTGATTCAGCCTTTGCAGGAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGGTGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTATCCAGCATTCACTGCTTAc                                                                                           <  2:605625/149‑1 (MQ=255)
                                                  agCCGTATTCACACTTGATGAAATCCATCTGATTCATCCTTTGCAGGAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCTTCATTCAGTCCTTACCCCATTTTCGCTTTATCCAGCATTCACTGCTTACGt                                                                                         <  2:738586/148‑1 (MQ=255)
                                                    ccGTATTCACACTTGATGAAATCCATCTGATTCAGCCTTTGCAGGAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTATCCAGCATTCACTGCTTACGttt                                                                                       <  2:874594/149‑1 (MQ=255)
                                                     cGTATTCACACTTGATGAAATCCATCTGATTCAGGCTTTGCAGGAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTATCCAGCATTCACTGCTTACGTTTg                                                                                      <  2:194111/149‑1 (MQ=255)
                                                          tCACACTTGATGAAATCCATCTGATTCAGCCTTTGCAGGAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCTTCATTCAGGTCTTACCCCATTTTCGCTTTATCCAGCATTCACTGCTTACGTTTGTTCCg                                                                                 >  1:473910/1‑149 (MQ=255)
                                                            acacTTGATGAAATCCAACTGATTCAGCCATTGCAGGAGACTTTTCGCCGCAATGACAGTCTGTGACGCCTGGTAGAAGACCCGGTCTTCCTTCAGGCCATACCCAATTTTCGCTTTATCGATCATTCACTGATAACATTTATTcatgc                                                                               >  2:57608/1‑145 (MQ=255)
                                                               cTTTATGAAATCCATCTGATTCAGCCTTTGCAGGAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTATCCAGCATTCACTGCTTACGTTTGTTCCGGCCGa                                                                            <  2:114475/149‑1 (MQ=255)
                                                                     gaAATCCATCTGATTCAGCCTTTGCAGGAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTATCCAGCATTCACTGCTTACGTTTGTTCCGGCCGACAGccc                                                                      >  2:1023577/1‑149 (MQ=255)
                                                                         tCCATCTGATTCAGCCTTTGCAGGAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTATCCAGCATTCACTGCTTACGTTTGTTCCGGCCGACAGCCCACCg                                                                  <  1:830992/149‑1 (MQ=255)
                                                                         tCCATCTGATTCAGCCTTTGCAGGAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTATCCAGCATTCACTGCTTACGTTTGTTCCGGCCGACAGCCCACCg                                                                  <  1:156559/149‑1 (MQ=255)
                                                                           cATCTGATTCAGCCTTTGCAGGAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTATCCAGCATTCACTGCTTACGTTTGTTCCGGCCGACAGCCCACCGGa                                                                <  2:823752/149‑1 (MQ=255)
                                                                              gtGATTCAGCCTTTGCAGGAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGCCCCGGTCTTCGTTCAGGCCTTACGCCTTTTTAGCATTATCCATCATTCACTGCTTACGTTTGTGCCGGCCGACAGCCCACCGGAATa                                                             <  1:67310/148‑1 (MQ=255)
                                                                              cTGATTCAGCCTTTGCAGGAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTATCCAGCATTCACTGCTTACGTTTGTTCCGGCCGACAGCCCACCGGAATa                                                             <  1:912626/149‑1 (MQ=255)
                                                                                  ttCAGCCTTTGCAGGAGACTCTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCTGTATCCAGCATTCACTGCTGTCGTTTGGTCCGGCCGCCAGCCGCCCGGaagtgatc                                                         >  2:208097/1‑143 (MQ=255)
                                                                                      gCCTTTACAGGAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTATCCAGCATTCACTGCTTACGTTTGTTCCGGCCGACAGCCCACCGGAATAGCGCAGGc                                                     <  2:115214/149‑1 (MQ=255)
                                                                                          ttGCAGTAGACTTTTCGCCTCAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCTTCATTCTGGCCTTACCCCATTTTCGCTTTATCCAGCATTCACTGCTTTCGTTTGTTCCGGCCGACAGCCCACCGGAATAGCGCAGGCgttg                                                 <  1:575866/149‑1 (MQ=255)
                                                                                              aGGAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCTTCATTCAGGCCTTACACCATTTTCGCTTTATCCAGCATTAACTGCTTACGTTTGTTCCGGCCGCCAGCCCACCGGAATAGCGCAGGCGTTGTTAt                                             >  1:56373/1‑149 (MQ=255)
                                                                                              aGGAGACTTATCGCCGCAATGATAGTCGGTGTCGCCTGGTGTAAGAGCGGGTCTTCATTCAGGCAATACCCGATTTTAGCTTTCTCCAGCATTCACTGCTTACGTTTGTTACGGCCGACAGCCCACCGGAATAGCGCAGTCGTGGTTAt                                             <  2:767529/149‑1 (MQ=255)
                                                                                               ggAGACTTTTCGCCGCAATGACAGTCTGTGTCTCCTGGTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTATCCAGCATTCACTGCTTACGTTTGTTCCGGCCGACAGCCCACCGGAATAGCGCAGGCGTTGTTATc                                            <  2:937414/149‑1 (MQ=255)
                                                                                               ggAGACATACCGCCGCATATACAGTCTGTCTCGCCTGGGTGAACACCCGCGCATCAGTCAGGCCTTACCCCATTTTCGCTTACTCCACCATACACTGCTTCCGTTTGTTCCGGCCGACAGCCCACCGGAATAGCGCAGGCGTTGTAATc                                            <  2:667556/149‑1 (MQ=255)
                                                                                                gagaGTTTTCGCCGCAATGACAGTCTGAGTCGCCTGGATGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTATCCAGCATTCACTGCTTACGTTTGTTCCGGCCGACAGCCCACCGGAATAGCGCAGGCGTTGTTATCa                                           >  1:827615/1‑149 (MQ=255)
                                                                                                gagaCTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCGGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTATCCAGCATTCACTGCTTACGTTTGTTCCGGCCGACAGCCCACCGGAATAGCGCAGGCGTTGTTATCa                                           >  1:338416/1‑149 (MQ=255)
                                                                                                       ttCGCCGCAATGACATTCTGTGTCGCCTGGTTGAAGACCCGGTCTTCTTTCAGGCCTTACCCCATTTTCGCTTTATCCAGCAATCACTGCTTACGTTTGTTCCGGCCGACAGCCCACCGGAATAGCGCAGGCGTTGTTATCAACCCGGc                                    <  2:997075/149‑1 (MQ=255)
                                                                                                             gcAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCTTCATTCAGGGCTTACCCCATTTTCGCTTTATCCAGCATTCACTGCTTACGTTTGTTCCGGCCGACAGCCCACCGGAATAGCGCAGGCGTTGTTATCAACCCGGCAACATg                              <  2:473910/149‑1 (MQ=255)
                                                                                                             gcAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTATCCAGCATTCACTGCTTACGTTTGTTCCGGCCGACAGCCCACCGGAATAGCGCAGGCGTTGTTATCAACCCGGCAACATg                              <  2:30264/149‑1 (MQ=255)
                                                                                                              cAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCAATTTCGCTTTATCCAGCATTCACTGCTTACGTTTGTTCCGGCCGACAGCCCACCGGAATAGCGCAGGCGTTGTTATCAACCCGGAAACATgg                             >  2:318389/1‑149 (MQ=255)
                                                                                                               aaTGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTATCCAGCATTCACTGCTTACGTTTGTTCCGGCCGACAGCCCACCGGAATAGCGCAGGCGTTGTTATCAACCCGGCAACATgga                            >  1:776918/1‑149 (MQ=255)
                                                                                                                 tGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTATCCAGCATTCACTGCTTACGTTTGTTCCGGCCGACAGCCCACCGGAATAGCGCAGGCGTTGTTATCAACCCGGCAACATggagg                          <  1:322619/149‑1 (MQ=255)
                                                                                                                     aGTCTGTGTCGCCTGGTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTATCCAGCATTCACTGCTTACGTTTGTTCCGGCCGACAGCCCACCGGAATAGCGCAGGCGTTGTTATCAACCCGGCAACATGGAGGAGcc                      >  2:360674/1‑149 (MQ=255)
                                                                                                                         ggtaaggCCTGAATGAAGACCGGGTCGTCAACCAGGCCTTACCCCAGTTTCGCTTTATCCAGCATTCACTGCTTACGTTTGTTCCGGCCGACAGCCCACCGGAATAGCGCAGGCGTTGTTATCAACCCGGCAACATGGAGGAGCCCAGt                  >  1:64735/7‑149 (MQ=255)
                                                                                                                             tCGCCTGGTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTATCCAGCATTCAGTGCTTACGTTTGTTCCGGCCGACAGCCCACCGGAATAGCGCAGGCGTTGTTATCAACCCGGCAACATGGAGGAGCCCAGTGccc              >  1:136739/1‑149 (MQ=255)
                                                                                                                               gCCTCGTTGAAGACACTGTCTTAAGTCAGGCCTGCCCCCATTTTCGCTTTATCCAGCATACACTGCTAACGTTTGTACCCGCCGAAAGCCCACCGGAATATCGCAGGCCTTGATATCAACCCGGAACCCTGACGGAGCCCAGTGCCCGt            >  2:463912/1‑149 (MQ=255)
                                                                                                                                ccTGGTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTATCCAGCATTCACTGCTTACGTTTGTTCCGGCCGACAGCCCACCGGAATAGCGCAGGCGTTGTTAGCAACCCGGCAACATGGAGGAGCCCAGTGCCCGtt           >  2:422925/1‑149 (MQ=255)
                                                                                                                                ccTGCTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCTTTTTCGCGTTATCCAGCATTCACTGCTTACGTTTGTTCCGGCGGACAGTCCACAGGAATAGCGCAGGCGTAGTTATAGACCCGCAATCATGGAGGAGCCCAGTGCCAGtt           >  2:59879/1‑149 (MQ=255)
                                                                                                                                 ctgaatGAAGACCGGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTATCCAGCATTCACTGCTTACGTTTGTTCCGGCCGACAGCCCACCGGAATAGCGCAGGCGTTGTTATCAACCCGGCAACATGGAGGAGCCCAGTGCCCGttt          >  1:715607/6‑148 (MQ=255)
                                                                                                                                  tGGTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTATCCAGCATTCACTGCTTACGTTTGTTCCGGCCGACAGCCCACCGGAATAGCGCAGGCGTTGTTATCAACCCGGCAACATGGAGGAGCCCAGTGCCCGTTAt         <  1:360674/149‑1 (MQ=255)
                                                                                                                                       taaGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTATCCAGCATTCACTGCTTACGTTTGTTCCGGCCGACAGCCCACCGGAATAGCGCAGGCGTTGTTATCAACCCGGCAACATGGAGGAGCCCAGTGCCCGTTATTGCCa    <  2:139727/148‑1 (MQ=255)
                                                                                                                                       gAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTATCCAGCATTCACTGCTTACGTTTGTTCCGGCCGACAGCCCACCGGAATAGCGCAGGCGTTGTTATCAACCCGGCAACATGGAGGCGCCCAGTGCCCGTTTTTGCCa    >  1:963355/1‑149 (MQ=255)
                                                                                                                                         aGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTATCCAGCATTCACTCCTTACGTTTGTTCCGGCCGACAGCCCACCGGAATAGCGCAGGCGTTGTGATCAACCCGGCAACATGGAGGAGCCCAGTGCCCGTGATTGCCAtt  >  1:478582/1‑149 (MQ=255)
                                                                                                                                              |                                                                                                                                               
CTCATCCACCAGCAGACCTGGTATCTGGTGGCTGAACACGAAGGGCATATCGCCGTATTCACACTTGATGAAATCCATCTGATTCAGCCTTTGCAGGAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTATCCAGCATTCACTGCTTACGTTTGTTCCGGCCGACAGCCCACCGGAATAGCGCAGGCGTTGTTATCAACCCGGCAACATGGAGGAGCCCAGTGCCCGTTATTGCCATT  >  NC_000913/2770177‑2770462

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