Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,705,971 C→A 26.2% intergenic (‑181/+127) dppB ← / ← dppA dipeptide/heme ABC transporter permease/dipeptide/heme ABC transporter periplasmic binding protein; dipeptide chemotaxis receptor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,705,9710CA26.2% 492.9 / 173.7 217intergenic (‑181/+127)dppB/dppAdipeptide/heme ABC transporter permease/dipeptide/heme ABC transporter periplasmic binding protein; dipeptide chemotaxis receptor
Reads supporting (aligned to +/- strand):  ref base C (85/75);  new base A (49/8);  total (134/83)
Fisher's exact test for biased strand distribution p-value = 6.73e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.63e-01

Only 100 of 226 total aligned reads displayed.

GGGAATGAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCT  >  NC_000913/3705873‑3706070
                                                                                                  |                                                                                                           
gggAATGAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATgcg                                                                                                            >  1:1343904/1‑100 (MQ=255)
  gAATGAAGTCTGGCTGCCGGATGGCGCTGCACGCAGTGGCACAGCGCCACCAGGTACACAGATTGTTGCTAATGTGCGCTTTGTTTATGCCGGATgcg                                                                                                            <  1:2110935/98‑1 (MQ=37)
  gAATGAAGTATTGCTGCCGGATGGCGCTGCGCGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATgcggc                                                                                                          <  2:2419228/100‑1 (MQ=255)
         gTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACg                                                                                                   >  1:1028113/1‑100 (MQ=255)
         gTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACg                                                                                                   <  2:2344140/100‑1 (MQ=255)
           aTTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGcc                                                                                                 <  2:3489775/100‑1 (MQ=255)
            ttGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCt                                                                                                <  2:1789875/100‑1 (MQ=255)
             tgctgcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCtt                                                                                               <  1:3238809/100‑1 (MQ=255)
               ctgcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTAt                                                                                             <  2:593103/100‑1 (MQ=255)
               ctgcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTAt                                                                                             >  2:2568847/1‑100 (MQ=255)
                 gcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATcc                                                                                           >  2:1485087/1‑100 (MQ=255)
                    ggATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATAAGGCGTTAACGCCGCATCCGGc                                                                                        <  2:247075/100‑1 (MQ=37)
                    ggATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGGATACGGATTGTTGACAAATTGCGCGTTGTTTATGCCGGATGCGGCGTTAACGCTTTATCCgg                                                                                         >  1:2311048/1‑99 (MQ=25)
                     gATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGcc                                                                                       <  2:3165294/100‑1 (MQ=255)
                       tGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTa                                                                                     <  2:1193575/100‑1 (MQ=255)
                       tGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGCTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCTTa                                                                                     >  1:1599127/1‑100 (MQ=255)
                        ggCGCTGCACGCAGAGGCACAGCGCCACAGGGTATACAGATTGTTGACAAAGTGCGCCTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTa                                                                                     <  1:2411340/99‑1 (MQ=255)
                         gcgcTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑a                                                                                  <  1:1800330/100‑1 (MQ=255)
                          cgcTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑aa                                                                                 >  1:1218937/1‑100 (MQ=255)
                          cgcTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑aa                                                                                 >  1:3870817/1‑100 (MQ=255)
                              gCACGCAGAGGCGCAGCGCCACCGGGTATACAGATTGTGGGCAACGTGCGCTCCGTTTATGCCGGATGCGGCGTTAACGCCTTCTCCGGCCTAC‑AAAATc                                                                             >  1:1184230/1‑100 (MQ=21)
                              gCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATAAGGCGTTAACGCCGCATCCGGcataaacaaag                                                                              <  2:415097/100‑11 (MQ=21)
                                aCGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGt                                                                           >  1:2882100/1‑100 (MQ=255)
                                 cGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTg                                                                          >  1:2087856/1‑100 (MQ=255)
                                 cGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTg                                                                          <  1:3354362/100‑1 (MQ=255)
                                 cGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTg                                                                          <  1:242753/100‑1 (MQ=255)
                                     gagGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCaaa                                                                      >  2:1742836/1‑100 (MQ=255)
                                       ggCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCaaa                                                                      >  2:624076/1‑98 (MQ=255)
                                       ggCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAAtt                                                                    >  1:3179000/1‑100 (MQ=255)
                                        gCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCGGGATGCGGCGTTAACGCCTTATCCAGCCTAC‑AAAATCGTGCAAATTc                                                                   >  1:2618647/1‑100 (MQ=255)
                                        gCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTc                                                                   <  2:1085894/100‑1 (MQ=255)
                                        gCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTc                                                                   <  2:3398772/100‑1 (MQ=255)
                                         cacaGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCa                                                                  >  2:947064/1‑100 (MQ=255)
                                           caGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAAt                                                                >  1:2222715/1‑100 (MQ=255)
                                            aGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATa                                                               >  2:887035/1‑100 (MQ=255)
                                             gcgcCACCGGGTATACAGATTGTTGACAAAGTGTGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATAc                                                              >  1:1660535/1‑100 (MQ=255)
                                             gcgcCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTACAAATTCAATa                                                               <  1:2250118/99‑1 (MQ=255)
                                               gcCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACAt                                                            >  2:238435/1‑100 (MQ=255)
                                               gcCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACAt                                                            >  2:693932/1‑100 (MQ=255)
                                                ccACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACAtt                                                           >  1:574814/1‑100 (MQ=255)
                                                    cGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑g                                                      <  1:3087587/100‑1 (MQ=255)
                                                     gggTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑Ga                                                     >  2:3861235/1‑100 (MQ=255)
                                                           aCAGATTGTTGACAACGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAACTCGTGCAAATTCAATACATTGCA‑GAAT‑T‑‑‑TTc                                           >  1:2583587/1‑100 (MQ=255)
                                                            tccggcataaaCAAAGCGCACTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑T‑‑‑TTC‑‑g                                        >  2:3114116/11‑100 (MQ=17)
                                                            cAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑T‑‑‑TTC‑‑g                                        >  2:3746948/1‑100 (MQ=255)
                                                            cAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑T‑‑‑TTC‑‑g                                        >  2:2095549/1‑100 (MQ=255)
                                                            cAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑T‑‑‑TTC‑‑g                                        <  1:1987050/100‑1 (MQ=255)
                                                                  gttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑T‑‑‑TTC‑‑GTAGGcc                                  <  2:2802649/100‑1 (MQ=255)
                                                                  gttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑T‑‑‑TTC‑‑GTAGGcc                                  <  1:3731034/100‑1 (MQ=255)
                                                                  gttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑T‑‑‑TTC‑‑GTAGGcc                                  <  2:3608908/100‑1 (MQ=255)
                                                                  gttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑T‑‑‑TTC‑‑GTAGGcc                                  >  2:2157340/1‑100 (MQ=255)
                                                                   ttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑T‑‑‑TTC‑‑GTAGGCCt                                 >  2:3259422/1‑100 (MQ=255)
                                                                    tgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑T‑‑‑TTC‑‑GTAGGCCTg                                >  2:1398521/1‑100 (MQ=255)
                                                                    tgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑T‑‑‑TTC‑‑GTAGGCCTg                                >  2:3683279/1‑100 (MQ=255)
                                                                     gACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑T‑‑‑TTC‑‑GTAGGCCTGa                               >  1:2017776/1‑100 (MQ=255)
                                                                      aCAAAGTGCGCGTTGTTTATGCCGGATGCGCCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTCCCAATTCAATACATTGCA‑GAAT‑T‑‑‑TTC‑‑TTCGGCCTGAt                              >  1:839425/1‑100 (MQ=17)
                                                                       cAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACACCTTATCCGTCCTAC‑AACATCGTGCAAATTCAATCCATTGCA‑GCAT‑T‑‑‑TTC‑‑GTAGTCCTGATa                             >  1:300540/1‑100 (MQ=17)
                                                                        aaaGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑T‑‑‑TTC‑‑GTAGGCCTGATaa                            >  1:2923977/1‑100 (MQ=255)
                                                                          aGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGc                          >  1:2996595/1‑100 (MQ=255)
                                                                            tGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCgg                        <  2:185404/100‑1 (MQ=255)
                                                                            tGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCgg                        <  1:2907061/100‑1 (MQ=255)
                                                                             gcgcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGa                       >  1:3663284/1‑100 (MQ=255)
                                                                              cgcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAg                      <  2:1914966/100‑1 (MQ=255)
                                                                              cgcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAg                      <  2:2897866/100‑1 (MQ=255)
                                                                              cgcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAg                      >  2:4818/1‑100 (MQ=255)
                                                                                cTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAgcg                    >  1:939811/1‑100 (MQ=255)
                                                                                cTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAgcg                    >  2:757969/1‑100 (MQ=255)
                                                                                 tttGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAgcgc                   >  1:1466270/1‑100 (MQ=255)
                                                                                   tGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATGCAATACATTGCA‑GAAT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCAt                 <  1:1359931/100‑1 (MQ=255)
                                                                                    ttgtaggCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAATCGTGCTAATTCAATATATTGCA‑GAAA‑C‑‑‑CTT‑‑GTAGGCCTGATAAGCGTAGCGCATc                >  1:650746/7‑100 (MQ=18)
                                                                                    gTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATc                >  2:834986/1‑100 (MQ=255)
                                                                                     tgtaggCCGGATAAGGCGTTCACGCCGCATCCGGCATAC‑‑AAATCATGCAAATCCAAAATATTGCA‑GGAG‑T‑‑‑TAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAg              >  2:3032693/6‑100 (MQ=18)
                                                                                      ttATGCCGGATGCGGCGTTAACGCCTTACCCGGCCTAT‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAg              <  1:3432207/100‑1 (MQ=255)
                                                                                        aggCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAATCGTGCAAATTCAACATATTGCA‑ATTC‑T‑‑‑CTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGc            >  1:3112356/3‑100 (MQ=18)
                                                                                        aggCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAAATGTGCAAATTCAATAAATTGC‑‑‑AAT‑TCAACTT‑‑GTAGGCCTGATAAGCGCAGCGCATCAgg             >  2:2656550/3‑100 (MQ=21)
                                                                                        aggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAATCGCTTAAATTCAATATATTGCA‑CCAT‑C‑‑‑CTC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGc            >  2:1603765/3‑100 (MQ=21)
                                                                                        aggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAAACCTGCAAATTCAATATATTGCA‑GGAT‑A‑‑‑CGT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGc            >  2:2740506/3‑100 (MQ=18)
                                                                                        aggCCGGATAAGGCGTTCACGCCGCATCCGGCATACAAAAATCTTGCCAATTCAATATATTGCA‑GGAC‑C‑‑‑AT‑‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGc            >  1:1100467/3‑100 (MQ=18)
                                                                                        aggCCGGATAAGGCGTTAACGCCGCATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGc            >  2:2069391/3‑100 (MQ=18)
                                                                                         ggCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAATTC‑TGCAAATTCAACAAATTGC‑‑GAAT‑CC‑‑CTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAt          >  1:434921/2‑100 (MQ=18)
                                                                                         ggCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAATCGTGCAAATTCAACATATTGC‑‑‑AAT‑T‑‑‑CTCTTGTAGGCCTGATAAGCGTAGCGCATCAGGct           >  1:206745/2‑99 (MQ=18)
                                                                                         ggCCGGATAAGGCGTTCACTCCGCATCCGGCCTAC‑AAAAGCACGCAAATTCAATATATTGCA‑GAGA‑T‑‑‑TGC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa           >  2:3296174/2‑100 (MQ=14)
                                                                                         ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAATCTTGCCAATTCAATATATTGCA‑GGAC‑C‑‑‑AT‑‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa           >  2:2067020/2‑100 (MQ=18)
                                                                                         ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAGTCGTATAAATTCAATGTATTGCA‑GTTT‑A‑‑‑TTG‑‑GTAGGCCTGATAAGCGCAGCGCATCAGGCa           >  2:2794303/2‑100 (MQ=25)
                                                                                         ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAATTGTGTAAATTCAATAAATTGAT‑CGAT‑C‑‑‑CAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa           >  1:2770704/2‑100 (MQ=21)
                                                                                         ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAACCTTGCTAATTCAATATATTGCAGGGAC‑T‑‑‑AT‑‑‑GTAGGCCTGATAAGCATAGCGCATCAGGCa           >  2:1209031/2‑100 (MQ=18)
                                                                                          gCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCAt          <  2:154352/100‑1 (MQ=255)
                                                                                          gCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAATTC‑TGCAAATTCAACAAATTGC‑‑GAATCC‑‑‑CTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAtt         >  2:2177810/1‑100 (MQ=18)
                                                                                          gCCGGATAAGGCGTTCACTCCGCATCCGGCCTAC‑AAAAGCACGCAAATTCAATATATTGCA‑GAGA‑T‑‑‑TGC‑‑GTAGGCCTGATAAGCGTAGCGCATCAgg             >  1:880512/1‑97 (MQ=11)
                                                                                          gCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAATCATGAAAATTCAATACATTGCA‑AGAT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCac          >  1:3042095/1‑99 (MQ=17)
                                                                                          gCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAACTTTGCAAATTCAATATATTGC‑‑‑AAT‑C‑‑‑TCCGTGTAGGCCTGATAAGCGTAGCGCATCAGGCaa          >  2:1222891/1‑99 (MQ=18)
                                                                                          gCCGGATAAGGCGTTAACGCCGCATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCAt          >  1:1761181/1‑100 (MQ=18)
                                                                                          gCCGGATAAGGCGTTAACGCCGCATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCAt          >  1:3835333/1‑100 (MQ=17)
                                                                                               aTGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttg     <  1:2451459/100‑1 (MQ=255)
                                                                                               aTGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttg     <  1:1268471/100‑1 (MQ=255)
                                                                                               aTGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttg     <  1:2886082/100‑1 (MQ=255)
                                                                                                tGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgc    >  2:2250515/1‑100 (MQ=255)
                                                                                                 gcggcgTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcc   >  1:2742838/1‑100 (MQ=255)
                                                                                                 gcggcgTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcc   <  2:232836/100‑1 (MQ=255)
                                                                                                  cggcgTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcct  >  2:2692477/1‑100 (MQ=255)
                                                                                                  |                                                                                                           
GGGAATGAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCT  >  NC_000913/3705873‑3706070

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