Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,822,445 T→A 6.5% E65D (GAA→GAT) ‡ recX ← regulatory protein for RecA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,822,4450TA 6.5% 270.7 / 7.1 77E65D (GAA→GAT) ‡recXregulatory protein for RecA
Reads supporting (aligned to +/- strand):  ref base T (37/35);  new base A (2/3);  total (39/38)
Fisher's exact test for biased strand distribution p-value = 6.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.37e-01

CACATTCACGCATCGCTTTTTCTGTCGCTTCGCGGGAAATACCTTTCTGATTCAGTTCCTGGCGAATACGCGCAGGTCCATAACCTTTGCGGCTACGGCTGGCGATAAAGCGCGCAACAAATCGGCTGTCATCGAGATAGCCATGTTCATGGCACCAGGCAATAACGCGCTCGTAATCTTCTGCCGTAGCATCAATCTCTTCTGGGCCATTTTTGCCCATAATCGGTGCCGCGAGTTTACGTCGCAGTTCTTGCTCACTGTGATCGCGCACCGCCAGAATGCGTACCGCA  >  NC_000913/2822300‑2822589
                                                                                                                                                 |                                                                                                                                                
cacaTTCACGCATCGCTTTTTCTGTCGCTTCGCGGGAAATACCTTTCTGATTCAGTTCCTGGCGAATACGCGCAGGTCCATAACCTTTGCGGCTACGGCTGGCGATAAAGCGCGCAACAAATCGGCTGTCATCGAGATAGCCATGTTCa                                                                                                                                               <  3:4023/149‑1 (MQ=255)
 acaTTCACGCATCGCTTTTTCTGTCGCTTCGCGGGAAATACCTTTCTGATTCAGTTCCTGGCGAATACGCGCAGGTCCATAACCTTTGCGGCTACGGCTGGCGATAAAGCGCGCAACAAATCGGCTGTCATCGAGATAGCCATGTTCAt                                                                                                                                              <  3:433009/149‑1 (MQ=255)
  caTTCACGCATCGCTTTTTCTGTCGCTTCGCGGGAAATACCTTTCTGATTCAGTTCCTGGCGAATACGCGCAGGTCCATAACCTTTGCGGCTACGGCTGGCGATAAAGCGCGCAACAAATCGGCTGTCATCGAGATAGCCATGTTCATg                                                                                                                                             <  1:123107/149‑1 (MQ=255)
   aTTCACGCATCGCTTTTTCTGTCGCTTCGCGGGAAATACCTTTCTGATTCAGTTCCTGGCGAATACGCGCAGGTCCATAACCTTTGCGGCTACGGCTGGCGATAAAGCGCGCAACAAATCGGCTGTCATCGAGATAGCCATGTTCATgg                                                                                                                                            >  3:528588/1‑149 (MQ=255)
   aTTCACGCATCGCTTTTTCTGTCGCTTCGCGGGAAATACCTTTCTGATTCAGTTCCTGGCGAATACGCGCAGGTCCATAACCTTTGCGGCTACGGCTGGCGATAAAGCGCGCAACAAATCGGCTGTCATCGAGATAGCCATGTTCATgg                                                                                                                                            >  3:545133/1‑149 (MQ=255)
     tCACGCATCGCTTTTTCTGTCGCTTCGCGGGAAATACCTTTCTGATTCAGTTCCTGGCGAATACGCGCAGGTCCATAACCTTTGCGGCTACGGCTGGCGATAAAGCGCGCAACAAATCGGCTGTCATCGAGATAGCCATGAACATGGct                                                                                                                                          <  2:510639/149‑2 (MQ=255)
         gCATCGCTTTTTCTGTCGCTTCGCGGGAAATACCTTTCTGATTCAGTTCCTGGCGAATACGCGCAGGTCCATAACCTTTGCGGCTACGGCTGGCGATAAAGCGCGCAACAAATCGGCTGTCATCGAGATAGCCATGTTCATGGCACCAg                                                                                                                                      <  1:599666/149‑1 (MQ=255)
             cGCTTTTTCTGTCGCTTCGCGGGAAATACCTTTCTGATTCAGTTCCTGGCGAATACGCGCAGGTCCATAACCTTTGCGGCTACGGCTGGCGATAAAGCGCGCAACAAATCGGCTGTCATCGAGATAGCCATGTTCATGGCACCAGGCaa                                                                                                                                  <  1:167332/149‑1 (MQ=255)
                tttttCTGTCGCTTCGCGGGAAATACCTTTCTGATTCAGTTCCTGGCGAATACGCGCAGGTCCATAACCTTTGCGGCTACGGCTGGCGATAAAGCGCGCAACAAATCGGCTGTCATCGAGATAGCCATGTTCATGGCACCAGGCAATaa                                                                                                                               <  2:383257/149‑1 (MQ=255)
                   ttCTGTCGCTTCGCGGGAAATACCTTTCTGATTCAGTTCCTGGCGAATACGCGCAGGTCCATAACCTTTGCGGCTACGGCTGGCGATAAAGCGCGCAACAAATCGGCTGTCATCGAGATAGCCATGTTCATGGCACCAGGCAATAAc                                                                                                                              <  2:156906/147‑1 (MQ=255)
                   ttCTGTCGCTTCGCGGGAAATACCTTTCTGATTCAGTTCCTGGCGAATACGCGCAGGTCCATAACCTTTGCGGCTACGGCTGGCGATAAAGCGCGCAACAAATCGGCTGTCATCGAGATAGCCATGTTCATGGCACCAGGCAATAAc                                                                                                                              >  1:156906/1‑147 (MQ=255)
                       gTCGCTTCGCGGGAAATACCTTTCTGATTCAGTTCCTGGCGAATACGCGCAGGTCCATAACCTTTGCGGCTACGGCTGGCGATAAAGCGCGCAACAAATCGGCTGTCATCGAGATAGCCATGTTCATGGCACCAGGCAATAACGCGCTc                                                                                                                        <  4:575511/149‑1 (MQ=255)
                        tCGCTTCGCGGGAAATACCTTTCTGATTCAGTTCCTGGCGAATACGCGCAGGTCCATAACCTTTGCGGCTACGGCTGGCGATAAAGCGCGCAACAAATCGGCTGTCATCGAGATAGCCATGTTCATGGCACCAGGCAATAACGCGCTCg                                                                                                                       >  1:281418/1‑149 (MQ=255)
                                cgGGAAATACCTTTCTGATTCAGTTCCTGGCGAATACGCGCAGGTCCATAACCTTTGCGGCTACGGCTGGCGATAAAGCGCGCAACAAATCGGCTGTCATCGAGATAGCCATGTTCATGCCACCAGGCAATAACGCGCTCGTAAtcttc                                                                                                               >  1:566518/1‑149 (MQ=255)
                                    aaaTACCTTTCTGATTCAGTTCCTGGCGAATACGCGCAGGTCCATAACCTTTGCGGCTACGGCTGGCGATAAAGCGCGCAACAAATCGGCTGTCATCGAGATAGCCATGTTCATGGCACCAGGCAATAACGCGCTCGTAATCTTCTGcc                                                                                                           <  3:350094/149‑1 (MQ=255)
                                       tACCTTTCTGATTCAGTTCCTGGCGAATACGCGCAGGTCCATAACCTTTGCGGCTACGGCTGGCGATAAAGCGCGCAACAAATCGGCTGTCATCGAGATAGCCATGTTCATGGCACCAGGCAATAACGCGCTCGTAATCTTCTGCCGTa                                                                                                        >  4:19313/1‑149 (MQ=255)
                                       tACCTTTCTGATTCAGTTCCTGGCGAATACGCGCAGGTCCATAACCTTTGCGGCTACGGCTGGCGATAAAGCGCGCAACAAATCGGCTGTCATCGAGATAGCCATGTTCATGGCACCAGGCAATAACGCGCTCGTAATCTTCTGCCGTa                                                                                                        <  2:353223/149‑1 (MQ=255)
                                       tACCTTTCTGATTCAGTTCCTGGCGAATACGCGCAGGTCCATAACCTTTGCGGCTACGGCTGGCGATAAAGCGCGCAACAAATCGGCTGTCATCGAGATAGCCATGTTCATGGCACCAGGCAATAACGCGCTCGTAATCTTCTGCCGTa                                                                                                        <  4:437760/149‑1 (MQ=255)
                                       tACCTTTCTGATTCAGTTCCTGGCGAATACGCGCAGGTCCATAACCTTTGCGGCTACGGCTGGCGATAAAGCGCGCAACAAATCGGCTGTCATCGAGATAGCCATGTTCATGGCACCAGGCAATAACGCGCTCGTAATCTTCTGCCGTa                                                                                                        >  2:349713/1‑149 (MQ=255)
                                           tttCTGATTCAGTTCCTGGCGAATACGCGCAGGTCCATAACCTTTGCGGCTACGGCTGGCGATAAAGCGCGCAACAAATCGGCTGTCATCGAGATAGCCATGTTCATGGCACCAGGCAATAACGCGCTCGTAATCTTCTGCCGTAGCAt                                                                                                    >  1:56286/1‑149 (MQ=255)
                                           tttCTGATTCAGTTCCTGGCGAATACGCGCAGGTCCATAACCTTTGCGGCTACAGCTGGCGATAAAGCGCGCAACAAATCGGCTGTCATCGAGATAGCCATGTTCATGGCACCAGGCAATAACGCGCTCGTAATCTTCTGCCGTAGCAt                                                                                                    <  4:423723/149‑1 (MQ=255)
                                               tGATTCAGTTCCTGGCGAATACGCGCAGGTCCATAACCTTTGCGGCTACGGCTGGCGATAAAGCGCGCAACAAATCGGCTGTCATCGAGATAGCCATGTTCATGGCACCAGGCAATAACGCGCTCGTAATCTTCTGCCGTAGCATCAAt                                                                                                >  2:169980/1‑149 (MQ=255)
                                                gATTCAGTTCCTGGCGAATACGCGCAGGTCCATAACCTTTGCGGCTACGGCTGGCGATAAAGCGCGCAACAAATCGGCTGTCATCGAGATAGCCATGTTCATGGCACCAGGCAATAACGCGCTCGTAATCTTCTGCCGTAGCATCAAtc                                                                                               >  1:267366/1‑149 (MQ=255)
                                                 aTTCAGTTCCTGGCGAATACGCGCAGGTCCATAACCTTTGCGGCTACGGCTGGCGATAAAGCGCGCAACAAATCGGCTGTCATCGAGATAGCCATGTTCATGGCACCAGGCAATAACGCGCTCGTAATCTTCTGCCGTAGCATCAAtct                                                                                              >  2:479311/1‑149 (MQ=255)
                                                   tCAGTTCCTGGCGAATACGCGCAGGTCCATAACCTTTGCGGCTACGGCTGGCGATAAAGCGCGCAACAAATCGGCTGTCATCGAGATAGCCATGTTCATGGCACCAGGCAATAACGCGCTCGTAATCTTCTGCCGTAGCATCAATctct                                                                                            >  1:487897/1‑149 (MQ=255)
                                                   tCAGTTCCTGGCGAATACGCGCAGGTCCATAACCTTTGCGGCTACGGCTGGCGATAAAGCGCGCAACAAATCGGCTGTCATCGAGATAGCCATGTTCATGGCACCAGGCAATAACGCGCTCGTAATCTTCTGCCGTAGCATCAATctct                                                                                            >  3:545247/1‑149 (MQ=255)
                                                    cAGTTCCTGGCGAATACGCGCAGGTCCATAACCTTTGCGGCTACGGCTGGCGATAAAGCGCGCAACAAATCGGCTGTCATCGAGATAGCCATGTTCATGGCACCAGGCAATAACGCGCTCGTAATCTTCTGCCGTAGCATCAATCtctt                                                                                           >  2:69625/1‑149 (MQ=255)
                                                           tGGCGAATACGCGCAGGTCCATAACCTTTGCGGCTACGGCTGGCGATAAAGCGCGCAACAAATCGGCTGTCATCGAGATAGCCCTGTTCATGGCACCAGGCAATAACGCGCTCGTAATCTTCTGCCGTAGCATCAATCTCTTCTGGGcc                                                                                    >  4:392018/1‑149 (MQ=255)
                                                               gAATACGCGCAGGTCCATAACCTTTGCGGCTACGGCTGGCGATAAAGCGCGCAACAAATCGGCTGTCATCGAGATAGCCATGTTCATGGCACCAGGCAATAACGCGCTCGTAATCTTCTGCCGTAGCATCAATCTCTTCTGGGCCAttt                                                                                >  1:11846/1‑149 (MQ=255)
                                                                aaTACGCGCAGGTCCATAACCTTTGCGGCTACGGCTGGCGATAAAGCGCGCAACAAATCGGCTGTCATCGAGATAGCCATGTTCATGGCACCAGGCAATAACGCGCTCGTAATCTTCTGCCGTAGCATCAATCTCTTCTGGGCCAtttt                                                                               >  2:489517/1‑149 (MQ=255)
                                                                 aTACGCGCAGGTCCATAACCTTTGCGGCTACGGCTGGCGATAAAGCGCGCAACAAATCGGCTGTCATCGAGATAGCCATGTTCATGGCACCAGGCAATAACGCGCTCGTAATCTTCTGCCGTAGCATCAATCTCTTCTGGGCCAttttt                                                                              >  4:257735/1‑149 (MQ=255)
                                                                       gcAGGTCCATAACCTTTGCGGCTACGGCTGGCGATAAAGCGCGCAACAAATCGGCTGTCATCGAGATAGCCATGTTCATGGCACCAGGTAATAACGCGCTCGTAATCTTCTGCCATAGCATCAATCTCTTCTGGGCCATTTTTGCCCAt                                                                        >  4:475737/1‑149 (MQ=255)
                                                                            tCCATAACCTTTGCGGCTACGGCTGGCGATAAAGCGCGCAACAAATCGGCTGTCATCGAGATAGCCATGTTCATGGCACCAGGCAATAACGCGCTCGTAATCTTCTGCCGTAGCATCAATCTCTTCTGGGCCATTTTTGCCCATAATCg                                                                   >  3:559987/1‑149 (MQ=255)
                                                                              cATAACCTTTGCGGCTACGGCTGGCGATAAAGCGCGCAACAAATCGGCTGTCATCGAGATAGCCATGTTCATGGCACCAGGCAATAACGCGCTCGTAATCTTCTGCCGTAGCATCAATCTCTTCTGGGCCATTTTTGCCCATAATCGGt                                                                 <  1:167644/149‑1 (MQ=255)
                                                                              cATAACCTTTGCGGCTACGGCTGGCGATAAAGCGCGCAACAAATCGGCTGTCATCGAGATAGCCATGTTCATGGCACCAGGCAATAACGCGCTCGTAATCTTCTGCCGTAGCATCAATCTCTTCTGGGCCATTTTTGCCCATAATCGGt                                                                 >  3:281851/1‑149 (MQ=255)
                                                                              cATAACCTTTGCGGCTACGGCTGGCGATAAAGCGCGCAACAAATCGGCTGTCATCGAGATAGCCATGTTCATGGCACCAGGCAATAACGCGCTCGTAATCTTCTGCCGTAGCATCAATCTCTTCTGGGCCATTTTTGCCCATAATCGGt                                                                 <  3:257735/149‑1 (MQ=255)
                                                                              cATAACCTTTGCGGCTACGGCTGGCGATAAAGCGCGCAACAAATCGGCTGTCATCGAGATAGCCATGTTCATGGCACCAGGCAATAACGCGCTCGTAATCTTCTGCCGTAGCATCAATCTCTTCTGGGCCATTTTTGCCCATAATCGGt                                                                 <  4:24063/149‑1 (MQ=255)
                                                                               aTAACCTTTGCGGCTACGGCTGGCGATAAAGCGCGCAACAAATCGGCTGTCATCGAGATAGCCATGTTCATGGCACCAGGCAATAACGCGCTCGTAATCTTCTGCCGTAGCATCAATCTCTTCTGGGCCATTTTTGCCCATAATCGGTg                                                                >  4:363352/1‑149 (MQ=255)
                                                                               aTAACCTTTGCGGCTACGGCTGGCGATAAAGCGCGCAACAAATCGGCTGTCATCGAGATAGCCATGTTCATGGCACCAGGCAATAACGCGCTCGTAATCTTCTGCCGTAGCATCAATCTCTTCTGGGCCATTTTTGCCCATAATCGGTg                                                                >  3:222629/1‑149 (MQ=255)
                                                                               aTAACCTTTGCGGCTACGGCTGGCGATAAAGCGCGCAACAAATCGGCTGTCATCGAGATAGCCATGTTCATGGCACCAGGCAATAACGCGCTCGTAATCTTCTGCCGTAGCATCAATCTCTTCTGGGCCATTTTTGCCCATAATCGGTg                                                                >  2:597312/1‑149 (MQ=255)
                                                                                tAACCTTTGCGGCTACGGCTGGCGATAAAGCGCGCAACAAATCGGCTGTCATCGAGATAGCCATGTTCATGGCACCAGGCAATAACGCGCTCGTAATCTTCTGCCGTAGCATCAATCTCTTCTGGGCCATTTTTGCCCATAATCGGTGc                                                               >  3:629839/1‑149 (MQ=255)
                                                                                     tttGCGGCTACGGCTGGCGATAAAGCGCGCAACAAATCGGCTGTCATCGAGATAGCCATGTTCATGGCACCAGGCAATAACGCGCTCGTAATCTTCTGCCGTAGCATCAATCTCTTCTGGGCCATTTTTGCCCATAATCGGTGCCGCGa                                                          <  4:528588/149‑1 (MQ=255)
                                                                                     tttGCGGCTACGGCTGGCGATAAAGCGCGCAACAAATCGGCTGTCATCGAGATAGCCATGTTCATGGCACCAGGCAATAACGCGCTCGTAATCTTCTGCCGTAGCATCAATCTCTTCTGGGCCATTTTTGCCCATAATCGGTGCCGCGa                                                          <  3:363352/149‑1 (MQ=255)
                                                                                      ttGCGGCTACGGCTGGCGATAAAGCGCGCAACAAATCGGCTGTCATCGAGATAGCCATGTTCATGGCACCAGGCAATAACGCGCTCGTAATCTTCTGCCGTAGCATCAATCTCTTCTGGGCCATTTTTGCCCATAATCGGTGCCGCGAg                                                         <  2:267366/149‑1 (MQ=255)
                                                                                          ggCTACGGCTGGCGATAAAGCGCGCAACAAATCGGCTGTCATCGAGATAGCCATGTTCATGGCACCAGGCAATAACGCGCTCGTAATCTTCTGCCGTAGCATCAATCTCTTCTGGGCCATTTTTGCCCATAATCGGTGCCGCGAGTTTa                                                     >  2:465713/1‑149 (MQ=255)
                                                                                            cTACGGCTGGCGATAAAGCGCGCAACAAATCGGCTGTCATCGAGATAGCCATGTTCATGGCACCAGGCAATAACGCGCTCGTAATCTTCTGCCGTAGCATCAATCTCTTCTGGGCCATTTTTGCCCATAATCGGTGCCGCGAGTTTACg                                                   >  3:250389/1‑149 (MQ=255)
                                                                                              aCGGCTGGCGATAAAGCGCGCAACAAATCGGCTGTCATCGAGATAGCCATGTTCATGGCACCAGGCAATAACGCGCTCGTAATCTTCTGCCGTAGCATCAATCTCTTCTGGGCCATTTTTGCCCATAATCGGTGCCGCGAGTTTACGTc                                                 >  4:531834/1‑149 (MQ=255)
                                                                                                    ggCGATAAAGCGCGCAACAAATCGGCTGTCATCGAGATAGCCATGTTCATGGCACCAGGCAATAACGCGCTCGTAATCTTCTGCCGTAGCATCAATCTCTTCTGGGCCATTTTTGCCCATAATCGGTGCCGCGAGTTTACGTCGCAGtt                                           <  1:169980/149‑1 (MQ=255)
                                                                                                     gCGATAAAGCGCGCAACAAATCGGCTGTCATCGAGATAGCCATGTTCATGGCACCAGGCAATAACGCGCTCGTAATCTTCTGCCGTAGCATCAATCTCTTCTGGGCCATTTTTGCCCATAATCGGTGCCGCGAGTTTACGTCGCAGTTc                                          <  3:125176/149‑1 (MQ=255)
                                                                                                       gATAAAGCGCGCAACAAATCGGCTGTCATCGAGATAGCCATGTTCATGGCACCAGGCAATAACGCGCTCGTAATCTTCTGCCGTAGCATCAATCTCTTCTGGGCCATTTTTGCCCATAATCGGTGCCGCGAGTTTACGTCGCAGTTCtt                                        >  3:679410/1‑149 (MQ=255)
                                                                                                        aTAAAGCGCGCAACAAATCGGCTGTCATCGAGATAGCCATGTTCATGGCACCAGGCAATAACGCGCTCGTAATCTTCTGCCGTAGCATCAATCTCTTCTGGGCCATTTTTGCCCATAATCGGTGCCGCGAGTTTACGTCGCAGTTCTTg                                       <  2:56286/149‑1 (MQ=255)
                                                                                                           aaGCGCGCAACAAATCGGCTGTCATCGAGATAGCCATGTTCATGGCACCAGGCAATAACGCGCTCGTAATCTTCTGCCGTAGCATCAATCTCTTCTGGGCCATTTTTGCCCATAATCGGTGCCGCGAGTTTACGTCGCAGTTCTTGCTc                                    >  3:197745/1‑149 (MQ=255)
                                                                                                           aaGCGCGCAACAAATCGGCTGTCATCGAGATAGCCATGTTCATGGCACCAGGCAATAACGCGCTCGTAATCTTCTGCCGTAGCATCAATCTCTTCTGGGCCATTTTTGCCCATAATCGGTGCCGCGAGTTTACGTCGCAGTTCTTGCTc                                    <  2:165725/149‑1 (MQ=255)
                                                                                                                  caacaaATCGGCTGTCGTCGAGATAGCCATGTTCATGGCACCAGGCTATAACGCGATCGTAATCGTCTGCCGTAGCATCAATCTCTTCTGGGCCATTTTTGCCCATAATCGGTGCCGCGAGTTTACGTCGCAGTTCTTGCTCACTGTGa                             <  4:413140/149‑1 (MQ=255)
                                                                                                                  caacaaATCGGCTGTCATCGAGATAGCCATGTTCATGGCACCAGGCAATAACGCGCTCGTAATCTTCTGCCGTAGCATCAATCTCTTCTGGGCCATTTTTGCCCATAATCGGTGCCGCGAGTTTACGTCGCAGt                                            <  1:86860/134‑1 (MQ=255)
                                                                                                                  caacaaATCGGCTGTCATCGAGATAGCCATGTTCATGGCACCAGGCAATAACGCGCTCGTAATCTTCTGCCGTAGCATCAATCTCTTCTGGGCCATTTTTGCCCATAATCGGTGCCGCGAGTTTACGTCGCAGt                                            >  2:86860/1‑134 (MQ=255)
                                                                                                                   aacaaATCGGCTGTCATCGAGATAGCCATGTTCATGGCACCAGGCAATAACGCGCTCGTAATCTTCTGCCGTAGCATCAATCTCTTCTGGGCCATTTTTGCCCATAATCGGTGCCGCGAGTTTACGTCGCAGTTCTTGCTCACTGTGAt                            >  4:604078/1‑149 (MQ=255)
                                                                                                                      aaaTCGGCTGTCATCGAGATAGCCATGTTCATGGCACCAGGCAATAACGCGCTCGTAATCTTCTGCCGTAGCATCAATCTCTTCTGGGCCATTTTTGCCCATAATCGGTGCCGCGAGtt                                                       >  1:624335/1‑119 (MQ=255)
                                                                                                                      aaaTCGGCTGTCATCGAGATAGCCATGTTCATGGCACCAGGCAATAACGCGCTCGTAATCTTCTGCCGTAGCATCAATCTCTTCTGGGCCATTTTTGCCCATAATCGGTGCCGCGAGtt                                                       <  2:624335/119‑1 (MQ=255)
                                                                                                                      aaaTCGGCTGTCATCGAGATAGCCATGTTCATGGCACCAGGCAATAACGCGCTCGTAATCTTCTGCCGTAGCATCAATCTCTTCTGGGCCATTTTTGCCCATAATCGGTGCCGCGAGTTTACGTCGCAGTTCTTGCTCACTGTGATcgc                         <  2:562507/149‑1 (MQ=255)
                                                                                                                      aaaTCGGCTGTCATCGAGATAGCCATGTTCATGGCACCAGGCAATAACGCGCTCGTAATCTTCTGCCGTAGCATCAATCTCTTCTGGGCCATTTTTGCCCATAATCGGTGCCGCGAGTTTACGTCGCAGTTCTTGCTCACTGTGATcgc                         <  1:69625/149‑1 (MQ=255)
                                                                                                                         tCGGCTGTCATCGAGATAGCCATGTTCATGGCACCAGGCAATAACGCGCTCGTAATCTTCTGCCGTAGCATCAATCTCTTCTGGGCCATTTTTGCCCATAATCGGTGCCGCGAGTTTACGTCGCAGTTCTTGCTCACTGTGATCGCGCa                      <  1:330807/149‑1 (MQ=255)
                                                                                                                           ggCTGTCATCGAGATAGCCATGTTCATGGCACCAGGCAATAACGCGCTCGTAATCTTCTGCCGTAGCATCAATCTCTTCTGGGCCATTTTTGCCCATAATCGGTGCCGCGAGTTTACGTCGCAGTTCTTGCTCACTGTGATCGCGCAcc                    <  2:568088/149‑1 (MQ=255)
                                                                                                                              tGTCATCGAGATAGCCATGTTCATGGCACCAGGCGATAACGCGCTCGTAATCTTCTGCCGTAGCATCAATCTCTTCTGGGCCATTTTTGCCCATAATCGGTGCCGCGAGTTTACGTCGCAGTTCTTGCTCACTGTGATCGCGCACCGcc                 <  1:479311/149‑1 (MQ=255)
                                                                                                                              tGTCATCGAGATAGCCATGTTCATGGCACCAGGCAATAACGCGCTCGTAATCTTCTGCCGTAGCATCAATCTCTTCTGGGCCATTTTTGCCCATAATCGGTGCCGCGAGTTTACGTCGCAGTTCTTGCTCACTGTGATCGCGCACCGcc                 <  4:545247/149‑1 (MQ=255)
                                                                                                                               gTCATCGAGATAGCCATGTTCATGGCACCAGGCAATAACGCGCTCGTAATCTTCTGCCGTAGCATCAATCTCTTCTGGGCCATTTTTGCCCATAATCGGTGCCGCGAGTTTACGTCGCAGTTCTTGCTCACTGTGATCGCGCACCGCCa                <  2:281418/149‑1 (MQ=255)
                                                                                                                                   tCGAGATAGCCATGTTCATGGCACCAGGCAATAACGCGCTCGTAATCTTCTGCCGTAGCATCAATCTCTTCTGGGCCATTTTTGCCCATAATCGGTGCCGCGAGTTTACGTCGCAGTTCTTGCTCACTGTGATCGCGCACCGCCAGAAt            >  3:211427/1‑149 (MQ=255)
                                                                                                                                   tCGAGATAGCCATGTTCATGGCACCAGGCAATAACGCGCTCGTAATCTTCTGCCGTAGCATCAATCTCTTCTGGGCCATTTTTGCCCATAATCGGTGCCGCGAGTTTACGTCGCAGTTCTTGCTCACTGTGATCGCGCACCGCCAGAAt            >  3:7001/1‑149 (MQ=255)
                                                                                                                                     gagaTAGCCATGTTCATGGCACCAGGCAATAACGCGCTCGTAATCTTCTGCCGTAGCATCAATCTCTTCTGGGCCATTTTTGCCCATAATCGGTGCCGCGAGTTTACGTCGCAGTTCTTGCTCACTGTGATCGCGCACCGCCAGAATGc          <  4:222629/149‑1 (MQ=255)
                                                                                                                                       tggtgCCATGAACATGGCACCAGGCAATAACGCGCTCGTAATCTTCTGCCGTAGCATCAATCTCTTCTGGGCCATTTTTGCCCATAATCGGTGCCGCGAGTTTACGTCGCAGTTCTTGCTCACtgt                               <  1:623795/122‑1 (MQ=255)
                                                                                                                                       tggtgCCATGAACATGGCACCAGGCAATAACGCGCTCGTAATCTTCTGCCGTAGCATCAATCTCTTCTGGGCCATTTTTGCCCATAATCGGTGCCGCGAGTTTACGTCGCAGTTCTTGCTCACtgt                               >  2:623795/5‑126 (MQ=255)
                                                                                                                                       tggtgCCATGAACATGGCACCAGGCAATAACGCGCTCGTAATCTTCTGCCGTAGCATCAATCTCTTCTGGGCCATTTTTGCCCATAATCGGTGCCGCGAGTTTACGTCGCAGTTCTTGCTCACTGTGATCGCGCACCGCCAGAATGCGt        >  1:577900/5‑149 (MQ=255)
                                                                                                                                       tggtgCCATGAACATGGCACCAGGCAATAACGCGCTCGTAATCTTCTGCCGTAGCATCAATCTCTTCTGGGCCATTTTTGCCCATAATCGGTGCCGCGAGTTTACGTCGCAGTTCTTGCTCACTGTGATCGCGCACCGCCAGAATGCGt        <  2:577900/145‑1 (MQ=255)
                                                                                                                                       gaTAGCCATGTTCATGGCACCAGGCAATAACGCGCTCGTAATCTTCTGCCGTAGCATCAATCTCTTCTGGGCCATTTTTGCCCATAATCGGTGCCGCGAGTTTACGTCGCAGTTCTTGCTCACTGTGATCGCGCACCGCCAGAATGCGt        <  4:211427/149‑1 (MQ=255)
                                                                                                                                           gCCATGTTCATGGCACCAGGCAATAACGCGCTCGTAATCTTCTGCCGTAGCATCAATCTCTTCTGGGCCATTTTTGCCCATAATCGGTGCCGCGAGTTTACGTCGCAGTTCTTGCTCACTGTGATCGCGCACCGCCAGAATGCGTACCg    <  3:604078/149‑1 (MQ=255)
                                                                                                                                            ccATGTTCATGGCACCAGGCAATAACGCGCTCGTAATCTTCTGCCGTAGCATCAATCTCTTCTGGGCCATTTTTGCCCATAATCGGTGCCGCGAGTTTACGTCGCAGTTCTTGCTCACTGTGATCGCGCACCGCCAGAATGCGTACCGc   >  3:573063/1‑149 (MQ=255)
                                                                                                                                             cATGTTCATGGCACCAGGCAATAACGCGCTCGTAATCTTCTGCCGTAGCATCAATCTCTTCTGGGCCATTTTTGCCCATAATCGGTGCCGCGAGTTTACGTCGCAGTTCTTGCTCACTGTGATCGCGCACCGCCAGAATGCGTACCGCa  >  3:117047/1‑149 (MQ=255)
                                                                                                                                                 |                                                                                                                                                
CACATTCACGCATCGCTTTTTCTGTCGCTTCGCGGGAAATACCTTTCTGATTCAGTTCCTGGCGAATACGCGCAGGTCCATAACCTTTGCGGCTACGGCTGGCGATAAAGCGCGCAACAAATCGGCTGTCATCGAGATAGCCATGTTCATGGCACCAGGCAATAACGCGCTCGTAATCTTCTGCCGTAGCATCAATCTCTTCTGGGCCATTTTTGCCCATAATCGGTGCCGCGAGTTTACGTCGCAGTTCTTGCTCACTGTGATCGCGCACCGCCAGAATGCGTACCGCA  >  NC_000913/2822300‑2822589

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