Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,592,095 A→T 5.1% K147* (AAA→TAA)  tsr → methyl‑accepting chemotaxis protein I, serine sensor receptor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,592,0950AT 5.1% 259.8 / 3.5 79K147* (AAA→TAA) tsrmethyl‑accepting chemotaxis protein I, serine sensor receptor
Reads supporting (aligned to +/- strand):  ref base A (35/40);  new base T (2/2);  total (37/42)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.08e-01

TGAAACAGGCGGAAAAAAACTGGGCGGATTACGAAGCGTTGCCGCGTGACCCGCGTCAGAGCACCGCCGCAGCGGCAGAGATCAAACGTAATTACGATATTTATCACAATGCGCTGGCGGAGCTGATCCAACTGTTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCCATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCT  >  NC_000913/4591949‑4592243
                                                                                                                                                  |                                                                                                                                                    
tGAAACAGGCGGAAAAAAACTGGGCGGATTACGAAGCGTTGCCGCGTGACCCGCGTCAGAGCACCGCCGCAGCGGCAGAGATCAAACGTAATTACGATATTTATCACAATGCGCTGGCGGAGCTGATCCAACTGTTAGGTGCAGGCaa                                                                                                                                                     >  7:78443/1‑148 (MQ=255)
tGAAACAGGCGGAAAAAAACTGGGCGGATTACGAAGCGTTGCCGCGTGACCCGCGTCAGAGCACCGCCGCAGCGGCAGAGATCAAACGTAATTACGATATTTATCACAATGCGCTGGCGGAGCTGATCCAACTGTTAGGTGCAGGCaa                                                                                                                                                     >  5:158835/1‑148 (MQ=255)
     cAGGCGGAAAAAAACTGGGCGGATTACGAAGCGTTGCCGCGTGACCCGCGTCAGAGCACCGCCGCAGCGGCAGAGATCAAACGTAATTACGATATTTATCACAATGCGCTGGCGGAGCTGATCCAACTGTTAGGTGCAGGCAAAATCaa                                                                                                                                               <  1:25217/149‑1 (MQ=255)
      aGGCGGAAAAAAACTGGGCGGATTACGAAGCGTTGCCGCGTGACCCGCGTCAGAGCACCGCCGCAGCGGCAGAGATCAAACGTAATTACGATATTTATCACAATGCGCTGGCGGAGCTGATCCAACTGTTAGGTGCAGGCAAAATCAAc                                                                                                                                              <  3:147526/149‑1 (MQ=255)
        gCGGAAAAAAACTGGGCGGATTACGAAGCGTTGCCGCGTGACCCGCGTCAGAGCACCGCCGCAGCGGCAGAGATCAAACGTAATTACGATATTTATCACAATGCGCTGGCGGAGCTGATCCAACTGTTAGGTGCAGGCAAAATCAACGa                                                                                                                                            <  2:126284/149‑1 (MQ=255)
          ggAAAAAAACTGGGCGGATTACGAAGCGTTGCCGCGTGACCCGCGTCAGAGCACCGCCGCAGCGGCAGAGATCACACGTAATTACGATATTTATCACAATGCGCTGGCGGAGCTGATCCACCTGATAGGTGCAGGCAAAATCAACTAGt                                                                                                                                          >  6:135746/1‑149 (MQ=255)
          acaaaaaaaCTGGGCGGATAACGAAGCGTTGCCGCGTGACCCGCGTCAGAGCACCGCCGCCGCCGCAGACAACAAACGAAAATACGATATTTATCACCATGCGCTGCCCGAGAAGATCCAACTCTTAGGTCCAGCCAACATCAACGAgg                                                                                                                                          >  2:351981/3‑148 (MQ=255)
            aaaaaaaCTGGGCGGATTACGAAGCGTTGCCGCGTGACCCGCGTCAGAGCACCGCCGCAGCGGCAGAGATCAAACGTAATTACGATATTTATCACAATGCGCTGGCGGAGCTGATCCAACTGTTAGGTGCAGGCAAAATCAACGAGTTc                                                                                                                                        >  3:16701/1‑149 (MQ=255)
             aaaaaaCTGGGCGGATTACGAAGCGTTGCCGCGTGACCCGCGTCAGAGCACCGCCGCAGCGGCAGAGATCAAACGTAATTACGATATTTATCACAATGCGCTGGCGGAGCTGATCCAACTGTTAGGTGCAGGCAAAATCAACGAGTTCt                                                                                                                                       <  5:282410/149‑1 (MQ=255)
             aaaaaaCTGGGCGGATTACGAAGCGTTGCCGCGTGACCCGCGTCAGAGCACCGCCGCAGCGGCAGAGATCAAACGTAATTACGATATTTATCACAATGCGCTGGCGGAGCTGATCCAACTGTTAGGTGCAGGCAAAATCAACGAGTTCt                                                                                                                                       <  4:250749/149‑1 (MQ=255)
             aaaaaaCTGGGCGGATTACGAAGCGTTGCCGCGTGACCCGCGTCAGAGCACCGCCGCAGCGGCAGAGATCAAACGTAATTACGATATTTATCACAATGCGCTGGCGGAGCTGATCCAACTGTTAGGTGCAGGCAAAATCAACGAGTTCt                                                                                                                                       >  4:383077/1‑149 (MQ=255)
               aaaaCTGGGCGGATTACGAAGCGTTGCCGCGTGACCCGCGTCAGAGCACCGCCGCAGCGGCAGAGATCAAACGTAATTACGATATTTATCCCAATGCGCTGGCGGAGCTGATCCAACTGTTACGTGCAGGCAAAATCAACGAGTTCttt                                                                                                                                     >  5:36175/1‑149 (MQ=255)
               aaaaCTGGGCGGATTACGAAGCGTTGCCGCGTGACCCGCGTCAGAGCACCGCCGCAGCGGCAGAGATCAAACGTAATTACGATATTTATCACAATGCGCTGGCGGAGCTGATCCAACTGTTAGGTGCAGGCAAAATCAACGAGTTCttt                                                                                                                                     <  3:194182/149‑1 (MQ=255)
                 aaCTGGGCGGATTACGAAGCGTTGCCGCGTGACCCGCGTCAGAGCACCGCCGCAGCGGCAGAGATCAAACGTAATTACGATATTTATCACAATGCGCTGGCGGAGCTGATCCAACTGTTAGGTGCAGGCAAAATCAACGAGTTCTTTGa                                                                                                                                   <  8:66822/149‑1 (MQ=255)
                                gAAGCGTTGCCGCGTGACCCGCGTCAGAGCACCGCCGCAGCGGCAGAGATCAAACGTAATTACGATATTTATCACAATGCGCTGGCGGAGCTGATCCAACTGTTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAgg                                                                                                                     >  8:286380/1‑148 (MQ=255)
                                 aaGCGTTGCCGCGTGACCCGCGTCAGAGCACCGCCGCAGCGGCAGAGATCAAACGTAATTACGATATTTATCACAATGCGCTGGCGGAGCTGATCCAACTGTTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGa                                                                                                                   <  5:250961/149‑1 (MQ=255)
                                      ttGCCGCGTGACCCGCGTCAGAGCACCGCCGCAGCGGCAGAGATCAAACGTAATTACGATATTTATCACAATGCGCTGGCGGAGCTGATCCAACTGTTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCa                                                                                                              <  1:141665/149‑1 (MQ=255)
                                      ttGCCGCGTGACCCGCGTCAGAGCACCGCCGCAGCGGCAGAGATCAAACGTAATTACGATATTTATCACAATGCGCTGGCGGAGCTGATCCAACTGTTAGGTGCAGGCAAAATCAACGAGTGCTTTGATCAGCCGACCCAGGGatat                                                                                                                >  5:326170/1‑147 (MQ=255)
                                              tGACCCGCGTCAGAGCACCGCCGCAGCGGCAGAGATCAAACGTAATTACGATATTTATCACAATGCGCTGGCGGAGCTGATCCAACTGTTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGtt                                                                                                      <  8:78443/149‑1 (MQ=255)
                                               gaaccGCGTCAGAGCACCGCCGCAGCGGCAGAGATCAAACGTAATTACGATATTTATCACAATGCGCTTGCGGAGCTGTTCCAACTGTTAGGTGCTGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGatt                                                                                                     <  6:384741/146‑1 (MQ=255)
                                                aCCCGCGTCAGAGCCCCGGCGCAGCGGCAGAGATCAAACGTACTAACGATAATTATCACAATGCGCTGGCGGAGCTGATCCAACTGTTAGGTGCCGGCAAAATCAACGAGTTCTTTGATCAGCCGACCAAGGGATATCAGGACGGttt                                                                                                     <  8:149485/148‑1 (MQ=255)
                                                  ccGCGTCAGAGCACCGCCGCAGCGGCAGAGATCAAACGTAATTACGATATTTATCACAATGCGCTGGCGGAGCTGATCCAACTGTTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGCCCCAGGGATAGCAGGACGGTTTc                                                                                                    >  2:375571/1‑147 (MQ=255)
                                                  ccGCGTCAGAGCACCGCCGCAGCGGCAGAGATCAAACGTAATTACGATATTTATCACAATGCGCTGGCGGAGCTGATCCAACTGTTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTc                                                                                                    >  8:376076/1‑147 (MQ=255)
                                                    gcgTCAGAGCACCGCCGCAGCGGCAGAGATCAAACGTAATTACGATATTTATCACAATGCGCTGGCGGAGCTGATCCAACTGTTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCgaga                                                                                                <  3:158751/149‑1 (MQ=255)
                                                    gcgTCAGAGCACCGCCGCAGCGGCAGAGATCAAACGTAATTACGATATTTATAACAATGCGCTGGCGGAGCTGATCCAACTGTTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGCCCCAGGGATCTCAGGACGGTTTCgaga                                                                                                >  6:341738/1‑149 (MQ=255)
                                                      gTCAGAGCACCGCCGCAGCGGCAGAGATCAAACGTAATTACGATATTTATCACAATGCGCTGGCGGAGCTGATCCAACTGTTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAg                                                                                              <  1:82841/149‑1 (MQ=255)
                                                         agagCACCGCCGCAGCGGCAGAGATCAAACGTAATTACGATATTTATCACAATGCGCTGGCGGAGCTGATCCAACTGTTAGGTGCAGGCTAAATCaa                                                                                                                                               >  2:409824/1‑97 (MQ=255)
                                                         agagCACCGCCGCAGCGGCAGAGATCAAACGTAATTACGATATTTATCACAATGCGCTGGCGGAGCTGATCCAACTGTTAGGTGCAGGCTAAATCaa                                                                                                                                               >  2:415710/1‑97 (MQ=255)
                                                         agagCACCGCCGCAGCGGCAGAGATCAAACGTAATTACGATATTTATCACAATGCGCTGGCGGAGCTGATCCAACTGTTAGGTGCAGGCTAAATCaa                                                                                                                                               <  1:409824/97‑1 (MQ=255)
                                                         agagCACCGCCGCAGCGGCAGAGATCAAACGTAATTACGATATTTATCACAATGCGCTGGCGGAGCTGATCCAACTGTTAGGTGCAGGCTAAATCaa                                                                                                                                               <  1:415710/97‑1 (MQ=255)
                                                         agagCACCGCCGCAGCGGCAGAGATCAAACGTAATTACGATATTTATCACAATGCGCTGGCGGAGCTGATCCAACTGTTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGa                                                                                                                   >  1:11685/1‑125 (MQ=255)
                                                         agagCACCGCCGCAGCGGCAGAGATCAAACGTAATTACGATATTTATCACAATGCGCTGGCGGAGCTGATCCAACTGTTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGa                                                                                                                   <  2:11685/125‑1 (MQ=255)
                                                         agagCACCGCCGCAGCGGCAGAGATCAAACGTAATTACGATATTTATCACAATGCGCTGGCGGAGCTGATCCAACTGTTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCa                                                                                            <  5:234962/148‑1 (MQ=255)
                                                         agagCACCGCCGCAGCGGCAGAGATCAAACGTAATTACGATATTTATCACAATGCGCTGGCGGAGCTGATCCAACTGTTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAg                                                                                           <  7:376076/149‑1 (MQ=255)
                                                            gCACCGCCGCAGCGGCAGAGATCAAACGTAATTACGATATTTATCACAATGCGCTGGCGGAGCTGATCCAACTGTTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTAt                                                                                        >  6:94156/1‑149 (MQ=255)
                                                                 gccgcAGCGGCAGAGATCAAACGTAATTACGATATTTATCACAATGCGCTGGCGGAGCTGATCCAACTGTTAGGTGCTGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGc                                                                                   <  1:144833/149‑1 (MQ=255)
                                                                     cAGCGGCAGAGATCAAACGTAATTACGATATTTATCACAATGCGCTGGCGGAGCTGATCCAACTGTTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTAc                                                                               <  3:383077/149‑1 (MQ=255)
                                                                      aGCGGCAGAGATCAAACGTAATTACGATATTTATCACAATGCGCTGGCGGAGCTTATACAACTGTTAGGTGCATGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACTGTGTCGATAAGCAGTATGTGGCTTACa                                                                              <  6:22462/149‑1 (MQ=255)
                                                                      aGCGGCAGAGATCAAACGTAATTACGATATTTATCACAATGCGCTGGCGGAGCTGATCCAACTGTTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACa                                                                              <  4:131310/149‑1 (MQ=255)
                                                                            agagATCAAACGTAATTACGATATTTATCACAATGCGCTGGCGGAGCTGATCCAACTGTTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAgg                                                                                                                     <  4:125563/104‑1 (MQ=255)
                                                                            agagATCAAACGTAATTACGATATTTATCACAATGCGCTGGCGGAGCTGATCCAACTGTTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAgg                                                                                                                     >  3:125563/1‑104 (MQ=255)
                                                                                aTCATACGTTATTACGATATTTATCACAATGCGCTGGCGGAGCTGATCCAACTGTTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCaaaa                                                                    <  6:36175/149‑1 (MQ=255)
                                                                                 tCAAACGTAATTACGATATTTATCACAATGCGCTGGCGGAGCTGATCCAACTGTTAGGTGCAGGCAAAATCAACGAGTACTTTGATCAGCCGACCCAGGGATATCTGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAAc                                                                   >  1:98163/1‑149 (MQ=255)
                                                                                  cAAACGTAATTACGATATTTATCACAATGCGCTGGCGGAGCTGATCCAACTGATAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGACGCAGTATGTGGCTTACATGGAGCAAAACg                                                                  <  2:337250/149‑1 (MQ=255)
                                                                                   aaaCGTAATTACGATATTTATCACAATGCGCTGGCGGAGCTGATCCAACTGTTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGa                                                                 >  2:379027/1‑149 (MQ=255)
                                                                                       gTAATTACGATATTTATCACAATGCGCTGGCGGAGCTGATCCAACTGTTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCa                                                                                                              <  8:340760/100‑1 (MQ=255)
                                                                                       gTAATTACGATATTTATCACAATGCGCTGGCGGAGCTGATCCAACTGTTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCa                                                                                                              >  7:340760/1‑100 (MQ=255)
                                                                                            tACGATATTTATCACAATGCGCTGGCGGAGCTGATCCAACTGTTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGa                                                                 >  4:64553/1‑140 (MQ=255)
                                                                                            tACGATATTTATCACAATGCGCTGGCGGAGCTGATCCAACTGTTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGa                                                                 <  3:64553/140‑1 (MQ=255)
                                                                                                atatTTATCACAATGCGCTGGCGGAGCTGATCCAACTGTTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCTGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCCATGat                                                    <  4:16701/149‑1 (MQ=255)
                                                                                                 tatTTATCACAATGCGCTGGCGGAGCTGATCCAACTGTTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGAt                                                                <  5:298416/136‑1 (MQ=255)
                                                                                                 tatTTATCACAATGCGCTGGCGGAGCTGATCCAACTGTTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGAt                                                                >  6:298416/1‑136 (MQ=255)
                                                                                                       tcgcaATGCGCTGGCGGAGCTGATCCAACTGTTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCCATGATATCGCCg                                             >  8:407212/4‑149 (MQ=255)
                                                                                                        cacaATGCGCTGGCGGAGCTGATCCAACTGTTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCCATGATATCGCCGt                                            >  5:288797/1‑149 (MQ=255)
                                                                                                           aaTGCGCTGGCGGAGCTGATCCAACTGTTAGGTGCAGGCAAAATCAACGAGTTCTGTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCCATGATATCGCCGTCAg                                         >  1:143957/1‑149 (MQ=255)
                                                                                                            aTGCGCTGGCGGAGCTGATCCAACTGTTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCCATGATATCGCCGTCAGc                                        <  8:356306/149‑1 (MQ=255)
                                                                                                               cgcTGGCGGAGCTGATCCAACTGGTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCCGCCGACCCAGGGATATCAGGCCGGTTTCGAGACGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCCATGATATCGCCGTCAGCGa                                      <  8:389186/148‑1 (MQ=255)
                                                                                                                gcTGGCGGAGCTGATCCAACTGTTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTc                                                                                                    <  4:84824/85‑1 (MQ=255)
                                                                                                                gcTGGCGGAGCTGATCCAACTGTTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTc                                                                                                    >  3:84824/1‑85 (MQ=255)
                                                                                                                gcTGGCGGAGCTGATCCAACTGTTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCCATGATATCGCCGTCAGCGAt                                     >  7:368432/1‑148 (MQ=255)
                                                                                                                    gcggagcggaTCCAACTGTTAGGGGCAGGCAAAATCAACGAGTGCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCCATGATATCGCCGTCAGCGATAACaa                                <  8:69075/149‑1 (MQ=255)
                                                                                                                    gCGGAGCTGATCCAACTGTTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCCATGATATCGCCGTCAGCGATAACaa                                <  7:286380/149‑1 (MQ=255)
                                                                                                                     cGGAGCTGATCCAACTGTTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCCATGATATCGCCGTCAGCGATAACAAt                               <  6:158835/149‑1 (MQ=255)
                                                                                                                     cGGAGCTGATCCAACTGTTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCCATGATATCGCCGTCAGCGATAACAAt                               <  5:94156/149‑1 (MQ=255)
                                                                                                                       gAGCTGATCCAACTGTTAGGTGCAGGCAAAATCAACGAGTTCTTTAATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCCATGATATCGCCGTCAGCGATAACAATGc                             <  8:264769/149‑1 (MQ=255)
                                                                                                                           tGATCCAACTGTTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCCATGATATCGCCGTCAGCGATAACAATGcctcc                         >  3:156302/1‑149 (MQ=255)
                                                                                                                            gATCCAACTGTTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCCATGATATCGCCGTCAGCGATAACAATGcctcc                         >  5:74313/1‑148 (MQ=255)
                                                                                                                             aTCCAACTGTTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCCATGATATCGCCGTCAGCGATAACAATGcctcct                        >  2:407742/1‑148 (MQ=255)
                                                                                                                              tCCAACTGTTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCCATGATATCGCCGTCAGCGATAACAATGCCTCCTa                       >  7:235331/1‑148 (MQ=255)
                                                                                                                              tCCAACTGTTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCCATGATATCGCCGTCAGCGATAACAATGCCTCCTa                       >  2:252535/1‑148 (MQ=255)
                                                                                                                               ccAACTGTTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCCATGATATCGCCGTCAGCGATAACAATGCCTCCTa                       >  3:257435/1‑147 (MQ=255)
                                                                                                                                 aaCTGTTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCCATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGc                   >  8:371391/1‑149 (MQ=255)
                                                                                                                                      ttAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCCATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGc              >  5:275574/1‑149 (MQ=255)
                                                                                                                                         ggTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCCATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGAt           <  7:407212/149‑1 (MQ=255)
                                                                                                                                          gTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTAc                                                                               <  6:196772/80‑1 (MQ=255)
                                                                                                                                          gTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTAc                                                                               >  5:196772/1‑80 (MQ=255)
                                                                                                                                          gTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCCATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGc                   <  6:365285/140‑1 (MQ=255)
                                                                                                                                          gTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCCATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGc                   >  5:365285/1‑140 (MQ=255)
                                                                                                                                           tGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCCATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGAtgt         <  8:235331/149‑1 (MQ=255)
                                                                                                                                               ggCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCCATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTggct     >  1:105832/1‑147 (MQ=255)
                                                                                                                                                gCAGAAACAACGAGTACTTTGCTCCGCCGACCCAGGGATATCAGGACGGTTTCGAGAATCAGTATGTGGCTTACATGTAGCAAAACGATAGGCTCCATGATATCGACTTCAGCGATAACAATGACTCCTACAGCCAGGAGATGTGGAtt    >  6:104136/1‑149 (MQ=255)
                                                                                                                                                  aaaaTCAACGAGTTCTTTGATCAGCCGCCCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCCATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCt  >  8:214027/1‑149 (MQ=255)
                                                                                                                                                  aaaaTCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCCATGATTTCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCt  >  8:277817/1‑149 (MQ=255)
                                                                                                                                                  aaaaTCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCCATGATATCGCCGTCAGCGATa                                    >  8:131832/1‑115 (MQ=255)
                                                                                                                                                  aaaaTCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCCATGATATCGCCGTCAGCGATa                                    <  7:131832/115‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
TGAAACAGGCGGAAAAAAACTGGGCGGATTACGAAGCGTTGCCGCGTGACCCGCGTCAGAGCACCGCCGCAGCGGCAGAGATCAAACGTAATTACGATATTTATCACAATGCGCTGGCGGAGCTGATCCAACTGTTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCCATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCT  >  NC_000913/4591949‑4592243

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