Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,553,210 A→G K287E (AAA→GAA)  amiA → N‑acetylmuramoyl‑l‑alanine amidase I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,553,2100AG99.2% 441.0 / NA 121K287E (AAA→GAA) amiAN‑acetylmuramoyl‑l‑alanine amidase I
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base G (62/58);  minor base C (0/1);  total (62/59)
Fisher's exact test for biased strand distribution p-value = 4.88e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 122 total aligned reads displayed.

TAGGCACGGCGGCGTTTCGTCAGAAAATCGCCACAGCGATTGCTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGAAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACGATTTGTCAGCAGCTGACCGCCGTCGAT  >  NC_000913/2553115‑2553305
                                                                                               |                                                                                               
tAGGCACGGCGGCGTTTCGTCAGAAAATCGCCACAGCGATTGCTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGgaaaa                                                                                             >  1:1814543/1‑100 (MQ=255)
 aGGCACGGCGGCGTTTCGTCAGAAAATCGCCACAGCGATTGCTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGGAAAAg                                                                                            <  2:1322122/100‑1 (MQ=255)
  ggCACGGCGGCGTTTCGTCAGAAAATCGCCACAGCGATTGCTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGGAAAAGc                                                                                           <  2:1995576/100‑1 (MQ=255)
    cACGGCGGCGTTTCGTCAGAAAATCGCCACAGCGATTGCTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGGAAAAGCGa                                                                                         <  2:2541648/100‑1 (MQ=255)
    cACGGCGGCGTTTCGTCAGAAAATCGCCACAGCGATTGCTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGGAAAAGCGa                                                                                         >  1:3125785/1‑100 (MQ=255)
         cggcgTTTCGTCAGAAAATCGCCACAGCGATTGCTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGGAAAAGCGATAAGt                                                                                    <  2:2668388/100‑1 (MQ=255)
           gcgTGTCGTCAGAAAATCGCCACAGCGATTGCTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTa                                                                                  <  2:1440332/100‑1 (MQ=255)
            cgTTTCGTCAGAAAATCGCCACAGCGATTGCTGAAGGCGCGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGGAAAAGCGACAAGTTAt                                                                                 >  2:286847/1‑100 (MQ=255)
             ttccgatctgAAAATCGCCACAGCGATTGCTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATg                                                                                <  2:1188974/91‑1 (MQ=255)
             gTTTCGTCAGAAAATCGCCACAGCGATTGCTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATg                                                                                <  1:2503655/100‑1 (MQ=255)
             gTTTCGTCAGAAAATCGCCACAGCGATTGCTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATg                                                                                <  2:1703284/100‑1 (MQ=255)
             gTTTCGTCAGAAAATCGCCACAGCGATTGCTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATg                                                                                <  1:780307/100‑1 (MQ=255)
                  gTCAGAAAATCGCCACAGCGATTGCTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAAcc                                                                           <  1:3034385/100‑1 (MQ=255)
                  gTCAGAAAATCGCCACAGCGATTGCTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAAcc                                                                           >  1:2793040/1‑100 (MQ=255)
                  gTCAGAAAATCGCCACAGCGATTGCTGAAGGCGTGATCAGTTATTTCCACTGGTTCGAAAACAAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAAcc                                                                           >  1:1986691/1‑100 (MQ=255)
                  gTCAGAAAATCGCCACAGCGATTGCTGAAGGCGTGATCAGTTATTTACACTGGTTCGACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAAcc                                                                           >  2:1894615/1‑100 (MQ=255)
                     aGAAAATCGCCACAGCGATTGCTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGa                                                                        >  2:1660454/1‑100 (MQ=255)
                      gAAAATCGCCACAGCGATTGCTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGagatcggaa                                                                       >  1:1188977/1‑92 (MQ=255)
                       aaaaTCGCCACAGCGATTGCTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACg                                                                      <  2:104293/100‑1 (MQ=255)
                       aaaaTCGCCACAGCGATTGCTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACg                                                                      <  2:3075794/100‑1 (MQ=255)
                       aaaaTCGCCACAGCGATTGCTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACg                                                                      >  2:2485526/1‑100 (MQ=255)
                       aaaaTCGCCACAGCGATTGCTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACg                                                                      <  1:1599884/100‑1 (MQ=255)
                          aTCGCCACAGCGATTGCTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGcac                                                                   <  1:694415/100‑1 (MQ=255)
                           tCGCCACAGCGATTGCTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGcaca                                                                  <  2:2707371/100‑1 (MQ=255)
                             gCCACAGCGATTGCTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGCACAcc                                                                >  1:1526643/1‑100 (MQ=255)
                             gCCACAGCGATTGCTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGCACAcc                                                                >  1:3073412/1‑100 (MQ=255)
                              ccACAGCGATTGCTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCa                                                               <  1:3099445/100‑1 (MQ=255)
                               cacaGCGATTGCTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGGAACCGCGATAAGTTATGAAACCCGACGCACACCAg                                                              >  2:2966720/1‑100 (MQ=255)
                               cacaGCGATTGCTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCAg                                                              >  1:720150/1‑100 (MQ=255)
                                   gCGATTGCTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTa                                                          >  2:2780281/1‑100 (MQ=255)
                                    cGATTGCTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTaa                                                         >  2:1732125/1‑100 (MQ=255)
                                     gATTGCTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTaaa                                                        >  2:1218671/1‑100 (MQ=255)
                                     gATTGCTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTaaa                                                        >  1:1762257/1‑100 (MQ=255)
                                     gATTGCTGAAGGCGTGATCAGTTATTTCCACTGGTCCGACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTaa                                                         <  1:1732128/99‑1 (MQ=255)
                                      aTTGCTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAAc                                                       <  1:1670178/100‑1 (MQ=255)
                                       ttGCTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACa                                                      <  1:1113437/100‑1 (MQ=255)
                                         gCTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGt                                                    >  2:2705161/1‑100 (MQ=255)
                                         gCTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGt                                                    >  1:932102/1‑100 (MQ=255)
                                         gCTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGt                                                    >  1:1633715/1‑100 (MQ=255)
                                         gCTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGt                                                    >  1:1550715/1‑100 (MQ=255)
                                         gCTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGt                                                    >  1:1311911/1‑100 (MQ=255)
                                           tGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGttt                                                  <  2:764895/100‑1 (MQ=255)
                                           tGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGttt                                                  >  1:1472170/1‑100 (MQ=255)
                                            gAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTc                                                 >  2:1930036/1‑100 (MQ=255)
                                            gAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTc                                                 <  1:661606/100‑1 (MQ=255)
                                              aGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTg                                               <  1:2764934/100‑1 (MQ=255)
                                               ggCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGc                                              <  1:2904886/100‑1 (MQ=255)
                                               gatctGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGc                                              <  1:825353/96‑1 (MQ=255)
                                                gCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCt                                             <  1:1660456/100‑1 (MQ=255)
                                                gCGTGATCAGTTATGTCCACTGGTTCGACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCt                                             <  1:1842322/100‑1 (MQ=255)
                                                 cGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTc                                            <  2:932101/100‑1 (MQ=255)
                                                  gTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCa                                           <  1:78257/100‑1 (MQ=255)
                                                  gTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCa                                           >  1:2694928/1‑100 (MQ=255)
                                                    gATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGcagatc                                         >  2:825353/1‑95 (MQ=255)
                                                    gATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAAc                                         >  2:2986796/1‑100 (MQ=255)
                                                    gATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAAc                                         >  2:12950/1‑100 (MQ=255)
                                                    gATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAAc                                         >  1:1672457/1‑100 (MQ=255)
                                                     aTCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAAcc                                        >  1:2876098/1‑100 (MQ=255)
                                                         gTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCa                                    >  1:1022852/1‑100 (MQ=255)
                                                         gTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCa                                    >  2:290230/1‑100 (MQ=255)
                                                         gTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCa                                    <  1:639054/100‑1 (MQ=255)
                                                         gTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCa                                    >  2:2518782/1‑100 (MQ=255)
                                                         gTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGCAAACCAGGTTAAACAGTTTCTGCTCAACCtt                                      >  1:299570/1‑98 (MQ=255)
                                                          ttATTTCCACTGGTTCGACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAg                                   <  2:942886/100‑1 (MQ=255)
                                                          ttATTTCCACTGGTTCGACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAg                                   <  1:2012417/100‑1 (MQ=255)
                                                           tATTTCCACTGGTTCGACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAgg                                  <  1:79045/100‑1 (MQ=255)
                                                             tttCCACTGGTTCGACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGAt                                >  2:2835453/1‑100 (MQ=255)
                                                               tCCACTGGTTCGACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATAc                              >  2:3329572/1‑100 (MQ=255)
                                                               tCCACTGGTTCGACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATAc                              <  2:2129742/100‑1 (MQ=255)
                                                                ccACTGGTTCGACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACg                             >  2:2211225/1‑100 (MQ=255)
                                                                 cACTGGTTCGACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACGa                            >  2:3392230/1‑100 (MQ=255)
                                                                   cTGGTTCGACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACGAtt                          >  2:1100823/1‑100 (MQ=255)
                                                                     ggTTCGACAACCAGAAGGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACGATTTg                        >  1:310500/1‑100 (MQ=255)
                                                                     gatctgACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACGATTTg                        <  2:697830/95‑1 (MQ=255)
                                                                      gTTCGATAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACGATTTGt                       >  2:532553/1‑100 (MQ=255)
                                                                      gTTCGACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACGATTTGt                       >  1:2536047/1‑100 (MQ=255)
                                                                       ttCGACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACGATTTGTc                      <  1:2518790/100‑1 (MQ=255)
                                                                       ttCGACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACGATTTGTc                      <  2:3125773/100‑1 (MQ=255)
                                                                       ttCGACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACGATTTGTc                      <  1:3168648/100‑1 (MQ=255)
                                                                       ttCGACAAACAGAAAGCATATTCGGAAAAGAGATAAGTTATGACAACAGAAGAACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGTTACGTTTTGTc                      <  2:1986687/100‑1 (MQ=255)
                                                                        tCGACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACGATTTGTca                     >  2:1246858/1‑100 (MQ=255)
                                                                          gACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCACCCTTCAGGATACGGTTTgagatc                   >  1:697830/1‑95 (MQ=255)
                                                                          gACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACGATTTGTcagc                   >  2:1440293/1‑100 (MQ=255)
                                                                          gACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACGATTTGTcagc                   <  2:1848334/100‑1 (MQ=255)
                                                                           aCAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACGATTTGTcagca                  <  1:1930040/100‑1 (MQ=255)
                                                                            ctaCCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACGATTTGTcagcag                 <  2:1279949/98‑1 (MQ=255)
                                                                            cgatctgAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACGATTTGTcagcag                 <  1:1126628/94‑1 (MQ=255)
                                                                            cAACCAGAACGCCAATACTCAAAGTCGATAAGTGATGAAGCCCTACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACGATTTGTcagca                  <  2:3074088/99‑1 (MQ=255)
                                                                             aaCCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACGATTTGTCagcagc                <  2:422786/100‑1 (MQ=255)
                                                                              aCCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACGATTTGTcagcagag               >  1:1279952/1‑98 (MQ=255)
                                                                              aCCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACGATTTGTCAGCAGCt               >  1:1703381/1‑100 (MQ=255)
                                                                                cAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACGATTTGTCagcagccga             >  2:2311357/1‑97 (MQ=255)
                                                                                  gAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACGATTTGTcagcagagatcg           >  2:1126626/1‑94 (MQ=255)
                                                                                     aGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACGATTTGTCAGCAGCTGAccgcc        >  1:2564613/1‑100 (MQ=255)
                                                                                      gCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACGATTTGTCAGCAGCTGAccgccg       <  2:1022851/100‑1 (MQ=255)
                                                                                         caTTCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACGATTTGTCAGCAGCTGACCGCCGTCg    >  1:2714508/1‑100 (MQ=255)
                                                                                         caTTCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACGATTTGTCAGCAGCTGACCGCCGTCg    >  1:295944/1‑100 (MQ=255)
                                                                                         caTTCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACGATTTGTCAGCAGCTGACCGCCGTCg    >  1:229203/1‑100 (MQ=255)
                                                                                           ttCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACGATTTGTCAGCAGCTGACCGCCGTCGAt  <  2:1672455/100‑1 (MQ=255)
                                                                                           ttCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACGATTTGTCAGCAGCTGACCGCCGTCGAt  <  2:1472167/100‑1 (MQ=255)
                                                                                               |                                                                                               
TAGGCACGGCGGCGTTTCGTCAGAAAATCGCCACAGCGATTGCTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGAAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACGATTTGTCAGCAGCTGACCGCCGTCGAT  >  NC_000913/2553115‑2553305

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