Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 885,494 A→C F213C (TTT→TGT)  ybjI ← 5‑amino‑6‑(5‑phospho‑D‑ribitylamino)uracil phosphatase; pyrimidine phosphatase; riboflavin synthesis

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913885,4940AC99.2% 453.0 / ‑4.3 127F213C (TTT→TGT) ybjI5‑amino‑6‑(5‑phospho‑D‑ribitylamino)uracil phosphatase; pyrimidine phosphatase; riboflavin synthesis
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base C (58/68);  total (58/69)
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.82e-01

Only 100 of 127 total aligned reads displayed.

GTATTTTGCCGCTGCGACGACCGCGCTACCGGCATTTTCCATTGCAAAACTAAAGCCTGCCTGACGCAGCATCTCAATATCGTTACCGCCATCGCCAAAGACCACCACTTCGCTGTCGTCTATTCCCCATAATTTCTGCAGTTGGCGAAGGCCATTGGCTTTATGTACGCCGGGGATAATCAGATCGATGCTGC  >  NC_000913/885397‑885590
                                                                                                 |                                                                                                
gTATTTTGCCGCTGCGACGACCGCGCTACCGGCATTTTCCATTGCAAAACTAAAGCCTGCCTGACGCAGCATCTCAATATCGTTACCGCCATCGCCAAAg                                                                                                <  1:440970/100‑1 (MQ=255)
 tATTTTGCCGCTGCGACGACCGCGCTACCGGCATTTTCCATTGCAAAACTAAAGCCTGCCTGACGCAGCATCTCAATATCGTTACCGCCATCGCCACAGa                                                                                               <  2:3891941/100‑1 (MQ=255)
 tATTTTGCCGCTGCGACGACCGCGCTACCGGCATTTTCCATTGCAAAACTAAAGCCTGCCTGACGCAGCATCTCAATATCGTTACCGCCATCGCCACAGa                                                                                               >  2:23114/1‑100 (MQ=255)
  aTTTTGCCGCTGCGACGACCGCGCTACCGGCATTTTCCATTGCAAAACTAAAGCCTGCCTGACGCAGCATCTCAATATCGTTACCGCCATCGCCACAGac                                                                                              >  2:1974502/1‑100 (MQ=255)
  aTTTTGCCGCTGCGACGACCGCGCTACCGGCATTTTCCATTGCAAAACTAAAGCCTGCCTGACGCAGCATCTCAATATCGTTACCGCCATCGCCACAGac                                                                                              >  2:2300539/1‑100 (MQ=255)
   tgttGCCGCTGCGACGACCGCGCTACCGGCATTTTCCAGTGCAAATCTAAGGCCTGCCTGACGGAGCATCTCAAGATCGTTACCGCCATCGCCACAGacc                                                                                             <  2:1557079/98‑1 (MQ=255)
    tttGCCGCTGCGACGACCGCGCTACCGGCATTTTCCATTGCAAAACTAAAGCCTGCCTGACGCAGCATCTCAATATCGTTACCGCCATCGCCACAGacca                                                                                            <  1:2622281/100‑1 (MQ=255)
     ttGCCGCTGCGACGACCGCGCTACCGGCATTTTCCATTGCAAAACTAAAGCCTGCCTGACGCAGCATCTCAATATCGTTACCGCCATCGCCACAGaccac                                                                                           <  1:910980/100‑1 (MQ=255)
        ccGCTGCGACGACCGCGCTACCGGCATTTTCCATTGCAAAACTAAAGCCTGCCTGACGCAGCATCTCAATATCGTTACCGCCATCGCCACAGACcaccac                                                                                        <  1:1248579/100‑1 (MQ=255)
        ccGCTGCGACGACCGCGCTACCGGCATTTTCCATTGCAAAACTAAAGCCTGCCTGACGCAGCATCTCAATATCGTTACCGCCATCGCCACAGACcaccac                                                                                        <  2:876572/100‑1 (MQ=255)
          gCTGCGACGACCGCGCTACCGGCATTTTCCATTGCAAAGCTAAAGCCTGCCTGACGCAGCATCTCAATATCGTTACCGCCATCGCCACAGACCACCACtt                                                                                      <  2:2637475/100‑1 (MQ=255)
          gCTGCGACGACCGCGCTACCGGCATTTTCCATTGCAAAACTAAAGCCTGCCTGACGCAGCATCTCAATATCGTTACCGCCATCGCCACAGACCACCACtt                                                                                      >  2:886196/1‑100 (MQ=255)
          gCTGCGACGACCGCGCTACCGGCATTTTCCATTGCAAAACTAAAGCCTGCCTGACGCAGCATCTCAATATCGTTACCGCCATCGCCACAGACCACCACtt                                                                                      >  2:3506808/1‑100 (MQ=255)
           cTGCGACGACCGCGCTACCGGCATTTTCCATTGCAAAACTAAAGCCTGCCTGACGCAGCATCTCAATATCGTTACCGCCATCGCCACAGACCACCACTTc                                                                                     <  2:3083305/100‑1 (MQ=255)
            tGCGACGACCGCGCTACCGGCATTTTCCATTGCAAAACTAAAGCCTGCCTGACGCAGCATCTCAATATCGTTACCGCCATCGCCACAGACCACCACTTCg                                                                                    <  1:3546756/100‑1 (MQ=255)
             gCGACGACCGCGCTACCGGCATTTTCCATTGCAAAACTAAAGCCTGCCTGACGCAGCATCTCAATATCGTTACCGCCATCGCCACAGACCACCACTTCGc                                                                                   <  2:1406705/100‑1 (MQ=255)
               gacgacCGCGCTACCGGCATTTTCCATTGCAAAACTAAAGCCTGCCTGACGCAGCATCTCACTATCGTTACCGCCATCGCCACCGACCACCACTTTGCTg                                                                                 <  2:3801580/100‑1 (MQ=255)
               gacgacCGCGCTACCGGCATTTTCCATTGCAAAACTAAAGCCTGCCTGACGCAGCATCTCAATATCGTTACCGCCATCGCCACAGACCACCACTTCGCTg                                                                                 <  2:3000044/100‑1 (MQ=255)
               gacgacCGCGCTACCGGCATTTTCCATTGCAAAACTAAAGCCTGCCTGACGCAGCATCTCAATATCGTTACCGCCATCGCCACAGACCACCACTTCGCTg                                                                                 >  2:1788611/1‑100 (MQ=255)
                 cgacCGCGCTACCGGCATTTTCCATTGCAAAACTAAAGCCTGCCTGACGCAGCATCTCAATATCGTTACCGCCATCGCCACAGACCACCACTTCGCTgtc                                                                               >  1:1701027/1‑100 (MQ=255)
                  gacCGCGCTACCGGCATTTTCCATTGCAAAACTAAAGCCTGCCTGACGCAGCATCTCAATATCGTTACCGCCATCGCCACAGACCACCACTTCGCTgtcg                                                                              <  1:3710593/100‑1 (MQ=255)
                  gacCGCGCTACCGGCATTTTCCATTGCAAAACTAAAGCCTGCCTGACGCAGCATCTCAATATCGTTACCGCCATCGCCACAGACCACCACTTCGCTgtcg                                                                              >  1:3622447/1‑100 (MQ=255)
                    ccGCGCTACCGGCATTTTCCATTGCAAAACTAAAGCCTGCCTGACGCAGCATCTCAATATCGTTACCGCCATCGCCACAGACCACCACTTCGCTgtcgtc                                                                            <  2:3455205/100‑1 (MQ=255)
                     cgcgCTACCGGCATTTTCCATTGCAAAACTAAAGCCTGCCTGACGCAGCATCTCAATATCGTTACCGCCATCGCCACAGACCACCACTTCGCTGTCGTCt                                                                           <  1:175379/100‑1 (MQ=255)
                     cgcgCTACCGGCATTTTCCATTGCAAAACTAAAGCCTGCCTGACGCAGCATCTCAATATCGTTACCGCCATCGCCACAGACCACCACTTCGCTGTCGTCt                                                                           >  2:3403319/1‑100 (MQ=255)
                          cagcGGCATTTTCCATTGCAAAACTAAAGCCTGCCTGACGCAGCATCTCAATATCGTTACCGCCATCGCCACAGACCACCACTTCGCTGTCGTCTATTcc                                                                      >  1:1597525/4‑100 (MQ=255)
                           aCCGGCATTTTCCATTGCAAATCTAAAGCCTGCCTGACGCAGCATCTCAATAGCGTTACCGCCATCGCCACAGACCACCACTTCGCTGTCGTCTATTccc                                                                     <  1:1783676/100‑1 (MQ=255)
                           aCCGGCATTTTCCATTGCAAAACTAAAGCCTGCCTGACGCAGCATCTCAATATCGTTACCGCCATCGCCACAGACCACCACTTCGCTGTCGTCTATTccc                                                                     <  2:111728/100‑1 (MQ=255)
                           aCCGGCATTTTCCATTGCAAAACTAAAGCCTGCCTGACGCAGCATCTCAATATCGTTACCGCCATCGCCACAGACCACCACTTCGCTGTCGTCTATTccc                                                                     <  2:3494602/100‑1 (MQ=255)
                            ccGGCATTTTCCATTGCAAAACTAAAGCCTGCCTGACGCAGCATCTCAATATCGTTACCGCCATCGCCACAGACCACCACTTCGCTGTCGTCTATTcccc                                                                    >  1:3064926/1‑100 (MQ=255)
                            ccGGCATTTTCCATTGCAAAACTAAAGCCTGCCTGACGCAGCATCTCAATATCGTTACCGCCATCGCCACAGACCACCACTTCGCTGTCGTCTATTcccc                                                                    >  1:1544624/1‑100 (MQ=255)
                            ccGGCATTTTCCATTGCAAAACTAAAGCCTGCCTGACGCAGCATCTCAATATCGTTACCGCCATCGCCACAGACCACCACTTCGCTATCGTCTATTcccc                                                                    >  2:3130410/1‑100 (MQ=255)
                              ggCATTTTCCATTGCAAAACTAAAGCCTGCCTGACGCAGCATCTCAATATCGTTACCGCCATCGCCACAGACCACCACTTCGCTGTCGTCTATTCCCCAt                                                                  <  1:1330015/100‑1 (MQ=255)
                              ggCATTTTCCATTGCAAAACTAAAGCCTGCCTGACGCAGCATCTCAATATCGTTACCGCCATCGCCACAGACCACCACTTCGCTGTCGTCTATTCCCCAt                                                                  <  1:1702432/100‑1 (MQ=255)
                                cATTTTCCATTGCAAAACTAAAGCCTGCCTGACGCAGCATCTCAATATCGTTACCGCCATCGCCACAGACCACCACTTCGCTGTCGTCTATTCCCCATaa                                                                >  2:1479224/1‑100 (MQ=255)
                                 aTTTTCCATTGCAAAACTAAAGCCTGCCTGACGCAGCATCTCAATATCGTTACCGCCATCGCCACAGACCACCACTTCGCTGTCGTCTATTCCCCATAAt                                                               >  1:2464819/1‑100 (MQ=255)
                                    ttCCATTGCAAAACTAAAGCCTGCCTGACGCAGCATCTCAATATCGTTACCGCCATCGCCACAGACCACCACTTCGCTGTCGTCTATTCCCCATAATTTc                                                            <  1:107899/100‑1 (MQ=255)
                                    ttCCATTGCAAAACTAAAGCCTGCCTGACGCAGCATCTCAATATCGTTACCGCCATCGCCACAGACCACCACTTCGCTGTCGTCTATTCCCCATAATTTc                                                            <  2:3206721/100‑1 (MQ=255)
                                    ttCCATTGCAAAACTAAAGCCTGCCTGACGCAGCATCTCAATATCGTTACCGCCATCGCCACAGACCACCACTTCGCTGTCGTCTATTCCCCATAATTTc                                                            >  2:2601783/1‑100 (MQ=255)
                                    ttCCATTGCAAAACTAAAGCCTGCCTGACGCAGCATCTCAATATCGTTACCGCCATCGCCACAGACCACCACTTCGCTGTCGTCTATTCCCCATAATTTc                                                            >  1:1906584/1‑100 (MQ=255)
                                      ccATTGCAAAACTAAAGCCTGCCTGACGCAGCATCTCAATATCGTTACCGCCATCGCCACAGACCACCACTTCGCTGTCGTCTATTCCCCATAATTTCTg                                                          >  2:528881/1‑100 (MQ=255)
                                       cATTGCAAAACTAAAGCCTGCCTGACGCAGCATCTCAATATCGTTACCGCCATCGCCACAGACCACCACTTCGCTGTCGTCTATTCCCCATAATTTCTGc                                                         <  1:2715944/100‑1 (MQ=255)
                                        atcgCAAAACTAAAGCCTGCCTGACGCAGCATCTCAATATCGTTACCGCCATCGCCACAGACCACCACTTCGCTGTCGTCTATTCCCCATAATTTCTGCa                                                        <  2:1296732/97‑1 (MQ=255)
                                        aTTGCAAAACTAAAGCCTGCCTGACGCAGCATCTCAATATCGTTACCGCCATCGCCACAGACCACCACTTCGCTGTCGTCTATTCCCCATAATTTCTGCa                                                        <  1:3635072/100‑1 (MQ=255)
                                        aTTGCAAAACTAAAGCCTGCCTGACGCAGCATCTCAATATCGTTACCGCCATCGCCACAGACCACCACTTCGCTGTCGTCTATTCCCCATAATTTCTGCa                                                        <  2:3776791/100‑1 (MQ=255)
                                        aTTGCAAAACTAAAGCCTGCCTGACGCAGCATCTCAATATCGTTACCGCCATCGCCACAGACCACCACTTCGCTGTCGTCTATTCCCCATAATTTCTGCa                                                        <  2:60916/100‑1 (MQ=255)
                                           gCAAAACTAAAGCCTGCCTGACGCAGCATCTCAATATCGTTACCGCCATCGCCACAGACCACCACTTCGCTGTCGTCTATTCCCCATAATTTCTGCAGtt                                                     >  1:966396/1‑100 (MQ=255)
                                           gCAAAACTAAAGCCTGCCTGACGCAGCATCTCAATATCGTTACCGCCATCGCCACAGACCACCACTTCGCTGTCGTCTATTCCCCATAATTTCTGCAGtt                                                     >  2:1495067/1‑100 (MQ=255)
                                            cAAAACTAAAGCCTGCCTGACGCAGTATCTCAATATCGTTACCGCCATCGCCACAGACCACCACTTCGCTGTCGTCTATTCCCCATAATTTCTGCAGTTg                                                    <  2:2964477/100‑1 (MQ=255)
                                            cAAAACTAAAGCCTGCCTGACGCAGCATCTCAATATGGTTCCCGCCATCGCCACGGACCACCACTTCGCTGTCGTCTATTCCCCATAATTTCTGCAGTTg                                                    >  2:3749835/1‑100 (MQ=255)
                                            cAAAACTAAAGCCTGCCTGACGCAGCATCTCAATATCGTTACCGCCATCGCCACAGACCACCACTTCGCTGTCGTCTATTCCCCATAATTTCTGCAGTTg                                                    <  2:2764567/100‑1 (MQ=255)
                                            cAAAACTAAAGCCTGCCTGACGCAGCATCTCAATATCGTTACCGCCATCGCCACAGACCACCACTTCGCTGTCGTCTATTCCCCATAATTTCTGCAGTTg                                                    >  2:1158174/1‑100 (MQ=255)
                                            cAAAACTAAAGCCTGCCTGACGCAGCATCTCAATATCGTTACCGCCATCGCCACAGACCACCACTTCGCTGTCGTCTATTCCCCATAATTTCTGCAGTTg                                                    >  1:398366/1‑100 (MQ=255)
                                             aaaaCTAAAGCCTGCCTGACGCAGCATCTCAATATCGTTACCGCCATCGCCACAGACCACCACTTCGCTGTCGTCTATTCCCCATAATTTCTGCAGTTgg                                                   >  2:738041/1‑100 (MQ=255)
                                             aaaaCTAAAGCCTGCCTGACGCAGCATCTCAATATCGTTACCGCCATCGCCACAGACCACCACTTCGCTGTCGTCTATTCCCCATAATTTCTGCAGTTgg                                                   >  1:2516481/1‑100 (MQ=255)
                                               aaCTAAAGCCTGCCTGACGCAGCATCTCAATATCGTTACCGCCATCGCCACAGACCACCACTTCGCTGTCGTCTATTCCCCATAATTTCTGCAGTTGGCg                                                 >  1:834300/1‑100 (MQ=255)
                                               aaCTAAAGCCTGCCTGACGCAGCATCTCAATATCGTTACCGCCATCGCCACAGACCACCACTTCGCTGTCGTCTATTCCCCATAATTTCTGCAGTTGGCg                                                 >  2:3985767/1‑100 (MQ=255)
                                               aaCTAAAGCCTGCCTGACGCAGCAGCTTAGGTTCGTTACTTCCATCGCCACAGACCACCACTTCGCTGTCGTCTATTCCCCATAATTTCTGCAGTTGGCg                                                 <  2:3807994/100‑1 (MQ=255)
                                                aCTAAAGCCTGCCTGACGCAGCATCTCAATATCGTTACCGCCATCGCCACAGACCACCACTTCGCTGTCGTCTATTCCCCATAATTTCTGCAGTTGGCGa                                                >  2:1014018/1‑100 (MQ=255)
                                                 cTAAAGCCTGCCTGACGCAGCATCTCAATATCGTTACCGCCATCGCCACAGACCACCACTTCGCTGTCGTCTATTCCCCATAATTTCTGCAGTTGGCGaa                                               <  2:497480/100‑1 (MQ=255)
                                                   aaaGCCTGCCTGACGCAGCATCTCAATATCGTTACCGCCATCGCCACAGACCACCACTTCGCTGTCGTCTATTCCCCATAATTTCTGCAGTTGGCGAAgg                                             <  2:56379/100‑1 (MQ=255)
                                                    aaGCCTGCCTGACGCAGCATCTCAATATCGTTACCGCCATCGCCACAGACCACCACTTCGCTGTCGTCTATTCCCCATAATTTCTGCAGTTGGCGAAGGc                                            <  1:2520318/100‑1 (MQ=255)
                                                         tgcctgACGCAGCATCTCAATATCGTTACCGCCATCGCCACAGACCACCACTTCGCTGTCGTCTATTCCCCATAATTTCTGCAGTTGGCGAAGGCCATTg                                       >  2:2110731/1‑100 (MQ=255)
                                                         tgcctgACGCAGCATCTCAATATCGTTACCGCCATCGCCACAGACCACCACTTCGCTGTCGTCTATTCCCCATAATTTCTGCAGTTGGCGAAGGCCATTg                                       <  1:2755987/100‑1 (MQ=255)
                                                         tgcctgACGCAGCATCTCAATATCGTTACCGCCATCGCCACAGACCACCACTTCGCTGTCGTCTATTCCCCATAATTTCTGCAGTTGGCGAAGGCCATTg                                       <  1:2336206/100‑1 (MQ=255)
                                                         tgcctgACGCAGCATCTCAATATCGTTACCGCCATCGCCACAGACCACCACTTCGCTGTCGTCTATTCCCCATAATTTCTGCAGTTGGCGAAGGCCATTg                                       <  1:1339622/100‑1 (MQ=255)
                                                           cctgACGCAGCATCTCAATATCGTTACCGCCATCGCCACAGACCACCACTTCGCTGTCGTCTATTCCCCATAATTTCTGCAGTTGGCGAAGGCCATTGGc                                     <  2:925077/100‑1 (MQ=255)
                                                             tgACGCAGCATCTCAATATCGTTACCGCCATCGCCACAGACCACCACTTCGCTGTCGTCTATTCCCCATAATTTCTGCAGTTGGCGAAGGCCATTGGCtt                                   >  2:2137389/1‑100 (MQ=255)
                                                             tgACGCAGCATCTCAATATCGTTACCGCCATCGCCACAGACCACCACTTCGCTGTCGTCTATTCCCCATAATTTCTGCAGTTGGCGAAGGCCATTGGCtt                                   <  2:1992133/100‑1 (MQ=255)
                                                             tgACGCAGCATCTCAATATCGTTACCGCCATCGCCACAGACCACCACTTCGCTGTCGTCTATTCCCCATAATTTCTGCAGTTGGCGAAGGCCATTGGCtt                                   >  2:574091/1‑100 (MQ=255)
                                                             tgACGCAGCATCTCAATATCGTTACCGCCATCGCCACAGACCACCACTTCGCTGTCGTCTATTCCCCATAATTTCTGCAGTTGGCGAAGGCCATTGGCtt                                   <  2:191777/100‑1 (MQ=255)
                                                               aCGCAGCATCTCAATATCGTTACCGCCATCGCCACAGACCACCACTTCGCTGTCGTCTATTCCCCATAATTTCTGCAGTTGGCGAAGGCCATTGGCTTTa                                 <  1:3475355/100‑1 (MQ=255)
                                                                cGCAGCATCTCAATATCGTTACCGCCATCGCCACAGACCACCACTTCGCTGTCGTCTATTCCCCATAATTTCTGCAGTTGGCGAAGGCCATTGGCTTTAt                                >  2:3839075/1‑100 (MQ=255)
                                                                    gcaTCTCAATATCGTTACCGCCATCGCCACAGACCACCACTTCGCTGTCGTCTATTCCCCATAATTTCTGCAGTTGGCGAAGGCCATTGGCTTTATGTAc                            <  2:1664859/100‑1 (MQ=255)
                                                                     caTCTCAATATCGTTACCGCCATCGCCACAGACCACCACTTCGCTGTCGTCTATTCCCCATAATTTCTGCAGTTGGCGAAGGCCATTGGCTTTATGTACg                           <  2:2952530/100‑1 (MQ=255)
                                                                      aTCTCAATATCGTTACCGCCATCGCCACAGACCACCACTTCGCTGTCGTCTATTCCCCATAATTTCTGCAGTTGGCGAAGGCCATTGGCTTTATGTACGc                          >  2:1018926/1‑100 (MQ=255)
                                                                       tctcAATATCGTTACCGCCATCGCCACAGACCACCACTTCGCTGTCGTCTATTCCCCATAATTTCTGCAGTTGGCGAAGGCCATTGGCTTTATGTACGcc                         <  2:807524/100‑1 (MQ=255)
                                                                       tctcAATATCGTTACCGCCATCGCCACAGACCACCACTTCGCTGTCGTCTATTCCCCATAATTTCTGCAGTTGGCGAAGGCCATTGGCTTTATGTACGcc                         <  1:3282247/100‑1 (MQ=255)
                                                                       tctcAATATCGTTACCGCCATCGCCACAGACCACCACTTCGCTGTCGTCTATTCCCCATAATTTCTGCAGTTGGCGAAGGCCATTGGCTTTATGTACGcc                         <  1:3932907/100‑1 (MQ=255)
                                                                        ctcAATATCGTTACCGCCATCGCCACAGACCACCACTTCGCTGTCGTCTATTCCCCATAATTTCTGCAGTTGGCGAAGGCCATTGGCTTTATGTACGCCg                        >  2:2500537/1‑100 (MQ=255)
                                                                          cAATATCGTTACCGCCATCGCCACAGACCACCACTTCGCTGTCGTCTATTCCCCATAATTTCTGCAGTTGGCGAAGGCCATTGGCTTTATGTACGCCggg                      >  2:2445508/1‑100 (MQ=255)
                                                                           aaTATCGTTACCGCCATCGCCACAGACCACCACTTCGCTGTCGTCTATTCCCCATAATTTCTGCAGTTGGCGAAGGCCATTGGCTTTATGTACGCCgggg                     <  1:1180908/100‑1 (MQ=255)
                                                                             tatCGTTACCGCCATCGCCACAGACCACCACTTCGCTGTCGTCTATTCCCCATAATTTCTGCAGTTGGCGAAGGCCATTGGCTTTATGTACGCCGGGGAt                   <  2:1889438/100‑1 (MQ=255)
                                                                                cGTTACCGCCATCGCCACAGACCACCACTTCGCTGTCGTCTATTCCCCATAATTTCTGCAGTTGGCGAAGGCCATTGGCTTTATGTACGCCGGGGATAAt                >  2:75178/1‑100 (MQ=255)
                                                                                  ttACCGCCATCGCCACAGACCACCACTTCGCTGTCGTCTATTCCCCATAATTTCTGCAGTTGGCGAAGGCCATTGGCTTTATGTACGCCGGGGATAATCa              >  1:563754/1‑100 (MQ=255)
                                                                                   tACCGCCATCGCCACAGACCACCACTTCGCTGTCGTCTATTCCCCATAATTTCTGCAGTTGGCGAAGGCCATTGGCTTTATGTACGCCGGGGATAAGCAg             >  2:724674/1‑100 (MQ=255)
                                                                                     ccGCCATCGCCACAGACCACCACTTCGCTGTCGTCTATTCCCCATAATTTCTGCAGTTGGCGAAGGCCATTGGCTTTATGTACGCCGGGGATAATCAGAt           >  2:1947695/1‑100 (MQ=255)
                                                                                      cGCCATCGCCACAGACCACCACTTCGCTGTCGTCTATTCCCCTTAATTTCTGCAGTTGGCGAAGGCCATTGGTTTTATGTACGCCGGGGATAATCAGATc          <  1:3449187/100‑1 (MQ=255)
                                                                                      cGCCATCGCCACAGACCACCACTTCGCTGTCGTCTATTCCCCATAATTTCTGCAGTTGGCGAAGGCCATTGGCTTTATGTACGCCGGGGATAATCAGATc          <  1:2214642/100‑1 (MQ=255)
                                                                                      cGCCATCGCCACAGACCACCACTTCGCTGTCGTCTATTCCCCATAATTTCTGCAGTTGGCGAAGGCCATTGGCTTTATGTACGCCGGGGATAATCAGATc          <  1:2539534/100‑1 (MQ=255)
                                                                                       gCCATCGCCACAGACCACCACTTCGCTGTCGTCTATTCCCCATAATTTCTGCAGTTGGCGAAGGCCATTGGCTTTATGTACGCCGGGGATAATCAGCTCg         <  2:3672328/100‑1 (MQ=255)
                                                                                       gCCATCGCCACAGACCACCACTTCGCTGTCGTCTATTCCCCATAATTTCTGCAGTTGGCGAAGGCCATTGGCTTTATGTACGCCGGGGATAATCAGATCg         <  2:1524444/100‑1 (MQ=255)
                                                                                       gCCATCGCCACAGACCACCACTTCGCTGTCGTCTATTCCCCATAATTTCTGCAGTTGGCGAAGGCCATTGGCTTTATGTACGCCGGGGATAATCAGATCg         <  2:445006/100‑1 (MQ=255)
                                                                                        ccATCGCCACAGACCACCACTTCGCTGTCGTCTATTCCCCATAATTTCTGCAGTTGGCGAAGGCCATTGGCTTTATGTACGCCGGGGATAATCAGATCGa        <  1:294922/100‑1 (MQ=255)
                                                                                          aTCGCCACAGACCACCACTTCGCTGTCGTCTATTCCCCATAATTTCTGCAGTTGGCGAAGGCCATTGGCTTTATGTACGCCGGGGATAATCAGATCGAtg      >  1:3443283/1‑100 (MQ=255)
                                                                                          aTCGCCACAGACCACCACTTCGCTGTCGTCTATTCCCCATAATTTCTGCAGTTGGCGAAGGCCATTGGCTTTATGTACGCCGGGGATAATCAGATCGAtg      >  1:994358/1‑100 (MQ=255)
                                                                                          aTCGCCACAGACCACCACTCCGCTGTCGTCTATTCCCCATAATTTCTGCAGTTGGCGAAGGCCATTGGCTTTATGTACGCCGGGGATAATCAGATCGAtg      >  1:1143902/1‑100 (MQ=255)
                                                                                             gCCACAGACCACCACTTCGCTGTCGTCTATTCCCCATAATTTCTGCAGTTGGCGAAGGCCATTGGCTTTATGTACGCCGGGGATAATCAGATCGAtgctg   >  1:314539/1‑100 (MQ=255)
                                                                                             gCCACAGACCACCACTTCGCTGTCGTCTATTCCCCATAATTTCTGCAGCTGGCGACGGCCACTGCACTTAAGTACTCCTGGGAGACTCAGACCGAtggtg   >  2:3369073/1‑100 (MQ=255)
                                                                                              ccACAGACCACCACTTCGCTGTCGTCTATTCCCCATAATTTCTGCAGTTGGCGAAGGCCATTTGCTTTATGTACGCCGGGGATAATCAGATCGAtgctgc  >  2:3635302/1‑100 (MQ=255)
                                                                                                 |                                                                                                
GTATTTTGCCGCTGCGACGACCGCGCTACCGGCATTTTCCATTGCAAAACTAAAGCCTGCCTGACGCAGCATCTCAATATCGTTACCGCCATCGCCAAAGACCACCACTTCGCTGTCGTCTATTCCCCATAATTTCTGCAGTTGGCGAAGGCCATTGGCTTTATGTACGCCGGGGATAATCAGATCGATGCTGC  >  NC_000913/885397‑885590

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