Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 934,731 A→G Q762R (CAG→CGG)  ftsK → cell division protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913934,7310AG84.1% 375.7 / 62.5 182Q762R (CAG→CGG) ftsKcell division protein
Reads supporting (aligned to +/- strand):  ref base A (18/11);  new base G (89/64);  total (107/75)
Fisher's exact test for biased strand distribution p-value = 8.37e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.76e-01

Only 100 of 205 total aligned reads displayed.

TTTACGCCAATTGTTGAACCTGTACAGCAGCCGCAACAACCGGTTGCACCGCAGCAGCAATATCAGCAGCCGCAACAACCAGTTCCGCCGCAGCCGCAGTATCAGCAGCCACAACAGCCGGTTGCGCCGCAGCCACAATATCAGCAGCCGCAACAACCGGTTGCGCCACAGCAGCAATATCAGC  >  NC_000913/934640‑934823
                                                                                            |                                                                                            
 ttACGCCAATTGTTGAACCTGTACAGCAGCCGCAACAACCGGTTGCACCGCAGCAGCAATATCAGCAGCCGCAACAACCAGTTCCgccgcggc                                                                                            <  2:6410918/93‑1 (MQ=255)
 ttACGCCAATTGTTGAACCTGTACAGCAGCCGCAACAACCGGTTGCACCGCAGCAGCAATATCAGCAGCCGCAACAACCAGTTCCgccgcggc                                                                                            >  2:2416645/1‑93 (MQ=255)
    cGCCAATTGTTGAACCTGTACAGCAGCCGCAACAACCGGTTGCACCGCAGCAGCAATATCAGCAGCCGCAACAACCAGTTCCgccgcggccgc                                                                                         >  1:771750/1‑93 (MQ=255)
     gCCAATTGTTGAACCTGTACAGCAGCCGCAACAACCGGTTGCACCGCAGCAGCAATATCAGCAGCCGCAACAACCAGTTCCgccgcggccgca                                                                                        <  2:5395802/93‑1 (MQ=255)
     ccaaaTTTTTAAACCTGTACAGCAGCCGCAACAACCGGTTGCACCGCAGCAGCAATATCAGCAGCCGCAACAACCAGTTCCgccgcggccgca                                                                                        <  1:4795335/90‑1 (MQ=18)
      ccAATTGTTGAACCTGTACAGCAGCCGCAACAACCGGTTGCACCGCAGCAGCAATATCAGCAGCCGCAACAACCAGTTCCGccgcggccgcag                                                                                       >  1:6091647/1‑93 (MQ=255)
       cAATTGTTGAACCTGTACAGCAGCCGCAACAACCGGTTGCACCGCAGCAGCAATATCAGCAGCCGCAACAACCAGTTCCGCCGCGGCCGCAGt                                                                                      >  2:2786063/1‑93 (MQ=255)
        aaTTGTTGAACCTGTACAGCAGCCGCAACAACCGGTTGCACCGCAGCAGCAATATCAGCAGCCGCAACAACCAGTTCCGCCGCGGCCGCAGTa                                                                                     >  2:7156529/1‑93 (MQ=255)
        aaTTGTTGAACCTGTACAGCAGCCGCAACAACCGGTTGCACCGCAGCAGCAATATCAGCAGCCGCAACAACCAGTTCCGCCGCGGCCGCAGTa                                                                                     <  2:4599924/93‑1 (MQ=255)
        aaTTGTTGAACCTGTACAGCAGCCGCAACAACCGGTTGCACCGCAGCAGCAATATCAGCAGCCGCAACAACCAGTTCCGCCGCAGCCGCAGTa                                                                                     >  1:402982/1‑93 (MQ=255)
        aaTTGTTGAACCTGTACAGCAGCCGCAACAACCGGTTGCACCGCAGCAGCAATATCAGCAGCCGCAACAACCAGTTCCGCCGCAGCCGCAGTa                                                                                     <  2:4439069/93‑1 (MQ=255)
        aaTTGTTAAACCTGTACAGCAGCCGCAACAACCGGTTGCACCGCAGCAGCAATATCAGCAGCCGCAACAACCAGTTCCGCCGCGGCCGCAGTa                                                                                     <  1:2571064/93‑1 (MQ=255)
         aTTGTTGAACCTGTACAGCAGCCGCAACAACCGGTTGCACCGCAGCAGCAATATCAGCAGCCGCAACAACCAGTTCCGCCGCGGCCGCAGTAt                                                                                    <  2:219814/93‑1 (MQ=255)
          ttgttgAACCTGTACAGCAGCCGCAACAACCGGTTGCACCGCAGCAGCAATATCAGCAGCCGCAACAACCAGTTGCGCCGCAACCGCAGTATc                                                                                   >  2:3782503/1‑93 (MQ=255)
          ttgttgAACCTGTACAGCAGCCGCAACAACCGGTTGCACCGCAGCAGCAATATCAGCAGCCGCAACAACCAGTTCCGCCGCGGCCGCAGTATc                                                                                   >  1:196171/1‑93 (MQ=255)
          ttgttgAACCTGTACAGCAGCCGCAACAACCGGTTGCACCGCAGCAGCAATATCAGCAGCCGCAACAACCAGTTCCGCCGCAGCCGCAGTATc                                                                                   <  2:4682206/93‑1 (MQ=255)
          ttgttgAACCTGTACAGCAGCCGCAACAACCGGTTGCACCGCAGCAGCAATATCAGCAGCCGCAACAACCAGTTCCGCCGCAGCCGCAGTATc                                                                                   >  1:6467124/1‑93 (MQ=255)
          ttgttaAACCTGAACAGCAGCCGCAACAACCGGTTGCACCGCAGCAGCAATATCAGCAGCCGCAACAACCAGTTCCGCCGCGGCCGCAGTATc                                                                                   <  1:1292440/93‑1 (MQ=18)
           tgttgACCCTGTACAGCAGCCGCAACAACCGGTTGCACCGCAGCAGCAATATCAGCAGCCGCAACAACCAGTTCCGCCGCGGCCGCAGTATca                                                                                  <  1:6831826/93‑1 (MQ=255)
            gttgAACCTGTACAGCAGCCGCAACAACCGGTTGCACCGCAGCAGCAATATCAGCAGCCGCAACAACCAGTTCCGCCGCGGCCGCAGTATcag                                                                                 >  1:4351674/1‑93 (MQ=255)
            gttgAACCTGTACAGCAGCCGCAACAACCGGTTGCACCGCAGCAGCAATATCAGCAGCCGCAACAACCAGTTCCGCCGCGGCCGCAGTATcag                                                                                 <  2:5017485/93‑1 (MQ=255)
             ttgAACCTGTACAGCAGCCGCAACAACCGGTTGCACCGCAGCAGCAATATCAGCAGCCGCAACAACCAGTTCCGCCGCGGCCGCAGTATcagc                                                                                >  2:1138079/1‑93 (MQ=255)
               gAACCTGTACAGCAGCCGCAACAACCGGTTGCACCGCAGCAGCAATATCAGCAGCCGCAACAACCAGTTCCGCCGCGGCCGCAGTATcagcag                                                                              >  1:265112/1‑93 (MQ=255)
                aaCCTGTACAGCAGCCGCAACAACCGGTTGCACCGCAGCAGCAATATCAGCAGCCGCAACAACCAGTTCCGCCGCGGCCGCAGTATCagcagc                                                                             >  2:5530001/1‑93 (MQ=255)
                 aCCTGTACAGCAGCCGCAACAACCGGTTGCACCGCAGCAGCAATATCAGCAGCCGCAACAACCAGTTCCGCCGCGGCCGCAGTATCAGCAGcc                                                                            >  2:4049256/1‑93 (MQ=255)
                  ccTGTACAGCAGCCGCAACAACCGGTTGCACCGCAGCAGCAATATCAGCAGCCGCAACAACCAGTTCCGCCGCGGCCGCAGTATCAGCAGCca                                                                           <  2:1386366/93‑1 (MQ=255)
                      tACAGCAGCCGCAACAACCGGTTGCACCGCAGCAGCAATATCAGCAGCCGCAACAACCAGTTCCGCCGCGGCCGCAGTATCAGCAGCCacaac                                                                       <  2:3034403/93‑1 (MQ=255)
                        cagcagCCGCAACAACCGGTTGCACCGCAGCAGCAATATCAGCAGCCGCAACAACCAGTTCCGCCGCGGCCGCAGTATCAGCAGCCACAACAg                                                                     <  2:1433317/93‑1 (MQ=255)
                         agcagcCGCAACAACCGGTTGCACCGCAGCAGCAATATCAGCAGCCGCAACAACCAGTTCCGCCGCGGCCGCAGTATCAGCAGCCACAACAGc                                                                    >  1:355241/1‑93 (MQ=255)
                         agcagcCGCAACAACCGGTTGCACCGCAGCAGCAATATCAGCAGCCGCAACAACCAGTTCCGCCGCGGCCGCAGTATCAGCAGCCACAACAGc                                                                    >  1:1670189/1‑93 (MQ=255)
                          gcagcCGCAACAACCGGTTGCACCGCAGCAGCAATATCAGCAGCCGCAACAACCAGTTCCGCCGCGGCCGCAGTATCAGCAGCCACAACAGcc                                                                   >  2:1961442/1‑93 (MQ=255)
                             cccGCAACACCGGGTTGCCCCGCAGCACAAATACCAGCACCCGCAACAACCAGTCCCGCCGCGGCCGCAGTATCAGCAGCCACAACAGCCGGt                                                                <  1:7026933/92‑1 (MQ=17)
                              ccGCAACACCGGGTTGCACCGCAGCAGAAATACCAGCAGCCGCAACAACCAGTCCCGCCGCGGCCGCAGTATCAGCAGCCACAACAGCCGGtt                                                               <  1:3841401/93‑1 (MQ=17)
                               caaaactacCGGTTGCACCGCAGCAGCAATATCGGCAGCCGCAACAACCAGTTCCGCCGCAGCCGCAGTATCAGCAGCCACAACAGCCGGTTg                                                              <  2:2052853/90‑1 (MQ=255)
                               cGCAACAACCGGTTGCACCGCAGCAGCAATATCAGCAGCCGCAACAACCAGTTCCGCCGCAGCCGCAGTATCAGCAGCCACAACAGCCGGTTg                                                              <  2:5069019/93‑1 (MQ=255)
                                gCAACAACCGGTTGCACCGCAGCAGCAATATCAGCAGCCGCAACAACCAGTTCCGCCGCGGCCGCAGTATCAGCAGCCACAACAGCCGGTTgc                                                             >  2:2885060/1‑93 (MQ=255)
                                gCAACAACCGGTTGCACCGCAGCAGCAATATCAGCAGCCGCAACAACCAGTTCCGCCGCGGCCGCAGTATCAGCAGCCACAACAGCCGGTTgc                                                             <  2:4071747/93‑1 (MQ=255)
                                 caacaaCCGGTTGCACCGCAGCAGCAATATCAGCAGCCGCAACAACCAGTTCCGCCGCAGCCGCAGTATCAGCAGCCACAACAGCCGGTTgcg                                                            >  2:4427555/1‑93 (MQ=255)
                                  aacaacCGGTTGCACCGCAGCAGCAATATCAGCAGCCGCAACAACCAGTTCCGCCGCGGCCGCAGTATCAGCAGCCACAACAGCCGGTTgcgc                                                           >  1:6951493/1‑93 (MQ=255)
                                    caacGGGTTGCACCGCAGCAGCAATATCAGCAGCCGCAACAACCAGTTCCGCCGCGGCCGCAGTATCAGCAGCCACAACAGCCGGTTGcgccg                                                         <  2:5908252/93‑1 (MQ=255)
                                     aacCGGTTGCACCGCAGCAGCAATATCAGCAGCCGCAACAACCAGTTCCGCCGCGGCCGCAGTATCAGCAGCCACAACAGCCGGTTGcgccgc                                                        >  1:6330680/1‑93 (MQ=255)
                                      acCGGTTGCACCGCAGCAGCAATATCAGCAGCCGCAACAACCAGTTCCGCCGCGGCCGCAGTATCAGCAGCCACAACAGCCGGTTGCGCCGCa                                                       >  1:2759496/1‑93 (MQ=255)
                                        cGGTTGCACCGCAGCAGCAATATCAGCAGCCGCAACAACCAGTTCCGCCGCGGCCGCAGTATCAGCAGCCACAACAGCCGGTTGCGCCGCaac                                                     >  1:6235888/1‑91 (MQ=255)
                                         ggTTGCACCGCAGCAGCAATATCAGCAGCCGCAACAACCAGTTCCGCCGCGGCCGCAGTATCAGCAGCCACAACAGCCGGTTGCGCCGCAGcc                                                    >  1:5040740/1‑93 (MQ=255)
                                          gTTGCACCGCAGCAGCAATATCAGCAGCCGCAACAACCAGTTCCGCCGCGGCCGCAGTATCAGCAGCCACAACAGCCGGTTGCGCCGCAGCca                                                   >  1:3606608/1‑93 (MQ=255)
                                          gTTGCACCGCAGCAGCAATATCAGCAGCCGCAACAACCAGTTCCGCCGCGGCCGCAGTATCAGCAGCCACAACAGCCGGTTGCGCCGCAGCca                                                   <  2:3902500/93‑1 (MQ=255)
                                          gTTGCACCGCAGCAGCAATATCAGCAGCCGCAACAACCAGTTCCGCCGCGGCCGCAGTATCAGCAGCCACAACAGCCGGTTGCGCCGCAGCca                                                   >  2:2545168/1‑93 (MQ=255)
                                             gCACCGCAGCAGCAATATCAGCAGCCGCAACAACCAGTTCCGCCGCAGCCGCAGTATCAGCAGCCACAACAGCCGGTTGCGCCGCAGCCACaa                                                >  2:3259328/1‑93 (MQ=255)
                                             gCACCGCAGCAAAAATACCAGCAGCCCCAACACCCAGTTCCGCCGCGGCCGAAGTACCACCAGCCACAACAGCCGGTTGCGCCGCAGCCACaa                                                <  1:3401638/93‑1 (MQ=21)
                                               cccGCAGCACAAATACCACCACCCGCAACAACCAGTCCCGCCGCGGCCGCAGTATCAGCAGCCACAACAGCCGGTTGCGCCGCAGCCACAata                                              <  1:6516985/92‑1 (MQ=21)
                                               aCCGCAGCAGCAATATCAGCAGCCGCAACAACCAGTTCCGCCGCGGCCGCAGTATCAGCAGCCACAACAGCCGGTTGCGCCGCAGCCACAata                                              >  1:4949784/1‑93 (MQ=255)
                                               aCCGCAGCAGCAATATCAGCAGCCGCAACAACCAGTTCCGCCGCGGCCGCAGTATCAGCAGCCACAACAGCCGGTTGCGCCGCAGCCACAata                                              <  2:2905803/93‑1 (MQ=255)
                                               aCCGCAGCAGCAATATCAGCAGCCGCAACAACCAGTTCCGCCGCAGCCGCAGTATCAGCAGCCACAACAGCCGGTTGCGCCGCAGCCACAata                                              <  2:2336976/93‑1 (MQ=255)
                                                ccGCAGCAGCAATATCAGCAGCCGCAACAACCAGTTCCGCCGCGGCCGCAGTATCAGCAGCCACAACAGCCGGTTGCGCCGCAGCCACAatat                                             <  2:1606513/93‑1 (MQ=255)
                                                ccGCAGCAGCAATATCAGCAGCCGCAACAACCAGTTCCGCCGCGGCCGCAGTATCAGCAGCCACAACAGCCGGTTGCGCCGCAGCCACAatat                                             >  1:5745582/1‑93 (MQ=255)
                                                ccGCAGCAGCAATATCAGCAGCCGCAACAACCAGTTCCGCCGCGGCCGCAGTATCAGCAGCCACAACAGCCGGTTGCGCCGCAGCCACAatat                                             >  2:6515950/1‑93 (MQ=255)
                                                  gcagcagcaATATCAGCAGCCGCAACAACCAGTTCCGCCGCGGCCGCAGTATCAGCAGCCACAACAGCCGGTTGCGCCGCAGCCACAATATca                                           <  2:124072/93‑1 (MQ=255)
                                                  gcagcagcaATATCAGCAGCCGCAACAACCAGTTCCGCCGCGGCCGCAGTATCAGCAGCCACAACAGCCGGTTGCGCCGCAGCCACAATATca                                           >  1:7150835/1‑93 (MQ=255)
                                                  gcagcagcaATATCAGCAGCCGCAACAACCAGTTCCGCCGCAGCCGCAGTATCAGCAGCCACAACAGCCGGTTGCGCCGCAGCCACAATATca                                           <  2:3489612/93‑1 (MQ=255)
                                                   cagcagcaATATCAGCAGCCGCAACAACCAGTTCCGCCGCGGCCGCAGTATCAGCAGCCACAACAGCCGGTTGCGCCGCAGCCACAATATcag                                          >  1:4151920/1‑93 (MQ=255)
                                                     gcagcaATATCAGCAGCCGCAACAACCAGTTCCGCCGCGGCCGCAGTATCAGCAGCCACAACAGCCGGTTGCGCCGCAGCCACAATATcagca                                        >  2:6278278/1‑93 (MQ=255)
                                                     gcagcaATATCAGCAGCCGCAACAACCAGTTCCGCCGCGGCCGCAGTATCAGCAGCCACAACAGCCGGTTGCGCCGCAGCCACAATATcagca                                        >  1:2173935/1‑93 (MQ=255)
                                                     gcagcaATATCAGCAGCCGCAACAACCAGTTCCGCCGCGGCCGCAGTATCAGCAGCCACAACAGCCGGTTGCGCCGCAGCCACAATATcagca                                        >  2:4962338/1‑93 (MQ=255)
                                                     gcagcaATATCAGCAGCCGCAACAACCAGTTCCGCCGCGGCCGCAGTATCAGCAGCCACAACAGCCGGTTGCGCCGCAGCCACAATATcagca                                        >  2:4335950/1‑93 (MQ=255)
                                                      cagcaATATCAGCAGCCGCAACAACCAGTTCCGCCGCGGCCGCAGTATCAGCAGCCACAACAGCCGGTTGCGCCGCAGCCACAATATcagcag                                       >  2:3630871/1‑93 (MQ=255)
                                                         caATATCAGCAGCCGCAACAACCAGTTCCGCCGCGGCCGCAGTATCAGCAGCCACAACAGCCGGTTGCGCCGCAGCCACAATATCAGCAGCCg                                    >  2:5721759/1‑93 (MQ=255)
                                                         caATATCAGCAGCCGCAACAACCAGTTCCGCCGCGGCCGCAGTATCAGCAGCCACAACAGCCGGTTGCGCCGCAGCCACAATATCAGCAGCCg                                    >  1:977616/1‑93 (MQ=255)
                                                         caATATCAGAGCCGCAACAACCAGTTCCGCCGCGGCCGCAGTATCAGCAGCCACAACAGCCGGGTGCGCCGCAGCCACAATATCAGCAGCCGc                                   >  2:3664063/1‑93 (MQ=25)
                                                           atatCAGCAGCCGCAACAACCAGTTCCGCCGCGGCCGCAGTATCAGCAGCCACAACAGCCGGTTGCGCCGCAGCCACAATATCAGCAGCCGca                                  >  2:5464432/1‑93 (MQ=255)
                                                           atatCAGCAGCCGCAACAACCAGTTCCGCCGCGGCCGCAGTATCAGCAGCCACAACAGCCGGTTGCGCCGCAGCCACAATATCAGCAGCCGca                                  >  1:5259305/1‑93 (MQ=255)
                                                              tCAGCAGCCGCAACAACCAGTTCCGCCGCGGCCGCAGTATCAGCAGCCACAACAGCCGGTTGCGCCGCAGCCACAATATCAGCAGCCGcaaca                               >  2:3154501/1‑93 (MQ=255)
                                                              tCAGCAGCCGCAACAACCAGTTCCGCCGCGGCCGCAGTATCAGCAGCCACAACAGCCGGTTGCGCCGCAGCCACAATATCAGCAGCCGcaaca                               >  1:5531163/1‑93 (MQ=255)
                                                               cagcagCCGCAACAACCAGTTCCGCCGCGGCCGCAGTATCAGCAGCCACAACAGCCGGTTGCGCCGCAGCCACAATATCAGCAGCCGcaacaa                              <  2:5252995/93‑1 (MQ=255)
                                                                 gcagcCGCAACAACCAGTTCCGCCGCGGCCGCAGTATCAGCAGCACAACAGCCGGTTGCGCCGCAGCCACAATATCAGCAGCCGCAACAACCg                           >  1:1108886/1‑93 (MQ=255)
                                                                  cagcCGCAACAACCAGTCCCGCCGCGGCCGCAGTATCGGCAGCCACAACAGCCGGTTGCGCCGCAGCCACAATATCAGCAGCCGCAACAACCg                           <  2:5129027/93‑1 (MQ=255)
                                                                   agcCGCAACAACCAGTTCCGCCGCAGCCGCAGTATCAGCAGCCACAACAGCCGGTTGCGCCGCAGCCACAATATCAGCAGCCGCAACAACCgg                          >  2:1861492/1‑93 (MQ=255)
                                                                    gcCGCAACAACCAGTTCCGCCGCGGCCGCAGTATCAGCAGCCACAACAGCCGGTTGCGCCGCAGCCACAATATCAGCAGCCGCAACAACCGGt                         <  2:6676024/93‑1 (MQ=255)
                                                                    gcCGCAACAACCAGTTCCGCCGCGGCCGCAGTATCAGCAGCCACAACAGCCGGTTGCGCCGCAGCCACAATATCAGCAGCCGCAACAACCGGt                         >  1:607377/1‑93 (MQ=255)
                                                                     cccgaaccacCAGGTCCCGCCGCGGCCGCATTATCACCAGCCACAACAGCCGGTTGCGCCGCAGCCACAATATCAGCAGCCGCAACAACCGGt                         <  1:5340748/89‑1 (MQ=17)
                                                                     ccGCAACAACCAGTTCCGCCGCGGCCGCAGTATCAGCAGCCACAACAGCCGGTTGCGCCGCAGCCACAATATCAGCAGCCGCAACAACCGGtt                        <  2:2962236/93‑1 (MQ=255)
                                                                      cGCAACACCCAGTTCCGCCGCGGCCGCAGTATCACCAGCCACAACAGCCGGTTGCGCCGCAGCCACAATATCAGCAGCCGCAACAACCGGTTg                       <  2:166070/93‑1 (MQ=255)
                                                                       gCAACACCCAGTCCCGCCGCGGCCGCAGTATCAGCAGCCACAACAGCCGGTTGCGCCGCAGCCACAATATCAGCAGCCGCAACAACCGGTTgc                      <  2:3289028/93‑1 (MQ=255)
                                                                            aacCAGTTCCGCCGCGGCCGCAGTATCAGCAGCCACAACAGCCGGTTGCGCCGCAGCCACAATATCAGCAGCCGCAACAACCGGTTGCGCcac                 >  1:1763496/1‑93 (MQ=255)
                                                                             acCAGTTCCGCCGCGGCCGCAGTATCAGCAGCCACAACAGCCGGTTGCGCCGCAGCCACAATATCAGCAGCCGCAACAACCGGTTGCGCcaca                >  2:7492510/1‑93 (MQ=255)
                                                                             acCAGTTCCGCCGCGGCCGCAGTATCAGCAGCCACAACAGCCGGTTGCGCCGCAGCCACAATATCAGCAGCCGCAACAACCGGTTGCGCcaca                >  2:6966335/1‑93 (MQ=255)
                                                                              ccAGTTCCGCCGCGGCCGCAGTATCAGCAGCCACAACAGCCGGTTGCGCCGCAGCCACAATATCAGCAGCCGCAACAACCGGTTGCGCCAcag               >  2:6949110/1‑93 (MQ=255)
                                                                               cAGTTCCGCCGCAGCCGCAGTATCAGCAGCCACAACAGCCGGTTGCGCCGCAGCCACAATATCAGCAGCCGCAACAACCGGTTGCGCCAcagc              <  2:6994082/93‑1 (MQ=255)
                                                                               cAGTCCCCCCGGGGCCGAATTACCACCAGCCACAACGGCCGGTTGCGCGGCGGCCACAATATCAGCAGCCGCAACAACCGGTTGCGCCAcagc              <  1:3735960/93‑1 (MQ=17)
                                                                                aGTTCCGCCGCGGCCGCAGTATCAGCAGCCACAACAGCCGGTTGCGCCGCAGCCACAATATCAGCAGCCGCAACAACCGGTTGCGCCAcagca             >  2:7120407/1‑93 (MQ=255)
                                                                                aGTCCCGCCGCAGCCGCAGTATCAGCAGCCACAACAGCCGGTTGCGCCGCAGCCACAATATCAGCAGCCGCAACAACCGGTTGCGCCAcagca             <  2:5960649/93‑1 (MQ=255)
                                                                                 gTTCCGCCGCGGCCGCAGTATCAGCAGCCACAACAGCCGGTTGCGCCGCAGCCACAATATCAGCAGCCGCAACAACCGGTTGCGCCAcagcag            >  1:3018677/1‑93 (MQ=255)
                                                                                 gTTCCGCCGCAGCCGCAGTATCAGCAGCCACAACAGCCGGTTGCGCCGCAGCCACAATATCAGCAGCCGCAACAACCGGTTGCGCCAcagcag            >  1:7415827/1‑93 (MQ=255)
                                                                                  ttCCGCCGCAGCCGCAGTATCAGCAGCCACAACAGCCGGTTGCGCCGCAGCCACAATATCAGCAGCCGCAACAACCGGTTGCGCCACagcagc           >  1:2128913/1‑93 (MQ=255)
                                                                                     ccgccgCGGCCGCAGTATCAGCAGCCACAACAGCCGGTTGCGCCGCAGCCACAATATCAGCAGCCGCAACAACCGGTTGCGCCACAGCAGCaa         <  2:6786347/93‑1 (MQ=255)
                                                                                     cccccgcggccgaagTCTCACCACCCACAACAGCCGGTTGCGCCGCAGCCACAATATCAGCAGCCGCAACAACCGGTTGCGCCACAGCAGCaa         <  1:5389557/90‑1 (MQ=18)
                                                                                       gccgcggccgcaGTATCAGCAGCCACAACAGCCGGTTGCGCCGCAGCCACAATATCAGCAGCCGCAACAACCGGTTGCGCCACAGCAGCAata       <  2:5699897/93‑1 (MQ=255)
                                                                                       gccgcggccgcaGTATCAGCAGCCACAACAGCCGGTTGCGCCGCAGCCACAATATCAGCAGCCGCAACAACCGGTTGCGCCACAGCAGCAata       >  2:6076482/1‑93 (MQ=255)
                                                                                        ccgcggccgcagTATCAGCAGCCACAACAGCCGGTTGCGCCGCAGCCACAATATCAGCAGCCGCAACAACCGGTTGCGCCACAGCAGCAatat      >  1:4674624/1‑93 (MQ=255)
                                                                                         cgcggccgcagTATCAGCAGCCACAACAGCCGGTTGCGCCGCAGCCACAATATCAGCAGCCGCAACAACCGGTTGCGCCACAGCAGCAATATc     >  1:7113685/1‑93 (MQ=255)
                                                                                          gcagccgcatTACCACCAGCCACAACAGCCGGTTGCGCCGCAGCCACAATATCAGCAGCCGCAACAACCGGTTGCGCCACAGCAGCAATATca    <  1:6745158/93‑1 (MQ=255)
                                                                                            |                                                                                            
TTTACGCCAATTGTTGAACCTGTACAGCAGCCGCAACAACCGGTTGCACCGCAGCAGCAATATCAGCAGCCGCAACAACCAGTTCCGCCGCAGCCGCAGTATCAGCAGCCACAACAGCCGGTTGCGCCGCAGCCACAATATCAGCAGCCGCAACAACCGGTTGCGCCACAGCAGCAATATCAGC  >  NC_000913/934640‑934823

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