Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,145,923 T→C intergenic (‑66/‑94) yceF ← / → yceD hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,145,9230TC99.7% 1045.8 / ‑5.3 326intergenic (‑66/‑94)yceF/yceDhypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base C (156/169);  total (157/169)
Fisher's exact test for biased strand distribution p-value = 4.82e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.86e-01

Only 100 of 342 total aligned reads displayed.

AGATATAGCGTATTGATGATAGCCATTTTAAACTATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGTAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTAATATGCGCG  >  NC_000913/1145832‑1146012
                                                                                           |                                                                                         
aGATATAGCGTATTGTTGATAGCCATTTTAAACTATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCATCACTGAACGAAAATAAAGCAGTa                                                                                          >  2:9259281/1‑93 (MQ=255)
 gATATAGCGTATTGATGATAGCCATTTTAAACTATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGCaa                                                                                         >  2:7134920/1‑93 (MQ=255)
 gATATAGCGTATTGATGATAGCCATTTTAAACTATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATAAGCACTGAACGAAAATAAAGCAGCaa                                                                                         >  2:3100746/1‑93 (MQ=255)
  atatAGCGTATTGATGATAGCCATTTTAAACTATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGCAAc                                                                                        >  2:9844163/1‑93 (MQ=255)
  atatAGCGTATTGATGATAGCCATTTTAAACTATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGCAAc                                                                                        >  2:653714/1‑93 (MQ=255)
  atatAGCGTATTGATGATAGCCATTTTAAACTATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGCAAc                                                                                        >  2:8383077/1‑93 (MQ=255)
     taGCGTATTGATGATAGCCATTTTAAACTATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGCAACCCg                                                                                     >  1:6795687/1‑93 (MQ=255)
        cGTATTGATGATAGCCATTTTAAACTATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGCAACCCGCaa                                                                                  <  2:7932525/93‑1 (MQ=255)
         gTATTGATGATAGCCATTTTAAACTATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGCAACCCGCAAt                                                                                 <  2:5393267/93‑1 (MQ=255)
         gTATTGATGATAGCCATTTTAAACTATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGCAACCCGCAAt                                                                                 >  2:4129855/1‑93 (MQ=255)
          tATTGATGATAGCCATTTTAAACTATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGCAACCCGCAAtg                                                                                >  2:528586/1‑93 (MQ=255)
           aTTGATGATAGCCATTTTAAACTATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGCAACCCGCAAtgt                                                                               <  1:9408708/93‑1 (MQ=255)
             tgatgaTAGCCATTTTAAACTATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGCAACCCGCAAtgtgt                                                                             >  2:1610151/1‑93 (MQ=255)
              gatgatAGCCATTTTAAACTATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGCAACCCGCAAtgtgtg                                                                            <  2:8231399/93‑1 (MQ=255)
                tgatAGCCATTTTAAACTATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGCAACCCGCAATGTGTGCg                                                                          >  2:4306605/1‑93 (MQ=255)
                tgatAGCCATTTTAAACTATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGCAACCCGCAATGTGTGCg                                                                          >  1:3213401/1‑93 (MQ=255)
                tgatAGCCATTTTAAACTATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAACCAGCAACCCGCAATGTGTGCg                                                                          >  2:6549799/1‑93 (MQ=255)
                  atAGCCATTTTAAACTATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGCAACCCGCAATGTGTGCGaa                                                                        >  1:8787613/1‑93 (MQ=255)
                  atAGCCATTTTAAACTATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGCAACCCGCAATGTGTGCGaa                                                                        >  1:8424725/1‑93 (MQ=255)
                  atAGCCATTTTAAACTATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGCAACCCGCAATGTGTGCGaa                                                                        >  1:5770407/1‑93 (MQ=255)
                   tAGCCATTTTAAACTATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGCAACCCGCAATGTGTGCGAat                                                                       >  1:9228424/1‑93 (MQ=255)
                   tAGCCATTTTAAACTATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGCAACCCGCAATGTGTGCGAat                                                                       >  2:7993805/1‑93 (MQ=255)
                    aGCCATTTTAAACTATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGCAACCCGCAATGTGTGCGAatt                                                                      <  1:4399184/93‑1 (MQ=255)
                    aGCCATTTTAAACTATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGCAACCCGCAATGTGTGCGAatt                                                                      >  2:5208301/1‑93 (MQ=255)
                    aGCCATTTTAAACTATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGCAACCCGCAATGTGTGCGAatt                                                                      <  1:695236/93‑1 (MQ=255)
                     gCCATTTTAAACTATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGCAACCCGCAATGTGTGCGAatta                                                                     <  1:1409368/93‑1 (MQ=255)
                     gCCATTTTAAACTATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGCAACCCGCAATGTGTGCGAatta                                                                     >  1:3750730/1‑93 (MQ=255)
                      ccATTTTAAACTATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGCAACCCGCAATGTGTGCGAattat                                                                    <  2:8161773/93‑1 (MQ=255)
                       cATTTTAAACTATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGCAACCCGCAATGTGTGCGAattatt                                                                   <  1:1886717/93‑1 (MQ=255)
                       cATTTTAAACTATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGCAACCCGCAATGTGTGCGAattatt                                                                   <  1:5099589/93‑1 (MQ=255)
                        aTTTTAAACTATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGCAACCCGCAATGTGTGCGAATTATTg                                                                  <  2:10039423/93‑1 (MQ=255)
                         ttttAAACTATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGCAACCCGCAATGTGTGCGAATTATTgg                                                                 >  2:3520588/1‑93 (MQ=255)
                          tttAAACTATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGCAACCCGCAATGTGTGCGAATTATTGGc                                                                <  1:6313783/93‑1 (MQ=255)
                          tttAAACTATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGCAACCCGCAATGTGTGCGAATTATTGGc                                                                >  2:4159790/1‑93 (MQ=255)
                           ttAAACTATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGCAACCCGCAATGTGTGCGAATTATTGGCa                                                               <  1:2763695/93‑1 (MQ=255)
                            tAAACTATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGCAACCCGCAATGTGTGCGAATTATTGGCaa                                                              <  2:9301189/93‑1 (MQ=255)
                             aaaCTATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGCAACCCGCAATGTGTGCGAATTATTGGCaaa                                                             >  1:2143743/1‑93 (MQ=255)
                              aaCTATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGCAACCCGCAATGTGTGCGAATTATTGGCaaaa                                                            <  2:1505528/93‑1 (MQ=255)
                              aaCTATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGCAACCCGCAATGTGTGCGAATTATTGGCaaaa                                                            >  2:4047643/1‑93 (MQ=255)
                              aaCTATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGCAACCCGCAATGTGTGCGAATTATTGGCaaaa                                                            >  1:7796946/1‑93 (MQ=255)
                               aCTATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGCAAATCATGTAGCAACCCGCAATGTGTGCGGATTATTGGCgaaat                                                           >  2:1734496/1‑92 (MQ=255)
                                cTATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGCAACCCGCAATGTGTGCGAATTATTGGCAAAAgg                                                          <  2:3033123/93‑1 (MQ=255)
                                cTATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGCAACCCGCAATGTGTGCGAATTATTGGCAAAAgg                                                          >  2:9004714/1‑93 (MQ=255)
                                 tATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGCAACCCGCAATGTGTGCGAATTATTGGCAAAAGGc                                                         >  2:3685200/1‑93 (MQ=255)
                                 tATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGCAACCCGCAATGTGTGCGAATTATTGGCAAAAGGc                                                         >  1:6857859/1‑93 (MQ=255)
                                  aTGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAATCAGCAACCCGCAATGTGTGCGAATTATTGGCAAAAggga                                                        >  1:6058786/1‑91 (MQ=255)
                                   tGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGCAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCaa                                                       >  2:4821404/1‑93 (MQ=255)
                                    gcgcTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGCAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAAc                                                      >  1:9907814/1‑93 (MQ=255)
                                     cgcTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGCAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAAcc                                                     <  2:577880/93‑1 (MQ=255)
                                      gcTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGCAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACca                                                    <  1:756116/93‑1 (MQ=255)
                                      gcTTCGTTTTGCAAGTTGATGTTTATTATCAGCACTGAACGCAAATAAAGCAGCAACCCGCAATGTGTGCGAAATATTGGCAAAAGGAAACca                                                    <  1:8962462/93‑1 (MQ=255)
                                       cTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGCAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACcac                                                   <  1:4514040/93‑1 (MQ=255)
                                       cTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGCAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACcac                                                   <  1:818699/93‑1 (MQ=255)
                                       cTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGCAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACcac                                                   >  1:5909727/1‑93 (MQ=255)
                                        ttCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGCAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACcaca                                                  >  2:5987886/1‑93 (MQ=255)
                                        ttCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGCAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACcaca                                                  <  2:2451394/93‑1 (MQ=255)
                                         tCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGCAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAg                                                 <  1:4157462/93‑1 (MQ=255)
                                         tCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGCAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAg                                                 >  2:4616055/1‑93 (MQ=255)
                                         tCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGCAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAg                                                 >  1:1968740/1‑93 (MQ=255)
                                           gTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGCAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGc                                               <  1:1016307/93‑1 (MQ=255)
                                             tttGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGCAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTg                                             <  1:4983523/93‑1 (MQ=255)
                                             tttGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGCAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTg                                             <  2:1815574/93‑1 (MQ=255)
                                                gCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGCAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCt                                          >  1:3633764/1‑93 (MQ=255)
                                                 cAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGCAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCtt                                         <  1:3440708/93‑1 (MQ=255)
                                                   ggTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGCAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCtttt                                       >  2:659074/1‑93 (MQ=255)
                                                    gTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGCAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCttttt                                      <  1:2801305/93‑1 (MQ=255)
                                                     ttGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGCAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTc                                     <  2:3852643/93‑1 (MQ=255)
                                                      tGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGCAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCt                                    <  2:7328423/93‑1 (MQ=255)
                                                      tGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGCAACCAGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCt                                    >  2:8677303/1‑93 (MQ=255)
                                                       gATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGCAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCtt                                   >  1:6161976/1‑93 (MQ=255)
                                                       gATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGCAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCtt                                   >  1:7449806/1‑93 (MQ=255)
                                                        aTGTTTGTTATCAGCACTGAACGAAAATAAAGCAGCAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCttt                                  <  1:2278121/93‑1 (MQ=255)
                                                        aTGTTTGTTATCAGCACTGAACGAAAATAAAGCAGCAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCttt                                  >  2:6281996/1‑93 (MQ=255)
                                                        aTGTTTGTTATCAGCACTGAACGAAAATAAAGCAGCAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCttt                                  <  2:4673632/93‑1 (MQ=255)
                                                         tgtttgttATCAGCACTGAACGAAAATAAAGCAGCAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTg                                 <  2:3779961/93‑1 (MQ=255)
                                                           tttgttATCAGCACTGAACGAAAATAAAGCAGCAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGAc                               >  2:7488662/1‑93 (MQ=255)
                                                               ttATCAGCACTGAACGAAAATAAAGCAGCAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTa                           <  1:524534/93‑1 (MQ=255)
                                                               ttATCAGCACTGAACGAAAATAAAGCAGCAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTa                           >  2:2708717/1‑93 (MQ=255)
                                                                tATCAGCACTGAACGAAAATAAAGCAGCAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTAt                          <  1:7956338/93‑1 (MQ=255)
                                                                   cAGCACTGAACGAAAATAAAGCAGCAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGAc                       <  2:6775266/93‑1 (MQ=255)
                                                                      cACTGAACGAAAATAAAGCAGCAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGtt                    <  2:7730985/93‑1 (MQ=255)
                                                                      cACTGAACGAAAATAAAGCAGCAACCCGCAA‑‑TGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTAc                  >  2:112900/1‑93 (MQ=255)
                                                                        cTGAACGAAAATAAAGCAGCAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGCTAc                  >  2:1334805/1‑93 (MQ=255)
                                                                          gAACGAAAATAAAGCAGCAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACaa                >  2:6307750/1‑93 (MQ=255)
                                                                          gAACGAAAATAAAGCAGCAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACaa                <  2:2137935/93‑1 (MQ=255)
                                                                           aaCGAAAATAAAGCAGCAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACaaa               >  1:5595482/1‑93 (MQ=255)
                                                                            aCGAAAATAAAGCAGCAACCCGCAATGTGTTCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAg              <  1:1564214/93‑1 (MQ=255)
                                                                               aaaaTAAAGCAGCAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTa           >  1:1166380/1‑93 (MQ=255)
                                                                                aaaTAAAGCAGCAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACACAGTTaa          >  2:9713739/1‑93 (MQ=255)
                                                                                aaaTAAAGCAGCAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTaa          >  2:3239211/1‑93 (MQ=255)
                                                                                aaaTAAAGCAGCAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTaa          <  1:8524344/93‑1 (MQ=255)
                                                                                 aaTAAAGCAGCAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTAat         <  1:7385465/93‑1 (MQ=255)
                                                                                 aaTAAAGCAGCAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTAat         <  1:9519646/93‑1 (MQ=255)
                                                                                  aTAAAGCAGCAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTAata        <  1:7901296/93‑1 (MQ=255)
                                                                                  aTAAAGCAGCAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTAata        <  1:2562626/93‑1 (MQ=255)
                                                                                  aTAAAGCAGCAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTAata        <  1:2513850/93‑1 (MQ=255)
                                                                                  aTAAAGCAGCAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTAata        <  2:2400848/93‑1 (MQ=255)
                                                                                    aaaGCAGCAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTAATATg      <  1:3304073/93‑1 (MQ=255)
                                                                                     aaGCAGCAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTAATATgc     <  2:284180/93‑1 (MQ=255)
                                                                                       gcagcaACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTAATATgcgc   >  1:7692812/1‑93 (MQ=255)
                                                                                           |                                                                                         
AGATATAGCGTATTGATGATAGCCATTTTAAACTATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGTAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTAATATGCGCG  >  NC_000913/1145832‑1146012

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