Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,043,899 A→T 5.1% I104I (ATT→ATA etk ← tyrosine‑protein kinase, role in O‑antigen capsule formation

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,043,8990AT 5.1% 253.2 / 2.8 78I104I (ATT→ATAetktyrosine‑protein kinase, role in O‑antigen capsule formation
Reads supporting (aligned to +/- strand):  ref base A (35/39);  new base T (2/2);  total (37/41)
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.34e-01

AGTTGCTGATCCTGACCATTCAGTTGTGGAATATGCATCCAGCTGATCGCCAGCTCACCTGGTTTTTCTTTGGTTAATCTCGCCCAGCCGCGACCCACAATCGGAAAATACTTCTGCTCAACTATGTCGCGCAGATTCAGTTCAGCAATGGTTTTACCGAGAATCATGCGCGATTGCAGCAGTTGGATCTCCGGTGCAGACTCGGGCGATGAGTTAGGGATCATATCGCTCAGGCCGCTGAGAATGGCGTTGCCCTGTTTTTGCTCAACCTGGACCAGAGTATCTGCCTGATA  >  NC_000913/1043753‑1044045
                                                                                                                                                  |                                                                                                                                                   
aGTTGCTGATCCTGACCATTCAGTTGTGGAATATGCATCCAGCTGATCGCCAGCTCACCTGGTTTTTCTTTGGTTAATCTCGCCCAGCCGCGACCCACAATCGGAAAATACTTCTGCTCAACTATGTCGCGCAGATTCAGTTCAGCAAt                                                                                                                                                   <  7:117083/149‑1 (MQ=255)
 gTTGCTGATCCTGACCATTCAGTTGTGGAATATGCATCCAGCTGATCGCCAGCTCACCTGGTTTTTCTTTGGTTAATCTCGCCCAGCCGCGACCCACAATCGGAAAATACTTCTGCTCAACTATGTCGCGCAGATTCAGTTCAGCAATg                                                                                                                                                  <  3:189059/149‑1 (MQ=255)
    gCTGATCCTGACCATTCAGTTGTGGAATATGCATCCAGCTGATCGCCAGCTCACCTGGTTTTTCTTTGGTTAATCTCGCCCAGCCGCGACCCACAATCGGAAAATACTTCTGCTCAACTATGTCGCGCAGATTCAGTTCAGCAATGGt                                                                                                                                                <  1:67378/148‑1 (MQ=255)
    gCTGATCCTGACCATTCAGTTGTGGAATATGCATCCAGCTGATCGCCAGCTCACCTGGTTTTTCTTTGGTTAATCTCGCCCAGCCGCGACCCACAATCGGAAAATACTTCTGCTCAACTATGACGCGCAGATTCAGTTCAGCAATGGtt                                                                                                                                               >  5:63895/1‑149 (MQ=255)
     cTGATCCTGACCATTCAGTTGTGGAATATGCATCCAGCTGATCGCCAGCTCACCTGGTTTTTCTTTGGTTAATCTCGCCCAGCCGCGACCCACAATCGGAAAATACTTCTGCTCAACTATGTCGCGCAGATTCAGTTCAGCAATGGttt                                                                                                                                              <  1:200453/149‑1 (MQ=255)
       gATCCTGACCATTCAGTTGTGGAATATGCATCCAGCTGATCGCCAGCTCACCTGGTTTTTCTTTGGTTAATCTCGCCCAGCCGCGACCCACAATCGGAAAATACTTCTGCTCAACTATGTCGCGCAGATTCAGTTCAGCAATGGTTTTa                                                                                                                                            >  7:2307/1‑149 (MQ=255)
        aTCCTGCCCATTCAGTTGTGGAATATGCATCCAGATGATCGCCAGATCACCTGGTTTTTATTTGGTTAATCTCGCCCAGCCGCGACCCACAATCGGAAAATACTTCTGCTCAACTATGTCGCGCAGATTCAGTTCAGCAATGGGTTTAc                                                                                                                                           >  5:157835/1‑149 (MQ=255)
        aTCCTGACCATTCAGTTGTGGAATATGCATCCAGCTGATCGCCAGCTCACCTGGTTTTTCTTTGGTTAATCTCGCCCAGCCGCGACCCACAATCGGAAAATACTTCTGCTCAACTATGTCGCGCAGATTCAGTTCAGCAATGGTTTTAc                                                                                                                                           <  6:291662/149‑1 (MQ=255)
          ccTGACCATTCAGTTGTGGAATATGCATCCAGCTGATCGCCAGCTCACCTGGTTTTTCTTTGGTTAATCTCGCCCAGCCGCGACCCACAATCGGAAAATACTTCTGCTCAACTATGTCGCGCAGATTCAGTTCAGCAATGGTTTTAcc                                                                                                                                          >  5:225245/1‑148 (MQ=255)
            tGACCATTCAGTTGTGGAATATGCATCCAGCTGATCGCAAGCTCACCTGGTTTTTCTTTGGTTAATCTCGCCCAGCCGCGCCCCACAATCGGAAAATACTTCTGCTCAACTATGTCGCGCAGATTCAGTTCAGCAATGGTTTTAcc                                                                                                                                          >  1:77195/1‑146 (MQ=255)
                  ttCAGTTGTGGAATATGCATCCAGCTGATCGCCAGCTCACCTGGTTTTTCTTTGGTTAATCTCGCCCAGCCGCGACCCACAATCGGAAAATACTTCTGCTCAACTATGTCGCGCAGATTCAGTTCAGCAATGGTTTTACCGAGAATCAt                                                                                                                                 <  6:160522/149‑1 (MQ=255)
                  ttCAGTTGTGGAATATGCATCCAGCTGATCGCCAGCTCACCTGGTTTTTCTTTGGTTAATCTCGCCCAGCCGCGACCCACAATCGGAAAATACTTCTGCTCAACTATGTCGCGCAGATTCAGTTCAGCAATGGTTTTACCGAGAATCAt                                                                                                                                 <  8:2307/149‑1 (MQ=255)
                  ttCAGTTGTGGAATATGCATCCAGCTGATCGCCAGCTCACCTGGTTTTTCTATGGTTAATCTCGCCCAGCCGCGACCCCAATTCGGAAAATACTTCTTCTCAACTATGTCGCGCAGATTGAGTTCAGCAATGGTTTTACCGAGAATc                                                                                                                                   >  1:84636/1‑147 (MQ=255)
                     aGTTGTGGAATATGCATCCAGCTGATCGCCAGCTCACCTGGTTTTTCTTTGGTTAATCTCGCCCAGCCGCGACCCACAATCGGAAAATACTTCTGCTCAACTATGTCGCGCATATTCAGTTCAGCAATGGTTTTACCGAGAATCATgcg                                                                                                                              >  3:44290/1‑149 (MQ=255)
                     aGTTGTGGAATATGCATCCAGCTGATCGCCAGCTCACCTGGTTTTTCTTTGGTTAATCTCGCCCAGCCGCGACCCACAATCGGAAAATACTTCTGCTCAACTATGTCGCGCAGATTCAGTTCAGCAATGGTTTTACCGAGAATCATgcg                                                                                                                              >  8:195347/1‑149 (MQ=255)
                            gAATATGCATCCAGCTGATCGCCAGCTCACCTGGTTTTTCTTTGGTTAATCTCGCCCAGCCGCGACCCACAATCGGAAAATACTTCTGCTCAACTATGTCGCGCAGATTCAGTTCAGCAATGGTTTTACCGAGAATCATGCGCGATTg                                                                                                                        <  8:18702/148‑1 (MQ=255)
                            gAATATGCATCCAGCTGATCGCCAGCTCACCTGGTTTTTCTTTGGTTAATCTCGCCCAGCCGCGACCCACAATCGGAAAATACTTCTGCTCAACTATGTCGCGCAGATTCAGTTCAGCAATGGTTTTACCGAGAATCATGCGCGATTg                                                                                                                        <  1:121950/148‑1 (MQ=255)
                                  gCATCCAGCTGATCGCCAGCTCACCTGGTTTTTCTTTGGTTAATCTCGCCCAGCCGCGACCCACAATCGGAAAATACTTCTGCTCAACTATGTCGCGCAGATTCAGTTCAGCAATGGTTTTACCg                                                                                                                                         >  3:277483/1‑125 (MQ=255)
                                  gCATCCAGCTGATCGCCAGCTCACCTGGTTTTTCTTTGGTTAATCTCGCCCAGCCGCGACCCACAATCGGAAAATACTTCTGCTCAACTATGTCGCGCAGATTCAGTTCAGCAATGGTTTTACCg                                                                                                                                         <  4:277483/125‑1 (MQ=255)
                                      ccAGCTGATCGCCAGCTCACCTGGTTTTTCTTTGGTTAATCTCGCCCAGCCGCGACCCACAATCGGAAAATACTTCTGCTCAACTATGTCGCGCAGATTCAGTTCAGCAATGGTTTTACCGAGAATCATGCGCGATTGCAGCAGTTGGa                                                                                                            >  4:50985/1‑149 (MQ=255)
                                          cTGATCGCCAGCTCACCTGGTTTTTCTTTGGTTAATCTCGCCCAGCCGCGACCCACAATCGGAAAATACTTCTGCTCAACTATGTCGCGCAGATTCAGTTCAGCAATGGTTTTACCGAGAATCATGCGCGATTGCAGCAGTTGGAtctc                                                                                                        >  5:16991/1‑149 (MQ=255)
                                             aTCTCCATCTCACCTGGTTTTTCTTTGGTTTATCTCGCCCAGCCGCGACCCACAATCGGAAAATCCTTCTGCTCAACTATGTCTAGCAGATTCAGTTCAGCAATGGTTTTAGCGAGTATCATGCGCGATTGCAGCAGTTGGATCTcc                                                                                                       <  7:214981/147‑1 (MQ=255)
                                                gCCAGCTCACCTGGTTTTTCTTTGGTTAATCTCGCCCAGCCGCGACCCACAATCGGAAAATACTTCTACTCAACTATGTCGCGCAGATTCAGTTCAGCAATGGTTTTACCGAGAATCATGCGCGATTGCAGCAGTTGGATCTCCGGTg                                                                                                   <  5:266368/148‑1 (MQ=255)
                                                 ccAGCTCACCTGGTTTTTCTTTGGTTAATCTCGCCCAGCCGCGACCCACAATCGGAAAATACTTCTGCTCAACTATGTCGCGCAGATTCAGTTCAGCAATGGTTTTACCGAGAATCATGCGCGATTGCAGCAGTTGGATCTCCGGTg                                                                                                   >  8:309474/1‑147 (MQ=255)
                                                 ccAGCTCACCTGGTTTTTCTTTGGTTAATCTCGCCCAGCCGCGACCCACAATCGGAAAATACTTCTGCTCAACTATGTCGCGCAGATTCAGTTCAGCAATGGTTTTACCGAGAATCATGCGCGATTGCAGCAGTTGGATCTCCGGTGCa                                                                                                 <  4:277205/149‑1 (MQ=255)
                                                   aGCTCACCTGGTTTTTCTTTGGTTAATCTCGCCCAGCCGCGACCCACAATCGGAAAATACTTCTGCTCAACTATGTCGCGCAGATTCAGTTCAGCAATGGTTTTACCGAGAATCATGCGCGATTGCAGCAGTTGGATCTCCGGTGCAGa                                                                                               <  4:44290/149‑1 (MQ=255)
                                                   aGCTCACCTGGTTTTTCTTTGGTTAATCTCGCCCAGCCGCGACCCACAATCGGAAAATACTTCTGCTCAACTATGTCGCGCAGATTCAGTTCAGCAATGGTTTTACCGAGAATCATGCGCGATTGCAGCAGTTGGATCTCCGGTGCAGa                                                                                               <  7:309474/149‑1 (MQ=255)
                                                     cTCCCCTGGTTTTTCTTTGGTTAATCTCGCCCAGCCGCGACCCACAATCGGAAAATACTTCTGCTCAACTATGTCGCGCAGATTCAGTTCAGCAATGGTTTTACCGGGATTCTTGCGGTTTGCCGGCAGTTGGATGTCGGTTGctgcgg                                                                                             >  5:272356/1‑144 (MQ=255)
                                                     cTCACCTGGTTTTTCTTTGGTTAATCTCGCCCAGCCGCGACCCACAATCGGAAAATACTTCTGCTCAACTATGTCGCGCAGATTCAGTTCAGCAATGGTTTTACCTCGAATCATGCGCGATTGCAGCAGTTGGATCTCCGGTGccgcat                                                                                             >  5:267769/1‑144 (MQ=255)
                                                        aCCTGGTTTTTCTTTGGTTAATCTCGCCCAGCCGCGACCCACAATCGGAAAATACTTCTGCTCAACTATGTCGCGCAGATTCAGTTCAGCAATGGTTTTACCGAGAATCATGCGCGATTGCAGCAGTTGGATCTCCGGTGCAGACTCgg                                                                                          <  8:41000/149‑1 (MQ=255)
                                                         ccTGGTTTTTCTTTGGTTAATCTCGCCCAGCCGCGACCCACAATCGGAAAATACTTCTGCTCAACTATGACGCGCAGATTCAGTTCAGCAATGGTTTTACCGAGAATCATGCGCGATTGCAGCAGTTGGATCTCCGGTGCAGACTCggg                                                                                         <  6:63895/149‑1 (MQ=255)
                                                          cTGGTTTTTCTTTGGTTAATCTCGCCCAGCCGCGACCCACAATCGGAAAATACTTCTGCTCAACTATGTCGCGCAGATTCAGTTCAGCAATGGTTTTACCGAGAATCATGCGCGATTGCAGCAGTTGGATCTCCGGTGCAGACTCGGGc                                                                                        >  8:170187/1‑149 (MQ=255)
                                                          cTGGTTTTTCTTTGGTTAATCTCGCCCAGCCGCGACCCACAATCGGAAAATACTTCTGCTCAACTATGTCGCGCAGATTCAGTTCAGCAATGGTTTTACCGAGAATCATGCGCGATTGCAGCAGTTGGATCTCCGGTGCAGACTCGGGc                                                                                        >  1:177516/1‑149 (MQ=255)
                                                             gTTTTTCTTTGGTTAATCTCGCCCAGCCGCGACCCACAATCGGAAAATACTTCTGCTCAACTATGTCGCGCAGATTCAGTTCAGCAATGGTTTTACCGAGAATCATGCGCGATTGCAGCAGTTGGATCTCCGGTGCGGACTCGGGCGAt                                                                                     <  7:84823/149‑1 (MQ=255)
                                                             gTTTTTCTTTGGTTAATCTCGCCCAGCCGCGACCCACAATCGGAAAATACTTCTGCTCAACTATGTCGCGCAGATTCAGTTCAGCAATGGTTTTACCGAGAATCATGCGCGATTGCAGCAGTTGGATCTCCGGTGCAGACTCGGGCGa                                                                                      <  5:50370/148‑1 (MQ=255)
                                                              tttttCTTTGGTTAATCTCGCCCAGCCGCGACCCACAATCGGAAAATACTTCTGCTCAACTATGTCGCGCAGATTCAGTTCAGCAATGGTTTTACCGAGAATCATGCGCGATTGCAGCAGTTGGATCTCCGGTGCAGACTCGGGCGATg                                                                                    <  6:225245/149‑1 (MQ=255)
                                                                    tttGGTTAATCTAGCCCAGCCGCGACCCACAATGGGAAAATACTTCTGCTCAACTATGTGGCGCAGATTCAGTTCAGCAATGGTTTTACCGAGAATCATGCGCGATTGCCGCAGTTGGATCTCCGGTGCAGACTCGGGCGATGAGttcg                                                                              <  4:180204/149‑3 (MQ=255)
                                                                        gTTAATCTCGCCCAGCCGCGACCCACAATCGGAAAATACTTCTGCTCAACTATGTCGCGCAGATTCAGTTCAGCAATGGTTTTACCGAGAATCATGCGCGATTGCAGCAGTTGGATCTCCGGTGCAGACTCGGGCGATGAGTTAGGGa                                                                           <  5:257793/148‑1 (MQ=255)
                                                                             tctcGCCCAGCCGCGACCCACAATCGGAAAATACTTCTGCTCAACTATGTCGCGCAGATTCAGTTCAGCAATGGTTTTACCGAGAATCATGCGCGATTGCAGCAGTTGGATCTCCgg                                                                                                     <  6:169372/117‑1 (MQ=255)
                                                                             tctcGCCCAGCCGCGACCCACAATCGGAAAATACTTCTGCTCAACTATGTCGCGCAGATTCAGTTCAGCAATGGTTTTACCGAGAATCATGCGCGATTGCAGCAGTTGGATCTCCgg                                                                                                     >  5:169372/1‑117 (MQ=255)
                                                                             tctcGCCCAGCCGCGACCCACAATCGGAAAATACTTCTGCTCAACTATGTCGCGCAGATTCAGTTCAGCAATGGTTTTACCGAGAATCATGCGCGATTGCAGCAGTTGGATCTCCGGTGCAGACTCGGGCGATGAGTTAGGGATCatat                                                                     >  7:345781/1‑149 (MQ=255)
                                                                               tcGCCCAGCCGCGACCCACAATCGGAAAATACTTCTGCTCAACTATGTCGCGCAGATTCAGTTCAGCAATGGTTTTACCGAGAATCATGCGCGATTGCAGCAGTTGGATCTCCGGTGCAGACTCGGGCGATGAGTTAGGGATCATATCg                                                                   >  6:227274/1‑149 (MQ=255)
                                                                                 gCCCAGCCGCGACCCACAATCGGAAAATACTTCTGCTCAACTATGTCGCGCAGATTCAGTTCAGCAATGGTTTTACCGAGAATCATGCGCGATTGCAGCAGTTGGATCTCCGGTGCAGACTCGGGCGATGAGTTAGGGATCatat                                                                     <  8:345781/145‑1 (MQ=255)
                                                                                  cccAGCCGCGACCCACAATCGGAAAATACTTCTGCTCAACTATGTCGCGCAGATTCAGTTCAGCTATGGTTTTACCGAGAATCATGCGCGATTGCAGCAGTTGGATCTCCGGTGCAGACTCGGGCGATGAGTTAGGGATCATATCGCt                                                                 >  5:211516/1‑148 (MQ=255)
                                                                                  cccAGCCGCGACCCACAATCGGAAAATACTTCTGCTCAACTATGTCGCGCAGATTCAGTTCAGCTATGGTTTTACCGAGAATCATGCGCGATTGCAGCAGTTGGATCTCCGGTGCAGACTCGGGCGATGAGTTAGGGATCATATCGCt                                                                 >  5:206625/1‑148 (MQ=255)
                                                                                       ccGCGACCCACAATCGGAAAATACTTCTGCTCAACTATGTCGCGCAGATTCAGTTCAGCAATGGTTTTACCGAGAATCATGCGCGATTGCAGCAGt                                                                                                                >  5:186773/1‑96 (MQ=255)
                                                                                       ccGCGACCCACAATCGGAAAATACTTCTGCTCAACTATGTCGCGCAGATTCAGTTCAGCAATGGTTTTACCGAGAATCATGCGCGATTGCAGCAGt                                                                                                                <  6:186773/96‑1 (MQ=255)
                                                                                             cccACAATCGGAAAATACTTCTGCTCAACTATGTCGCGCAGATTCAGTTCAGCAATGGTTTTACCGAGAATCATGCGCGATTGCAGCAGTTGGATCTCCGGTGCAGACTCGGGCGATGAGTTAGGGATCATATCGCTCAGGCCGCt                                                        >  3:16712/1‑146 (MQ=255)
                                                                                               cacaATCGGAAAATACTTCTGCTCAACTATGTCGCGCAGATTCAGTTCAGCAATGGTTTTACCGAGAATCATGCGCGATTGCAGCAGTTGGATCTCCGGTGCAGACTCGGGCGATGAGTTAg                                                                              >  2:129142/1‑122 (MQ=255)
                                                                                               cacaATCGGAAAATACTTCTGCTCAACTATGTCGCGCAGATTCAGTTCAGCAATGGTTTTACCGAGAATCATGCGCGATTGCAGCAGTTGGATCTCCGGTGCAGACTCGGGCGATGAGTTAg                                                                              <  1:129142/122‑1 (MQ=255)
                                                                                                acaATCGGAAAATACTTCTGCTCAACTATGTCGCGCAGATTCAGTTCAGCAATGGGTTTACCGAGAATCATGCGCGATTGCAGCAGTTGGATCTCCGGTGCAGACTCGGGCGATGAGTTAGGGATCATATCGCTCAGGCCGCTGAGAAt                                                  <  6:157835/149‑1 (MQ=255)
                                                                                                  aaTCGGAAAATACTTCTGCTCAACTATGTCGCGCAGATTCAGTTCAGCAATGGTTTTACCGAGAATCATGCGCGATTGCAGCAGTTGGATCTCCGGTGCAGACTCGGGCGATGAGTTAGGGATCATATCGCTCAGGCCGCTGAGAATgg                                                <  6:16991/149‑1 (MQ=255)
                                                                                                    tCGGAAAATACTTCTGATCAACTATGTCGCGCAGATTCAGTTCAGCAATGTTTTTACCGAGAATCATGCGCGATTGCAGCAGTTGGATCTCCTGTGCAGACTCTGGCGATGAGTTAGGGATCATATCGCTCAGGCCGCTGAGAATGGCg                                              <  5:325013/149‑1 (MQ=255)
                                                                                                     cGGAAAATACTTCTGCTCAACTATGTCGCGGAGATTCAGTTCAGCAATGGTTTTACGGAGAATCATGCGCGATTGGAGCTGGTGGATCTCAGGTGCAGACTCGGGCGATGAGTTAGGGATAATATCGCTCAGGCCGCTGAGCATGGCGt                                             >  4:12898/1‑149 (MQ=255)
                                                                                                          aaTACTTCTGCTCAACTATGTCGCGCAGATTCAGTTCAGCAATGGTTTTACCGAGAATCATGCGCGATTGCAGCAGTTGGATCTCCGGTTAAGACTCGGGCGATGAGTTAGGGATCATATCGCTCAGGCCGAAGAGAATGGCGTTGccc                                        >  3:17629/1‑149 (MQ=255)
                                                                                                           aTACTTCTGCTCAACTATGTCGCGCAGATGCAGTTCAGCAATGGTTTTACCGAGAATCATGCGCGATTGCAGCAGTTGGATCTCCGGTGCAGACTCGGGGGATGAGTTAGGGATCATATCGCTCAGGCCGCTGAGAATGGCGTTGCCCt                                       <  7:304958/149‑1 (MQ=255)
                                                                                                                tCTGCTCAACTATGTCGCGCAGATTCAGTTCAGGAATGGTTTTTCCGAGAATCATGCGCGATTGCAGCAGTTGGATCTCCGGTGAAGAGTCGTGCTCTGAGTTAGGGAGCTTATCGCTCAGGCCTCTTAGAATGGCGGTGCCCTGttgt                                  >  4:14353/1‑147 (MQ=255)
                                                                                                                tCTGCTCAACTATGTCGCGCAGATTCAGTTCAGCAATGGTTTTACCGAGAATCATGCGCGAGTGCAGCAGTTGGTTCTCCGGAGCCGCCTCGGGCGAGGCGTGAGGCAGCCTATCGCTGAGGCCGCTGCGCATGGCGTTGCCCTggtgt                                  >  4:9856/1‑145 (MQ=255)
                                                                                                                 cTGCTCAACTATGTCGCGCAGATTCAGTTCAGCAATGGTTTTACCGAGAATCATGCGCGATTGCAGCAGTTGGATCTCCGGTGCAGACTCGGGCGATGAGTTAGGGATCATATCGCTCAGGCCGCTGAGAATGGCGTTGCCCTGttttt                                 >  5:128009/1‑149 (MQ=255)
                                                                                                                    ctctACTATGTCGCGCAGATTCAGTTCAGCAATGGTTTTACCGAGAATCATGCGCGATTGCAGCAGTTGGATCTCCGGTGCAGACTCGGGCGATGAGTTAGGGATCATATCGCTCAGGCCGCTGAGAATGGCGTTGCCCTGTTTTTGCt                              >  2:259445/1‑149 (MQ=255)
                                                                                                                    cgcAACTATGTCGGGCCGCTTCAGGTCATGAATGGGTTTACCGAGAATCATGCGCGATTGCAGCAGTTGTATCTCCGGTGCAGACTCGGGCGATGAGTTAGGGATCATATCGCTCAGGCCGCTGAGAATGGCGTTGCCCTGTTTTTGCt                              <  6:267769/147‑1 (MQ=255)
                                                                                                                    cTCAACTATGTCGCGCAGATTCAGTTCAGCTATGGTTTTACCGAGAATCATGCGCGATTGCAGCAGTTGGATCTCCGGTGCAGACTCGGGCGATGAGTTAGGGATCATATCGCTCAGGCCGCTGAGAATGGCGTTGCCCTGTTTTTg                                <  6:206625/147‑1 (MQ=255)
                                                                                                                    cTCAACTATGTCGCGCAGATTCAGTTCAGCTATGGTTTTACCGAGAATCATGCGCGATTGCAGCAGTTGGATCTCCGGTGCAGACTCGGGCGATGAGTTAGGGATCATATCGCTCAGGCCGCTGAGAATGGCGTTGCCCTGTTTTTg                                <  6:211516/147‑1 (MQ=255)
                                                                                                                    cTCAACTATGTCGCGCAGATTCAGTTCAGCAATGGTTTTACCGAGAATCATGCGCGATTGCAGCAGTTGGATCTCCGGTGCAGACTCGGGCGATGAGTTAGGGATCATATCGCTCAGGCCGCTGAGAATGGCGTTGCCCTGTTTTTGc                               <  7:170187/148‑1 (MQ=255)
                                                                                                                    cTCAACTATGTCGCGCAGATTCAGTTCAGCAATGGTTTTACCGAGAATCATGCGCGATTGCAGCAGTTGGATCTCCGGTGCAGACTCGGGCGATGAGTTAGGGATCATATCGCTCAGGCCGCTGAGAATGGCGTTGCCCTGTTTTTGCt                              >  5:287359/1‑149 (MQ=255)
                                                                                                                       aaCTATGTCGCGCAGATTCAGTTCAGCAATGGTTTTCCCGAGAATCATGCGCGATTGCAGCAGGTGGATCTCCGGTGCAGACTCGGGCGATGAGTTAGGGATCATATCGCTCAGGCCGCTGAGAATGGCGTTGCCCTGTTTTTGCTc                             <  8:148740/147‑1 (MQ=255)
                                                                                                                       aaCTATGTCGCGCAGATTCAGTTCAGCAATGGTTTTACCGAGAATCATGCGCGATTGCAGCAGTTGGATCTCCGGTGCAGACTCGGGCGATGAGTTAGGGATCATATCTCTCAGGCCGCTGAGAATGGCGTTGCCCTGTTTTTGCTCac                           >  1:261939/1‑148 (MQ=255)
                                                                                                                       aaCTATGTCGCGCAGATTCAGTTCAGCAATGGTTTTACCGAGAATCATGCGCGATTGCAGCAGTTGGATCTCCGGTGCAGACTCGGGCGATGAGTTAGGGATCATATCGCTCAGGCCGCTGAGAATGGCGTTGCCCTGTTTTTGCTCaa                           >  2:176895/1‑149 (MQ=255)
                                                                                                                            tGTCGCGCAGATTCAGTTCAGCAATGGTTTTACCGAGAATCATGCGCGATTGCAGCAGTTGGATCTCCGGTGCAGACTCGGGCGATGAGTTAGGGATCATATCGCTCAGGCCGCTGAGAATGGCGTTGCCCTGTTTTTGCTCAACCTg                       <  7:223831/148‑1 (MQ=255)
                                                                                                                                  gcAGATTCAGTTCAGCAATGGTTTTACCGAGAATCATGCGCGATTGCAGCAGTTGGATCTCCGGTGCAGACTCGGGCGATGAGTTAGGGATCATATCGCTCAGGCCGCTGAGAATGGCGTTGCCCTGTTTTTGCTCAACCTGGACCaga                <  1:259445/149‑1 (MQ=255)
                                                                                                                                         cagttcagCAATGGTTTTACCGAGAATCATGCGCGATTGCAGCAGTTGGATCTCCGGTGCAGACTCGGGCGATGAGTTAGGGATCATATCGCTCAGGCCGCTGAGAATGGCGTTGCCCTGTTTTTGCTCAACCTGGACCAGAGTATCt          <  3:50985/148‑1 (MQ=255)
                                                                                                                                            ttcagCAATGGTTTTACCGAGAATCATGCGCGATTGCAGCAGTTGGATCTCCGGTGCAGACTCGGGCGATGAGTTAGGGATCATATCGCTCAGGCCGCTGAGAATGGCGTAGCCCTGTTTTTGCTCAACCAGGAACAGAGTATCTGCCt      >  4:214866/1‑149 (MQ=255)
                                                                                                                                               agCAATGGTTTTACCGAGAATCATGCGCGATTGCAGCAGTTGGATCTCCGGTGCAGACTCGGGCGATGAGTTAGGGATCATATCGCTCAGGCCGCTGAGAATGGCGTTGCCCTGTTTTTGCTCAACCTGGACCAGAGTATCTGCCTGAt   <  6:292346/149‑1 (MQ=255)
                                                                                                                                               agCAATGGTTTTACCGAGAATCATGCGCGATTGCAGCAGTTGGATCTCCGGTGCAGACTCGGGCGATGAGTTAGGGATCATATCGCTCAGGCCGCTGAGAATGGCGTTGCCCTGTTTTTGCTCAACCTGGACCAGAGTATCTGCCTGAt   <  1:167865/149‑1 (MQ=255)
                                                                                                                                                gCAATGGTTTTACCGAGAATCATGCGCGATTGCAGCAGTTGGATCTCCGGTGCAGACTCGGGCGATGAGTTAGGGATCATATCGCTCAGGCCGCTGAGAATGGCGTTGCCCTGTTTTTGCTCAACCTGGACCAGAGTATCTGCCTGAt   >  5:218812/1‑148 (MQ=255)
                                                                                                                                                 cAATGGTTTTACCGAGAATCATGCGCGATTGCAGCAGTTGGATCTCCGGTGCAGACTCGGGCGATGAGTTAGGGATCATATCGCTCAGGCCGCTGAGAATGGCGTTGCCCTGTTTTTGCTCACCCTGGACCAGAGTATCTGCCTGATa  >  8:223297/1‑148 (MQ=255)
                                                                                                                                                 cAATGGTTTTACCGAGAATCATGCGCGATTGCAGCAGTTGGATCTCCGGTGCAGACTCGGGCGATGAGTTAGGGATCATATCGCTCAGGCCGCTGAGAATGGCGTTGCCCTGTTTTTGCTCAACCTGGACCAGAGTATc           <  7:29826/139‑1 (MQ=255)
                                                                                                                                                 cAATGGTTTTACCGAGAATCATGCGCGATTGCAGCAGTTGGATCTCCGGTGCAGACTCGGGCGATGAGTTAGGGATCATATCGCTCAGGCCGCTGAGAATGGCGTTGCCCTGTTTTTGCTCAACCTGGACCAGAGTATc           >  8:29826/1‑139 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
AGTTGCTGATCCTGACCATTCAGTTGTGGAATATGCATCCAGCTGATCGCCAGCTCACCTGGTTTTTCTTTGGTTAATCTCGCCCAGCCGCGACCCACAATCGGAAAATACTTCTGCTCAACTATGTCGCGCAGATTCAGTTCAGCAATGGTTTTACCGAGAATCATGCGCGATTGCAGCAGTTGGATCTCCGGTGCAGACTCGGGCGATGAGTTAGGGATCATATCGCTCAGGCCGCTGAGAATGGCGTTGCCCTGTTTTTGCTCAACCTGGACCAGAGTATCTGCCTGATA  >  NC_000913/1043753‑1044045

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: