Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,867,181 T→C N124S (AAC→AGC)  rpoS ← RNA polymerase, sigma S (sigma 38) factor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,867,1810TC98.5% 245.2 / ‑4.8 66N124S (AAC→AGC) rpoSRNA polymerase, sigma S (sigma 38) factor
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base C (27/38);  total (28/38)
Fisher's exact test for biased strand distribution p-value = 4.24e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GAATCGGCAAACGAATAGTACGGGTTTGGTTCATAATCGCCCGTTCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCG  >  NC_000913/2867044‑2867313
                                                                                                                                         |                                                                                                                                    
gAATCGGCAAACGAATAGTACGGGTTTGGTTCATAATCGCCCGTTCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGCGCGGATCAGCCCCAGGCTg                                                                                                                                    >  1:286672/1‑140 (MQ=255)
  aTCGGCAAACGAATAGTACGGGTTTGGTTCATCATCGCCCGTTCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGCGCGGATCAGCCCCAGGCTGcc                                                                                                                                  <  1:1200951/140‑1 (MQ=255)
  aTCGGCAAACGAATAGTACGGGTTTGGTTCATAATCGCCCGTTCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGCAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGCGCGGATCAGCCCCAGGCTGcc                                                                                                                                  <  2:1187959/140‑1 (MQ=255)
       cAAACGAATAGTACGGGTTTGGTTCATAATCGCCCGTTCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGCGCGGATCAGCCCCAGGCTg                                                                                                                                    >  1:483805/1‑133 (MQ=255)
       cAAACGAATAGTACGGGTTTGGTTCATAATCGCCCGTTCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGCGCGGATCAGCCCCAGGCTg                                                                                                                                    <  2:483805/133‑1 (MQ=255)
       cAAACGAATAGTACGGGTTTGGTTCATAATCGCCCGTTCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGCGCGGATCAGCCCCAGGCTGCCCTCtt                                                                                                                             <  2:932660/140‑1 (MQ=255)
       cAAACGAATAGTACGGGTTTGGTTCATAATCGCCCGTTCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGCGCGGATCAGCCCCAGGCTGCCCTCtt                                                                                                                             >  1:1376044/1‑140 (MQ=255)
           cGAATAGTACGGGTTTGGTTCATAATCGCCCGTTCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGCGCGGATCAGCCCCAGGCTGCCCTCTTCGAt                                                                                                                         <  2:742059/140‑1 (MQ=255)
            gAATAGTACGGGTTTGGTTCATAATCGCCCGTTCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGCGCGGATCAGCCCCAGGCTGCCCTCTTCGATa                                                                                                                        <  2:66576/140‑1 (MQ=255)
                    cGGGTTTGGTTCATAATCGCCCGTTCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGCGCGGATCAGCCCCAGGCTGCCCTCTTCGATAAGGTCCAg                                                                                                                <  2:482219/140‑1 (MQ=255)
                     gggTTTGGTTCATAATCGCCCGTTCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGCGCGGATCAGCCCCAGGCTGCCCTCTTCGATAAGGTCCAGc                                                                                                               <  1:1582692/140‑1 (MQ=255)
                     gggTTTGGTTCATAATCGCCCGTTCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGCGCGGATCAGCCCCAGGCTGCCCTCTTCGATAAGGTCCAGc                                                                                                               <  2:85706/140‑1 (MQ=255)
                      ggTTTGGTTCATAATCGCCCGTTCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGCGCGGATCAGCCCCAGGCTGCCCTCTTCGATAAGGTCCAGCa                                                                                                              >  2:1204652/1‑140 (MQ=255)
                       gTTTGGTTCATAATCGCCCGTTCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGCGCGGATCAGCCCCAGGCTGCCCTCTTCGATAAGGTCCAGCaa                                                                                                             >  1:472814/1‑140 (MQ=255)
                          tGGTTCATAATCGCCCGTTCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGCGCGGATCAGCCCCAGGCTGCCCTCTTCGATAAGGTCCAGCAACGc                                                                                                          >  1:898523/1‑140 (MQ=255)
                          tGGTTCATAATCGCCCGTTCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGCGCGGATCAGCCCCAGGCTGCCCTCTTCGATAAGGTCCAGCAACGc                                                                                                          >  1:740497/1‑140 (MQ=255)
                            gTTCATAATCGCCCGTTCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGCGCGGATCAGCCCCAGGCTGCCCTCTTCGATAAGGTCCAGCAACGCCa                                                                                                        >  2:356337/1‑140 (MQ=255)
                            gTTCATAATCGCCCGTTCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGCGCGGATCAGCCCCAGGCTGCCCTCTTCGATAAGGTCCAGCAACGCCa                                                                                                        >  2:338379/1‑140 (MQ=255)
                            gTTCATAATCGCCCGTTCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGCGCGGATCAGCCCCAGGCTGCCCTCTTCGATAAGGTCCAGCAACGCCa                                                                                                        >  1:1176092/1‑140 (MQ=255)
                               cATAATCGCCCGTTCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGCGCGGATCAGCCCCAGGCTGCCCTCTTCGATAAGGACCAGCAACGCCAGAc                                                                                                     <  1:779804/140‑1 (MQ=255)
                                   aTCGCCCGTTCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGCGCGGATCAGCCCCAGGCTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACg                                                                                                 <  2:521828/140‑1 (MQ=255)
                                   aTCGCCCGTTCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGCGCGGATCAGCCCCAGGCTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACg                                                                                                 <  2:969640/140‑1 (MQ=255)
                                    tCGCCCGTTCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGCGCGGATCAGCCCCAGGCTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGa                                                                                                >  1:1280126/1‑140 (MQ=255)
                                      gCCCGTTCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGAtt                                                                                              >  2:1511069/1‑140 (MQ=255)
                                      gCCCGTTCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGCGCGGATCAGCCCCAGGCTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGAtt                                                                                              >  1:160924/1‑140 (MQ=255)
                                       cccGTTCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGCGCGGATCAGCCCCAGGCTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTg                                                                                             >  2:1306280/1‑140 (MQ=255)
                                          gTTCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGCGCGGATCAGCCCCAGGCTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCa                                                                                          >  1:205257/1‑140 (MQ=255)
                                            tCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGCGCGGATCAGCCCCAGGCTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATa                                                                                        >  2:741192/1‑140 (MQ=255)
                                                    cTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGCGCGGATCAGCCCCAGGCTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCggg                                                                                <  2:447900/140‑1 (MQ=255)
                                                    cTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGCGCGGATCAGCCCCAGGCTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCggg                                                                                <  1:1318294/140‑1 (MQ=255)
                                                     tGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGCGCGGATCAGCCCCAGGCTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGc                                                                               >  1:289646/1‑140 (MQ=255)
                                                             ccaccaGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGCGCGGATCAGCCCCAGGCTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTa                                                                       <  2:1509364/140‑1 (MQ=255)
                                                             ccaccaGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGCGCGGATCAGCCCCAGGCTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTa                                                                       <  2:1193586/140‑1 (MQ=255)
                                                                ccaGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGCGCGGATCAGCCCCAGGCTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTacca                                                                    <  1:638537/140‑1 (MQ=255)
                                                                 caGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGCGCGGATCAGCCCCAGGCTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTaccac                                                                   >  1:1216654/1‑140 (MQ=255)
                                                                 caGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGCGCGGATCAGCCCCAGGCTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTaccac                                                                   <  2:289646/140‑1 (MQ=255)
                                                                 caGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGCGCGGATCAGCCCCAGGCTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAAATTTTacccc                                                                   >  2:948440/1‑138 (MQ=255)
                                                                   ggTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGCGCGGATCAGCCCCAGGCTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTAccacca                                                                 <  2:1383650/140‑1 (MQ=255)
                                                                      tGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGCGCGGATCAGCCCCAGGCTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGAc                                                              <  1:953817/140‑1 (MQ=255)
                                                                       gCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGCGCGGATCAGCCCCAGGCTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACg                                                             <  2:1376044/140‑1 (MQ=255)
                                                                           aTGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGCGCGGATCAGCCCCAGGCTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAg                                                         <  2:740497/140‑1 (MQ=255)
                                                                           aTGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGCGCGGATCAGCCCCAGGCTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAg                                                         <  2:745968/140‑1 (MQ=255)
                                                                           aTGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGCGCGGATCAGCCCCAGGCTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAg                                                         <  1:1065885/140‑1 (MQ=255)
                                                                           aTGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGCGCGGATCAGCCCCAGGCTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAg                                                         <  2:205257/140‑1 (MQ=255)
                                                                           aTGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGCGCGGATCAGCCCCAGGCTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAg                                                         <  1:356337/140‑1 (MQ=255)
                                                                              ttGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGCGCGGATCAGCCCCAGGCTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTa                                                      >  1:1372763/1‑140 (MQ=255)
                                                                                gagaAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGCGCGGATCAGCCCCAGGCTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGGTAct                                                    >  1:1371986/1‑140 (MQ=255)
                                                                                   aaGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGCGCGGATCAGCCCCAGGCTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACtctc                                                 >  1:123333/1‑140 (MQ=255)
                                                                                    aGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGCGCGGATCAGCCCCAGGCTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCg                                                <  1:512493/140‑1 (MQ=255)
                                                                                     gCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGCGCGGATCAGCCCCAGGCTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGa                                               <  1:1277628/140‑1 (MQ=255)
                                                                                        gAAACCACGTTCCGGGTCAAACTTCTCTACCGCGCGGATCAGCCCCAGGCTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTTCCACCAGACGCAAGGTACTCTCCatca                                            >  1:1040531/1‑140 (MQ=255)
                                                                                               cGTTCCGGGTCAAACTTCTCTACCGCGCGGATCAGCCCCAGGCTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTTCCACCAGACGCAAGTTACTCTCGATCATCcggcg                                     >  1:8025/1‑140 (MQ=255)
                                                                                                gTTCCGGGTCAAACTTCTCTACCGCGCGGATCAGCCCCAGGCTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCCCCCGACGCAAGTTACTCTCCATCATCcggcgg                                    >  2:600963/1‑140 (MQ=255)
                                                                                                          aaaCTTCTCTACCGCGCGGATCAGCCCCAGGCTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGa                          <  2:286672/140‑1 (MQ=255)
                                                                                                          aaaCTTCTCTACCGCGCGGATCAGCCCCAGGCTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGa                          <  1:338379/140‑1 (MQ=255)
                                                                                                          aaaCTTCTCTACCGCGCGGATCAGCCCCAGGCTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGa                          <  1:125288/140‑1 (MQ=255)
                                                                                                          aaaCTTCTCTACCGCGCGGATCAGCCCCAGGCTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGa                          <  2:1280126/140‑1 (MQ=255)
                                                                                                           aaCTTCTCTACCGCGCGGATCAGCCCCAGGCTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGGCATATGTTACTCCGAACCCAAATAATTCTCGATCATACGGCGGCGGGAGGCGAc                         >  2:647560/1‑140 (MQ=255)
                                                                                                            aCTTCTCTACCGCGCGGATCAGCCCCAGGCTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACa                        <  2:1216654/140‑1 (MQ=255)
                                                                                                            aCTTCTCTACCGCGCGGATCAGCCCCAGGCTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACa                        <  2:898523/140‑1 (MQ=255)
                                                                                                                   tACCGCGCGGATCAGCCCCAGGCTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCAc                 <  2:312103/140‑1 (MQ=255)
                                                                                                                    aCCGCGCGGATCAGCCCCAGGCTGCCCTCTTCGATAAGGCCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACg                <  2:1255333/140‑1 (MQ=255)
                                                                                                                     ccGCGCGGATCAGCCCCAGGCTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTTCTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGc               <  2:1427164/140‑1 (MQ=255)
                                                                                                                     ccGCGCGGATCAGCCCCAGGCTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGc               <  2:160924/140‑1 (MQ=255)
                                                                                                                        cgcgGATCAGCCCCAGGCTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGt            >  1:785501/1‑140 (MQ=255)
                                                                                                                                  ccccAGGCTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGGGAGGGGACATCTCCACCCCGggcgcggcgcc  >  2:100824/1‑139 (MQ=255)
                                                                                                                                         |                                                                                                                                    
GAATCGGCAAACGAATAGTACGGGTTTGGTTCATAATCGCCCGTTCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCG  >  NC_000913/2867044‑2867313

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