Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,158,881 C→T S338L (TCG→TTG)  ptsG → fused glucose‑specific PTS enzymes: IIB component/IIC component

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,158,8810CT98.9% 321.2 / ‑4.5 91S338L (TCG→TTG) ptsGfused glucose‑specific PTS enzymes: IIB component/IIC component
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base T (48/42);  total (49/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 = 9.39e-01

CCTCGTTCCTGACCGGTATCACCGAGCCGATCGAGTTCTCCTTCATGTTCGTTGCGCCGATCCTGTACATCATCCACGCGATTCTGGCAGGCCTGGCATTCCCAATCTGTATTCTTCTGGGGATGCGTGACGGTACGTCGTTCTCGCACGGTCTGATCGACTTCATCGTTCTGTCTGGTAACAGCAGCAAACTGTGGCTGTTCCCGATCGTCGGTATCGGTTATGCGATTGTTTACTACACCATCTTCCGCGTGCTGATTAAAGCACTGGATCTGAAAACGCCGG  >  NC_000913/1158737‑1159021
                                                                                                                                                |                                                                                                                                            
ccTCGTTCCTGACCGGTATCACCGAGCCGATCGAGTTCTCCTTCATGTTCGTTGCGCCGATCCTGTACATCATCCACGCGATTCTGGCAGGACTGGCATTCCCAATCTGTATTCTTCTGGGGATGCGTGACGGTACGTCGTTCTTGCAc                                                                                                                                          >  5:425949/1‑149 (MQ=255)
 cTCGTTCCTGACCGGTATCACCGAGCCGATCGAGTTCTCCTTCATGTTCGTTGCGCCGATCCTGTACATCATCCACGCGATTCTGGCAGGCCTGGCATTCCCAATCTGTATTCTTCTGGGGATGCGTGACGGTACGTCGTTCTTGCACg                                                                                                                                         >  3:197682/1‑149 (MQ=255)
 cTCGTTCCTGACCGGTATCACCGAGCCGATCGAGTTCTCCTTCATGTTCGTTGCGCCGATCCTGTACATCATCCACGCGATTCTGGCAGGCCTGGCATTCCCAATCTGTATTCTTCTGGGGATGCGTGACGGTACGTCGTTCTTGCACg                                                                                                                                         >  7:42208/1‑149 (MQ=255)
 cTCGTTCCTGACCGGTATCACCGAGCCGATCGAGTTCTCCTTCATGTTCGTTGCGCCGATCCTGTACATCATCCACGCGATTCTGGCAGGCCTGGCATTCCCAATCTGTATTCTTCTGGGGATGCGTGACGGTACGTCGTTCTTGCACg                                                                                                                                         <  5:41141/149‑1 (MQ=255)
    gTTCCTGACCGGTATCACCGAGCCGATCGAGTTCTCCTTCATGTTCGTTGCGCCGATCCTGTACATCATCCACGCGATTCTGGCAGGCCTGGCATTCCCAATCTGTATTCTTCTGGGGATGCGTGACGGTACGTCGTTCTTGCACGGTc                                                                                                                                      <  3:166114/149‑1 (MQ=255)
     ttCCTGACCGGTATCACCGAGCCGATCGAGTTCTCCTTCATGTTCGTTGCGCCGATCCTGTACATCATCCACGCGATTCTGGCAGGCCTGGCATTCCCAATCTGTATTCTTCTGGGGATGCGTGACGGTACGTCGTTCTTGCACGGTCt                                                                                                                                     >  5:303944/1‑149 (MQ=255)
        cTGACCGGTATCACCGAGCCGATCGAGTTCTCCTTCATGTTCGTTGCGCCGATCCTGTACATCATCCACGCGATTCTGGCAGGCCTGGCATTCCCAATCTGTATTCTTCTGGGGATGCGTGCGGGTACGTCGTTCTTGCACGGTCTggt                                                                                                                                  >  6:189823/1‑147 (MQ=255)
        cTGACCGGTATCACCGAGCCGATCGAGTTCTCCTTCATGTTCGTTGCGCCGATCCTGTACATCATCCACGCGATTCTGGCAGGCCTGGCATTCCCAATCTGTAGTCTGCTGGGGATGCGTGACGGTACGTCGTTCTTGCACGGTCTGAt                                                                                                                                  >  6:333945/1‑149 (MQ=255)
         tgcccGGTATCACCGAGCCGATCGAGTTCTCCTTCATGTTCGTTGCGCCGATCCTGTACATCATCCACGCGATTCTGGCAGGCCTGGCATTCCCAATCTGTATTCTTCTGGGGATGCGTGACGGTACGTCGTTCTTGCCCGGTCTGATc                                                                                                                                 >  5:381507/4‑149 (MQ=255)
         tGACCGGTATCACCGAGCCGAACGAGATCTCCTTCATGATCGTTGCGCCGATCCTGTACATCATCCACTCGAATCTGGCAGGCCTCGCATTCCCAAAATTTATTATTCTTGTGATGCGTGACGATACTTCCTTCTTGCACCGACTGAta                                                                                                                                 >  2:144949/1‑148 (MQ=255)
           aCCGGTATCACCGAGCCGATCGAGTTCTCCTTCATGTTCGTTGCGCCGATCCTGTACATCATCCACGCGATTCTGGCAGGCCTGGCATTCCCAATCTGTATTCTTCTGGGGATGCGTGACGGTACGTCGTTCTTGCACGGTCTGATc                                                                                                                                 <  6:232966/147‑1 (MQ=255)
            ccGGTATCACCGAGCCGATCGAGTTCTCCTTCATGTTCGTTGCGCCGATCCTGTACATCATCCACGCGATTCTGGCAGGCCTGGCATTCCCAATCTGTATTCTTCTGGGGATGCGTGACGGTACGTCGTTCTTGCACGGTCTGATCGAc                                                                                                                              <  2:389375/149‑1 (MQ=255)
                tATCACCGAGCCGATCGAGTTCTCCTTCATGTTCGTTGCGCCGATCCTGTACATCATCCACGCGATTCTGGCAGGCCTGGCATTCCCAATCTGTATTCTTCTGGGGATGCGTGACGGTACGTCGTTCTTGCACGGTCTGATCGACTTCa                                                                                                                          >  7:126915/1‑149 (MQ=255)
                tATCACCGAGCCGATCGAGTTCTCCTTCATGTTCGTTGCGCCGATCCTGTACATCATCCACGCGATTCTGGCAGGCCTGGCATTCCCAATCTGTATTCTTCTGGGGATGCGTGACGGTACGTCGTTCTTGCACGGTCTGATCGACTTCa                                                                                                                          >  1:412285/1‑149 (MQ=255)
                   cACCGAGCCGATCGAGTTCTCCTTCATGTTCGTTGCGCCGATCCTGTACATCATCCACGCGATTCTGGCAGGCCTGGCATTCCCAATCTGTATTCTTCTGGGGATGCGTGACGGTACGTCGTTCTTGCACGGTCTGATCGACTTCATCg                                                                                                                       >  7:24767/1‑149 (MQ=255)
                    aCCGAGCCGATCGAGTTCTCCTTCATGTTCGTTGCGCCGATCCTGTACATCATCCACGCGATTCTGGCAGGCCTGGCATTCCCAATCTGTATTCTTCTGGGGATGCGTGACGGTACGTCGTTCTTGCACGGTCTGATCGACTTCATCGt                                                                                                                      >  2:202305/1‑149 (MQ=255)
                    aCCGAGCCGATCGAGTTCTCCTTCATGTTCGTTGCGCCGATCCTGTACATCATCCACGCGATTCTGGCAGGCCTGGCATTCCCAATCTGTATTCTTCTGGGGATGCGTGACGGTACGTCGTTCTTGCACGGTCTGATCGACTTCATCGt                                                                                                                      <  4:197682/149‑1 (MQ=255)
                    aCCGAGCCGATCGAGGTCTCCTTCATGTTCGTTGCGCCGATCCTGTACATCATCCACGCGATTCTGGCAGGCCTGGCATTCCCAATCTGTATTCTTCTGGGGATGCGTGACGGTACGTCGTTCTTGCACGGTCTGATCGACTTCATCGt                                                                                                                      >  8:43984/1‑149 (MQ=255)
                     ccGAGCCGATCGAGTTCTCCTTCATGTTCGTTGCGCCGATCCTGTACATCATCCACGCGATTCTGGCAGGCCTGGCATTCCCAATCTGTATTCTTCTGGGGATGCGTGACGGTACGTCGTTCTTGCACGGTCTGATCGACTTCATCGtt                                                                                                                     <  2:166296/149‑1 (MQ=255)
                       gAGCCGATCGAGTTCTCCTTCATGTTCGTTGCGCCGATCCTGTACATCATCCACGCGATTCTGGCAGGCCTGGCATTCCCAATCTGTATTCTACTGGGGATGCGTGACGGTACGTCGTTCTCGCACGGTCTGATCGACTTCATCGAtct                                                                                                                   >  2:375880/1‑149 (MQ=255)
                          ccGATCGAGTCCCCCTTCATGTTCGGTGCGCCGATCCTGTACATCATCCACGCGATTCTGGCACGCCTTCCATTCCCCCTCTCTATTCTTCTGGGGATGCGTGACGGCACGTCGTTCTTGCACGGTCTGATCGACTTCATCGTtctgtc                                                                                                                <  2:2357/149‑1 (MQ=255)
                            gATCGAGTTCTCCTTCATGTTCGTTGCGCCGATCCTGTACATCATCCACGCGATTCTGGCAGGCCTGGCATTCCCAATCTGTATTCTTCTGGGGATGCGTGACGGTACGTCGTTCTTGCACGGTCTGATCGACTTCATCGTtctgtctg                                                                                                              >  5:158185/1‑149 (MQ=255)
                            gATCGAGTTCTCCTTCATGTTCGTTGCGCCGATCCTGTACATCATCCACGCGATTCTGGCAGGCCTGGCATTCCCAATCTGTATTCTTCTGGGGATGCGTGACGGTACGTCGTTCTTGCACGGTCTGATCGACTTCATCGTtctgtctg                                                                                                              >  1:429664/1‑149 (MQ=255)
                            gATCGAGTTCTCCTTCATGTTCGTTGCGCCGATCCTGTACATCATCCACGCGATTCTGGCAGGCCTGGCATTCCCAATCTGTATTCTTCTGGGGATGCGTGACGGTACGTCGTTCTTGCACGGTCTGATCGACTTCATCGTtctgtctg                                                                                                              >  8:179724/1‑149 (MQ=255)
                            gATCGAGTTCTCCTTCATGTTCGTTGCGCCGATCCTGTACATCATCCACGCGATTCTGGCAGGCCTGGCATTCCCAATCTGTATTCTTCTGGGGATGCGTGACGGTACGTCGTTCTTGCACGGTCTGATCGACTTCATCGTtctgtctg                                                                                                              >  3:213535/1‑149 (MQ=255)
                            gATCGAGGTCTCCTTCATGTTCGTTGCGCCGATCCTGTACATCATCCACGCGATTCTGGCAGGCCTGGCATTCCCAATCTGTCTTCTTCTGGGGATGCGTGACGGTACGTCGTTCTTGCACGGTCTGATCGACTTCATCGTtctgtctg                                                                                                              <  6:34017/149‑1 (MQ=255)
                              tataGTTTTTCTTGATTTTCGTTGCGACGATCCTGTACATCATCCACGCGATTATTGCTGGTGTGGCATTCCCAATCTGTATTCTTCTGGGGATGCGTGACGTTACGTCGTTCTTGCACGGTCTGATCGACTTCATCGTTCTGTCTgg                                                                                                             <  6:425949/145‑1 (MQ=255)
                              tCGAGTTCTCCTTCATGTTCGTTGCGCCGATCCTGTACATCATCCACGCGATTCTGGCAGGCCTGGCATTCACAATCTGTATTCTTCTGGGGATGCGTGACGGTACGTCGTTCTTGCACGGTCTGATCGACTTCATCGTTCTGTCTgg                                                                                                             >  6:191845/1‑148 (MQ=255)
                                    tctcCTTCATGTTCGTTGCGCCGATCCTGTACATCATCCACGCGATTCTGGCAGGCCTGGCATTCCCAATCTGTATTCTTCTGGGGATGCGTGACGGTACGTCGTTCTTGCACGGTCTGATCGACTTCATCGTTCTGTCTGGTAAcag                                                                                                       >  3:57917/1‑148 (MQ=255)
                                        cTTCATGTTCGTTGCGCCGATCCTGTACATCATCCACGCGATTCTGGCAGGCCTGGCATTCCCAATCTGTATTCTTCTGGGGATGCGTGACGGTACGTCGTTCTTGCACGGTCTGATCGACTTCATCGTTCTGTCTGGTAAcagca                                                                                                     <  8:54784/146‑1 (MQ=255)
                                           cATGTTCGTTGCGCCGATCCTGTACATCATCCACGCGATTCTGGCAGGCCTGGCATTCCCAATCTGTATTCTTCTGGGGATGCGTGACGGTACGTCGTTCTTGCACGGTCTGATCGACTTCATCGTTCTGTCTGGTAACAGCAGCAAAc                                                                                               >  5:139038/1‑149 (MQ=255)
                                                 cGTTGCGCCGATCCTGTACATCATCCACGCGATTCTGGCAGGCCTGGCATTCCCAATCTGTATTCTTCTGGGGATGCGTGACGGTAGTTCGTTCTTGCACGGTCTGATCGACTTCATCGTTCTGTCTGGTAACAGCAGCAAACTGTgg                                                                                          <  6:162113/148‑1 (MQ=255)
                                                   ttGCGCCGATCCTGTACATCATCCACGCGATTCTGGCAGGCCTGGCATTCCCAATCTGTATTCTTCTGGGGATGCGTGACGGTACGTCGTTCTTGCACGGTCTGATCGACTTCATCGTTCTGTCTGGTAACAGCAGCAAACTGTgg                                                                                          >  3:95954/1‑146 (MQ=255)
                                                   ttGCGCCGATCCTGTACATCATCCACGCGATTCTGGCAGGCCTGGCATTCCCAATCTGTATTCTTCTGGGGATGCGTGACGGTACGTCGTTCTTGCACGGTCTGATCGACTTCATCGTTCTGTCTGGTAACAGCAGCAAACTGTgg                                                                                          <  4:95954/146‑1 (MQ=255)
                                                     gcgcCGATCCTGTACATCATCCACGCGATTCTGGCAGGCCTGGCATTCGCAATCTGTATTCTTCTGGGGATGCGTGACGGTACGTCGTTCTTGCACGGTCTGATCGACTTCATCGTTCTGTCTGGTAACAGCAGCAAACTGTGGCTGt                                                                                      <  2:66054/148‑1 (MQ=255)
                                                      cgcCGATCCTGTACATCATCCACGCGATTCTGGCAGGCCTGGCATTCCCAATCTGTATTCTTCTGGGGATGCGTGACGGTACGTCGTTCTTGCACGGTCTGATCGACTTCATCGTTCTGTCTGGTAACAGCAGCAAACTGTGGCTGtt                                                                                     <  7:179724/148‑1 (MQ=255)
                                                      cgcCGATCCTGTACATCATCCACGCGATTCTGGCAGGCCTGGCATTCCCAATCTGTATTCTTCTGGGGATGCGTGACGGTACGTCGTTCTTGCACGGTCTGATCGACTTCATCGTTCTGTCTGGTAACAGCAGCAAACTGTGGCTGTTc                                                                                    >  7:431165/1‑149 (MQ=255)
                                                               tGTACATCATCCACGCGATTCTGGCAGGCCTGGCATTCCCAATCTGTATTCTTCTGGGGATGCGTGACGGTACGTCGTTCTTGCACGGTCTGATCGACTTCATCGTTCTGTCTGGt                                                                                                            >  8:179633/1‑116 (MQ=255)
                                                               tGTACATCATCCACGCGATTCTGGCAGGCCTGGCATTCCCAATCTGTATTCTTCTGGGGATGCGTGACGGTACGTCGTTCTTGCACGGTCTGATCGACTTCATCGTTCTGTCTGGt                                                                                                            <  7:179633/116‑1 (MQ=255)
                                                                 tACATCATCCCCGCGATTCTGGCAGGCCTGGCATTCCCAAGCTGTATTCTTCTGGGGATGCGGGACGGTACGTCGTTCTTGCACGGTCTGATCGACTTCATCGTTCTGTCTGGTAACAGCAGCAAACTGTGGCTGTTCCCGAtcgtcg                                                                          <  1:338121/148‑1 (MQ=255)
                                                                 tACATCATCCACGCGATTCTGGCAGGCCTGGCATTCCCAATCTGTATTCTTCTGGGGATGCGTGACGGTACGTCGTTCTTGCACGGTCTGATCGACTTCATCGTTCTGTCTGGTAACAGCAGCAAACTGTGGCTGTTCCCGAtcgtcg                                                                          <  8:23339/148‑1 (MQ=255)
                                                                 tACATCATCCACGCGATTCTGGCAGGCCTGGCATTCCCAATCTGTATTCTTCTGGGGATGCGTGACGGTACGTCGTTCTTGCACGGTCTGATCGACTTCATCGTTCTGTCTGGTAACAGCAGCAAACTGTGGCTGTTCCCGAtcgtcg                                                                          <  4:138823/148‑1 (MQ=255)
                                                                    atcatcCACGCGATTCTGGCAGGCCTGGCATTCCCAATCTGTATTCTTCTGGGGATGCGTGACGGTACGTCGTTCTTGCACGGTCTGATCGACTTCATCGTTCTGTCTGGTAACAGCAGCAAACTGTGGCTGTTCCCGATCGTCGGTAt                                                                      >  1:88890/1‑149 (MQ=255)
                                                                    atcatcCACGCGATTCTGGCAGGCCTGGCATTCCCAATCTGTATTCTTCTGGGGATGCGTGACGGTACGTCGTTCTTGCACGGTCTGATCGACTTCATCGTTCTGTCTGGTAACAGCAGCAAACTGTGGCTGATCCCGATCGTCGGTAt                                                                      >  3:417966/1‑149 (MQ=255)
                                                                      catcCACGCGATTCTGGCAGGCCTGGCATTCCCAATCTGTATTCTTCTGGGGATGCGTGACGGTACGTCGTTCTTGCACGGTCTGATCGACTTCATCGTTCTGTCTGGTAACAGCAGCAAACTGTGGCTGTTCCCGATCGTCGGTATCg                                                                    >  7:189911/1‑149 (MQ=255)
                                                                       atcCACGCGATTCTGGCAGGCCTGGCATTCCCAATCTGTATTCTTCTGGGGATGCGTGACGGTACGTCGTTCTTGCACGGTCTGATCGACTTCATCGTTCTGTCTGGTAACAGCAGCAAACTGTGGCTGTTCCCGATCGTCGGTATCg                                                                    <  8:42208/148‑1 (MQ=255)
                                                                       atcCACGCGATTCTGGCAGGCCTGGCATTCCCAATCTGTATTCTTCTGGGGATGCGTGACGGTACGTCGTTCTTGCACGGTCTGATCGACTTCATCGTTCTGTCTGGTAACAGCAGCAAACTGTGGCTGTTCCCGATCGTCGGTATCg                                                                    <  1:144949/148‑1 (MQ=255)
                                                                          cACGCGATTCTGGCAGGCCTGGCATTCCCAATCTGTATTCTTCTGGGGATGCGTGACGGTACGGCGTTCTTGCACGGGCTGATCGACTTCATCGTTCTGTCTGGTAACAGCAGCAAACTGTGGCTGTTCCCGATCGTCGGTAGCGGTTa                                                                <  4:98752/149‑1 (MQ=255)
                                                                            cgcgATTCTGGCAGGCCTGGCATTCCCAATCTGTATTCTTCTGGGGATGCGTGACGGTACGTCGTTCTTGCACGGTCTGATCGACTTCATCGTTCTGTCTGGTAACAGCAGCAAACtgt                                                                                            <  4:38962/119‑1 (MQ=255)
                                                                            cgcgATTCTGGCAGGCCTGGCATTCCCAATCTGTATTCTTCTGGGGATGCGTGACGGTACGTCGTTCTTGCACGGTCTGATCGACTTCATCGTTCTGTCTGGTAACAGCAGCAAACtgt                                                                                            >  3:38962/1‑119 (MQ=255)
                                                                            cgcgATTCTGGCAGGCCTGGCATTCCCAATCTGTATTCTTCTGGGGATGCGTGACGGTACGTCGTTCTTGCACGGTCTGATCGACTTCATCGTTCTGTCTGGTAACAGCAGCAAACTGTGGCTGTTCCCGATCGTCGGTATCGGTTATg                                                              <  8:126915/149‑1 (MQ=255)
                                                                            cgcgATTCTGGCAGGCCTGGCATTCCCAATCTGTATTCTTCTGGGGATGCGTGACGGTACGTCGTTCTTGCACGGTCTGATCGACTTCATCGTTCTGTCTGGTAACAGCAGCAAACTGTGGCTGTTCCCGATCGTCGGTATCGGTTATg                                                              <  8:431165/149‑1 (MQ=255)
                                                                             gcgATTCTGGCAGGCCTGGCATTCCCAATCTGTATTCTTCTGGGGATGCGTGACGGTACGTCGTTCTTGCACGGTCTGATCGACTTCATCGTTCTGTCTGGTAACAGCAGCAAACTGTGGCTGTTCCCGATCGTCGGTATCGGTTATGc                                                             >  6:27344/1‑149 (MQ=255)
                                                                             gcgATTCTGGCAGGCCTGGCATTCCCAATCTGTATTCTTCTGGGGATGCGTGACGGTACGTCGTTCTTGCACGGTCTGATCGACTTCATCGTTCTGTCTGGTAACAGCAGCAAACTGTGGCTGTTCCCGATCGTCGGTATCGGTTATGc                                                             >  4:102638/1‑149 (MQ=255)
                                                                                 ttCTGGCAGGCCTGGCATTCCCAATCTGTATTCTTCTGGGGATGCGTGACGGTACGTCGTTCTTGCACGGTCTGATCGACTTCATCGTTCTGTCTGGTAACAGCAGCAAACTGTGGCTGTTCCCGATCGTCGGTATCGGTTATGCGAtt                                                         >  1:442316/1‑149 (MQ=255)
                                                                                      gacggCCTGGCATTCACAATATCTATTATTCTGGTGACGCTGGACGGTACGTCTCTCTTGCACGGTCTGATCGACTTCCTCGTTCTGTCTGGTACCAGCAGCAAACTGTGGCTGTTCCCGATCGTCGGTCTCGGTTCTGCGATTGTTta                                                    <  2:328606/146‑1 (MQ=255)
                                                                                      gCAGGCCTGGCATTCCCAATCTGTATTCTTCTGGGGATGCGTGACGGTACGTCGTTCTTGCACGGTCTGATCGACTTCATCGTTCTGTCTGGTAACAGCAGCAAACTGTGGCTGTTCCCGATCGTCGGTATCGGTTATGCGATTGt                                                       <  8:24767/146‑1 (MQ=255)
                                                                                      gCAGGCCTGGCATTCCCAATCTGTATTCTTCTGGGGATGCGTGACGGTACGTCGTTCTTGCACGGTCTGATCGACTTCATCGTTCTGTCTGGTAACAGCAGCAAACTGTGGCTGTTCCCGATCGTCGGTATCGGTTATGCGATTGTTta                                                    >  8:262755/1‑149 (MQ=255)
                                                                                       cAGGCCTGGCATTCCCAATCTGTATTCTTCTGGGGATGCGTGACGGTACGTCGTTCTTGCACGGTCTGATCGACTTCATCGTTCTGTCTGGTAACAGCAGCAAACTGTGGCTGTTCCCGATCGTCGGTATCGGTTATGCGATTGTTtac                                                   <  7:262755/149‑1 (MQ=255)
                                                                                         ggCCTGGCATTCCCAATCTGTATTCTTCTGGGGATGCGTGACGGTACGTCGTTCTTGCACGGTCTGATCGACTTCATCGTTCTGTCTGGTAACAGCAGCAAACTGTGGCTGTTCCCGATCGTCGGTATCGGTTATGCGATTGTTtacta                                                 <  1:202305/149‑1 (MQ=255)
                                                                                         ggCCTGGCATTCCCAATCTGTATTCTTCTGGGGATGCGTGACGGTACGTCGTTCTTGCACGGTCTGATCGACTTCATCGTTCTGTCTGGTAACAGCAGCAAACTGTGGCTGTTCCCGATCGTCGGTATCGGTTATGCGATTGTTtacta                                                 >  5:225491/1‑149 (MQ=255)
                                                                                          gCCTGGCATTCCCAATCTGTATTCTTCTGGGGATGCGTGACGGTACGTCGTTCTTGCACGGTCTGATCGACTTCATCGTTCTGTCTGGTAACAGCAGCAAACTGTGGCTGTTCCCGATCGTCGGTATCGGTTATGCGATTGTTtactac                                                <  7:43984/149‑1 (MQ=255)
                                                                                            cTGGCATTCCCAATCTGTATTCTTCTGGGGATGCGTGACGGTACGTCGTTCTTGCACGGTCTGATCGACTTCATCGTTCTGTCTGGTAACAGCAGCAAACTGTGGCTGTTCCCGATCGTCGGTATCGGTTATGCGATTGTTTACTaca                                               >  7:11315/1‑148 (MQ=255)
                                                                                            cTGGCATTCCCAATCTGTATTCTTCTGGGGATGCGTGACGGTACGTCGTTCTTGCACGGTCTGATCGACTTCATCGTTCTGTCTGGTAACAACAGCAAACTGTGGCTGTTCCCGATCGTCGGTATCGGTTATGCGATTGTTTACTacac                                              >  7:60021/1‑149 (MQ=255)
                                                                                             tGGCATTCCCAATCTGTATTCTTCTGGGGATGCGTGACGGTACGTCGTTCTTGCACGGTCTGATCGACTTCATCGTTCTGTCTGGTAACAGCAGCAAACTGTGGCTGTTCCCGATCGTCGGTATCGGTTATGCGATTGTTTACTacac                                              >  1:449906/1‑148 (MQ=255)
                                                                                                cATTCCCAATCTGTATTCTTCTGGGGATGCGTGACGGTACGTCGTTCTTGCACGGTCTGATCGACTTCATCGTTCTGTCTGGTAACAGCAGCAAACTGTGGCTGTTCCCGATCGTCGGTATCGGTTATGCGATTGTTTACTACACCATc                                          >  6:413280/1‑149 (MQ=255)
                                                                                                  ctCCCCATCTGAATGCTTCTGGGGCTGCGGCACGGTCCGGCGATCCTGCCCGCTCTGCTCGCCCTCCTCGCTCAGCCTGGTCACAGCAGCAAACTCTGGCGGTTCCCGATCGGCGCTATCGGTTATGCGCTTGTTTACTACACCATCtt                                        <  6:381507/148‑1 (MQ=255)
                                                                                                   tCCCAATCTGTATTCTTCTGGGGATGCGTGACGGTACGTCGTTCTTGCACGGTCTGATCGACTTCATCGTTCTGTCTGGTAACAGCAGCAAACTGTGGCTGTTCCCGATCGTCGGTATCGGTTATGCGATTGTTTACTACACCATCTTc                                       >  6:305944/1‑149 (MQ=255)
                                                                                                     ccAATCTGTATTCTTCTGGGGATGCGTGACGGTACGTCGTTCTTGCACGGTCTGATCGACTTCATCGTTCTGTCTGGTAACCGCCGCAAACTGTGGCTGTTCCCGCTCGTCGGTATCGGTTATGCGATTGTTTACTACACCATCTTCcg                                     >  6:415167/1‑149 (MQ=255)
                                                                                                     ccAATCTGTATTCTTCTGGGGATGCGTGACGGTACGTCGTTCTTGCACGGTCTGATCGACTTCATCGTTCTGTCTGGTAACAGCAGCAAACTGTGGCTGTTCCCGATCGTCGGTATCGGTTATGCGATTGTTTACTACACCATCTTCcg                                     >  3:352483/1‑149 (MQ=255)
                                                                                                      cAATCTGTATTCTTCTGGGGATGCGTGACGGTACGTCGTTCTTGCACGGTCTGATCGACTTCATCGTTCTGTCTGGTAACAGCAGCAAACTGTGGCTGTTCCCGATCGTCGGTATCGGTTATTCGATTGTTTACTACACCATCTTCcgc                                    >  3:245138/1‑149 (MQ=255)
                                                                                                            gTATTCTTCTGGGGATGCGTGACGGTACGTCGTTCTTGCACGGTCTGATCGACTTCATCGTTCTGTCTGGTAACAGCAGCAAACTGTGGCTGTTCCCGATCGTCGGTATCGGTTATGc                                                             <  6:334477/118‑1 (MQ=255)
                                                                                                            gTATTCTTCTGGGGATGCGTGACGGTACGTCGTTCTTGCACGGTCTGATCGACTTCATCGTTCTGTCTGGTAACAGCAGCAAACTGTGGCTGTTCCCGATCGTCGGTATCGGTTATGc                                                             >  5:334477/1‑118 (MQ=255)
                                                                                                             tATTCTTCTGGGGATGCGTGACGGTACGTCGTTCTTGCACGGTCTGATCGACTTCATCGTTCTGTCTGGTAACAGCAGCAAACTGTGGCTGTTCCCGATCGTCGGTATCGGTTATGCGATTGTTTACTACACCATCTTCCGCGTGCTGa                             <  3:102638/149‑1 (MQ=255)
                                                                                                             tATTCTTCTGGGGATGCGTGACGGTACGTCGTTCTTGCACGGTCTGATCGACTTCATCGTTCTGTCTGGTAACAGCAGCAAACTGTGGCTGTTCCCGATCGTCGGTATCGGTTATGCGATTGTTTACTACACCATCTTCCGCGTGCTGa                             <  6:158185/149‑1 (MQ=255)
                                                                                                             gattattCTGGGGATGAGTGACGGTAAGTCGATCTTGCACGGTCTGATCGACTTCATCGTTCTGTCTGGTAACAGCAGCAAACTGTGGCTGTTCCCGATCGTCGGTATCTGTTATGCGATTGTTTACTACAACATCTTCCGCGTGCTGa                             >  1:424708/2‑149 (MQ=255)
                                                                                                              aTTCTTCTGGGGATGCGTGACGGTACGTCGTTCTTGCACGGTCTGATCGACTTCATCGTTCTGTCTGGTAACAGCAGCAAACTGTGGCTGTTCCCGATCGTCGGTATCGGTTATGCGATTGTTTACTACACCATCTTCCGCGTGCTGAt                            >  2:181209/1‑149 (MQ=255)
                                                                                                               ttcttcTGGGGATGCGTGACGGTACGTCGTTCTTGCACGGTCTGATCGACTTCATCGTTCTGTCTGGTAACAGCAGCAAACTGTGGCTGTTCCCGATCGTCGGTATCGGTTATGCGATTGTTTACTACACCATCTTCCGCGTGCTGAtt                           <  5:305944/149‑1 (MQ=255)
                                                                                                                 cttctGGGGATGCGTGACGGTACGTCGTTCTTGCACGGTCTGATCGACTTCATCGTTCTGTCTGGTAACAGCAGCAAACTGTGGCTGTTCCCGATCGTCGGTATCGGTTATGCGATTGTTTACTACACCATCTTCCGCGTGCTGATTa                          <  1:122764/148‑1 (MQ=255)
                                                                                                                    ctGGGGATGCGTGACGGTACGTCGTTCTTGCACGGTCTGATCGACTTCATCGTTCTGTCTGGTAACAGCAGCAAACTGTGGCTGTTCCCGATCGTCGGTATCGGTTATGCGATTGTTTACTACACCATCTTCCGCGTGCTGATTaaa                        <  4:57917/147‑1 (MQ=255)
                                                                                                                    ctGGGGATGCGTGACGGTACGTCGTTCTTGCACGGTCTGATCGACTTCATCGTTCTGTCTGGTAACAGCAGCAAACTGTGGCTGATCCCGATCGTCGGTATCGGTTATGCGATTGTTTACTACACCATCTTCCGCGTGCTGATTAAAg                       <  4:417966/148‑1 (MQ=255)
                                                                                                                       gggATGCGTGACGGTACGTCGTTCTTGCACGGTCTGATCGACTTCATCGTTCTGTCTGGTAACAGCAGCAAACTGTGGCTGTTCCCGATCGTCGGTATCGGTTATGCGATTGTTTACTACACCATCTTCCGCGTGCTGATTAAAGCACt                   <  5:27344/149‑1 (MQ=255)
                                                                                                                        ggATGCGTGACGGTACGTCGTTCTTGCACGGTCTGATCGACTTCATCGTTCTGTCTGGTAACAGCAGCAAACTGTGGCTGTTCCCGATCGTCGGTATCGGTTATGCGATTGTTTACTACACCATCTTCCGCGTGCTGATTAAAGCACTg                  <  2:412285/149‑1 (MQ=255)
                                                                                                                              gTGACGGTACGTCGTTCTTGCACGGTCTGATCGACTTCATCGTTCTGTCTGGTAACAGCAGCAAACTGTGGCTGTTCCCGATCGTCGGTATCGGTTATGCGATTGTTTACTACACCATCTTCCGCGTGCTGATTAAAGCACTGGATCt             >  4:250115/1‑148 (MQ=255)
                                                                                                                              gTGACGGTACGTCGTTCTTGCACGGTCTGATCGACTTCATCGTTCTGTCTGGTAACAGCAGCAAACTGTGGCTGTTCCCGATCGTCGGTATCGGTTATGCGATTGTTTACTACACCATCTTCCGCGTGCTGATTAAAGCACTGGATCt             >  7:263840/1‑148 (MQ=255)
                                                                                                                              gTGACGGTACGTCGTTCTTGCACGGTCTGATCGACTTCATCGTTCTGTCTGGTAACAGCAGCAAACTGTGGCTGTTCCCGATCGTCGGTATCGGTTATGCGATTGTTTACTACACCATCTTCCGCGTGCTGATTAAAGCACTGGATCt             >  1:445214/1‑148 (MQ=255)
                                                                                                                               tGACGGTACGTCGTTCTTGCACGGTCTGATCGACTTCATCGTTCTGTCTGGTAACAGCAGCAAACTGTGGCTGTTCCCGATCGTCGGTATCGGTTATGCGATTGTTTACTACACCATCTTCCGCGTGCTGATTAAAGCACTGGATCTg            >  2:427637/1‑148 (MQ=255)
                                                                                                                                   ggTACGTCGTTCTTGCACGGTCTGATCGACTTCATCGTTCTGTCTGGTAACAGCAGCAAACTGTGGCTGTTCCCGATCGTCGGTATCGGTTATGCGATTGTTTACTACACCATCTTCCGCGTGCTGATTAAAGCACTGGATCTGAAAAc       >  2:235794/1‑149 (MQ=255)
                                                                                                                                     tACGTCGTTCTTGCACGGTCTGATCGACTTCATCGTGCTGTCTGGTAACAGCAGCAAACTGTGGCTGTTCCCGATCGTCGGTATCGGTTATGCGATTGTTTACTACACCATCTTCCGCGTGCTGATTAAAGCACTGGATCTGAAAACGc     <  6:225491/149‑1 (MQ=255)
                                                                                                                                       cgtcgtTCTTGCACGGTCTGATCGACTTCATCGTTCTGTCTGGTAACAGCAGCAAACTGTGGCTGTTCCCGATCGTCGGTATCGGTTATGCGATTGTTTACTACACCATCTTCCGCGTGCTGATTAAAGCACTGGATCTGAAAACGCCg   <  1:235794/149‑1 (MQ=255)
                                                                                                                                        ggcgtTCTTGCACGGTCTGATCGACTTCATCGTTCTGTCTGGTAACAGCAGCAAACTGTGGCTGTTCCCGATCGTCGGTATCGGTTATGCGATTGTTTACTACACCATCTTCCGCGTGCTGATTAAAGCACTGGATCTGAAAACGCCgg  <  4:125664/147‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                            
CCTCGTTCCTGACCGGTATCACCGAGCCGATCGAGTTCTCCTTCATGTTCGTTGCGCCGATCCTGTACATCATCCACGCGATTCTGGCAGGCCTGGCATTCCCAATCTGTATTCTTCTGGGGATGCGTGACGGTACGTCGTTCTCGCACGGTCTGATCGACTTCATCGTTCTGTCTGGTAACAGCAGCAAACTGTGGCTGTTCCCGATCGTCGGTATCGGTTATGCGATTGTTTACTACACCATCTTCCGCGTGCTGATTAAAGCACTGGATCTGAAAACGCCGG  >  NC_000913/1158737‑1159021

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