Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,946,711 C→A 15.4% noncoding (12/120 nt) rrfC → 5S ribosomal RNA of rrnC operon

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,946,7110CA15.4% 474.0 / 75.0 193noncoding (12/120 nt)rrfC5S ribosomal RNA of rrnC operon
Reads supporting (aligned to +/- strand):  ref base C (107/56);  new base A (26/4);  total (133/60)
Fisher's exact test for biased strand distribution p-value = 3.02e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.57e-01

Only 100 of 194 total aligned reads displayed.

CGCCGAAGATGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATGAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCT  >  NC_000913/3946612‑3946788
                                                                                                   |                                                                             
  ccGAAGATGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATGAAACAGAATTTGCCTGGCGGCCGt                                                                             >  2:2801149/1‑100 (MQ=38)
    gAAGATGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATGAAACAGAATTTGCCTGGCGGCCGTAg                                                                           <  1:3185024/100‑1 (MQ=38)
      aGATGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATGAAACAGAATTTGCCTGGCGGCCGTAgcg                                                                         >  2:1515483/1‑100 (MQ=38)
       gATGTTTTGACGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATGAAACAGAATTTGCCTGGCGGCCGTAgcgc                                                                        >  2:1568934/1‑100 (MQ=25)
        aTGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATGAAACAGAATTTGCCTGGCGGCCGTAGcgcg                                                                       >  1:1024429/1‑100 (MQ=35)
        aTGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATGAAACAGAATTTGCCTGGCGGCCGTAGcgcg                                                                       >  2:1638221/1‑100 (MQ=35)
         tGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATGAAACAGAATTTGCCTGGCGGCCGTAGCGCgg                                                                      >  1:1592887/1‑100 (MQ=32)
         tGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATGAAACAGAATTTGCCTAGCGGCCGTAGCGCgg                                                                      <  1:521869/100‑1 (MQ=18)
          gTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATGAAACAGAATTTGCCTGGCGGCCGTAGCGCggt                                                                     <  1:2819235/100‑1 (MQ=33)
             ttGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATGAAACAGAATTTGCCTGGCGGCCGTAGCGCggtggt                                                                  <  1:1111408/100‑1 (MQ=33)
               ggCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATGAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTcc                                                                >  1:802344/1‑100 (MQ=33)
                 nggATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATGAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCa                                                              >  1:2054033/2‑100 (MQ=18)
                 cGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATGAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCa                                                              >  1:2752617/1‑100 (MQ=33)
                    aTGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATGAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCt                                                           <  1:633429/100‑1 (MQ=33)
                        gagaAGATTTTCAGCCTGATACAGATTAAATCAGAACTCAGAAGCGGTCTGATGAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGAcc                                                       <  2:2448428/100‑1 (MQ=18)
                        gagaAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATGAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGAcc                                                       >  2:1558515/1‑100 (MQ=33)
                        gagaAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATGAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGAcc                                                       <  1:1394351/100‑1 (MQ=33)
                        gagaAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATGAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGAcc                                                       >  1:337366/1‑100 (MQ=33)
                        gagaAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATGAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGAcc                                                       <  1:110538/100‑1 (MQ=33)
                        gagaAGATTGTCAGCCTGATACAGATGAAATCAGAACGCAGAAGCGGTCTGATGAAACAGAATTTTCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGAc                                                        <  2:2390399/99‑1 (MQ=255)
                           aagaTTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATGAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCa                                                    <  1:1134419/100‑1 (MQ=33)
                                tttCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATGAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGCCg                                               <  1:438410/100‑1 (MQ=33)
                                 ttCAGCCTGATACAGATTAAATCAGAACGCAGCAGCGGTCTGATGAAACAGAATTTGCTTGACGGCCGTAGCGCGGTGCTCCCACCGGCCCCCATgacta                                              >  1:3038110/1‑96 (MQ=11)
                                 ttCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATGAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGCCGa                                              >  2:2695288/1‑100 (MQ=33)
                                     gCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATGAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTc                                          >  2:603406/1‑100 (MQ=33)
                                      gatctattctgttttatcagaccgcttctgcgttctgatAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCa                                         <  2:974691‑M2/61‑1 (MQ=255)
                                      ccTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATGAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCa                                         <  1:612227/100‑1 (MQ=33)
                                      ccTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATGAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCa                                         >  2:176924/1‑100 (MQ=33)
                                        tGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATGAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGa                                       >  1:2509365/1‑100 (MQ=33)
                                          aTACAGATTAAATCAGAACGCAGAAGCGGTCTGATGAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAg                                     >  1:2284131/1‑100 (MQ=33)
                                          aTACAGATTAAATCAGAACGCAGAAGCGGTCTGATGAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAg                                     <  2:2997575/100‑1 (MQ=33)
                                           attctgttttatcagaccgcttctgcgttctgatAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGaactcaagatc                                    >  1:974698‑M2/35‑95 (MQ=255)
                                           attctgtttcatcagaccgcttctgcgttctgatGAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGt                                    >  2:1997369‑M2/35‑100 (MQ=255)
                                            aCAGATTAAATCAGAACGCAGAAGCGGTCTGATGAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTg                                   >  1:2777022/1‑100 (MQ=33)
                                             cAGATTAAATCAGAACGCAGAAGCGGTCTGATGAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGa                                  <  1:1817397/100‑1 (MQ=33)
                                              ctgttttatcagaccgcttctgcgttctgatAAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGaa                                 >  1:1493263‑M2/32‑100 (MQ=255)
                                              ctgttttatcagaccgcttctgcgttctgatAAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGaa                                 >  1:693096‑M2/32‑100 (MQ=255)
                                              ctgttttatcagaccgcttctgcgttctgatAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGaa                                 <  2:1493255‑M2/69‑1 (MQ=255)
                                              ctgttttatcagaccgcttctgcgttctgatAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGaa                                 >  2:1839242‑M2/32‑100 (MQ=255)
                                              ctgttttatcagaccgcttctgcgttctgatAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGaa                                 >  2:1754273‑M2/32‑100 (MQ=255)
                                              ctgttttatcagaccgcttctgcgttctgatAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGaa                                 >  1:1286954‑M2/32‑100 (MQ=255)
                                              aGATTAAATCAGAACGCAGAAGCGGTCTGATGAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGaa                                 <  1:555928/100‑1 (MQ=33)
                                               tgttttatcagaccgcttctgcgttctgatAAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGaaa                                >  1:1967175‑M2/31‑100 (MQ=255)
                                               tgttttatcagaccgcttctgcgttctgatAAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGaaa                                >  1:3205532‑M2/31‑100 (MQ=255)
                                               tgttttatcagaccgcttctgcgttctgatAAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGaaa                                >  1:1351857‑M2/31‑100 (MQ=255)
                                               tgttttatcagaccgcttctgcgttctgatAAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGaaa                                >  2:3388249‑M2/31‑100 (MQ=255)
                                               tgttttatcagaccgcttctgcgttctgatAAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGa                                  >  1:2625920‑M2/31‑98 (MQ=255)
                                               gATTAAATCAGAACGCAGAAGCGGTCTGATGAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGaca                                <  1:2461198/100‑3 (MQ=18)
                                               cttccgATCTGAACGCAGAAGCGGTCTGATGAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGaaa                                >  1:1365908‑M1/7‑100 (MQ=255)
                                                nttttatcagaccgcttctgcgttctgatAAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAAc                               >  1:1408253‑M2/30‑100 (MQ=255)
                                                gttttattagaccgcttctgcgttctgatAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAAc                               >  1:1918414‑M2/30‑100 (MQ=255)
                                                gttttatcagaccgcttctgcgttctgatAAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAAc                               >  2:2605364‑M2/30‑100 (MQ=255)
                                                gttttatcagaccgcttctgcgttctgatAAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAAc                               >  1:2964835‑M2/30‑100 (MQ=255)
                                                gttttatcagaccgcttctgcgttctgatAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGaaa                                >  2:1387407‑M2/30‑99 (MQ=255)
                                                gttttatcagaccgcttctgcgttctgatAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAAc                               >  2:1930940‑M2/30‑100 (MQ=255)
                                                gttttatcagaccgcttctgcgttctgatAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAAc                               >  1:1178110‑M2/30‑100 (MQ=255)
                                                gttttatcagaccgcttctgcgttctgatAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAAc                               >  1:341285‑M2/30‑100 (MQ=255)
                                                gttttatcagaccgcttctgcgttctgatAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAAc                               >  1:2172575‑M2/30‑100 (MQ=255)
                                                gttttatcagaccgcttctgcgttctgatAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAAc                               >  1:821533‑M2/30‑100 (MQ=255)
                                                gttttatcagaccgcttctgcgttctgatAAAAAAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACATGACCCCATGCCGAACTCAGAAGTGAAAc                               >  2:2188718‑M2/30‑100 (MQ=255)
                                                 ttttatcagaccgcttctgcgttctgatAAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACg                              >  1:2140425‑M2/29‑100 (MQ=255)
                                                 ttttatcagaccgcttctgcgttctgatAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACg                              >  2:189370‑M2/29‑100 (MQ=255)
                                                  tttatcagaccgcttctgcgttctgatAAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGc                             >  2:3055590‑M2/28‑100 (MQ=255)
                                                  tttatcagaccgcttctgcgttctgatAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGc                             >  2:1054023‑M2/28‑100 (MQ=255)
                                                  tttatcagaccgcttctgcgttctgatAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGc                             >  1:2375129‑M2/28‑100 (MQ=255)
                                                   ttatcagaccgcttctgcgttctgatAAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGcc                            >  1:846529‑M2/27‑100 (MQ=255)
                                                   tcatcagaccgcttctgcgttctgatGAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGcc                            >  1:3279899‑M2/27‑100 (MQ=255)
                                                   aaaTCAGAACGCAGAAGCGGTCTGATGAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGcc                            >  2:2157545/1‑100 (MQ=32)
                                                   aaaTCAGAACGCAGAAGCGGTCTGATGAAACAGAATTTGCCTGCCGGCCGTAGCGCGGTGGTCCCACCTGACCCCACGCCGAACTCAGAAGTGAAACGcc                            >  1:1579247/1‑100 (MQ=12)
                                                    aaTCAGAACGCAGAAGCGGTCTGATGAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCg                           <  2:2726230/100‑1 (MQ=32)
                                                     atcagaccgcttctgcgttctgatGAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGt                          >  2:3237733‑M2/25‑100 (MQ=255)
                                                     atcagaccgcttctgcgttctgatGAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGATCCCACCTGACCCCATGCCGCACTCAGAAGTGAAACGc                             >  2:3026171‑M2/25‑97 (MQ=255)
                                                      tcagaccgcttctgcgttctgatAAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTa                         >  2:2658651‑M2/24‑100 (MQ=255)
                                                      tcagaccgcttctgcgttctgatAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTa                         >  2:1244153‑M2/24‑100 (MQ=255)
                                                      TCAGAACGCAGAAGCGGTCTGATGAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTa                         >  1:1769157‑M1/1‑100 (MQ=255)
                                                       CAGAACGCAGAAGCGGTCTGATGAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAg                        >  2:1203702‑M1/1‑100 (MQ=255)
                                                        AGAACGCAGAAGCGGTCTGATGAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAgc                       <  2:1284862‑M1/100‑1 (MQ=255)
                                                        AGAACGCAGAAGCGGTCTGATGAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAgc                       <  1:2671721‑M1/100‑1 (MQ=255)
                                                         GAACGCAGAAGCGGTCTGATGAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAActcagaagtgaaaagatcggaag                      <  2:1365323‑M1/100‑11 (MQ=255)
                                                            cGCAGAAGCGGTCTGATGAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCACACCTGACCCCGTGCCGAACTCAGAAGTGAAACGCCGTAGCGCCg                   >  2:3411272/1‑100 (MQ=16)
                                                             gcttctgcgttctgatAAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGa                  >  1:3009177‑M2/17‑100 (MQ=255)
                                                              cAGAAGCGGTCTGATGAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGAt                 >  1:3334414/1‑100 (MQ=33)
                                                                gAAGCGGTCTGATGAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATgg               >  2:473914/1‑100 (MQ=33)
                                                                gAAGCGGTCTGATGAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATgg               >  2:2888956/1‑100 (MQ=32)
                                                                gAAGCGGTCTGATGAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATgg               <  2:3314780/100‑1 (MQ=33)
                                                                 aaGCGGTCTGATGAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCGGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGt              <  1:3487355/100‑1 (MQ=17)
                                                                  aGCGGTCTGATGAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTa             <  2:2024627/100‑1 (MQ=32)
                                                                     ggTCTGATGAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAgtg          >  1:289786/1‑100 (MQ=33)
                                                                     ggTCTGATGAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAgtg          <  1:2333331/100‑1 (MQ=33)
                                                                      gTCTGATGAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAgtgt         >  2:190867/1‑100 (MQ=33)
                                                                        ctgatGAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGTCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTgg       >  2:134213‑M2/6‑100 (MQ=255)
                                                                        ctgatGAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTgg       <  2:3073160‑M2/95‑1 (MQ=255)
                                                                         tgatGAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGGAGTGTgg       >  2:2779066‑M2/5‑99 (MQ=255)
                                                                         tgatGAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGGACCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGtgtg        >  2:3458528‑M2/5‑98 (MQ=255)
                                                                          gatGAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTATTGTgggg     <  1:351095‑M2/97‑1 (MQ=255)
                                                                          gatGAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTgggg     >  2:3126336‑M2/4‑100 (MQ=255)
                                                                          gatGAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTgggg     <  2:1116330‑M2/97‑1 (MQ=255)
                                                                          gatGAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTgggg     >  2:1807333‑M2/4‑100 (MQ=255)
                                                                          gatGAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGCCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTgggg     >  1:2331369‑M2/4‑100 (MQ=255)
                                                                            tGAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTc   >  2:426771‑M2/2‑100 (MQ=255)
                                                                                                   |                                                                             
CGCCGAAGATGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATGAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCT  >  NC_000913/3946612‑3946788

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