Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,553,210 A→G 100% 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,2100AG93.2% 284.9 / 4.5 95K287E (AAA→GAA) amiAN‑acetylmuramoyl‑l‑alanine amidase I
Reads supporting (aligned to +/- strand):  ref base A (6/1);  new base G (52/34);  total (59/36)
Fisher's exact test for biased strand distribution p-value = 2.49e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TAGGCACGGCGGCGTTTCGTCAGAAAATCGCCACAGCGATTGCTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGAAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACGATTTGTCAGCAGCTGACCGCCGTCGAT  >  NC_000913/2553115‑2553305
                                                                                               |                                                                                               
tAGGCACGGCGGCGTTTCGTCAGAAAATCGCCACAGCGATTGCTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGgaaaa                                                                                             >  2:276629/1‑100 (MQ=255)
 aGGCACGGCGGCGTTTCGTCAGAAAATCGCCACAGCGATTGCTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGGAAAAg                                                                                            >  2:1846053/1‑100 (MQ=255)
 aGGCACGGCGGCGTTTCGTCAGAAAATCGCCACAGCGATTGCTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGGAAAAg                                                                                            >  2:220808/1‑100 (MQ=255)
 aGGCACGGCGGCGTTTCGTCAGAAAATCGCCACAGCGATTGCTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGAAAAAg                                                                                            >  1:947888/1‑100 (MQ=255)
  ggCACGGCGGCGTTTCGTCAGAAAATCGCCACAGCGATTGCTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGagatcgg                                                                                           >  1:2062667/1‑94 (MQ=255)
    cACGGCGGCGTTTCGTCAGAAAATCGCCACAGCGATTGCTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGGAAAAGCGa                                                                                         >  1:59542/1‑100 (MQ=255)
    cACGGCGGCGTTTCGTCAGAAAATCGCCACAGCGATTGCTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGGAAAAGCGa                                                                                         >  2:707183/1‑100 (MQ=255)
    cACGGCGGCGTTTCGTCAGAAAATCGCCACAGCGATTGCTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGGAAAAGCGa                                                                                         >  1:933314/1‑100 (MQ=255)
    cACGGCGGCGTTTCGTCAGAAAATCGCCACAGCGATTGCTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGGAAAAGCGa                                                                                         >  2:140410/1‑100 (MQ=255)
    cACGGCGGCGTTTCGTCAGAAAATCGCCACAGCGATTGCTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGGAAAAGCGa                                                                                         <  2:1787317/100‑1 (MQ=255)
      cggcggCGTTTCGTCAGAAAATCGCCACAGCGATTGCTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGGAAAAGCGATa                                                                                       >  1:1091731/1‑100 (MQ=255)
       ggcggcGTTTCGTCAGAAAATCGCCACAGCGATTGCTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGGAAAAGCGATaa                                                                                      >  2:1284450/1‑100 (MQ=255)
        gcggcgTTTCGTCAGAAAATCGCCACAGCGATTGCTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGGAAAAGCGATAAg                                                                                     <  1:192892/100‑1 (MQ=255)
        gcggcgTTTCGTCAGAAAATCGCCACAGCGATTGCTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGGAAAAGCGATAAg                                                                                     >  2:638677/1‑100 (MQ=255)
            cgTTTCGTCAGAAAATCGCCACAGCGATTGCTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGAAAAAGCGATAAGTTAt                                                                                 >  2:1101085/1‑100 (MQ=255)
             gTTTCGTCAGAAAATCGTCACAGCGATTGCTGAAGGCGTGATCAGTTATTTCCAATGGTTCGACAACCAGAAAGCACATTCGGAAAAGCGATACGTTATg                                                                                <  1:1353526/100‑1 (MQ=255)
             gTTTCGTCAGAAAATCGCCACAGCGATTGCTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATg                                                                                >  1:2351290/1‑100 (MQ=255)
             gTTTCGTCAGAAAATCGCCACAGCGATTGCTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGAAAAAGCGATAAGTTATg                                                                                >  1:1466763/1‑100 (MQ=255)
              tttCGTCAGAAAATCGCCACAGCGATTGCTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGa                                                                               >  1:2258336/1‑100 (MQ=255)
                 cGTCAGAAAATCGCCACAGCGATTGCTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAAc                                                                            >  2:482030/1‑100 (MQ=255)
                 cGTCAGAAAATCGCCACAGCGATTGCTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAAc                                                                            >  2:717340/1‑100 (MQ=255)
                  gTCAGAAAATCGCCACAGCGATTGCTGAAGGCTTGTGCTGTTATTTCCGCTGGTTCGACAACCAGAAAGCACATTCGGAAAAGCGATTAGTTATGAAAc                                                                            <  1:202718/99‑1 (MQ=255)
                    cAGAAAATCGCCACAGCGATTGCTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAAccc                                                                          <  1:2144848/99‑1 (MQ=255)
                     aGAAAATCGCCACAGCGATTGCTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGa                                                                        >  2:2478196/1‑100 (MQ=255)
                     aGAAAATCGCCACAGCGATTGCTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGa                                                                        <  2:997510/100‑1 (MQ=255)
                     aGAAAATCGCCACAGCGATTGCAGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGa                                                                        >  1:1499448/1‑100 (MQ=255)
                      gAAAATCGCCACAGCGATTGCTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGAc                                                                       >  2:1587202/1‑100 (MQ=255)
                       aataTCGCCACAGCGATTGCTGAAGGAGTGATCAGTTATTTCCACTTGTTCGACAACCCGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGAc                                                                       <  1:1805555/96‑1 (MQ=255)
                        aaaTCGCCACAGCGATTGCTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGc                                                                     <  1:468750/100‑1 (MQ=255)
                           tCGCCACAGCGATTGCTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGcaca                                                                  <  1:276626/100‑1 (MQ=255)
                             gCCACAGCGATTGCTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACCACCAGAAAGCCCATTCGGAAAAGCGATAAGTTATGAAACACGACGcac                                                                   >  1:247552/1‑97 (MQ=255)
                             gCCACAGCGATTGCTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGCACAcc                                                                <  1:2102875/100‑1 (MQ=255)
                             gCCACAGCGATTGCTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGCACAcc                                                                >  2:2146505/1‑100 (MQ=255)
                             gCCACAGCGATTGCTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGCACAcc                                                                >  1:2593693/1‑100 (MQ=255)
                             gCCACAGCGATTGCTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGCACAcc                                                                <  2:530511/100‑1 (MQ=255)
                             gCCACAGCGATTGCTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGCACAcc                                                                >  1:1308983/1‑100 (MQ=255)
                             gCCACAGCGACTGCTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGAAAAAGCGATACGATATGAAACCCGACACACAcc                                                                >  1:573005/1‑100 (MQ=255)
                               cacaGCGATTGCTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGTAAAAGCGATAAGTTATGAAACCAGACGCACACCAg                                                              >  1:592689/1‑100 (MQ=255)
                               cacaGCGATTGCTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCan                                                              <  2:1921346/100‑2 (MQ=255)
                                   gCGATTGCTGCAGGCGAGATCAGTTATTTCCACTGCTTCGACAACCAGCAAGCACACCCGGAATAGCGATACGTTATCAAACCCGACGCacacacgttta                                                          >  1:2456775/1‑93 (MQ=255)
                                   gCGATTGCTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTa                                                          >  2:1932614/1‑100 (MQ=255)
                                       ttGCTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACa                                                      <  1:2084967/100‑1 (MQ=255)
                                         gCTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGt                                                    >  2:1434028/1‑100 (MQ=255)
                                          cTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGGAAAAGCGATCAGTTATGAAACCCGACGCACACCAGGTTAAACAGtt                                                   >  1:179504/1‑100 (MQ=255)
                                          cTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGtt                                                   >  2:1895661/1‑100 (MQ=255)
                                          cTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGAAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAg                                                     >  2:2411607/1‑98 (MQ=255)
                                           tGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGttt                                                  >  1:1447157/1‑100 (MQ=255)
                                            gAAGGCGTGATCAGTTATTTCCACTGGTTCGCCAACCAGAAAGCACATTCGGAAAAGCGATACGTTATGAGACCGGACGCACACCAGGGTAAACAGtt                                                   >  1:100517/1‑98 (MQ=255)
                                            gAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTc                                                 <  1:712880/100‑1 (MQ=255)
                                              aGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTg                                               <  2:1315137/100‑1 (MQ=255)
                                               ggCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGc                                              <  1:1463423/100‑1 (MQ=255)
                                                gCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCt                                             <  2:817571/100‑1 (MQ=255)
                                                gCGTGAGCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCt                                             <  2:2156727/100‑1 (MQ=255)
                                                  gTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCa                                           <  2:1028214/100‑1 (MQ=255)
                                                  gTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCa                                           >  2:2492388/1‑100 (MQ=255)
                                                  gTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGAAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCa                                           >  2:1638877/1‑100 (MQ=255)
                                                   tGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCaa                                          <  1:330167/100‑1 (MQ=255)
                                                    gATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAAc                                         <  2:1091735/100‑1 (MQ=255)
                                                    gATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAAc                                         <  1:758051/100‑1 (MQ=255)
                                                     aTCTGTTATTTCTACTGGTTCGACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAAcc                                        <  2:2160897/100‑1 (MQ=255)
                                                        aGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTc                                     >  2:2004864/1‑100 (MQ=255)
                                                         gTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCa                                    <  2:1559744/100‑1 (MQ=255)
                                                           tATTTCCACTGGTTCGACAACCAGAAAGCACATTCGGAGAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAgg                                  >  1:1600029/1‑100 (MQ=255)
                                                              ttCCACTGGTTCGACAACCAGAAAGCACATTCGGCAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATa                               >  2:1875177/1‑100 (MQ=255)
                                                              ttCCACTGGTTCGACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATa                               >  1:2170112/1‑100 (MQ=255)
                                                               tCCCCTGGTTCGACAAACAGAAAGCACATTAGGAAAAGCGATAAGTTATGACAACCGACGCACACCAGGCTAAACAGTTTCTCCTCAACCTTCAGAATAc                              >  1:755283/1‑100 (MQ=255)
                                                               tCCACTGGTTCGACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATAc                              <  1:2419679/100‑1 (MQ=255)
                                                               tCCACTGGTTCGACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATAc                              >  2:409676/1‑100 (MQ=255)
                                                               tCCACTGGTTCGACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATAc                              >  2:1556041/1‑100 (MQ=255)
                                                                ccACTGGTTCGACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACg                             >  2:2359056/1‑100 (MQ=255)
                                                                 cACTGGTTCGACAACCAGAAAGCACATTCGGAAACGCGTTAAGTTATGAAACCCGACGCACACCAGTCTAAACCGTTTCTGCTCAACCTCAAGGATACGa                            >  1:658797/1‑100 (MQ=255)
                                                                 cACTGGTTCGACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACGa                            >  2:1769543/1‑100 (MQ=255)
                                                                  aCTGGTTCGACAACCAGAAAGCACATTCGCAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGCTAAACAGTTTCTGCTCAACCTTCAGGATACGAt                           <  1:1556032/100‑1 (MQ=255)
                                                                     ggTTCGCCAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATCAAACACGACGCACACCAAGCTAAACAGCTTCTGCCCAACCTTCACGATCCg                             >  1:2014589/1‑95 (MQ=255)
                                                                     ggTTCGACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACGATTTg                        <  2:2279462/100‑1 (MQ=255)
                                                                     ggTTCGACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACGATTTg                        >  2:2443502/1‑100 (MQ=255)
                                                                     ggTTCGACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACGATTTg                        >  2:1829538/1‑100 (MQ=255)
                                                                       ttCGACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACGATTTGTc                      <  2:1912/100‑1 (MQ=255)
                                                                        tCGACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACGATTTGTca                     <  1:694871/100‑1 (MQ=255)
                                                                         cgatctccAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACGATTTGTcag                    <  2:525807/94‑1 (MQ=255)
                                                                         cGACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACGATTTGTcag                    <  2:1416071/100‑1 (MQ=255)
                                                                         cGACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACGATTTGTcag                    >  2:1076841/1‑100 (MQ=255)
                                                                          gACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACGATTTGTcagc                   >  1:1835432/1‑100 (MQ=255)
                                                                          gACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACGATTTGTcagc                   <  2:2287931/100‑1 (MQ=255)
                                                                            cAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACGATTTGTcagcag                 >  1:117544/1‑100 (MQ=255)
                                                                             aaCCAGAAAGCACATTCGAAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATATGATTTGTCagcagc                <  1:1638868/100‑1 (MQ=255)
                                                                              aCCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACGATTTGTCAGCAGCt               >  2:2585605/1‑100 (MQ=255)
                                                                              aCCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACGATTTGTCAGCAGCt               >  2:486933/1‑100 (MQ=255)
                                                                               ccAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACGATTTGTcagagatcg              >  1:525805/1‑94 (MQ=255)
                                                                               ccAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACGATTTGTCAGCAGCTg              >  2:392761/1‑100 (MQ=255)
                                                                                  gAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACGATTTGTCAGCAGCTGAcc           >  1:1099752/1‑100 (MQ=255)
                                                                                   aaaGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACGATTTGTCAGCAGCTGAccg          <  2:1766042/100‑1 (MQ=255)
                                                                                   aaaGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACGATTTGTCAGCAGCTGAccg          <  1:1788955/100‑1 (MQ=255)
                                                                                        acaTTCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACGATTTGTCAGCAGCTGACCGCCGTc     >  2:413676/1‑100 (MQ=255)
                                                                                         caTTCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACGATTTGTCAGCAGCTGACCGCCGTCg    >  1:2500380/1‑100 (MQ=255)
                                                                                         caTTCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACGATTTGTCAGCAGCTGACCGCCGTCg    >  2:287488/1‑100 (MQ=255)
                                                                                           ttCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACGATTTGTCAGCAGCTGACCGCCGTCGAt  <  1:2492369/100‑1 (MQ=255)
                                                                                               |                                                                                               
TAGGCACGGCGGCGTTTCGTCAGAAAATCGCCACAGCGATTGCTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGAAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACGATTTGTCAGCAGCTGACCGCCGTCGAT  >  NC_000913/2553115‑2553305

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