Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 922,316:1 +C 18.6% intergenic (+23/+50) macB → / ← cspD macrolide ABC transporter peremase/ATPase/inhibitor of DNA replication, cold shock protein homolog

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913922,3141.C18.6% 499.4 / 133.6 161intergenic (+21/+52)macB/cspDmacrolide ABC transporter peremase/ATPase/inhibitor of DNA replication, cold shock protein homolog
Reads supporting (aligned to +/- strand):  ref base . (65/66);  new base C (13/17);  total (78/83)
Fisher's exact test for biased strand distribution p-value = 5.51e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.70e-01

Only 100 of 170 total aligned reads displayed.

CACCGGGATTTTATTTGGCTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTG  >  NC_000913/922215‑922413
                                                                                                    |                                                                                                   
cACCGGGATTTTATTTGGCTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAg                                                                                                      >  1:1398193/1‑100 (MQ=255)
 aCCGGGATTTTATTTTGCTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑c                                                                                                    <  2:563457/100‑1 (MQ=255)
 aCCGGGATTTTATTTGGCTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑c                                                                                                    <  1:902143/100‑1 (MQ=255)
  ccGGGATTTTATTTGGCTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑cc                                                                                                   >  1:1281874/1‑100 (MQ=255)
   cGGGATTTTATTTGGCTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCg                                                                                                  <  2:517381/100‑1 (MQ=255)
     ggATTTTATTTGGCTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGAt                                                                                                <  2:501836/100‑1 (MQ=255)
       aTTTTATTTGGCTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGATATTTATATGCTAG‑CCGATc                                                                                               >  1:1808566/1‑99 (MQ=255)
       aTTTTATTTGGCTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCg                                                                                              <  2:348234/100‑1 (MQ=255)
        ttttATTTGGCTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCgg                                                                                             >  1:4211448/1‑100 (MQ=255)
          ttATTTGGCTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGc                                                                                           >  1:4382071/1‑100 (MQ=255)
          ttATTTGGCTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGc                                                                                           >  1:2171076/1‑100 (MQ=255)
            aTTTGGCTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTg                                                                                         <  2:2790550/100‑1 (MQ=255)
               tGGCTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAGCCCGATC‑GGCTGGCa                                                                                      <  2:4216404/100‑1 (MQ=255)
               tGGCTGGTTACCCGCACGAAATGCGGCACGACTGGAGCCAGTAGATGCTCTGGCACGAGAGTAGTTTTTGAGATAAAAATGCCAGCCCGATC‑GGCTGGCa                                                                                      <  2:3100680/100‑1 (MQ=255)
                ggCTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGTTAAAAATGCCAGCCCGATC‑GGCTGGCAt                                                                                     <  1:1150168/100‑1 (MQ=39)
                 gCTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAGCCCGATC‑GGCTGGCAtt                                                                                    <  1:2264898/100‑1 (MQ=39)
                  cTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAGCCCGAGC‑GGCTGGCAttt                                                                                   <  2:2897545/100‑1 (MQ=255)
                  cTGGGTACCCGCACGAAATGCGGCAAGACTGGATCCAGAAAATACTCGGGCACGAGAGTAATTTTTGAGATAAAAATGCCAGCCCGATC‑GGCTGGCa                                                                                      <  1:2179409/97‑1 (MQ=25)
                    ggTTCCCCGCACGAAATGCGGCACGACTGGATGCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAACATGCCAG‑CCGATCTGGCTGGCATTTTg                                                                                 <  1:1597505/100‑1 (MQ=255)
                    ggTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTg                                                                                 <  2:1552828/100‑1 (MQ=255)
                      ttACCCGCACGAAATGCGGCACGACTGTATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGcc                                                                               <  1:2524636/100‑1 (MQ=255)
                         nccGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCttt                                                                            >  2:4285028/2‑100 (MQ=255)
                         ccgtcACGAAATGCGGCACGTCTGGATCCAGTGGATGCTCTGGCACGATAGTAATTTTTGAGATAAAAATGCCAGCCCGATC‑GGCTGGCAtttttatct                                                                             <  2:2883512/95‑6 (MQ=21)
                         cccGCACGAAGTGCGGCAGGACTGGATCCAGTAGATGCGCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAGCCCGATC‑GGCTGGCAtttttatc                                                                              <  1:1915800/98‑5 (MQ=25)
                          ccGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTa                                                                           >  2:3362822/1‑100 (MQ=255)
                           cGCAGGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGTTAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAg                                                                          <  2:1446349/100‑1 (MQ=255)
                           cGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTa                                                                           >  2:229212/1‑99 (MQ=255)
                           cGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAg                                                                          >  2:2452398/1‑100 (MQ=255)
                           cGCACGAAATGCGGCACGACTGGATCCAGTAGATGCCCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAg                                                                          >  2:1153433/1‑100 (MQ=255)
                            gCACGAAGTGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGcc                                                                               <  1:1271438/94‑1 (MQ=255)
                            gCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAgg                                                                         <  2:3659787/100‑1 (MQ=255)
                            gCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAgg                                                                         <  2:2005234/100‑1 (MQ=255)
                             cACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGa                                                                        <  1:1201708/100‑1 (MQ=255)
                              aCGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGAt                                                                       >  1:87082/1‑100 (MQ=255)
                               cGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGTGATAAAATTGCCAG‑CCGATCGGGCTGCCATTTTGCCTCTAggc                                                                        >  1:3058293/1‑97 (MQ=255)
                                 aaaTGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTa                                                                    <  2:3042842/100‑1 (MQ=255)
                                     gCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTAcaca                                                                >  2:1538166/1‑100 (MQ=255)
                                       ggCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCGAGCTGGCATTTTGCCTTTAGGATGTACACAAt                                                              <  2:2977375/100‑1 (MQ=255)
                                         cACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATga                                                            <  2:2280356/100‑1 (MQ=255)
                                          aCGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATgag                                                           <  2:3752885/100‑1 (MQ=255)
                                            gACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGAc                                                         <  1:766113/100‑1 (MQ=255)
                                            gACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATGTTTGAGATGAAAATGCCAG‑CCGTTCAGGCTGGCCTTTTGCATTTAGGAGGTACACTATGAGAc                                                         >  1:3937492/1‑100 (MQ=255)
                                             aCTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACa                                                        >  1:3005069/1‑100 (MQ=255)
                                              cTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACag                                                       <  2:1634035/100‑1 (MQ=255)
                                               tGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACaga                                                      >  2:688641/1‑100 (MQ=255)
                                               tGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACaga                                                      >  2:692646/1‑100 (MQ=255)
                                                ggATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACagaa                                                     >  2:1761462/1‑100 (MQ=255)
                                                ggATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACagaa                                                     >  2:1211799/1‑100 (MQ=255)
                                                ggATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACagaa                                                     >  1:157348/1‑100 (MQ=255)
                                                 gATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACagaag                                                    >  2:373362/1‑100 (MQ=255)
                                                  aTCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACagaaga                                                   <  2:3016707/100‑1 (MQ=255)
                                                  aTCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACagaaga                                                   <  2:4132095/100‑1 (MQ=255)
                                                   tCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAagag                                                  >  2:3225320/1‑100 (MQ=255)
                                                    ccAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGACGTACACAATGAGACAGAAGAGc                                                 <  2:3310692/100‑1 (MQ=255)
                                                    ccAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CAGATCGGGCTGGTTTTTTGCCTTTAGGATGTACACAATGAGACAGAagag                                                  >  1:341608/1‑99 (MQ=255)
                                                     cAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCt                                                <  1:4301363/100‑1 (MQ=255)
                                                     cAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCATTTAGGATGTACACAATGAGACAGAAGAGCt                                                >  1:1079276/1‑100 (MQ=255)
                                                       gTAGATGCGCTGGCACGAGAGTAATTTTTGAGACAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTa                                               <  1:4162958/99‑1 (MQ=255)
                                                           aTGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGa                                          <  2:43106/100‑1 (MQ=255)
                                                              ctctGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACt                                        >  1:3723082/1‑99 (MQ=255)
                                                                ctGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGcc                                     <  1:2062250/100‑1 (MQ=255)
                                                                 tGGCACGAGAGTAATTTTTGGTAGAAGAATGCCAG‑CCGCTAGGGCTGGCATTTGGCCTTTAGGATGTACACAATGCGACAGAAGAGCTATGGGACTg                                       <  1:1809079/97‑1 (MQ=255)
                                                                  ggCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGAGTGCCGc                                   >  2:3593935/1‑100 (MQ=255)
                                                                   gCACGAGATTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCt                                  <  1:229211/100‑1 (MQ=255)
                                                                   gCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCt                                  <  1:2743135/100‑1 (MQ=255)
                                                                    cACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCtt                                 <  2:4039971/100‑1 (MQ=255)
                                                                    cACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCt                                  <  2:2219559/99‑1 (MQ=255)
                                                                      cGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCt                               >  1:1011376/1‑100 (MQ=255)
                                                                      cGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCt                               >  1:985002/1‑100 (MQ=255)
                                                                        agagTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTAc                             <  2:4341727/100‑1 (MQ=255)
                                                                        agagTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTAc                             >  1:2466391/1‑100 (MQ=255)
                                                                         gagTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACt                            <  2:1011380/100‑1 (MQ=255)
                                                                          agTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACtt                           >  1:2191085/1‑100 (MQ=255)
                                                                           gTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTc                          >  1:1643220/1‑100 (MQ=255)
                                                                           gTAATTTTTGAGATAAAAATGCCAG‑CCGACCGGGCTGTCATTTTGCCTTAAGGACGTCCACACTGAGCCCGAAGTGCTATGCGACGGCAGCTTCTACTTc                          >  1:539548/1‑100 (MQ=255)
                                                                               tttttGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACg                      >  2:2114963/1‑100 (MQ=255)
                                                                                ttttGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACg                      >  1:2023022/1‑99 (MQ=255)
                                                                                 tttGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACggg                    <  2:133667/100‑1 (MQ=255)
                                                                                 tttGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACggg                    <  1:2843502/100‑1 (MQ=255)
                                                                                   tGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGca                  <  2:1355274/100‑1 (MQ=255)
                                                                                      gaTAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACaa               <  1:3373494/100‑1 (MQ=255)
                                                                                       aTAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACAAt              <  1:1992727/100‑1 (MQ=255)
                                                                                        tAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACAATa             >  2:2638808/1‑100 (MQ=255)
                                                                                        tAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACAATa             >  2:1972565/1‑100 (MQ=255)
                                                                                        tAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACAATa             <  2:2105733/100‑1 (MQ=255)
                                                                                         caaaaTGCCAGCCCGATC‑GGCTGGCATTTTGCCTTTAGGATGTACACAATGGGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGAGGGGCACAATaa            >  1:3422452/2‑100 (MQ=38)
                                                                                         aaaaaTGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACAAt              <  2:1690618/98‑1 (MQ=255)
                                                                                          aaaaTGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAac           >  2:476375/1‑100 (MQ=255)
                                                                                            aaTGCCAGCCCGATC‑GGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAacac         >  1:1483242/1‑100 (MQ=39)
                                                                                            aaTGCCAGCCCGATC‑GGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAacac         >  2:773446/1‑100 (MQ=39)
                                                                                              tGCCAGCCCGATC‑GGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTg       >  1:4129032/1‑100 (MQ=255)
                                                                                              tGCCAGCCCGATC‑GGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTg       >  2:3040054/1‑100 (MQ=255)
                                                                                              tGCCAGCCCGATC‑GGCTGGCATTTTGCATTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTg       >  1:1887758/1‑100 (MQ=255)
                                                                                              tGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGTGACTGCCGCTTCTACTTCGACGGGCACAATAACACTg       <  1:2342471/100‑1 (MQ=255)
                                                                                               gCCAGCCCGATC‑GGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTgg      >  2:865158/1‑100 (MQ=255)
                                                                                               gCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTgg      >  2:4183319/1‑100 (MQ=255)
                                                                                                ccAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGc     >  2:3568196/1‑100 (MQ=255)
                                                                                                 cAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCg    <  2:2568615/100‑1 (MQ=255)
                                                                                                  aG‑CCGATCGGGCTGGCATCTTGCCTTTAGGATGTACACAATGAGACAGAAGGGCTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCg    <  1:1020832/99‑1 (MQ=255)
                                                                                                   g‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTg  <  1:1134962/100‑1 (MQ=255)
                                                                                                    |                                                                                                   
CACCGGGATTTTATTTGGCTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTG  >  NC_000913/922215‑922413

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