Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,592,188 C→T H178Y (CAT→TAT)  tsr → methyl‑accepting chemotaxis protein Tsr

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,592,1880CT97.7% 129.6 / NA 44H178Y (CAT→TAT) tsrmethyl‑accepting chemotaxis protein Tsr
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (19/23);  minor base G (0/1);  total (20/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.73e-01

AACTGGGCGGATTACGAAGCGTTGCCGCGTGACCCGCGTCAGAGCACCGCCGCAGCGGCAGAGATCAAACGTAATTACGATATTTATCACAATGCGCTGGCGGAGCTGATCCAACTGTTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCCATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTTCGCCGTCTGGTTCGGTATTAAAGCCTCGCTGGTAGCGCCAATGAATCGCCTGATTGACAGCATTCGTCATATTGCAGGCGGCGATCTGGTGAAACCGATTGAGGTGGATGGCTCTAATGAGATGGGGCAACTGGCAGAGAGTTTGCGCCATATG  >  NC_000913/4591966‑4592433
                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                      
aaCTGGGCAGATTACGAAGCGTTGCCGCGTGACCCGCGTCAGAGCACCGCCGCAGCGGCAGAGTTCAAACGTAATTACGATATTTATCACAATGCGCTGGCGGAGCTGATCCAACTATTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCTATGATATCGCCGTCAGCGATAACAATGc                                                                                                                                                                                                                            >  1:140855/1‑251 (MQ=255)
 aCTGGGCAGATTACGAAGCGTTGCCGCGTGACCCGCGTCAGAGCACCGACGCAGAGGCAGAGATCAAACGTAATTACGATATTTATCACAATGAGCCGGCGGAGCTGATCCAACTATAAGGTGCAGACAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCTATGATATCGCCGTCAGCGATAACAATGcc                                                                                                                                                                                                                           >  1:465826/1‑251 (MQ=255)
            tACGAAGCGTTGCCGCGTGACCCGCGTCAGAGCACCGCCGCAGCGGCAGAGATCAAACGTAATTACGATATTTATCACAATGCGCTGGCGGAGCTGATCCAACTATTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCTATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGcc                                                                                                                                                                                                                 <  1:453242/250‑1 (MQ=255)
               gAAGCGTTGCCGCGTGACCCGCGTCAGAGCACCGCCGCAGCGGCAGAGATCAAACGTAATTACGATATTTATCACAATGCGCTGGCGGAGCTGATCCAACTATTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCTATGATATCACCGTCAGCGATAACAATGCCTCCTACAGCCAGGc                                                                                                                                                                                                             >  1:92280/1‑251 (MQ=255)
                   cGTTGCCGCGTGACCCGCGTCAGAGCACCGCCGCAGCGGCAGAGATCAAACGTAATTACGATATTTATCACAATGCGCTGGCGGAGCTGATCCAACTATTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCTATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCgttg                                                                                                                                                                                                         >  1:214662/1‑248 (MQ=255)
                    gtggCCGCGTGCCCCGCGTCAGAGCACCGCCGCAGCGGCAGAGATCAAACGTAATTACGATATTTATCACAATGCGCTGGCGGAGCTGATCCAACTATTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCTATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGAtg                                                                                                                                                                                                         <  1:469759/247‑1 (MQ=255)
                    gTTGCCGCGTGACCCGCGTCAGAGCACCGCCGCAGCGGCAGAGATCAAACGTAATTACGATATTTATCACAATGCGCTGGCGGAGCTGATCCAACTATTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCTATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGAtgt                                                                                                                                                                                                        >  1:102127/1‑251 (MQ=255)
                          gcgTGACCCGCGTCAGAGCACCGCCGCAGCGGCAGAGATCAAACGTAATTACGATATTTATCACAATGCGCTGGCGGAGCTGATCCAACTATTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCTATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTc                                                                                                                                                                                                  >  1:149299/1‑251 (MQ=255)
                             tGACCCGCGTCAGAGCACCGCCGCAGCGGCAGAGATCAAACGTAATTACGATATTTATCACAATGCGCTGGCGGAGCTGATCCAACTATTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCTATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCtgg                                                                                                                                                                                               >  1:405953/1‑251 (MQ=255)
                               ccccGCGTCAGAGCACCGCCGCAGCGGCAGAGATCAAACGTAATTACGATATTTATCACAATGCGCTGGCGGAGCTGATCCAACTATTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCTATGATATCGACGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCtggtg                                                                                                                                                                                             <  2:340304/250‑1 (MQ=255)
                                   gcgacagacaaCCGCCGCAGCGGCAGAGATCAAACGTAATTACGATATTTATCACAATGCGCTGGCGAAGCTGATCCAACTATTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCTATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCg                                                                                                                                                                                         <  2:70547/241‑1 (MQ=255)
                                            cACCGCCGCAGCGGCAGAGATCAAACGTAATTACGATATTTATCACAATGCGCTGGCGGAGCTGATCCAACTATTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCTATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGAtcg                                                                                                                                                                                >  1:79411/1‑251 (MQ=255)
                                             aCCGCCGCAGCGGCAGAGATCAAACGTAATTACGATATTTATCCCAATGCGCGGGCGGAGCTGATCCAACTATAAGGTGAAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCGGGGATATAAGGCGGGTTTCGAGAAGCAGTATGGGGCTTACATGGAGAAAAACGATGGGCTCAATGATATGGCAGTCAGCGATAACAATGCATCCTAAAGCCAGGCGATGTGGATCTTGGTGGGCgggagaatggg                                                                                                                                                                               >  1:127860/1‑241 (MQ=255)
                                                            gagaTCAAACGTAATTACGATATTTATCACAATGCGCTGGCGGAGCTGATCCAACTATTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCTATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCAt                                                                                                                                                                <  2:102127/251‑1 (MQ=255)
                                                                        aaTTACGATATTTATCACAATGCGCTGGCGGAGCTGATCCAACTATTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAGAACGATCGGCTCTATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTTCGCCGTCt                                                                                                                                                     <  1:18625/250‑1 (MQ=255)
                                                                         aTTACGATATTTATAACAATGCGCTGGCGGAGCTGATCCAACTATTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCTATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTTCTCTTTGCCGTCTgg                                                                                                                                                   >  1:155532/1‑251 (MQ=255)
                                                                                     aTCACAATGCGCTGGCGGAGCTGATCCAACTATTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCTATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTAATCTTCGCCGTcctggtcggtgttaa                                                                                                                                        >  1:278934/1‑236 (MQ=255)
                                                                                                  ggCGGAGCTGATCCAACTATTAGGTGCAGGCAAAATCAACGAGTTCTGTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCTATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTTCGCCGTCTGGTTCGGTATTAAAGCCTCGCTGGTAg                                                                                                                          <  1:162056/251‑1 (MQ=255)
                                                                                                                   tattAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCTATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATATTTGCCCGCTGGTTCGGTGTTTAAGCCTCCCTGGGAGCGCCAATGAATCGCCTg                                                                                                         >  1:143667/3‑251 (MQ=255)
                                                                                                                    attAGGTGCAGGCAAAATCAACGAGTTCTTGGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCTATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCAGCTTCGCCGTCTGGTTCGGTATTAAAGCCTCGCTGGTAGCGCCAATGAATCGCCt                                                                                                          <  2:212250/248‑1 (MQ=255)
                                                                                                                         gTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCTATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCTTGATGATCGTCGTACTGGCGGTCATCTTCGCCGTCTGGTTCGGTATTAAAGACCCGCTGGTAGCGCCACTGAATCGCCCGATTggc                                                                                                  >  1:415220/1‑249 (MQ=255)
                                                                                                                                       aacgatgtCTTTGTTCATCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATTTGGCTTACATGTAGCTATACGATCGGCTCGATGATATCGCCGTCAGCGATAACAATGCCTCCTGCAGCCCGGCGATGTGGATTCTTTTGGGCGTTACGAGCGTCGTACTGGCGTTCATCTTCGCCGTCTGGTACGGTATGAAAGCCTCTCTTTTAGCGCCAATGAATCGCCTGATTCACAGCATTcattatat                                                                                    <  2:212207/244‑8 (MQ=255)
                                                                                                                                             ttCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCTATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTTCGCCGTCTGGTTCGGTATTAAAGCCTCGCTGGGAGCGCCCATGAATCGCCTGATTTACAGCATTCGTCATATTGccgg                                                                              >  1:150093/1‑248 (MQ=255)
                                                                                                                                                    aTCAGCCGTCCCAGGGATATAAGGAAGTGTTCGAGAAGCAGTATGTGGCTTACAGGGAGCAAAACGAACGGCTCTATGATATCGCCGGCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGAGGATCGTCGTACTGGCGGTCATCTTCGCCGTCTGGTGCGGTATTAAAGCCTCGCTGGTAGCGCCAATGAATCGCCTGATTGACAGCATTCGTCATATTGCAGGCGGCGAt                                                                       <  2:214662/251‑1 (MQ=255)
                                                                                                                                                    aTCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCTATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTTCGCCGTCTGGGTCGGGATTAAAGCCTCGCTGGTAGCGCCCATGAATCGCCTGATTGACAGCATTCGTCATATTTCAGggggcggt                                                                       >  1:94777/1‑249 (MQ=255)
                                                                                                                                                          cGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCTATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTTCGCCGTCTGGTTCGGTATTAAAGCCTCGCTGGTAGCGCCAATGAATCGCCTGATTGACAGCATTCGTCATATTGCAGGCGGCGATCTGGTg                                                                 <  2:405953/251‑1 (MQ=255)
                                                                                                                                                               cAGGGATATCAGGTCGGTTTCGAGAAGCAGTATGTGGCTTACAAGGAGCAAAACGATCGGCTCTATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCAGCTTCGCCGTCTGGTTCGGTATTAAAGCCTCGCTGGTAGCGCCAATGAATCGCCTGATTGACAGCATTCGTCATATTGCAGGCGGCGATCTGGTGAAAcc                                                            <  2:58504/251‑1 (MQ=255)
                                                                                                                                                                aGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCTATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTTCGCCGTCTGGTTCGGTATTAAAGCCTCGCTGGTAGCGCCAATGAATCGCCTGATTGACAGCGTTCGTCATATTGCAGGCGGCGATCTGGTGAAACCg                                                           <  1:460539/251‑1 (MQ=255)
                                                                                                                                                                    atatCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCTATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTTCGCCGTCTGGTTCGGTATTAAAGCCTCGCTGGTAGCGCCAATGAATCGCCTGATTGACAGCATTCGTCATATTGCAGGCGGCGATCTGGTGAAACCGATTg                                                       <  2:140855/251‑1 (MQ=255)
                                                                                                                                                                             cGGTTGCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCTATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTTCGCCGTCTGGTTCGGTATTAAAGCCTCGCTGGTAGCGCCAATGAATCGCCTGATTGACAGCATTCGTCATATTGCAGGCGGCGATCTGGTGAAACCGATTGAGGTGGATg                                              <  2:143667/251‑1 (MQ=255)
                                                                                                                                                                                          cAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCTATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTTTGCCGTCTGGGTCGGGATTAAAACCTCGCTGGGAGCGCCAATGAATCGCCGGATTGACAGCATTCGGCATATTGCAGGGGGCGGGCTGGGGGAAACGATTGAggggggggggcctaatgaga                                  >  1:139878/1‑232 (MQ=255)
                                                                                                                                                                                           aGTATGTGGCTTACATGGAGCAAAACGATCGGCTCTATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTTTGCCGTCTGGTTCGGTATTTAAGCCTCGCTGGTAGCGCCCATGAATCGCCTGATTTACAGCATTCGTCATAATGCAGGCGGCGATCTGGGGAAACCGACTGAGGGGGATGGGTCTAAttaagtg                                >  1:322114/1‑245 (MQ=255)
                                                                                                                                                                                                gtgGCTTACATGGAGCAAAACGATCGGCTCTATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTTCGCCGTCTGGTTCGGTATTAAAGCCTCGCTGGTAGCGCCAATGAATCGCCTGATTGACAGCATTCGTCATAGTGCAGGCGGCGATCTGGTGAAACCGATTGAGGTGGATGGCTCTAATGAGATGGGGc                            <  1:272646/250‑1 (MQ=255)
                                                                                                                                                                                                 tgGCTTACATGGAGCAAAACGATCGGCTCTATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGCTGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTTCGCCGTCTGGGTTGGGATTAAAACCTCGCTGGTAGCGCCAATGAATAGCCTGTTTGACAGCATTCGTCATATTTCAGGCGGCGGTCTGGTGAAAACGATTTGGGGGGAAGGCTCTAATTTGATggggga                           >  1:285964/1‑248 (MQ=255)
                                                                                                                                                                                                  ggCTTACATGGAGCAAAACGATCGGATCTATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTTTGCCGTCTGGGTCGGGATTAAAGCCTCGCTGGGAGAGCCAATGAATTGCCTGATTTGAAGAATTCGTAATATTTCAGGCGGCGGACCGGTGGAAACGAATTAGGTGGATGGCTCTAATGTGATGGGGGAAc                         >  1:379496/1‑251 (MQ=255)
                                                                                                                                                                                                   gCTTACATGGAGCAAAACGATCGGCTCTATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTTCGCCGTCTGGTTCGGTATTAAAGCCTCGCTGGTAGCGCCAATGAATCGCCTGATTGACAGCATTCGTCATATTGCAGGCGGCGATCTGGTGAAACCGATTGAGGTGGATGGCTCTAATGAGATGGGGCAACt                        <  1:75663/251‑1 (MQ=255)
                                                                                                                                                                                                           ggAGCAAAACGATCGGCTCTATGATAGCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTTCGCCGTCTGGTTCGGTATTAAAGCCTCGCTGGTAGCGCCAATGAATCGCCTTATTGACAGCATTCGTCATATTGCATGCGGTGATCTGGTGAAACCGATTGAGGTGGATGGCTCTAATGAGATGGGGCAACTGGCagag                 <  1:271322/250‑1 (MQ=255)
                                                                                                                                                                                                            gAGCAAAACGATCGGCTCTATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTTCGCCGTCTGGTTCGGTATTAAAGCCTCGCTGGTAGCGCCAATGAATCGCCTGATTGACAGCATTCGTCATATTGCAGGCGGCGATCt                                                                     <  1:10569/197‑1 (MQ=255)
                                                                                                                                                                                                             aGCAAAACGATCGGCTCTATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTTCGCCGTCt                                                                                                                                                     <  1:308365/117‑1 (MQ=255)
                                                                                                                                                                                                               aaaaaCGATCGGCTTTATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCTTACTGGCGGTCATCTTCGCCGTCTGGTTCGGTAGTAAAGCCTCGCTGGTAGCGCCAATGAATCGCCTGATTGACAGCATTCGTCATATTGCAGGCGGCGATCTGGTGAAACCGATTGAGGTGGATGGCTCTAATGAGATGGGGCAACTGGCAGAGAGttt            <  2:424487/250‑1 (MQ=255)
                                                                                                                                                                                                                  aaCGATCGGCTCTATGATAGCGGCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACGGGCGGTCATCTTCGCCGTCTGGTGCGGTATTAAAGCCTCGCTGGTAGCGCCAATGAATCGCCTGATTGACTGCATTCGTCATATTGCAGGCGGCGATCTGGTGAAACCGATTGAGGTGGATGGCTCTAATGAGATGGGGCAACTGGCAGAGAGTTTgcg         <  1:46008/251‑1 (MQ=255)
                                                                                                                                                                                                                    cGATCGGCTCTATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTTCGCCGTCTGGTTCGGTATTAAAGCCTCGCTGGTAGCGCCAATGAATCGCCTGATTGACAGCATTCGTCATATTGCAGGCGGCGATCTGGTGAAACCGATTGAGGTGGATGGCTCTAATGAGATGGGGCAACTGGCAGAGAGTTTGCGcc       <  1:217826/251‑1 (MQ=255)
                                                                                                                                                                                                                       tCGGCTCTATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATATTCGCCGGCTGGTTCGGGATTAAAAACCCGCTGGTAGCGCCCATTGATCCCCCGGGTGTCAGACCTCGTTATATTTAAGGGGGCGGTCTGGCGGAACCGGTTGAGGTGGGTGGCGCTAATGAGATgggggcaccgggcggaagattgcgccctt    >  1:75348/1‑226 (MQ=255)
                                                                                                                                                                                                                         ggCTCTATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGAGTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTTCGCCGTCTGGTTCGGTATTAAAGCCTCGCTGGTAGCGCCAATGAATCGCCTGATTGACAGCATTCGTCATATTGCAGGCGGCGATCTGGTGAAACCGATTGAGGTGGATGGCTCTAATGAGATGGGGCAACTGGCAGAGAGTTTGCGCCATATg  <  2:391541/251‑1 (MQ=255)
                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                      
AACTGGGCGGATTACGAAGCGTTGCCGCGTGACCCGCGTCAGAGCACCGCCGCAGCGGCAGAGATCAAACGTAATTACGATATTTATCACAATGCGCTGGCGGAGCTGATCCAACTGTTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCCATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTTCGCCGTCTGGTTCGGTATTAAAGCCTCGCTGGTAGCGCCAATGAATCGCCTGATTGACAGCATTCGTCATATTGCAGGCGGCGATCTGGTGAAACCGATTGAGGTGGATGGCTCTAATGAGATGGGGCAACTGGCAGAGAGTTTGCGCCATATG  >  NC_000913/4591966‑4592433

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