Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,196,120 G→A 100% intergenic (‑21/+212) thiC ← / ← rsd phosphomethylpyrimidine synthase/stationary phase protein, binds sigma 70 RNA polymerase subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,196,1200GA96.4% 284.3 / 0.3 84intergenic (‑21/+212)thiC/rsdphosphomethylpyrimidine synthase/stationary phase protein, binds sigma 70 RNA polymerase subunit
Reads supporting (aligned to +/- strand):  ref base G (0/3);  new base A (49/32);  total (49/35)
Fisher's exact test for biased strand distribution p-value = 6.87e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.69e-01

GCGGTGCCTTCCAGGGTGTCGATAAAATGTTGGGCCCGGGCGCGTTGTTCGCGGCGGGTCAGTTTTGTTGCAGACATAGCTCATTCCAAAAAGTTAAGGACGTGGCTTGTCAGACGACGGATGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCGCAGGTAT  >  NC_000913/4196023‑4196211
                                                                                                 |                                                                                           
gCGGTGCCTTCCAGGGTGTCGATAAAATGTTGGGCCCGGGCGCGTTGTTCGCGGCGGGTCAGTTTTGTTGCAGACATAGATAAATCAAAAAAGTTAAAGa                                                                                           <  2:58778/100‑1 (MQ=255)
 cGGTGCCTTCCAGGGTGTCGATAAAATGTTGGGCCCGGGTGCGTTGTTCGCGGCGGGGCACTTTTTTGGCAGACATAGCGCATTCCAAAAAGTTAAAGa                                                                                           >  2:1454826/1‑99 (MQ=255)
     gCCTTCCAGGGTGTCGATAAAATGTTGGGCCCGGGCGCGTTGTTCGCGGCGGGTCAGTTTTGTTGCAGACATAGCTCATTCCAAAAAGTTAAAGACGTgg                                                                                      >  1:308750/1‑100 (MQ=255)
     gCCTTCCAGGGTGTCGATAAAATGTTGGGCCCGGGCGCGTTGTTCGCGGCGGGTCAGTTTTGTTGCAGACATAGCTCATTCCAAAAAGTTAAAGACGTgg                                                                                      >  1:308761/1‑100 (MQ=255)
       cTTCCAGGGTGTCGATAAAATGTTGGGCCCGGGCGCGTTGTTCGCGGCGGGTCAGTTTTGTTGCAGACATAGCTCATTCCAAAAAGTTAAAGACGTGGCt                                                                                    >  1:1798612/1‑100 (MQ=255)
       cTTCCAGGGTGTCGATAAAATGTTGGGCCCGGGCGCGTTGTTCGCGGCGGGTCAGTTTTGTTGCAGACATAGCTCATTCCAAAAAGTTAAAGACGTGGCt                                                                                    >  1:1444823/1‑100 (MQ=255)
       cTTCCAGGGTGTCGATAAAACGTTGGGCCCGGGCGCGTTGTTCGCGGCGGGTCAGTTTTGTTGCAGACATAGCTCATTCCACAAAGTTAAAGACGTGGCt                                                                                    >  2:637939/1‑100 (MQ=255)
         tCCAGGGTGTCGATAAAATGTTGGGCCCGGGCGCGTTGTTCGCGGCGGGTCAGTTTTGTTGCAGACATAGCTCATTCCAAAAAGTTAAAGACGTGGCTTg                                                                                  <  1:1326128/100‑1 (MQ=255)
         tCCAGGGTGTCGATAAAATGTTGGGCCCGGGCGCGTTGTTCGCGGCGGGTCAGTTTTGTTGCAGACATAGCTCATTCCAAAAAGTTAAAGACGTGGCTTg                                                                                  >  1:818525/1‑100 (MQ=255)
         tCCAGGGTGTCGATAAAATGTTGGGCCCGGGCGCGTTGTTCGCGGCGGGTCAGTTTTGTTGCAGACATAGCTCATTCCAAAAAGTTAAAGACGTGGCTTg                                                                                  <  2:556479/100‑1 (MQ=255)
         tCCAGGGTGTCGATAAAATGTTGGGCCCGGGCGCGTTGTTCGCGGCGGGTCAGTTTTGTTGCAGACATAGCTCATTCCAAAAAGTTAAAGACGTGGCTTg                                                                                  >  2:464403/1‑100 (MQ=255)
             gggTGTCGATAAAATGTTGGGCCCGGGCGCGTTGTTCGCGGCGGGTCAGTTTTGTTGCAGACATAGCTCATTCCAAAAAGTTAAAGACGTGGCTTGTCAg                                                                              <  1:837111/100‑1 (MQ=255)
                tgtTGATTAAATGTTGGGCCCGGGCAAGTGATTCGCGTCGGGACAGTTTTGTTGCAGACATAGCTAAGTCCTAAATGTTAAGGTCGTGGCTTGTCAgacg                                                                           <  2:1862542/100‑1 (MQ=255)
                  tCGATAAAATGTTGGGCCCGGGCGCGTTGTTCGCGGCGGGTCAGTTTTGTTGCAGACATAGCTCATTCCAAAAAGTTAAAGACGTGGCTTGTCAgacgac                                                                         >  2:2004704/1‑100 (MQ=255)
                   cGATAAAATGTTGGGCCCGGGCGCGTTGTTCGCGGCGGGTCAGTTTTGTTGCAGACATAGCTCATTCCAAAAAGTTAAAGACGTGGCTTGTCAGacgacg                                                                        <  1:1002733/100‑1 (MQ=255)
                    gATAAAATGTTGGGCCCGGGCGCGTTGTTCGCGGCGGGTCAGTTTTGTTGCAGACATAGCTCATTCCAAAAAGTTAAAGACGTGGCTTGTCAGACGACgg                                                                       >  2:773671/1‑100 (MQ=255)
                      tAAAATGTTGGGCCCGGGCGCGTTGTTCGCGGCGGGTCAGTTTTGTTGCAGACATAGCTCATTCCAAAAAGTTAAAGACGTGGCTTGTCAGACGACGGAt                                                                     <  2:1106199/100‑1 (MQ=255)
                       aaaaTGTTGGGCCCGGGCGCGTTGTTCGCGGCGGGTCAGTTTTGTTGCAGACATAGCTCATTCCAAAAAGTTAAAGACGTGGCTTGTCAGACGACGGATg                                                                    >  2:1503586/1‑100 (MQ=255)
                        aaaTGTTGGGCCCGGGCGCGTTGTTCGCGGCGGGTCAGTTTTGTTGCAGACATAGCTCATTCCAAAAAGTTAAAGACGTGGCTTGTCAGACGACGGATGa                                                                   <  2:1705651/100‑1 (MQ=255)
                        aaaTGTTGGGCCCGGGCGCGTTGTTCGCGGCGGGTCAGTTTTGTTGCAGACATAGCTCATTCCAAAAAGTTAAAGACGTGGCTTGTCAGACGACGGATGa                                                                   <  1:1180407/100‑1 (MQ=255)
                           tGTTGGGCCCGGGCGCGTTGTTCGCGGCGGGTCAGTTTTGTTGCAGACATAGCTCATTCCAAAAAGTTAAAGACGTGGCTTGTCAGACGACGGATGAAGc                                                                <  1:1762733/100‑1 (MQ=255)
                             ttGGGCCCGGGCGCGTTGTTCGCGGCGGGTCAGTTTTGTTGCAGACATAGCTCATTCCAAAAAGTTAAAGACGTGGCTTGTCAGACGACGGATGAAGCaa                                                              <  2:818522/100‑1 (MQ=255)
                              tctgCCCGGGCGCGTTGTTCGCGGCGGGTCAGTTTTGTTGCAGACATAGCTCATTCCAAAAAGTTAAAGACGTGGCTTGTCAGACGACGGATGAAGCAag                                                             <  2:408064/97‑1 (MQ=255)
                              tGGGCCCGGGCGCGTTGTTCGCGGCGGGTCAGTTTTGTTGCAGACATAGCTCATTCCAAAAAGTTAAAGACGTGGCTTGTCAGACGACGGATGAAGCAag                                                             >  2:796640/1‑100 (MQ=255)
                               gggCCCGGGCGCGTTGTTCGCGGCGGGTCAGTTTTGTTGCAGACATAGCTCATTCCAAAAAGTTAAAGACGTGGCTTGTCAGACGACGGATGAAGCAaga                                                            >  2:405959/1‑100 (MQ=255)
                                ggCCCGGGCGCGTTGTTCGCGGCGGGTCAGTTTTGTTGCAGACATAGCTCATTCCAAAAAGTTAAAGACGTGGCTTGTCAGACGACGGATGAAGCAagag                                                           <  1:553380/100‑1 (MQ=255)
                                 gCCCGGGCGCGTTGTTCGCGGCGGGTCAGTTTTGTTGCAGACATAGCTCATTCCAAAAAGTTAAAGACGTGGCTTGTCAGACGACGGATGAAGCAAgaga                                                          >  2:1193121/1‑100 (MQ=255)
                                 gCCCGGGCGCGTTGTTCGCGGCGGGTCAGTTTTGTTGCAGACATAGCTCATTCCAAAAAGTTAAAGACGTGGCTTGTCAGACGACGGATGAAGCAAgaga                                                          >  1:1550579/1‑100 (MQ=255)
                                 gCCCGGGCGCGTTGTTCGCGGCGGGTCAGTTTTGTTGCAGACATAGCTCATTCCAAAAAGTTAAAGACGTGGCTTGTCAGACGACGGATGAAGCAAgaga                                                          <  2:719205/100‑1 (MQ=255)
                                 gCCCGGGCGCGTTGTTCGCGGCGGGTCAGTTTTGTTGCAGACATAGCTCATTCCAAAAAGTTAAAGACGTGGCTTGTCAGACGACGGATGAAGCAAgaga                                                          >  1:408067/1‑100 (MQ=255)
                                    cGGGCGCGTTGTTCGCGGCGGGTCAGTTTTGTTGCAGACATAGCTCATTCCAAAAAGTTAAAGACGTGGCTTGTCAGACGACGGATGAAGCAAGAGACGa                                                       >  2:998544/1‑100 (MQ=255)
                                        cgcgTTGTTCGCGGCGGGTCAGTTTTGTTGCAGACATAGCTCATTCCAAAAAGTTAAGGACGTGGCTTGTCAGACGACGGATGAAGCAAGAGACGATcgc                                                   <  2:1435172/100‑1 (MQ=255)
                                          cgTTGTTCGCGGCGGGTCAGTTTTGTTGCAGACATAGCTCATTCCAAAAAGTTAAAGACGTGGCTTGTCAGACGACGGATGAAGCAAGAGACGATcgccg                                                 <  1:1828775/100‑1 (MQ=255)
                                          cgTTGTTCGCGGCGGGTCAGTTTTGTTGCAGACATAGCTCATTCCAAAAAGTTAAAGACGTGGCTTGTCAGACGACGGATGAAGCAAGAGACGATcgccg                                                 >  1:240968/1‑100 (MQ=255)
                                          cgTTGTTCGCGGCGGGTCAGTTTTGTTGCAGACATAGCTCATTCCAAAAAGTTAAAGACGTGGCTTGTCAGACGACGGATGAAGCAAGAGACGATcgccg                                                 >  1:991077/1‑100 (MQ=255)
                                            ttgttCGCGGCGGGTCAGTTTTGTTGCAGACATAGCTCATTCCAAAAAGTTAAAGACGTGGCTTGTCAGACGACGGATGAAGCAAGAGACGATCGCCGCa                                               >  2:618131/1‑100 (MQ=255)
                                              gttCGCGGCGGGTCAGTTTTGTTGCAGACATAGCTCATTCCAAAAAGTTAAAGACGTGGCTTGTCAGACGACGGATGAAGCAAGAGACGATCGCCGCAgg                                             >  2:550261/1‑100 (MQ=255)
                                               ttCGCGGCGGGTCAGTTTTGTTGCAGACATAGCTCATTCCAAAAAGTTAAAGACGTGGCTTGTCAGACGACGGATGAAGCAAGAGACGATCGCCGCAggg                                            <  1:1045779/100‑1 (MQ=255)
                                               ttCGCGGCGGGTCAGTTTTGTTGCAGACATAGCTCATTCCAAAAAGTTAAAGACGTGGCTTGTCAGACGACGGATGAAGCAAGAGACGATCGCCGCAggg                                            <  1:1644405/100‑1 (MQ=255)
                                                tCGCGGCGGGTCAGTTTTGTTGCAGACATAGCTCATTCCAAAAAGTTAAAGACGTGGCTTGTCAGACGACGGATGAAGCAAGAGACGATCGCCGCAgggg                                           <  1:474284/100‑1 (MQ=255)
                                                 cgcgGCGGGTCAGTTTTGTTGCAGACATAGCTCATTCCAAAAAGTTAAAGACGTGGCTTGTCAGACGACGGATGAAGCAAGAGACGATCGCCGCAGGGgc                                          <  1:1740159/100‑1 (MQ=255)
                                                 cgcgGCGGGTCAGTTTTGTTGCAGACATAGCTCATTCCAAAAAGTTAAAGACGTGGCTTGTCAGACGACGGATGAAGCAAGAGACGATCGCCGCAGGGgc                                          <  2:1444813/100‑1 (MQ=255)
                                                 cgccgcGGGTCAGTTTTGTTGCAGACATAGCTCATTCCAAAAAGTTAAAGACGTGGCTTGTCAGACGACGGATGAAGCAAGAGACGATCGCCGCAGGGgc                                          <  1:187872/100‑1 (MQ=255)
                                                      cggGTCAGTTTTGTTGCAGACATAGCTCATTCCAAAAAGTTAAAGACGTGGCTTGTCAGACGACGGATGAAGCAAGAGACGATCGCCGCAGGGgcgatgc                                     >  2:21052/1‑100 (MQ=255)
                                                      cggGTCAGTTTTGTTGCAGACATAGCTCATTCCAAAAAGTTAAAGACGTGGCTTGTCAGACGACGGATGAAGCAAGAGACGATCGCCGCAGGGgcgat                                       >  2:965099/1‑98 (MQ=255)
                                                        ggTCAGTTTTGTTGCAGACATAGCTCATTCCAAAAAGTTAAAGACGTGGCTTGTCAGACGACGGATGAAGCAAGAGACGATCGCCGCAGGGgcgatgcga                                   <  1:173059/100‑1 (MQ=255)
                                                         gTCAGTTTTGTTGCAGACATAGCTCATTCCAAAAAGTTAAAGACGTGGCTTGTCAGACGACGGATGAAGCAAGAGACGATCGCCGCAGGGgcgctgc                                     >  2:1599608/1‑97 (MQ=255)
                                                            aGTTTTGTTGCAGACATAGCTCATTCCAAAAAGTTAAAGACGTGGCTTGTCAGACGACGGATGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGc                               >  1:637252/1‑100 (MQ=255)
                                                             gTTTTGTTGCAGACATAGCTCATTCCAAAAAGTTAAAGACGTGGCTTGTCAGACGACGGATGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCa                              >  1:948564/1‑100 (MQ=255)
                                                             gTTTTGTTGCAGACATAGCTCATTCCAAAAAGTTAAAGACGTGGCTTGTCAGACGACGGATGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCa                              <  1:1717757/100‑1 (MQ=255)
                                                             gTTTTGTTGCAGACATAGCTCATTCCAAAAAGTTAAAGACGTGGCTTGTCAGACGACGGATGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCa                              >  1:678910/1‑100 (MQ=255)
                                                              ttttGTTGCAGACATAGCTCATTCCAAAAAGTTAAAGACGTGGCTTGTCAGACGACGGATGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAg                             <  1:1480100/100‑1 (MQ=255)
                                                                 tgttgCAGACATAGCTCATTCCAAAAAGTTAAAGACGTGGCTTGTCAGACGACGGATGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGAtt                          >  2:251789/1‑100 (MQ=255)
                                                                 tgttgCAGACATAGCTCATTCCAAAAAGTTAAAGACGTGGCTTGTCAGACGACGGATGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGAtt                          >  1:1633208/1‑100 (MQ=255)
                                                                 tgttgCAGACATAGCTCATTACAAAAAGTTAAAGACGTGGCTTGTCAGACGACGGATGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGAtt                          <  1:1530305/100‑1 (MQ=255)
                                                                  gttgCAGACATAGCTCATTCCAAAAAGTTAAAGACGTGGCTTGTCAGACGACGGATGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTa                         <  1:672838/100‑1 (MQ=255)
                                                                  gttgCAGACATAGCTCATTCCAAAAAGTTAAAGACGTGGCTTGTCAGACGACGGATGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTa                         <  1:280445/100‑1 (MQ=255)
                                                                    tgCAGACATAGCTCATTCCAAAAAGTTAAAGACGTGGCTTGTCAGACGACGGATGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAAc                       >  1:1217492/1‑100 (MQ=255)
                                                                     gCAGACATAGCTCATTCCAAAAAGTTAAAGACGTGGCTTGTCAGACGACGGATGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAAct                      >  1:311216/1‑100 (MQ=255)
                                                                     gCAGACATAGCTCATTCCAAAAAGTTAAAGACGTGGCTTGTCAGACGACGGATGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAAct                      <  1:901150/100‑1 (MQ=255)
                                                                     gCAGACATAGCTCATTCCAAAAAGTTAAAGACGTGGCTTGTCAGACGACGGAAGAAGCAAGAGACGATCGCCGCAGGGGCGAAGCGATAGCAGATTAAc                       >  2:1453099/1‑99 (MQ=255)
                                                                      cAGACATAGCTCATTCCAAAAAGTTAAAGACGTGGCTTGTCAGACGACGGATGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAActc                     >  2:887537/1‑100 (MQ=255)
                                                                      cAGACATAGCTCATTCCAAAAAGTTAAAGACGTGGCTTGTCAGACGACGGATGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAActc                     <  1:745047/100‑1 (MQ=255)
                                                                       aGACATAGCTCATTCCAAAAAGTTAAAGACGTGGCTTGTCAGACGACGGATGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAActct                    >  1:328319/1‑100 (MQ=255)
                                                                       aGACATAGCTCATTCCAAAAAGTTAAAGACGTGGCTTGTCAGACGACGGATGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAActct                    <  2:540281/100‑1 (MQ=255)
                                                                        gACATAGCTCAGTCCAAAAAGTTAAAGACGTGGCTTGTCAGACGACGGATGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCtt                   <  1:744010/100‑1 (MQ=255)
                                                                            tAGCTCATTCCAAAAAGTTAAAGACGTGGCTTGTCAGACGACGGATGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTc               >  2:1961486/1‑100 (MQ=255)
                                                                             aGCTCATTCCAAAAAGTTAAAGACGTGGCTTGTCAGACGACGGATGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTcc              >  2:469218/1‑100 (MQ=255)
                                                                              gCTCATTCCAAAAAGTTAAAGACGTGGCTTGTCAGACGACGGATGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTccc             >  2:619462/1‑100 (MQ=255)
                                                                               cTCATTCCAAAAAGTTAAAGACGTGGCTTGTCAGACGACGGATGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCt            >  2:1036853/1‑100 (MQ=255)
                                                                                tCATTCCAAAAAGTTAAAGACGTGGCTTGTCAGACGACGGATGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCtt           >  2:1547662/1‑100 (MQ=255)
                                                                                tCATTCCAAAAAGTTAAAGACGTGGCTTGTCAGACGACGGATGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCtt           >  2:1982396/1‑100 (MQ=255)
                                                                                tCATTCCAAAAAGTTAAAGACGTGGCTTGTCAGACGACGGATGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCtt           >  1:1836911/1‑100 (MQ=255)
                                                                                 cATTCCAAAAAGTTAAAGACGTGGCTTGTCAGACGACGGATGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTc          >  1:700440/1‑100 (MQ=255)
                                                                                   ttCCAAAAAGTTAAAGACGTGGCTTGTCAGACGACGGATGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCGc        <  2:1511256/100‑1 (MQ=255)
                                                                                    tCCAAAAAGTTAAAGACGTGGCTTGTCAGACGACGGATGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCGCa       >  1:1250978/1‑100 (MQ=255)
                                                                                    tCCAAAAAGTTAAAGACGTGGCTTGTCAGACGACGGATGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCGCa       >  1:1238313/1‑100 (MQ=255)
                                                                                     ccAAAAAGTTAAAGACGTGGCTTGTCAGACGACGGATGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCGCAg      >  1:655840/1‑100 (MQ=255)
                                                                                     ccAAAAAGTTAAAGACGTGGCTTGTCAGACGACGGATGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCGCAg      >  2:844277/1‑100 (MQ=255)
                                                                                      cAAAAAGTTAAGGACGTGGCTTGTCAGACGACGGATGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCGCAgg     <  1:61638/100‑1 (MQ=255)
                                                                                      cAAAAAGTTAAAGACGTGGCTTGTCAGACGACGGATGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCGCAgg     >  2:744652/1‑100 (MQ=255)
                                                                                      cAAAAAGTTAAAGACGTGGCTTGTCAGACGACGGATGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCGCAgg     >  1:2062210/1‑100 (MQ=255)
                                                                                         aaaGTTAAAGACGTGGCTTGTCAGACGACGGATGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCGCAGGTAt  >  1:752650/1‑100 (MQ=255)
                                                                                         aaaGTTAAAGACGTGGCTTGTCAGACGACGGATGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCGCAGGTAt  <  1:229159/100‑1 (MQ=255)
                                                                                                 |                                                                                           
GCGGTGCCTTCCAGGGTGTCGATAAAATGTTGGGCCCGGGCGCGTTGTTCGCGGCGGGTCAGTTTTGTTGCAGACATAGCTCATTCCAAAAAGTTAAGGACGTGGCTTGTCAGACGACGGATGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCGCAGGTAT  >  NC_000913/4196023‑4196211

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