Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,158,920 G→A 7.0% M351I (ATG→ATA gtsA → mannose/glucose ABC transporter substrate‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,158,9200GA 7.0% 253.5 / 5.4 85M351I (ATG→ATAgtsAmannose/glucose ABC transporter substrate‑binding protein
Reads supporting (aligned to +/- strand):  ref base G (42/37);  new base A (4/2);  total (46/39)
Fisher's exact test for biased strand distribution p-value = 6.83e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

TACAACATCGACTCGTTGGTGGTGTTCAAGCAGAACAACGCTGGCACCTCTGCTGGTCAGCAGGACATCGCCCGCAAGGTGTTGGGTGAGGACTTCCAGAAGGTCTTCAGCATCAACAAGGGTTCGATCCCGGTACGCAACGACATGCTTGCCGACATGGGCAAGTATGGTTTCGATGCCTGTGCTCAAACCTCCGCCAAGGACTTCCTGGCTGACGCCAAAACGGGCGGCCTGCAACCGAGCATGGCGCACAACATGGCCACCACGCTGGCCGTGCAGGGCGCGTTCTTCGATG  >  NC_002947/1158774‑1159068
                                                                                                                                                  |                                                                                                                                                    
tACAACATCGACTCGTTGGTGGTGTTCAAGCAGAACAACGCTGGCACCTCTGCTGGTCAGCAGGACATCGCCCGCAAGGTGTTGGGTGAGGACTTCCAGAAGGTCTTCAGCATCAACAAGGGTTCGATCCCGGTACGCAACGACATGCt                                                                                                                                                    >  8:305407/1‑149 (MQ=255)
    acaTCGACTCGTTGGTGGTGTTCAAGCAGAACAACGCTGGCACCTCTGCTGGTCAGCAGGACATCGCCCGCAAGGTGTTGGGTGAGGACTTCCAGAAGGTCTTCAGCATCAACAAGGGTTCGATCCCGGTACGCAACGACATGCTTGcc                                                                                                                                                >  8:203278/1‑149 (MQ=255)
      aTCGACTCGTTGGTGGTGTTCAAGCAGAACAACGCTGGCACCTCTGCTGGTCAGCAGGACATCGCCCGCAAGGTGTTGGGTGAGGACTTCCAGAAGGTCTTCAGCATCAACAAGGGTTCGATCCCGGTACGCAACGACATGCTTGCCg                                                                                                                                               >  7:64634/1‑148 (MQ=255)
       tCGACTCGTTGGTGGTGTTCAAGCAGAACAACGCTGGCACCTCTGCTGGTCAGCAGGACATCGCCCGCAAGGTGTTGGGTGAGGACTTCCAGAAGGTCTTCAGCATCAACAAGGGTTCGATCCCGGTACGCAACGACATGCTTGCCGa                                                                                                                                              >  7:256700/1‑148 (MQ=255)
          aCTCGTTGGTGGTGTTCAAGCAGAACAACGCTGGCACCTCTGCTGGTCAGCAGGACATCGCCCGCAAGGTGTTGGGTGAGGACTTCCAGAAGGTCTTCAGCATCAACAAGGGTTCGATCCCGGTACGCAACGACATACTTGcccatgtc                                                                                                                                          >  6:353435/1‑143 (MQ=255)
                   tggtgTTCAAGCAGAACAACGCTGGCACCTCTGCTGGTCAGAAGGACATCGCCCGCAAGGTGTTGGGTGAGGACTTCCAGAAGGTGTTCAGCATCAACAAGGGGTCGATCCCGGTACGCAACGCCTTGCTTGCCGACATGGGCAAGTAt                                                                                                                                 >  2:351186/1‑149 (MQ=255)
                        ttCAAGCATAACAACGCTGCCACCTCTGCTGGTCAGAAGTAAATCGCCCGCAAGGTGTTGGGTTAGGAATTCCAGACGGTCTTCAGCATCAACAAGGTTACGATCCCGGTACGCAACGACATACTTGCCCATGTCGGCAAGCATgtcg                                                                                                                             <  1:286846/148‑4 (MQ=255)
                        ttCAAGCAGAACAACGCTGGCACCTCTGCTTGTCAGCAGGACATCGCCCGCAAGGTGTTGGGTGATGACTTCCAGAAGGTCTTCAGCATCAACAAGGGTTCGATCCCGGTACGCAACGACATGCTTGCCGACATTGGCAAGTATGGttt                                                                                                                            <  2:290868/149‑1 (MQ=255)
                        ttCAAGCAGAACAACGCTGGCACCTCTGCTGGTCAGCAGGACATCGCCCGCAAGGTGTTGGGTGAGGACTTCCAGAAGGTCTTCAGCATCAACAAGGGTTCGATCCCGGTACGCAACGACATGCTTGCCGACATGGGCAAGTATGGttt                                                                                                                            >  7:309325/1‑149 (MQ=255)
                        ttCAAGCAGAACAACGCTGGCACCTCTGCTGGTCAGCAGGACATCGCCCGCAAGGTGTTGGGTGAGGACTTCCAGAAGGTCTTCAGCATCAACAAGGGTTCGATCCCGGTACGCAACGACATGCTTGCCGACATGGGCAAGTATGGttt                                                                                                                            <  3:191775/149‑1 (MQ=255)
                        ttCAAGCAGAACAACGCTGGCACCTCTGCTGGTCAGCAGGACATCGCCCGCAAGGTGTTGGGTGAGGACTTCCAGAAGGTCTTCAGCATCAACAAGGGTTCGATCCCGGTACGCAACGACATACTTGCCCATGTCGGCAAGCATgtcg                                                                                                                             <  5:353435/148‑4 (MQ=255)
                        ttCAAGCAGAACAACGCTGGCACCTCTGCTGGTCAGCAGGACATCGCCCGCAACGTGTTGGGTGAGGACTTCCAGAAGGTCTTCAGCATCAACAAGCGTTCGATCCCGGTACGCAACGACATGCTTGCCGACATGCGCAAGTATCGttt                                                                                                                            >  8:195817/1‑149 (MQ=255)
                          cAAGCAGAACAACGCTGGCACCTCTGCTGGTCAGCAGGACATCGCCCGCAAGGTGGTGGGTTAGGACTTCCAGAAGGTCTTCAGCATCAACAAGGGTTCGATCGCGGTACGCAACGACATGCTTGCCGACATGGGCAAGTATGGTTTc                                                                                                                           <  7:195817/148‑1 (MQ=255)
                             gCAGAACAACGCTGGCACCTCTGCTGGTCAGCAGGACATCGCCCGCAAGGTGTTGGGTGAGGACTTCCAGAAGGTCTTCAGCATCAACAAGGGTTCGATCCCGGTACGCAACGACATGCTTGCCGACATGGGCAAGTATGGTTTCGATg                                                                                                                       >  8:59811/1‑149 (MQ=255)
                              cAGAACAACGCTGGCACCTCTGCTGGTCAGCAGGACATCGCCCGCAAGGTGTTGGGTGAGGACTTCCAGAAGGTCTTCAGCATCAACAAGGGTTCGATCCCGGTACGCAACGACATGCTTGCCGACATGGGCAAGTATGGTTTCGATGc                                                                                                                      <  1:206953/149‑1 (MQ=255)
                                gAACAACGCTGGCACCTCTGCTGGTCAGCAGGACATCGCCCGCAAGGTGTTGGGTGAGGACTTCCAGAAGGTCTTCAGCATCAACAAGGGTTCGATCCCGGTACGCAACGACATGCTTGCCGACATGGGCAAGTATGGTTTCGATGCCt                                                                                                                    <  8:64634/149‑1 (MQ=255)
                                  acaacGCTGGCACCTCTGCTGGTCAGCAGGACATCGCCCGCAAGGTGTTGGGTGAGGACTTCCAGAAGGTCTTCAGCATCAACAAGGGTTCGATCCCGGTACGCAACGACATGCTTGCCGACATGGGCAAGTATGGTTTCGATGCCtgt                                                                                                                  <  6:122602/149‑1 (MQ=255)
                                   caacGCTGGCACCTCTGCTGGTCAGCAGGACATCGCCCGCAAGGTGTTGGGTGAGGACTTCCAGAAGGTCTTCAGCATCAACAAGGGTTCGATCCCGGTACGCAACGACATGCTTGCCg                                                                                                                                               <  8:247703/119‑1 (MQ=255)
                                   caacGCTGGCACCTCTGCTGGTCAGCAGGACATCGCCCGCAAGGTGTTGGGTGAGGACTTCCAGAAGGTCTTCAGCATCAACAAGGGTTCGATCCCGGTACGCAACGACATGCTTGCCg                                                                                                                                               >  7:247703/1‑119 (MQ=255)
                                    aacGCTGGCACCTCTGCTGGTCAGCAGGACATCGCCCGCAAGGTGTTGGGTGAGGACTTCCAGAAGGTCTTCAGCATCAACAAGGGTTCGATCCCGGTACGCAACGACATGCTTGCCGACATGGGCAAGTATGGTTTCGATGCCTGTGc                                                                                                                <  3:132435/149‑1 (MQ=255)
                                    aacGCTGGCACCTCTGCTGGTCAGCAGGACATCGCCCGCAAGGTGTTGGGTGAGGACTTCCAGAAGGTCTTCAGCATCAACAAGGGTTCGATCCCGGTACGCAACGACATGCTTGCCGACATGGGCAAGTATGGTTTCGATGCCTGTGc                                                                                                                >  8:213902/1‑149 (MQ=255)
                                        cTGGCACCTCTGCTGGTCAGCAGGACATCGCCCGCAAGGTGTTGGGTGAGGACTTCCAGAAGGTCTTCAGCATCAACAAGGGTTCGATCCCGGTACGCAACGACATGCTTGCCGACATGGGCAAGTATGGTTTCGATGCCTGTGCTCaa                                                                                                            <  7:305407/149‑1 (MQ=255)
                                           gCACCTCTGCTGGTCAGCAGGACATCGCCCGCAAGGTGTTGGGTGAGGACTTCCAGAAGGTCTTCAGCATCAACAAGGGTTCGATCCCGGTACGCAACGACATGCTTGCCGACATGGGCAAGTATGGTTTCGATGCCTGTGCTCaaa                                                                                                           <  8:256700/147‑1 (MQ=255)
                                             aCCTCTGCTGGTCAGCAGGACATCGCCCGCAAGGTGTTGGGTGAGGACTTCCAGAAGGTCTTCAGCATCAACAAGGGTTCGATCCCGGTACGCAACGACATGCTTGCCGCCATGGGCAAGTATGGTTTCGATGCCTGTGCTCAAACCTc                                                                                                       >  2:167287/1‑149 (MQ=255)
                                              ccactgctgGTGAGCAGCACATCGCCCCCAAGGTGGTGGGTGAGGACTTCCAGAAGGTCTTCAGCATCAACTAGGGTTCGATCCCGGTACGCAACGACATGCTTGCCGACATTGGCAAGTATGGTTTCGATGCCTGTGCTCAAACCTcc                                                                                                      <  8:361255/146‑1 (MQ=255)
                                              ccTCTGCTGGTCAGCAGGACATCGCCCGCAAGGTGTTGGGTGAGGACTTCCAGAAGGTCTTCAGCATCAACAAGGGTTCGATCCCGGTACGCAACGACATGCTTGCCGACATGGGCAAGTATGGTTTCGATGCCTGTGCTCAAACCTcc                                                                                                      <  8:212050/149‑1 (MQ=255)
                                                     tgGTCAGCAGGACATCGCCCGCAAGGTGTTGGGTGAGGACTTCCAGAAGGTCTTCAGCATCAACAAGGGTTCGATCCCGGTACGCAACGACATGCTTGCCGACATGGGCAAGTATGGTTTCGATGCCTGTGCTCAAACCTCCGCCAAg                                                                                                >  2:347648/1‑148 (MQ=255)
                                                       gTCAGCAGGACATCGCCCGCAAGGTGTTGGGTGAGGACTTCCAGAAGGTCTTCAGCATCAACAAGGGTTCGATCCCGGTACGCAACGACATGCTTGCCGACATGGGCAAGTATGGTTTCGATGCCTGTGCTCAAACCTCCGCCAAGGAc                                                                                             >  2:381504/1‑149 (MQ=255)
                                                       gTCAGCAGGACATCGCCCGCAAGGTGTTGGGTGAGGACTTCCAGAAGGTCTTCAGCATCAACAAGGGTTCGATCCCGGTACGCAACGACATGCTTGCCGACATGGGCAAGTATGGTTTCGATGCCTGTGCTCAAACCTCCGCCAAGGAc                                                                                             >  7:151891/1‑149 (MQ=255)
                                                         cagcagGACATCGCCCGCAAGGTGTTGGGTGAGGACTTCCAGAAGGTCTTCAGCATCAACAAGGGTTCGATCCCGGTACGCAACGACATGCTTGCCGACATGGGCAAGTATGGTTTCGATGCCTGTGCTCAAACCTCCGCCAAGGACtt                                                                                           >  2:85874/1‑149 (MQ=255)
                                                          agcagGACATCGCCCGCAAGGTGTTGGGTGAGGACTTCCAGAAGGTCTTCAGCATCAACAAGGGTTCGATCCCGGTACGCAACGACATGCTTGCCGACATGGGCAAGTATGGTTTCGATGCCTGTGCTCAAACCTCCGCCAAGGACTTc                                                                                          <  6:329128/149‑1 (MQ=255)
                                                             agGACATCGCCCGCAAGGTGTTGGGTGAGGACTTCCAGAAGGTCTTCAGCATCAACAAGGGTTCGATCCCGGTACGCAACGACATGCTTGCCGACATGGGCAAGTATGGTTTCGATGCCTGTGCTCAAACCTCCGCCAAGGACTTCCTg                                                                                       <  8:309325/149‑1 (MQ=255)
                                                             agGACATCGCCAGCAAGGTGTTGGGTGAGGACTTCCAGAAGGTCTTCAGCATCAACAAGGGTTCGATCCCGGTACGCAACGACATGCTTGCCGACATGGGCAAGTATGGTTTCGATGCCTGTGCTCAAACCTCCGCCAAGGACTTCCTg                                                                                       <  3:282896/149‑1 (MQ=255)
                                                                   tCGCCCGCAAGGTGTTGGGTGAGGACTTCCAGAAGGTCTTCAGCATCAACAAGGGTTCGATCCCGGTACGCAACGACATGCTTGCCGACATGGGCAAGTATGGTTTCGATGCCTGTGCTCAAACCTCCGCCAAGGACTTCCTGGCTGAc                                                                                 >  8:311236/1‑149 (MQ=255)
                                                                   tCGCCCGCAAGGTGTTGGGTGAGGACTTCCAGAAGGTCTTCAGCATCAACAAGGGTTCGATCCCGGTACGCAACGACATGCTTGCCGACATGGGCAAGTATGGTTTCGATGCCTGTGCTCAAACCTCCGCCAAGGACTTCCTGGCTGAc                                                                                 >  5:392591/1‑149 (MQ=255)
                                                                       ccGCAAGGTGTTGGGTGAGGACTTCCAGAAGGTCTTCAGCATCAACAAGGGTTCGATCCCGGTACGCAACGACATGCTTGCCGACATGGGCAAGTATGGTTTCGATGCCTGTGCGCAAACCACCGCCAAGGACTTCCTGGCTGACGcc                                                                              >  7:55853/1‑148 (MQ=255)
                                                                                   gggTGAGGACTTCCAGAAGGTCTTCAGCATCAACAAGGGTTCGATCCCGGTACGCAACGACATGCTTGCCGACATGGGCAAGTATGGTTTCGATGCCTGTGCTCAAACCTCCGCCAAGGACTTCCTGGCTGACGCCAAAACGggcggc                                                                  <  5:214792/148‑1 (MQ=255)
                                                                                   gggTGAGGACTTCCAGAAGGTCTTCAGCATCAACAAGGGTTCGATCCCGGTACGCAACGACATGCTTGCCGACATGGGCAAGTATGGTTTCGATGCCTGTGCTCAAACCTCCGCCAAGGACTTCCTGGCTGACGCCAAAACGGGCGGcc                                                                 >  8:383/1‑149 (MQ=255)
                                                                                    ggTGAGGACTTCCAGAAGGTCTTCAGCATCAACAAGGGTTCGATCCCGGTACGCAACGACATGCTTGCCGACATGGGCAAGTATGGTTTCGATGCCTGTGCTCAAACCTCCGCCAAGGACTTCCTGGCTGACGCCAAAACGGGCGGCCt                                                                >  3:185216/1‑149 (MQ=255)
                                                                                    ggTGAGGACTTCCAGAAGGTCTTCAGCATCAACAAGGGTTCGATCCCGGTACGCAACGACATGCTTGCCGACATGGGCAAGTATGGTTTCGATGCCTGTGCTCAAACCTCCGCCAAGGACTTCCTGGCTGACGCAAAAACGGGCGGCCt                                                                <  8:151891/149‑1 (MQ=255)
                                                                                       gAGGACTTCCAGAAGGTCTTCAGCATCAACAAGGTTCGATCCCGGTACGCAACGACATGCTTGCCGACATGGGCAAGTATGGTTTCGATGCCTGTGCTCAAACCTCCGCCAAGGACTTCCTGGCTGACGCCAAAACGGGCGGCCTGCaa                                                            <  5:196560/149‑1 (MQ=255)
                                                                                         ggACTTCCAGAAGGTCTTCAGCATCAACAAGGGTTCGATCCCGGTACGCAACGACATGCTTGCCGACATGGGCAAGTATGGTTTCGATGCCTGTGCTCAAACCTCCGCCAAGGACTTCCTGGCTGACGCCAAAACGGGCGGCCTGCAAc                                                           >  1:233114/1‑149 (MQ=255)
                                                                                            cTTCCAGAAGGTCTTCAGCATCAACAAGGGTTCGATCCCGGTACGCAACGACATGCTTGCCGACATGGGCAAGTATGGTTTCGATGCCTGTGCTCAAACCTCCGCCAAGGACTTCCTGGCTGACGCCAAAAc                                                                         <  5:104545/132‑1 (MQ=255)
                                                                                            cTTCCAGAAGGTCTTCAGCATCAACAAGGGTTCGATCCCGGTACGCAACGACATGCTTGCCGACATGGGCAAGTATGGTTTCGATGCCTGTGCTCAAACCTCCGCCAAGGACTTCCTGGCTGACGCCAAAAc                                                                         >  6:104545/1‑132 (MQ=255)
                                                                                            cTTCCAGAAGGTCTTCAGCATCAACAAGGGTTCGATCCCGGTACGCAACGACATGCTTGCCGACATGGGCAAGTATGGTTTCGATGCCTGTGCTCAAACCTCCGCCAAGGACTTCCTGGCTGACGCCAAAACGGGCGGCCTGCAAcc                                                          >  1:115957/1‑147 (MQ=255)
                                                                                            cTTCCAGAAGGTCTTCAGCATCAACAAGGGTTCGATCCCGGTACGCAACGACATGCTTGCCGACATGGGCAAGTATGGTTTCGATGCCTGTGCTCAAACCTCCGCCAAGGACTTCCTGGCTGACGCCAAAACGGGCGGCCTGCAACCGa                                                        <  2:115957/149‑1 (MQ=255)
                                                                                             ttCCAGAAGGTCTTCAGCATCAACAAGGGTTCGATCCCGGTACGCAACGACATGCTTGCCGACATGGGCAAGTATGGTTTCGATGCCTGTGCTCAAACCTCCGCCAAGGACTTCCTGGCTGACGCCAAAACGGGCGGCCTGCAACCGAg                                                       >  6:332358/1‑149 (MQ=255)
                                                                                               ccAGAAGGTCTTCAGCATCAACAAGGGTTCGATCCCGGTACGCAACGACATGCTTGCCGACATGGGCAAGTATGGTTTCGATGCCTGTGCTCAAACCTCCGCCAAGGACTTCCTGGCTGACGCCAAAACGGGCGGCCTGCAACCGAGCa                                                     >  8:246452/1‑149 (MQ=255)
                                                                                                  gAAGGTCTTCAGCATCAACAAGGGTTCGATCCCGGTACGCAACGACATGCTTGCCGACATGGGCAAGTATGGTTTCGATGCCTGTGCTCAAACCTCCGCCAAGGACTTCCTGGCTGACGCCAAAACGGGCGGCCTGCAACCGAGCATgg                                                  <  1:347648/149‑1 (MQ=255)
                                                                                                     ggTCTTCAGCATCAACAAGGGTTCGATCCCGGTACGCAACGACATGCTTGCCGACATGGGCAAGTATGGTTTCGATGCCTGTGCTCAAACCTCCGCCAAGGACTTCCTGGCTGACGCCAAAACGGGCGGCCTGCAACCGAGCATGgcg                                                >  5:375056/1‑148 (MQ=255)
                                                                                                       tCTTCAGCATCAACAAGGGTTCGATCCCGGTACGCAACGACATGCTTGCCGACATGGGCAAGTATGGTTTCGATGCCTGTGCTCAAACCTCCGCCAAGGACTTCCTgg                                                                                      >  6:154159/1‑108 (MQ=255)
                                                                                                       tCTTCAGCATCAACAAGGGTTCGATCCCGGTACGCAACGACATGCTTGCCGACATGGGCAAGTATGGTTTCGATGCCTGTGCTCAAACCTCCGCCAAGGACTTCCTgg                                                                                      <  5:154159/108‑1 (MQ=255)
                                                                                                       tCTTCAGCATCAACAAGGGTTCGATCCCGGTACGCAACGACATGCTTGCCGACATGGGCAAGTATGGTTTCGATGCCTGTGCTCAAACCTCCGCCAAGGACTTCCTGGCTGACGCCAAAACGGGCGGCCTGCAACCGAGCATGGCGca                                              >  3:84092/1‑148 (MQ=255)
                                                                                                        cTTCAGCATCAACAAGGGTTCGATCCCGGTACGCAACGACATGCTTGCCGACATGGGCAAGTATGGTTTCGATGCCTGTGCTCAAACCTCCGCCAAGGACTTCCTGGCTGACGCCAAAACGGGCGTCCTGCAACCGAGCATGGCGcaca                                            <  7:59811/149‑1 (MQ=255)
                                                                                                         ttCAGCATCAACAAGGGTTCGATCCCGGTACGCAACGACATGCTTGCCGACATGGGCAAGTATGGTTTCGATGCCTGTGCTCAAACCTCCGCCAAGGACTTCCTGGCTGACGCCAAAACGGGCGGCCTGCAACCGAGCATGGCGCacaa                                           >  6:104693/1‑149 (MQ=255)
                                                                                                             gCATCAACAAGGGTTCGATCCCGGTACGCAACGACATGCTTGCCGACATGGGCAAGTATGGTTTCGATGCCTGTGCTCAAACCTCCGCCAAGGACTTCCTGGCTGACGCCAAAACGGGCGGCCTGCAACCGAGCATGGCGCACAACATg                                       >  1:92097/1‑149 (MQ=255)
                                                                                                               aTCAACAAGGGTTCGATCCCGGTACGCAACGACATGCTTGCCGACATGGGCAAGTATGGTTTCGATGCCTGTGCTCAAACCTCCGCCAAGGACTTCCTGACTGACGCCAAAACGGGCGGCCTGCAACCGAGCATGGCGCACAACATGGc                                     >  6:294125/1‑149 (MQ=255)
                                                                                                                tCAACAAGGGTTCGATCCCGGTACGCAACGACATGCTTGCCGACATGGGCAAGTATGGTTTCGATGCCTGTGCTCAAACCTCCGCCAAGGACTTCCTGGCTGACGCCAAAACGGGCGGCCTGCAACCGAGCATGGCGCACAACATGGcc                                    <  2:233114/149‑1 (MQ=255)
                                                                                                                tCAACAAGGGTTCGATCCCGGTACGCAACGACATGCTTGCCGACATGGGCAAGTATGGTTTCGATGCCTGTGCTCAAACCTCCGCCAAGGACTTCCTGGCTGACGCCAAAACGGGCGGCCTGCAACCGAGCATGGCGCACAACATGGcc                                    <  5:104693/149‑1 (MQ=255)
                                                                                                                tCAACAAGGGTTCGATCCCGGTACGCAACGACATGCTTGCCGACATGGGCAAGTATGGTTTCGATGCCTGTGCTCAAACCTCCGCCAAGGACTTCCTGGCTGACGCCAAAACGGGCGGCCTGCAACCGAGCATGGCGCACAACATGGc                                     >  3:138006/1‑148 (MQ=255)
                                                                                                                   gcaaGGGTGCGATCCCGGTAGGCAGCGACATGCTTGCCGACATGGGGATGTATGGTTTCGATGCCTGTGCTGAGACCTCCGCCAGGGAGTTCCTGGCTGACGCCAAGCCGGGCGGCCTGCAACCGAGCATGGCGGACAACATGGccac                                  >  4:54518/2‑148 (MQ=255)
                                                                                                                   acaaGGGTTCGATCCCGGTACGCAACGACATGCTTGCCGACATGGGCAAGTATGGTTTCGATGCCTGTGCTCAAACCTCCGCCAAGGACTTCCTGGCTGACGCCAAAACGGGCGGCCTGCAACCGAGCATGGCGCACAACATGGccacc                                 >  7:252544/1‑149 (MQ=255)
                                                                                                                   acaaGGGTTCGATCCCGGTACGCAACGACATGCTTGCCGACATGGGCAAGTATGGTTTCGATGCCTGTGCTCAAACCTCCGCCAAGGACTTCCTGGCTGACGCCAAAACGGGCGGCCTGCAACCGAGCATGGCGCACAACATGGccac                                  >  8:213903/1‑148 (MQ=255)
                                                                                                                     aaGGGTTCGATCCCGGTACGCAACGACATGCTTGCCGACATGGGCAAGTATGGTTTCGATGCCTGTGCTCAAACCTCCGCCAAGGACTTCCTGGCTGACGCCAAAACGGGCGGCCTGCAACCGAGCATGGCGCACAACATGGCcaccac                               >  3:337423/1‑149 (MQ=255)
                                                                                                                        ggTTCGATCCCGGTACGCAACGACATGCTTGCCGACATGGGCAAGTATGGTTTCGATGCCTGTGCTCAAACCTCCGCCAAGGACTTCCTGGCTGACGCCAAAACGGGCGGCCTGCAACCGAGCATGGCGCACAACATGGCCACCACGCt                            >  2:111149/1‑149 (MQ=255)
                                                                                                                            cGATCCCGGTACGCAACGACATGCTTGCCGACATGGGCAAGTATGGTTTCGATGCCTGTGCTCAAACCTCCGCCAAGGACTTCCTGGCTGACGCCAAAACGGGCGGCCTGCAACCGAGCATGGCGCACAACATGGCCACCACGCTGGc                         <  7:213902/148‑1 (MQ=255)
                                                                                                                             gATCCCGGTACGCAACGACATGCTTGCCGACATGGGCAAGTATGGTTTCGATGCCTGTGCTCAAACCTCCGCCAAGGACTTCCTGGCTGACGCCAAAACGGGCGGCCTGCAACCGAGCATGGCGCACAACATGGCCACCACGCTGGcc                        <  8:252544/148‑1 (MQ=255)
                                                                                                                             gATCCCGGTACGCAACGACATGCTTGCCGACATGGGCAAGTATGGTTTCGATGCCTGTGCTCAAACCTCCGCCAAGGACTTCCTGGCTGACGCCAAAACGGGCGGCCTGCAACCGAGCATGGCGCACAACATGGCCACCACGCTGGcc                        <  6:375056/148‑1 (MQ=255)
                                                                                                                              aTCCCGGTACGCAACGACATGCTTGCCGACATGGGCAAGTATGGTTTCGATGCCTGTGCTCAAACCTCCGCCAAGGACTTCCTGGCTGACGCCAAAACGGGCGGCCTGCAACCGAGCATGGCGCACAACATGGCCACCACGCTGGCCGt                      <  7:383/149‑1 (MQ=255)
                                                                                                                                 ccGGTACGCAACGACATGCTTGCCGACATGGGCAAGTATGGTTTCGATGCCTGTGCTCAAACCTCCGCCAAGGACTTCCTGGCTGACGCCAAAACGGGCGGCCTGCAACCGAGCATGGCGCACAACATGGCCACCACGCTGGCCGTGCa                   <  2:92097/149‑1 (MQ=255)
                                                                                                                                 ccGGTACGCAACGACATGCTTGCCGACATGGGCAAGTATGGTTTCGATGCCTGTGCTCAAACCTCCGCCAAGGACTTCCTGGCTGACGCCAAAACGGGCGGCCTGCAACCGAGCATGGCGCACAACATGGCCACCACGCTGGCCGTGCa                   <  4:185216/149‑1 (MQ=255)
                                                                                                                                  cGGTACGCAACGACATGCTTGCCGACATGGGCAAGTATGGTTTCGATGCCTGTGCTCAAACCTCCGCCAAGGACTTCCTGGCTGACGCCAAAACGGGCGGCCTGCAACCGAGCATGGCGCACAACATGGCCACCACGCTGGCCGTGCAg                  >  4:10139/1‑149 (MQ=255)
                                                                                                                                  cGGTACGCAACGACATGCTTGCCGACATGGGCAAGTATGGTTTCGATGCCTGTGCTCAAACCTCCGCCAAGGACTTCCTGGCTGACGCCAAAACGGGCGGCCTGCAACCGAGCATGGCGCACAACATGGCCACCACGCTGGCCGTGCAg                  >  2:181441/1‑149 (MQ=255)
                                                                                                                                     tACGCAACGACATGCTTGCCGACATGGGCAAGTATGGTTTCGATGCCTGTGCTCAAACCTCCGCCAAGGACTTCCTGGCTGACGCCAAAACGGGCGGCCTGCAACCGAGCATGGCGCACAACATGGCCACCACGCTGGCCGTGCAGGgc               <  5:332358/149‑1 (MQ=255)
                                                                                                                                        gCAACGACATGCTTGCCGACATGGGCAAGTATGGTTTCGATGCCTGTGCTCAAACCTCCGCCAAGGACTTCCTGGCTGACGCCAAAACGGGCGGCCTGCAACCGAGCATGGCGCACAACATGGCCACCACGCTGGCCGTGCAGGGcgcg            <  7:213903/149‑1 (MQ=255)
                                                                                                                                         cAACGACATGCTTGCCGACATGGGCAAGTATGGTTTCGATGCCTGTGCTCAAACCTCCGCCAAGGACTTCCTGGCTAACGCCAAAACGGGCGGCCTGCAACCGAGCATGGCGCACAACATGGCCACCACGCTGGCCGTGCAGGGCGCGt           >  4:54603/1‑149 (MQ=255)
                                                                                                                                          aaCGACATGCTTGCCGACATGGGCAAGTATGGTTTCGATGCCTGTGCTCAAACCTCCGCCAAGGACTTCCTGGCTGACGCCAAAACGGGCGGCCTGCAACCGAGCATGGCGCACAACATGGCCACCACGCTGGCCGTGCAGGGCGCGtt          >  7:205374/1‑149 (MQ=255)
                                                                                                                                           aCGACATGCTTGCCGACATGGGCAAGTATGGTTTCGATGCCTGTGCTCAAACCTCCGCCAAGGACTTCCTGGCTGACGCCAAAACGGGCGGCCTGCAACCGAGCATGGCGCACAACATGGCCACCACGCTGGCCGTGCAGGGCGCGttc         <  6:392591/149‑1 (MQ=255)
                                                                                                                                             gACATGCTTGCCGACAGGGGCAAGTATGGTTTCGATGCCTGTGCTCAAACCTCCGCCAAGGACTTCCTGGCTGACGCCAAAACGGGCGGCCTGCAACCGAGCATGGCGCACAACATGGCCACCACGCTGGCCGTGCAGGGCGCGttctt       >  7:331019/1‑149 (MQ=255)
                                                                                                                                             accATACTTGCCCATGTCGGCAAGCATGGTTTCGATGCCTGTGCTCAAACCTCCGCCAAGGACTTCCTGGCTGACGCCAAAACGGGCGGCCTGCAACCGAGCATGGCGCACAACATGGCCACCACGCTGGCCGTGCAGGGCGCGttctt       >  2:18878/3‑149 (MQ=255)
                                                                                                                                               cATGCTTGCCGACATGGGCAAGTATGGTTTCGATGCCTGTGCTCAAACCTCCGCCAAGGACTTCCTGGCTGACGCCAAAACGGGCGGCCTGCAACCGAGCATGGCGCACAACATGGCCACCACGCTGGCCGTGCAGGGCGCGTTCTTCg     <  1:381504/149‑1 (MQ=255)
                                                                                                                                               cATACTTGCCCATGTCGGCAAGTATGGTTTCGATGCCTGTGCTCAAACCTCCGCCAAGGACTTCCTGGCTGACGCCAAAACGGGCGGCCTGCAACCGAGCATGGCGCACACCATGGCCACCACGCTGGCCGTGCAGGGCGCGttcttc      >  7:265859/1‑148 (MQ=255)
                                                                                                                                               cATACTTGCCCATGTCGGCAAGTATGGTTTCGATGCCTGTGCTCAAACCTCCGCCAAGGACTTCCTGGCTGACGCCAAAACGGGCGGCCTGCAACCGAGCATGGCGCACAACATGGCCACCACGCTGGCCGTGCAGGGCGCGttcttc      >  7:270883/1‑148 (MQ=255)
                                                                                                                                                 tGCTTGCCGACATGGGCAAGTATGGTTTCGATGCCTGTGCTCAAACCTCCGCCAAGGACTTCCTGGCTGACGCCAAAACGGGCGGCCTGCAACCGAGCATGGCGCACAACATGGCCACCACGCTGGCCGTGCAGGGCGCGTTCTTCGAt   <  7:311236/149‑1 (MQ=255)
                                                                                                                                                 tGCTTGCCGACATGGGCAAGTATGGTTTCGATGCCTGTGCTCAAACCTCCGCCAAGGACTTCCTGACTGACGCCAAAACGGGCGGCCTGCAACCGAGCATGGCGCACAACATGGCCACCACGCTGGCCGTGCAGGGCGCGTTCTTCGAt   <  5:294125/149‑1 (MQ=255)
                                                                                                                                                  gCTTGCCGACATGGGCAAGTATGGTTTCGATGCCTGTGCTCAAACCTCCGCCAAGGACTTCCTGGCTGACGCCAAAACGGGCGGCCTGCAACCGAGCATGGCGCACAACATGGCCACCACGCTGGCCGTGCAGGGCGCGTTCTTCGAtg  >  3:140723/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
TACAACATCGACTCGTTGGTGGTGTTCAAGCAGAACAACGCTGGCACCTCTGCTGGTCAGCAGGACATCGCCCGCAAGGTGTTGGGTGAGGACTTCCAGAAGGTCTTCAGCATCAACAAGGGTTCGATCCCGGTACGCAACGACATGCTTGCCGACATGGGCAAGTATGGTTTCGATGCCTGTGCTCAAACCTCCGCCAAGGACTTCCTGGCTGACGCCAAAACGGGCGGCCTGCAACCGAGCATGGCGCACAACATGGCCACCACGCTGGCCGTGCAGGGCGCGTTCTTCGATG  >  NC_002947/1158774‑1159068

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: