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.5% 708.2 / ‑4.9 218intergenic (‑66/‑94)yceF/yceDhypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base C (108/109);  total (108/110)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.05e-01

Only 100 of 223 total aligned reads displayed.

GATATAGCGTATTGATGATAGCCATTTTAAACTATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGTAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTAATATGCGC  >  NC_000913/1145833‑1146011
                                                                                          |                                                                                        
gATATAGCGTATTGATGATAGCCATTTTAAACTATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCACCACTGAACGAAAATAAAGCAGCaa                                                                                        <  2:4554592/93‑1 (MQ=255)
 atatAGCGTATTGATGATAGCCATTTTAAACTATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGCAAc                                                                                       <  1:1129001/93‑1 (MQ=255)
  tataGCGTATTGATGATAGCCATTTTAAACTATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGCAAcc                                                                                      <  1:2833660/93‑1 (MQ=255)
   ataGCGTATTGATGATAGCCATTTTAAACTATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGCAAccc                                                                                     <  1:828439/93‑1 (MQ=255)
   ataGCGTATTGATGATAGCCATTTTAAACTATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGCAAccc                                                                                     >  1:3938668/1‑93 (MQ=255)
   ataGCGTATTGATGATAGCCATTTTAAACTATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGCAAccc                                                                                     >  2:3522060/1‑93 (MQ=255)
   ataGCGTATTGATGATAGCCATTTTAAACTATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGCAAccc                                                                                     >  1:1340463/1‑93 (MQ=255)
   ataGCGTATTGATGATAGCCATTTTAAACTATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGCAAccc                                                                                     <  1:1345207/93‑1 (MQ=255)
   ataGCGTATTGATGATAGCCATTTTAAACTATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGCAAccc                                                                                     <  1:1468211/93‑1 (MQ=255)
        gTATTGATGATAGCCATTTTAAACTATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGCAACCCGCAAt                                                                                >  2:3037209/1‑93 (MQ=255)
         tATTGATGATAGCCATTTTAAACTATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGCAACCCGCAAtg                                                                               >  1:161263/1‑93 (MQ=255)
         tATTGATGATAGCCATTTTAAACTATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGCAACCCGCAAtg                                                                               <  2:383434/93‑1 (MQ=255)
         tATTGATGATAGCCATTTTAAACTATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGCAACCCGCAAtg                                                                               >  2:4000423/1‑93 (MQ=255)
           ttGATGATAGCCATTTTAAACTATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGCAACCCGCAAtgtg                                                                             <  2:1232721/93‑1 (MQ=255)
           ttGATGATAGCCATTTTAAACTATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGCAACCCGCAAtgtg                                                                             <  2:723931/93‑1 (MQ=255)
           ttGATGATAGCCATTTTAAACTATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAACCAGCAAGCTGCAAtgtg                                                                             >  1:3319405/1‑93 (MQ=255)
             gatgatAGCCATTTTAAACTATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGCAACCCGCAAtgtgtg                                                                           >  2:1429073/1‑93 (MQ=255)
             gatgatAGCCATTTTAAACTATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGCAACCCGCAAtgtgtg                                                                           >  2:621030/1‑93 (MQ=255)
              atgatAGCCATTTTAAACTATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGCAACCCGCAATGTGTGc                                                                          <  1:2794766/93‑1 (MQ=255)
              atgatAGCCATTTTAAACTATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGCAACCCGCAATGTGTGc                                                                          >  1:3856511/1‑93 (MQ=255)
                  tAGCCATTTTAAACTATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGGAGCAACCCGCAATGTGTGCGAat                                                                      >  2:4277206/1‑93 (MQ=255)
                  tAGCCATTTTAAACTATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGCAACCCGCAATGTGTGCGAat                                                                      >  2:3905011/1‑93 (MQ=255)
                    gCCATTTTAAACTATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGCAACCCGCAATGTGTGCGAatta                                                                    >  1:2115261/1‑93 (MQ=255)
                    gCCATTTTAAACTATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGCAACCCGCAATGTGTGCGAatta                                                                    <  1:565301/93‑1 (MQ=255)
                     ccATTTTAAACTATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGCAACCCGCAATGTGTGCGAattat                                                                   >  2:1939304/1‑93 (MQ=255)
                      cATTTTAAACTATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGCAACCCGCAATGTGTGCGAattatt                                                                  <  2:2913234/93‑1 (MQ=255)
                      cATTTTAAACTATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGCAACCCGCAATGTGTGCGAattatt                                                                  >  1:2212922/1‑93 (MQ=255)
                       aTTTTAAACTATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGCAACCCGCAATGTGTGCGAATTATTg                                                                 <  1:2972689/93‑1 (MQ=255)
                       aTTTTAAACTATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGCAACCCGCAATGTGTGCGAATTATTg                                                                 <  2:4589151/93‑1 (MQ=255)
                       aTTTTAAACTATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGCAACCCGCAATGTGTGCGAATTATTg                                                                 <  1:2071576/93‑1 (MQ=255)
                         tttAAACTATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGCAACCCGCAATGTGTGCGAATTATTGGc                                                               <  2:3780400/93‑1 (MQ=255)
                            aaaCTATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGCAACCCGCAATGTGTGCGAATTATTGGCaaa                                                            <  1:3145723/93‑1 (MQ=255)
                            aaaCTATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGCAACCCGCAATGTGTGCGAATTATTGGCaaa                                                            >  2:3455311/1‑93 (MQ=255)
                             aaCTATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGCAACCCGCAATGTGTGCGAATTATTGGCaaaa                                                           <  1:3031809/93‑1 (MQ=255)
                             aaCTATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGCAACCCGCAATGTGTGCGAATTATTGGCaaaa                                                           <  2:1011511/93‑1 (MQ=255)
                              aCTATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGCAACCCGCAATGTGTGCGAATTATTGGCAAAAg                                                          >  2:2027155/1‑93 (MQ=255)
                              aCTATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGCAACCCGCAATGTGTGCGAATTATTGGCAAAAg                                                          <  2:4243085/93‑1 (MQ=255)
                              aCTATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGCAACCCGCAATGTGTGCGAATTATTGGCAAAAg                                                          <  1:4429331/93‑1 (MQ=255)
                                tATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGCAACCCGCAATGTGTGCGAATTATTGGCAAAAGGc                                                        >  2:809809/1‑93 (MQ=255)
                                 aTGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGCAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCa                                                       >  2:2539460/1‑93 (MQ=255)
                                 aTGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGCAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCa                                                       >  1:642192/1‑93 (MQ=255)
                                   gcgcTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGCAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAAc                                                     >  1:3186662/1‑93 (MQ=255)
                                     gcTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGCAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACca                                                   <  1:1228695/93‑1 (MQ=255)
                                      cTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGCAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACcac                                                  >  2:4397555/1‑93 (MQ=255)
                                      cTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGCAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACcac                                                  >  2:4715002/1‑93 (MQ=255)
                                       ttCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGCAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACcaca                                                 >  2:3932273/1‑93 (MQ=255)
                                        tCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGCAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAg                                                <  1:4747252/93‑1 (MQ=255)
                                        tCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGCAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAg                                                >  1:4077382/1‑93 (MQ=255)
                                        tCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGCAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAg                                                >  2:3274198/1‑93 (MQ=255)
                                        tCGTTTCGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGCAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAg                                                >  1:2939855/1‑93 (MQ=255)
                                          gTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGCAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGc                                              >  2:2838117/1‑93 (MQ=255)
                                           ttttGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGCAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAggat                                             <  1:1523966/93‑3 (MQ=255)
                                           ttttGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGCAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCt                                             <  2:1265079/93‑1 (MQ=255)
                                           ttttGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGCAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCt                                             >  2:403994/1‑93 (MQ=255)
                                           ttttGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGCAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCt                                             >  1:2417247/1‑93 (MQ=255)
                                               gCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGCAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCt                                         >  1:1749699/1‑93 (MQ=255)
                                                cAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGCAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCtt                                        <  2:359924/93‑1 (MQ=255)
                                                 aGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGCAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCttt                                       <  2:2714414/93‑1 (MQ=255)
                                                 aGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGCAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCttt                                       >  2:2529066/1‑93 (MQ=255)
                                                    ttGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGCAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTc                                    >  1:2387144/1‑93 (MQ=255)
                                                      gATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGCAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCtt                                  >  1:2747096/1‑93 (MQ=255)
                                                      gATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGCAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCtt                                  <  1:3599801/93‑1 (MQ=255)
                                                      gATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGCAACCCGCAATGCGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCtt                                  >  1:556858/1‑93 (MQ=255)
                                                        tgtttgttATCAGCACTGAACGAAAATAAAGCAGCAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTg                                <  1:2385852/93‑1 (MQ=255)
                                                         gtttgttATCAGCACTGAACGAAAATAAAGCCGCAACCCGCAATGTGTGCGAATTATTGGAAAAAGGCAACCACAGACTGCCTTTTTTTTTgc                               >  1:1653123/1‑92 (MQ=255)
                                                          tttgttATCAGCACTGAACGAAAATAAAGCAGCAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGAc                              <  1:3756664/93‑1 (MQ=255)
                                                           ttgttATCAGCACTGAACGAAAATAAAGCAGCAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGAct                             >  1:3174669/1‑93 (MQ=255)
                                                           ttgttATCAGCACTGAACGAAAATAAAGCAGCAACCCGCAATGTGTGCAAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGatt                             <  2:2243115/93‑3 (MQ=255)
                                                            tgttATCAGCACTGAACGAAAATAAAGCAGCAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGActc                            <  2:2184740/93‑1 (MQ=255)
                                                             gttATCAGCACTGAACGAAAATAAAGCAGCAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGActct                           <  2:648300/93‑1 (MQ=255)
                                                              ttATCAGCACTGAACGAAAATAAAGCAGCAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTa                          >  2:1298252/1‑93 (MQ=255)
                                                               tATCAGCACTGAACGAAAATAAAGCAGCAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTAt                         <  1:3013243/93‑1 (MQ=255)
                                                                 tCAGCACTGAACGAAAATAAAGCAGCAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGa                       <  1:4314715/93‑1 (MQ=255)
                                                                  cAGCACTGAACGAAAATAAAGCAGCAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGAc                      <  2:547915/93‑1 (MQ=255)
                                                                   aGCACTGAACGAAAATAAAGCAGCAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACg                     >  2:163519/1‑93 (MQ=255)
                                                                    gCACTGAACGAAAATAAAGCAGCAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGt                    <  2:1595349/93‑1 (MQ=255)
                                                                      aCTGAACGAAAATAAAGCAGCAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTa                  >  1:2851117/1‑93 (MQ=255)
                                                                      aCTGAACGAAAATAAAGCAGCAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTa                  >  1:1972627/1‑93 (MQ=255)
                                                                        tGAACGAAAATAAAGCAGCAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACa                >  1:4090887/1‑93 (MQ=255)
                                                                        tGAACGAAAATAAAGCAGCAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACa                >  1:415075/1‑93 (MQ=255)
                                                                        tGAACGAAAATAAAGCAGCAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTGTTCTTTGACTCTATGACGTTACa                >  2:3394219/1‑93 (MQ=255)
                                                                         gAACGAAAATAAAGCAGCAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACaa               >  1:3461007/1‑93 (MQ=255)
                                                                         gAACGAAAATAAAGCAGCAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACaa               >  2:1803658/1‑93 (MQ=255)
                                                                         gAACGAAAATAAAGCAGCAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACaa               <  1:681494/93‑1 (MQ=255)
                                                                         gAACGAAAATAAAGCAGCAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACaa               >  1:1398098/1‑93 (MQ=255)
                                                                          aaCGAAAATAAAGCAGCGACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCATTTTTCTTTGACTCTATGACGTTACaaa              <  1:4215814/93‑1 (MQ=255)
                                                                          aaCGAAAATAAAGCAGCAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACaaa              <  2:4457866/93‑1 (MQ=255)
                                                                           aCGAAAATAAAGCAGCAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAg             <  2:614018/93‑1 (MQ=255)
                                                                            cGAAAATAAAGCAGCAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGt            >  1:3979711/1‑93 (MQ=255)
                                                                             gAAAATAAAGCAGCAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGtt           >  2:1079217/1‑93 (MQ=255)
                                                                               aaaTAAAGCAGCAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTaa         <  2:2929786/93‑1 (MQ=255)
                                                                               aaaTAAAGCAGCAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTaa         <  2:1728418/93‑1 (MQ=255)
                                                                                aaTAAAGCAGCAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTAat        <  2:874222/93‑1 (MQ=255)
                                                                                aaTAAAGCAGCAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTAat        <  2:500671/93‑1 (MQ=255)
                                                                                aaTAAAGCAGCAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTAat        >  2:2436972/1‑93 (MQ=255)
                                                                                aaTAAAGCAGCAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTAat        >  2:2024019/1‑93 (MQ=255)
                                                                                 aTAAAGCAGCAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTAata       <  1:3371653/93‑1 (MQ=255)
                                                                                 aTAAAGCAGCAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTAata       <  2:1336657/93‑1 (MQ=255)
                                                                                     agcagcAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTAATATgcg   <  1:4604486/93‑1 (MQ=255)
                                                                                      gcagcaACCCGCAATGTGTGCGAATTATTGGCAAAAGGAAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTAATATgcgc  >  1:2450537/1‑93 (MQ=255)
                                                                                          |                                                                                        
GATATAGCGTATTGATGATAGCCATTTTAAACTATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGTAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTAATATGCGC  >  NC_000913/1145833‑1146011

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