Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,591,807 T→A 5.2% S51T (TCC→ACC)  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,591,8070TA 5.2% 252.1 / 3.4 77S51T (TCC→ACC) tsrmethyl‑accepting chemotaxis protein I, serine sensor receptor
Reads supporting (aligned to +/- strand):  ref base T (32/41);  new base A (2/2);  total (34/43)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.95e-01

TTAAAACGTATCAAAATTGTGACCAGCTTACTGCTGGTTTTGGCCGTTTTTGGCCTTTTACAACTGACATCAGGCGGTCTGTTCTTTAATGCCTTAAAGAATGACAAAGAAAATTTCACTGTTTTACAAACCATTCGCCAGCAGCAATCCACGCTGAATGGCAGCTGGGTCGCGTTGTTGCAGACGCGTAACACCCTCAACCGCGCGGGTATCCGCTACATGATGGATCAGAATAATATTGGTAGCGGTTCAACCGTTGCTGAGCTGATGGAGAGTGCCAGTATTTCGCTGAAA  >  NC_000913/4591660‑4591953
                                                                                                                                                   |                                                                                                                                                  
ttAAAACGTATCAAAATTGTGACCAGCTTACTGCTGGTTTTGGCCGTTTTTGGCCTTTTACAACTGACATCAGGCGGTCTGTTCTTTAATGCCTTAAAGAATGACAAAGAAAATTTCACTGTTTTACAAACCATTCGCCAGCAGCAATc                                                                                                                                                   >  7:278217/1‑149 (MQ=255)
 tAAAACGTATCAAAATTGTGACCAGCTTACTGCTGGTTTTGGCCGTTTTTGGCCTTTTACAACTGACATCAGGCGGTCTGTTCTTTAATGCCTTAAAGAATGACAAAGAAAATTTCACTGTTTTACAAACCATTAGCCAGCAGCAATc                                                                                                                                                   >  5:354881/1‑148 (MQ=255)
    aaCGTATCAAAATTGTGACCAGCTTACTGCTGGTTTTGGCCGTTTTTGGCCTTTTACAACTGACATCAGGCGGTCTGTTCTTTAATGCCTTAAAGAATGACAAAGAAAATTTCACTGTTTTACAAACCATTCGCCAGCAGCAATCCACg                                                                                                                                               >  5:376041/1‑149 (MQ=255)
    aaCGTATCAAAATTGTGACCAGCTTACTGCTGGTTTTGGCCGTTTTTGGCCTTTTACAACTGACATCAGGCGGTCTGTTCTTTAATGCCTTAAAGAATGACAAAGAAAATTTCACTGTTTTACAAACCATTCGCCAGCAGCAATCCACg                                                                                                                                               <  7:114697/149‑1 (MQ=255)
          tCAAAATTGTGACCAGCTTACTGCTGGTTTTGGCCGTTTTTGGCCTTTTACAACTGACATCAGGCGGTCTGTTCTTTAATGCCTTAAAGAATGACAAAGAAAATTTCACTGTTTTACAAACCATTCGCCAGCAGCAATCCACGCTGAAt                                                                                                                                         <  6:69812/149‑1 (MQ=255)
              aaTTGTGACCAGCTTACTGCTGGTTTTGGCCGTTTTTGGCCTTTTACAACTGACATCAGGCGGTCTGTTCTTTAATGCCTTAAAGAATGACAAAGAAAATTTCACTGTTTTACAAACCATTCGCCAGCAGCAATCCACGCTGAATgg                                                                                                                                       >  3:353958/1‑147 (MQ=255)
                ttGTGACCAGCTTACTGCTGGTTTTGGCAGTTTTTGGCCTTTTACAAATGACATCAGGCGGTCTGTTCTTTAATGCCTTAAAGAATGACAAAGAAAATTTCACTGTTTTACAAACCATTCGCCAGCAGCAATCCACGCTGAATGGCAGc                                                                                                                                   <  5:284992/149‑1 (MQ=255)
                     aCCAGCTTACTGCTGGTTTTGGCCGTTTTTGGCCTTTTACAACTGACATCAGGCGGTCTGTTCTTTAATGCCTTAAAGAATGACAAAGAAAATTTCACTGTTTTACAAACCATTCGCCAGCAGCAATCCACGCTGAATGGCAGCTGGGt                                                                                                                              >  7:312233/1‑149 (MQ=255)
                      ccAGCTTACTGCTGGTTTTGGCCGTTTTTGGCCTTTTACAACTGACATCAGGCGGTCTGTTCTTTAATGCCTTAAAGAATGACAAAGAAAATTTCACTGTTTTACAAACCAATCGCCAGCAGCAATCCACGCTGAATGGCAGCTGGGTc                                                                                                                             <  1:337859/149‑1 (MQ=255)
                       cAGCTTACTGCTGGTTTTGGCCGTTTTTGGCCTTTTACAACTGACATCAGGCGGTCTGTTCTTTAATGCCTTAAAGAATGACAAAGAAAATTTCACTGTTTTACAAACCATTCGCCAGCAGCAATCCACGCTGAATGGCAGCTGGGTcg                                                                                                                            >  5:202410/1‑149 (MQ=255)
                         gCTTACTGCTGGTTTTGGCCGTTTTTGGCCTTTTACAACTGACATCAGGCGGTCTGTTCTTTAATGCCTTAAAGAATGACAAAGAAAATTTCACTGTTTTACAAACCATTCGCCAGCAGCAATCCACGCTGAATGGCAGCTGGGTcgcg                                                                                                                          >  1:44927/1‑149 (MQ=255)
                          cTTACTGCTGGTTTTGGCCGTTTTTGGCCTTTTACAACTGACATCAGGCGGTCTGTTCTTTAATGCCTTAAAGAATGACAAAGAAAATTTCACTGTTTTACAAACCATTCGCCAGCAGCAATCCACGCTGAATGGCAGCTGGGTCGCgt                                                                                                                         <  3:225953/149‑1 (MQ=255)
                             aCTGCTGGTTTTGGCCGTTTTTGGCCTTTTACAACTGACATCAGGCGGTCTGTTCTTTAATGCCTTAAAGAATGACAAAGAAAATTTCACTGTTTTACAAACCATTCGCCAGCAGCAATCCACGCTGAATGGCAGCTGGGTCGCgttg                                                                                                                       <  2:44927/148‑1 (MQ=255)
                              ctgctgGTTTTGGCCGTTTTTGGCCTTTTACAACTGACATCAGGCGGTCTGTTCTTTAATGCCTTAAAGAATGACAAAGAAAATTTCACTGTTTTACAAACCATTCGCCAGCAGCAATCCACGCTGAATGGCAGCTGGGTCGCgttgtt                                                                                                                     <  5:190999/149‑1 (MQ=255)
                              ctgctgGTTTTGGCCGTTTTTGGCCTTTTACAACTGACATCAGGCGGTCTGTTCTTTAATGCCTTAAAGAATGACAAAGAAAATTTCACTGTTTTACAAACCATTCGCCAGCAGCAATCCACGCTGAATGGCAGCTGGGTCGCgttgt                                                                                                                      <  5:176569/148‑1 (MQ=255)
                                 ctgGTTTTGGCCGTTTTTGGCCTTTTACAACTGACATCAGGCGGTCTGTTCTTTAATGCCTTAAAGAATGACAAAGAAAATTTCACTGTTTTACAAACCATTCGCCAGCAGCAATCCACGCTGAATGGCAGCTGGGTCGCgtt                                                                                                                        >  3:214892/1‑143 (MQ=255)
                                 ctgGTTTTGGCCGTTTTTGGCCTTTTACAACTGACATCAGGCGGTCTGTTCTTTAATGCCTTAAAGAATGACAAAGAAAATTTCACTGTTTTACAAACCATTCGCCAGCAGCAATCCACGCTGAATGGCAGCTGGGTCGCgtt                                                                                                                        <  4:214892/143‑1 (MQ=255)
                                 ctgGTTTTGGCCGTTTTTGGCCTTTTACAACTGACATCAGGCGGTCTGTTCTTTAATGCCTTAAAGAATGACAAAGAAAATTTCACTGTTTTACAAACCATTCGCCAGCAGCAATCCACGCTGAATGGCAGCTGGGTCGCGttgttg                                                                                                                    <  6:376041/147‑1 (MQ=255)
                                 ctgGTTTTGGCCGTTTTTGGCCTTTTACAACTGACATCAGGCGGTCTGTTCTTTAATGCCTTAAAGAATGACAAAGAAAATTTCACTGTTTTACAAACCATTCGCCAGCAGCAATCCACGCTGAATGGCAGCTGGGTCGCGTTGTTGc                                                                                                                   <  4:119792/148‑1 (MQ=255)
                                 ctgGTTTTGGCCGTTTTTGGCCTTTTACAACTGACATCAGGCGGTCTGTTCTTTAATGCCTTAAAGAATGACAAAGAAAATTTCACTGCTTTACAAACCATTCGCCAGCAGCAATCCACGCTGAATGGc                                                                                                                                      <  1:83643/129‑1 (MQ=255)
                                 ctgGTTTTGGCCGTTTTTGGCCTTTTACAACTGACATCAGGCGGTCTGTTCTTTAATGCCTTAAAGAATGACAAAGAAAATTTCACTGCTTTACAAACCATTCGCCAGCAGCAATCCACGCTGAATGGc                                                                                                                                      >  2:83643/1‑129 (MQ=255)
                                    gTTTTGGCCGTTTTTGGCCTTTTACAACTGACATCAGGCGGTCTGTTCTTTAATGCCTTAAAGAATGACAAAGAAAATTTCACTGTTTTACAAACCATTCGCCAGCAGCAATCCACGCTGAATGGCAGCTGGGTCGCGTTGTTGCAGAc                                                                                                               <  8:278217/149‑1 (MQ=255)
                                    gTTTTGGCCGTTTTTGGCCTTTTACAACTGACATCAGGCGGTCTGTTCTTTAATGCCTTAAAGAATGACAAAGAAAATTTCACTGTTTTACAAACCATTCGCCAGCAGCAATCCACGCTGAATGGCAGCTGGGTCGCGTTGTTGCAGAc                                                                                                               <  4:107422/149‑1 (MQ=255)
                                       ttGGCCGTTTTTGGCCTTTTACAACTGACATCAGGCGGTCTGTTCTTTAATGCCTTAAAGAATGACAAAGAAAATTTCACTGTTTTACAAACCATTCGCCAGCAGCAATCCACGCTGAATGGCa                                                                                                                                     <  7:58736/124‑1 (MQ=255)
                                       ttGGCCGTTTTTGGCCTTTTACAACTGACATCAGGCGGTCTGTTCTTTAATGCCTTAAAGAATGACAAAGAAAATTTCACTGTTTTACAAACCATTCGCCAGCAGCAATCCACGCTGAATGGCa                                                                                                                                     >  8:58736/1‑124 (MQ=255)
                                       ttGGCAGTTTTTGGCCTTTTACAACTGACATCATGCGTTCTGTTCTTTAATGCCTTAAAGAATGACAAAGAAAATTTCACTGTTTTACAAACCATTCGCCTGCAGCAATCCAGGCTGAATGGCAGCTGGGTCGCGTTGTTGCAGAcgc                                                                                                             <  8:333912/148‑1 (MQ=255)
                                          gacGTTTTTGGCCTTTTACAACTGACATCACGCGGTCTGTTCTTTAATGCCTTAAAGAATGACAAAGAAAATTTCACTGTTTTACAACCCATTCGCCAGCAGCAATCCACGCTGAATGGCAGCTCGGTCGCGTTGTTGCAGACGCGTaa                                                                                                         >  2:252937/3‑149 (MQ=255)
                                          gCCGTTTTTGGCCTTTTACAACTGACATCAGGCGGTCTGTTCTTTAATGCCTTAAAGAATGACAAAGAAAATTTCACTGTTTTACAAACCATTCGCCAGCAGCAATCCACGCTGAATGGCAGCTGGGTCGCGTTGTTGCAGACGCGTaa                                                                                                         >  8:125053/1‑149 (MQ=255)
                                          gCCGTTTTTGGCCTTTTACAACTGACATCAGGCGGTCTGTTCTTTAATGCCTTAAAGAATGACAAAGAAAATTTCACTGTTTTACAAACCATTCGCCAGCAGCAATCCACGCTGAATGGCAGCTGGGTCGCGTTGTTGCAGACGCGTa                                                                                                          <  8:372391/148‑1 (MQ=255)
                                            cGTTTTTGGCCTTTTACAACTGACATCAGGCGGTCTGTTCTTTAATGCCTTAAAGAATGACAAAGAAAATTTCACTGTTTTACAAACCATTCGCCAGCAGCAATCCACGCTGAATGGCAGCTGGGTCGCGTTGTTGCAGACGCGTAaca                                                                                                       >  1:231239/1‑149 (MQ=255)
                                             gTTTTTGGCCTTTTACAACTGACATCAGGCGGTCTGTTCTTTAATGCCTTAAAGAATGACAAAGAAAATTTCACTGTTTTACAAACCATTCGCCAGCAGCAATCCACGCTGAATGGCAGCTGGGTCGCGTTGTTGCAGACGCGTAacac                                                                                                      >  7:118003/1‑149 (MQ=255)
                                              tttttGGCCTTTTACAACTGACATCAGGCGGTCTGTTCTTTAATGCCTTAAAGAATGACAAAGAAAATTTCACTGTTTTACAAACCATTCGCCAGCAGCAATCCACGCGGAATGGCAGCTGGGTCGCGTTGTTGCAGACGCGTAACAcc                                                                                                     <  1:270773/149‑1 (MQ=255)
                                                         ttACAACTGACATCAGGCGGTCTGTTCTTTAATGCCTTAAAGAATGACAAAGAAAATTTCACTGTTTTACAAACCATTCGCCAGCAGCAATCCACGCTGAATTGCATCTGGGTCGCGTTGTAGCAGACGCGTAACACCCTCAACcgcgc                                                                                          >  1:259625/1‑149 (MQ=255)
                                                             aaCTGACATCAGGCGGTCTGTTCTTTAATGCCTTAAAGAATGACAAAGAAAATTTCACTGTTTTACAAACCATTCGCCAGCAGCAATCCACGCTGAATGGCAGCTGGGTCGCGTTGTTGCAGACGCGTAACACCCTCAACCGCGCGGGt                                                                                      <  2:231239/149‑1 (MQ=255)
                                                                     tCAGGCGGTCTGTTCTTTAATGCCTTAAAGAATGACAAAGAAAATTTCACTGTTTTACAAACCATTCGCCAGCAGCAATCCACGCTGAATGGCAGCTGGGTCGCGTTGTTGCAGACGCGTAACACCCTCAACCGCGCGGGTATCCGCTa                                                                              <  8:90878/149‑1 (MQ=255)
                                                                        ggCGGTCTGTTCTTTAATGCCTTAAAGAATGACAAAGAAAATTTCACTGTTTTACAAACCATTCGCCAGCAGCAATCCACGCTGAATGGCAGCTGGGTCGCGTTGTTGCAGACGCGTAACACCCTCAACCGCGCgg                                                                                        >  3:343291/1‑136 (MQ=255)
                                                                        ggCGGTCTGTTCTTTAATGCCTTAAAGAATGACAAAGAAAATTTCACTGTTTTACAAACCATTCGCCAGCAGCAATCCACGCTGAATGGCAGCTGGGTCGCGTTGTTGCAGACGCGTAACACCCTCAACCGCGCgg                                                                                        <  4:343291/136‑1 (MQ=255)
                                                                            gTCTGTTCTTTAATGCCTTAAAGAATGACAAAGAAAATTTCACTGTTTTACAAACCATTCGCCAGCAGCAATCCACGCTGAATGGCAGCTGGGTCGCGTTGTTGCAGACGCGTAACACCCTCAACCGCGCGGGTATCCGCTACatgatg                                                                       >  1:136641/1‑149 (MQ=255)
                                                                               tGTTCTTTAATGCCTTAAAGAATGACAAAGAAAATTTCACTGTTTTACAAACCATTCGCCAGCAGCAATCCACGCTGAATGGCAGCTGGGTCGCGTTGTTGCAGACGCGTAACACCCTCAACCGCGCGGGTATCCGCTACATGATGGAt                                                                    >  6:224404/1‑149 (MQ=255)
                                                                               tGTTCTTTAATGCCTTAAAGAATGACAAAGAAAATTTCACTGTTTTACAAACCATTCGCCAGCAGCAATCCACGCTGAATGGCAGCTGGGTCGCGTTGTTGCAGACGCGTAACACCCTCAACCGCGCGGGTATCCGCTACATGATGGAt                                                                    >  6:5502/1‑149 (MQ=255)
                                                                                 ttCTTTAATGCCTTAAAGAATGACAAAGAAAATTTCACTGTTTTACAAACCATTCGCCAGCAGCAATCCACGCTGAATGGCAGCTGGGTCGCGTTGTTGCAGACGCGTAACACCCTCAACCGCGCGGGTATCCGCTACATGATGGATCa                                                                  <  4:353958/149‑1 (MQ=255)
                                                                                 ttCTTTAATGCCTTAAAGAATGACAAAGAAAATTTCACTGTTTTACAAACCATTCGCCAGCAGCAATCCACGCTGAATGGCAGCTGGGTCGCGTTGTTGCAGACGCGTAACACCCTCAACCGCGCGGGTATCCGCTACATGATGGATCa                                                                  <  2:136641/149‑1 (MQ=255)
                                                                                 ttCTATAATGCCTTAAAGAATGACAAAGAAAATTTCACTGTTTTACAAACCATTCGCCAGCAGCAATCCACGCTGAATGGCAGCTGGGTCGCGGTGTTGCAGACGCGTAACACCCTCAACCGCGCGGGTATCCGCTACATGATGGAt                                                                    <  4:107550/147‑1 (MQ=255)
                                                                                    tttAATGCCTTAAAGAATGACAAAGAAAATTTCACTGTTTTACAAACCATTCGCCAGCAGCAATCCACGCTGAATGGCAGCTGGGTCGCGTTGTTGCAGACGCGTAACACCCTCAACCGCGCgg                                                                                        >  1:73288/1‑124 (MQ=255)
                                                                                    tttAATGCCTTAAAGAATGACAAAGAAAATTTCACTGTTTTACAAACCATTCGCCAGCAGCAATCCACGCTGAATGGCAGCTGGGTCGCGTTGTTGCAGACGCGTAACACCCTCAACCGCGCgg                                                                                        <  2:73288/124‑1 (MQ=255)
                                                                                       aaTGCCTTAAAGAATGACAAAGAAAATTTCACTGTTTTACAAACCATTCGCCAGCAGCAATCCACGCTGAATGGCAGCTGGGTCGCGTTGTTGCAGACGCGTAACACCCTCAACCGCGCGGGTATCCGCTACATGATGGATCAGaataa                                                            >  6:62911/1‑149 (MQ=255)
                                                                                         tGCCTTAAAGAATGACAAAGAAAATTTCACTGTTTTACAAACCATTCGCCAGCAGCAATCCACGCTGAATGGCAGCTGGGTCGCGTTGTTGCAGACGCGTAACACCCTCAACCGCGCGGGTATCCGCTACATGATGGATCAGAataata                                                          <  6:202410/149‑1 (MQ=255)
                                                                                            cTTAAAGAATGACAAAGAAAATTTCACTGTTTTACAAACCATTCGCCAGCAGCAATCCACGCTGAATGGCAGCTGGGTCGCGTTGTTGCAGACGCGTAACACCCTCAACCGCGCGGGTATCCGCTACATGATGGATCAGAATAATATTg                                                       >  1:239921/1‑149 (MQ=255)
                                                                                             ttAAAGAATGACAAAGAATATTTGACTGTTTTACAAACCATTCGCCAGCAGCAATCCACGCTGAATGGCAGCTGGGTCGCGTTGTTGCAGACGCGTAACACCCTCaa                                                                                                <  5:245177/107‑1 (MQ=255)
                                                                                             ttAAAGAATGACAAAGAAAATTTCACTGTTTTACAAACCATTCGCCAGCAGCAATCCACGCTGAATGGCAGCTGGGTCGCGTTGTTGCAGACGCGTAACACCCTCaa                                                                                                >  6:245177/1‑107 (MQ=255)
                                                                                             ttAAAGAATGACAAAGAAAATTTCACTGTTTTACAAACCATTCGCCAGCAGCAATCCACGCTGAATGGCAGCTGGGTCGCGTTGCTGCAGACGCGTAACACCCTCAACCGCGCGGGTATCCGCTACATGATGGATc                                                                   >  1:126197/1‑136 (MQ=255)
                                                                                             ttAAAGAATGACAAAGAAAATTTCACTGTTTTACAAACCATTCGCCAGCAGCAATCCACGCTGAATGGCAGCTGGGTCGCGTTGCTGCAGACGCGTAACACCCTCAACCGCGCGGGTATCCGCTACATGATGGATc                                                                   <  2:126197/136‑1 (MQ=255)
                                                                                               aaaGAATGACAAAGAAAATTTCACTGTTTTACAAACCATTCGCCAGCAGCAATCCACGCTGAATGGCAGCTGGGTCGCGTTGTTGCAGACGCGTAACACCCTCAACCGCGCGGGTATCCGCTACATGATGGATCAGAATAATATTgg                                                      <  8:118003/147‑1 (MQ=255)
                                                                                                   aaTGACAAAGAAAATTTCACTGTTTTACAAACCATTCGCCAGCAGCAATCCACGCTGAATGGCAGCTGGGTCGCGTTGTTGCAGACGCGTAACACCCTCAACCGCGCGGGTATCCGCTACATGATGGATCAGAATAATATTGGTAGCgg                                                >  6:172779/1‑149 (MQ=255)
                                                                                                         aaaGAAAATTTCACTGTTTTACAAACCATTCGCCAGCAGCAATCCACGCTGAATGGCAGCTGGGTCGCGTTGTTGCAGACGCGTAACACCCTCAACCGCGCGGGTATCCGCTACATGATGGATCAGaata                                                             <  8:47251/130‑1 (MQ=255)
                                                                                                         aaaGAAAATTTCACTGTTTTACAAACCATTCGCCAGCAGCAATCCACGCTGAATGGCAGCTGGGTCGCGTTGTTGCAGACGCGTAACACCCTCAACCGCGCGGGTATCCGCTACATGATGGATCAGaata                                                             >  7:47251/1‑130 (MQ=255)
                                                                                                         aaaGAAAATTTCACTGTTTTACAAACCATTCGCCAGCAGCAATCCACGCTGAATGGCACCTGGGTCGCGTTGTTGCAGCCGCGTAACACCCTCAACCGCGCCGGTAACCGCAACATGATGGATCAGAAGAATATTGGTAGCGGTTCAAc                                          >  1:288900/1‑149 (MQ=255)
                                                                                                           aGAAAATTTCACTGTTTTACAAACCATTCGCCAGCAGCAAACCACGCTGAATgg                                                                                                                                       <  5:268885/54‑1 (MQ=255)
                                                                                                           aGAAAATTTCACTGTTTTACAAACCATTCGCCAGCAGCAAACCACGCTGAATgg                                                                                                                                       >  6:274075/1‑54 (MQ=255)
                                                                                                           aGAAAATTTCACTGTTTTACAAACCATTCGCCAGCAGCAAACCACGCTGAATgg                                                                                                                                       <  5:274075/54‑1 (MQ=255)
                                                                                                           aGAAAATTTCACTGTTTTACAAAACATTCGCCAGCAGCAAACCACGCTGAATgg                                                                                                                                       >  6:268885/1‑54 (MQ=255)
                                                                                                             aaaaTTTCACTGTTTTACAAACCATTCGCCAGCAGCAATCCACGCTGAATGGCAGCTGGGTCGCGTTGTTGCAGACGCGTAACACCCTCAACCGCGCGGGTATCCGCTACATGATGGATCAGAATAATATTGGTAGCGGTTCAACCGtt                                      <  1:252937/149‑1 (MQ=255)
                                                                                                               aaTTTCACTGTTTTACAAACCATTCGCCAGCAGCAATCCACGCTGAATGGCAGCTGGGTCGCGTTGTTGCAGACGCGTAACACCCTCAACCGCGCGGGTATCCGCTACATGATGGATCAGAATAATATTGGTAGCGGTTCAACCGTTgc                                    <  2:239921/149‑1 (MQ=255)
                                                                                                                   tCACTGTTTTACAAACCATTCGCCAGCAGCAATCCACGCTGAATGGCAGCTGGGTCGCGGTGTTGCAGACGCGTAACACCCTCAACCGCGCGGGTATCCGCTACATGATGGATCAGAATAATATTGGTAGCGGTTCAACCGTTgctga                                 <  5:5502/148‑1 (MQ=255)
                                                                                                                     aCTGTTTTACAAACCATTCGCCAGCAGCAATCCACGCTGAATGGCAGCTGGGTCGCGTTGTTGCAGACGCGTAACACCCTCAACCGCGCGGGTATCCGCTACATGATGGATCAGAATAATATTGGTAGCGGTTCAACCGTTgctgagc                               <  1:243438/148‑1 (MQ=255)
                                                                                                                         ttttACAAACCATTCGCCAGCAGCAATCCACGCTGAATGGCAGCTGGGTCGCGTTGTTGCAGACGCGTAACACCCTCAACCGCGCGGGTATCCGCTACATGATGGATCAGAATAATATTGGTAGCGGTTCAACCGTTGCTGAGCTGATg                          >  8:61388/1‑149 (MQ=255)
                                                                                                                           ttACAAACCATTCGCCAGCAGCAATCCACGCTGAATGGCAGCTGGGTCGCGTTGTTGCAGACGCGTAACACCCTCAACCGCGCGGGTATCCGCTACATGATGGATCAGAATAATATTGGTAGCGGTTCAACCGTTGCTGAGCTGATgg                         <  6:354881/148‑1 (MQ=255)
                                                                                                                            tACAAACCATTCGCCAGCAGCAATCCACGCTGAATGGCAGCTGGGTCGCGTTGTTGCAGACGCGTAACACCCTCAACCGCGCGGGTATCCGCTACATGATGGATCAGAATAATATTGGTAGCGGTTCAACCGTTGCTGAGCTGATGga                        <  7:125053/148‑1 (MQ=255)
                                                                                                                             aCAAACCATTCGCCAGCAGCAATCCACGCTGAATGGCAGCTGGGTCGCGTTGTTGCAGACGCGTAACACCCTCAACCGCGCGGGTATCCGCTACATGATGGATCAGAATAATATTGGTAGCGGTTCAACCGTTGCTGAGCTGATGga                        >  4:247804/1‑147 (MQ=255)
                                                                                                                               aaaCCATTCGCCAGCAGCAATCCACGCTGAATGGCAGCTGGGTCGCGTTGTTGCAGACGCGTAACACCCTCAACCGCGCGGGTATCCGCTACATGATGGATCAGAATAATATTGGTAGCGGTTCAACCGTTGCTGAGCTGATGGAGAGt                    <  7:238279/149‑1 (MQ=255)
                                                                                                                                            gcagcaATCCACGCTGAATGGCAGCTGGGTCGCGTTGTTGCAGACGCGTAACACCCTCAACCGCGCGGGTATCCGCTACATGATGGATCAGAATAATATTGGTAGCGGTTCAACCGTTGCTGAGCTGATGGAGAGTGCCa                >  4:301067/1‑140 (MQ=255)
                                                                                                                                            gcagcaATCCACGCTGAATGGCAGCTGGGTCGCGTTGTTGCAGACGCGTAACACCCTCAACCGCGCGGGTATCCGCTACATGATGGATCAGAATAATATTGGTAGCGGTTCAACCGTTGCTGAGCTGATGGAGAGTGCCa                <  3:301067/140‑1 (MQ=255)
                                                                                                                                            gcagcaATCCACGCTGAATGGCAGCTGGGTCGCGTTGTTGCAGACGCGTAACACCCTCAACCGCGCGGGTATCCGCTACATGATGGATCAGAATAATATTGGTAGCGGTTCAACCGTTGCTGAGCTGATGGAGAGTGCCAGTATTTCGc       <  5:224404/149‑1 (MQ=255)
                                                                                                                                                caATCCACGCTGAATGGCAGCTGGGTCGCGTTGTTGCAGACGCGTAACACCCTCAACCGCGCGGGTATCCGCTACATGATGGATCAGAATAATATTGGTAGCGGTTc                                             >  8:100013/1‑107 (MQ=255)
                                                                                                                                                caATCCACGCTGAATGGCAGCTGGGTCGCGTTGTTGCAGACGCGTAACACCCTCAACCGCGCGGGTATCCGCTACATGATGGATCAGAATAATATTGGTAGCGGTTc                                             <  7:100013/107‑1 (MQ=255)
                                                                                                                                                  aTCCACGCTGAATGGCAGCTGGGTCGCGTTGTTGCAGACGCGTAACACCCTCAACCGCGCGGGTATCCGCTACATGATGGATCAGAATAATATTGGTAGCGGTTCAACCGTTGCTGAGCTGATGGAGAGTGCCAGTATTTc         >  8:324151/1‑141 (MQ=255)
                                                                                                                                                  aTCCACGCTGAATGGCAGCTGGGTCGCGTTGTTGCAGACGCGTAACACCCTCAACCGCGCGGGTATCCGCTACATGATGGATCAGAATAATATTGGTAGCGGTTCAACCGTTGCTGAGCTGATGGAGAGTGCCAGTATTTc         <  7:324151/141‑1 (MQ=255)
                                                                                                                                                   tCCACGCTGAATGGCAGCTGGGTCGCGTTGTTGCAGACGCGTAACACCCTCAACCGCGCGGGTATCCGCTACATGATGGATCAGAATAATATTGGTAGCGGTTCAACCGTTGCTGAGCGGCTGGAGAGTGCCAGTATTTGGCTGaaa  >  8:307026/1‑147 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
TTAAAACGTATCAAAATTGTGACCAGCTTACTGCTGGTTTTGGCCGTTTTTGGCCTTTTACAACTGACATCAGGCGGTCTGTTCTTTAATGCCTTAAAGAATGACAAAGAAAATTTCACTGTTTTACAAACCATTCGCCAGCAGCAATCCACGCTGAATGGCAGCTGGGTCGCGTTGTTGCAGACGCGTAACACCCTCAACCGCGCGGGTATCCGCTACATGATGGATCAGAATAATATTGGTAGCGGTTCAACCGTTGCTGAGCTGATGGAGAGTGCCAGTATTTCGCTGAAA  >  NC_000913/4591660‑4591953

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