Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 101,342 C→T T193I (ACC→ATC)  murC → UDP‑N‑acetylmuramate:L‑alanine ligase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913101,3420CT99.0% 316.6 / ‑5.0 95T193I (ACC→ATC) murCUDP‑N‑acetylmuramate:L‑alanine ligase
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base T (48/46);  total (48/47)
Fisher's exact test for biased strand distribution p-value = 4.95e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.68e-01

TTCATGCGCGTTTGGGGCATGGTCGGTACCTGATTGCCGAAGCAGATGAGAGTGATGCATCGTTCCTGCATCTGCAACCGATGGTGGCGATTGTCACCAATATCGAAGCCGACCACATGGATACCTACCAGGGCGACTTTGAGAATTTAAAACAGACTTTTATTAATTTTCTGCACAACCTGCCGTTTTACG  >  NC_000913/101246‑101437
                                                                                                |                                                                                               
ttCATGCGCGTGTGGGGCATGGTCGGTACCTGATTGCCGAAGCACATGTGAGTGATGCATCGTTCCTGCATCTGCAACCGATGGTGGCGATTGTCATCaa                                                                                              >  1:2319480/1‑100 (MQ=255)
 tCATGCGCGTTTGGGGCATGGTCGGTACCTGATTGCCGAAGCAGATGAGAGTGATGCATCGTTCCTGCATCTGCAACCGATGGTGGCGATTGTCATCAat                                                                                             >  1:1806694/1‑100 (MQ=255)
  cATTCGCGTTTGGAGCATGGTCGGTACCTGATTGCCGAAGCAGATGAGAGTGATGCATCGTTCCTGCATCTGCAACCGATGCTGGCGATTGTCATCAata                                                                                            <  2:1654617/100‑1 (MQ=255)
  cATGCGCGTTTGGGGCATGGTCGGTACCTGATTGCCGAAGCAGATGCGAGTGATGCCTCGTTCCTGCATCTGCAACCGATGGTGGCGATTGTCATCAata                                                                                            >  1:1580646/1‑100 (MQ=255)
  cATGCGCGTTTGGGGCATGGTCGGTACCTGATTGCCGAAGCAGATGAGAGTGATGCATCGTTCCTGCATCTGCAACCGATGGTGGCGATTGTCATCAata                                                                                            >  1:1867857/1‑100 (MQ=255)
    tGCGCGTTTGGGGCATGGTCGGTACCTGATTGCCGAAGCAGATGAGAGTTATGCATCGTTCCTGCATCTGCAACCGATGGTGGCGATTGTCATCAATATc                                                                                          <  1:1168233/100‑1 (MQ=255)
    tGCGCGTTTGGGGCATGGTCGGTACCTGATTGCCGAAGCAGATGAGAGTGATGCATCGTTCCTGCATCTGCAACCGATGGTGGCGATTGTCATCAATATc                                                                                          <  2:1843116/100‑1 (MQ=255)
    tGCGCGTTTGGGGCATGGTCGGTACCTGATTGCCGAAGCAGATGAGAGTGATGCATCGTTCCTGCATCTGCAACCGATGGTGGCGATTGTCATCAATATc                                                                                          >  2:1028868/1‑100 (MQ=255)
     gcgcGTTTGGGGCATGGTCGGTACCTGATTGCCGAAGCAGATGAGAGTGATGCATCGTTCCTGCATCTGCAACCGATGGTGGCGATTGTCATCAATATCg                                                                                         >  1:1915672/1‑100 (MQ=255)
      cgcgTTTGGGGCATGGTCGGTACCTGATTGCCGAAGCAGATGAGAGTGATGCATCGTTCCTGCATCTGCAACCGATGGTGGCGATTGTCATCAATATCGa                                                                                        >  2:1605309/1‑100 (MQ=255)
        cgTTTGGGGCATGGTCGGTACCTGATTGCCGAAGCAGATGAGAGTGATGCATCGTTCCTGCATCTGCAACCGATGGTGGCGATTGTCATCAATATCGAAg                                                                                      <  2:1553389/100‑1 (MQ=255)
         gTTTGGGGCATGGTCGGTACCTGATTGCCGAAGCAGATGAGAGTGATGCATCGTTCCTGCATCTGCAACCGATGGTGGCGATTGTCATCAATATCGAAGc                                                                                     >  1:1072216/1‑100 (MQ=255)
          tttGGGGCATGGTCGGTACCTGATTGCCGAAGCAGATGAGAGTGATGCATCGTTCCTGCATCTGCAACCGATGGTGGCGATTGTCATCAATATCGAAGcc                                                                                    <  2:1681107/100‑1 (MQ=255)
           ttGGGGCATGGTCGGTACCTGATTGCCGAAGCAGATGAGAGTGATGCATCGTTCCTGCATCTGCAACCGATGGTGGCGATTGTCATCAATATCGAAGCCg                                                                                   <  2:1721644/100‑1 (MQ=255)
            tGGGGCATGGTCGGTACCTGATTGCCGAAGCAGATGAGAGTGATGCATCGTTCCTGCATCTGTAACCGATGGTGGCGATTGTCATCAATATCGAAGCCGa                                                                                  >  1:1787134/1‑100 (MQ=255)
            tGGGGCATGGTCGGTACCTGATTGCCGAAGCAGATGAGAGTGATGCATCGTTCCTGCATCTGCAACCGATGGTGGCGATTGTCATCAATATCGAAGCCGa                                                                                  <  1:2221920/100‑1 (MQ=255)
            tGGGGCATGGTCGGTACCTGATTGCCGAAGCAGATGAGAGTGATGCATCGTTCCTGCATCTGCAACCGATGGTGGCGATTGTCATCAATATCGAAGCCGa                                                                                  >  2:970069/1‑100 (MQ=255)
              gggCATGGTCGGTACCTGATTGCCGAAGCAGATGAGAGTGATGCATCGTTCCTGCATCTGCAACCGATGGTGGCGATTGTCATCAATATCGAAGCCGAcc                                                                                >  1:1927006/1‑100 (MQ=255)
              gggCATGGTCGGTACCTGATTGCCGAAGCAGATGAGAGTGATGCATCGATCCTGCATCTGCAACCGATGGTGGCGATTGTCATCAATATCGAAGCCGaac                                                                                >  2:1192912/1‑98 (MQ=255)
               ggCATGGTCGGTACCTGATTGCCGAAGCAGATGAGAGTGATGCATCGTTCCTGCATCTGCAACCGATGGTGGCGATTGTCATCAATATCGAAGCCGACca                                                                               <  1:2611569/100‑1 (MQ=255)
                ggaTGGTCGGTACCTGATTGCCGAAGCAGATGAGAGTGATGCATCGTTCCTGCATCTGCAACCGATGGTGGCGATTGTCATCAATATCGAAGCCGACcac                                                                              <  2:1946732/98‑1 (MQ=255)
                gCATGGTCGGTACCTGATTGCCGAAGCAGATGAGAGTGATGCATCGTTCCTGCATCTGCAACCGATGGTGGCGATTGTCATCAATATCGAAGCCGACcac                                                                              >  1:1953376/1‑100 (MQ=255)
                  aTGGTCGGTACCTGATTGCCGAAGCAGATGAGAGTGATGCATCGTTCCTGCATCTGCAACCGATGGTGGCGATTGTCATCAATATCGAAGCCGACCACAt                                                                            <  2:481302/100‑1 (MQ=255)
                   tGGTCGGTACCTGATTGCCGAAGCAGATGAGAGTGATGCATCGTTCCTGCATCTGCAACCGATGGTGGCGATTGTCATCAATATCGAAGCCGACCACATg                                                                           <  2:168756/100‑1 (MQ=255)
                     gTCGGTACCTGATTGCCGAAGCAGATGAGAGTGATGCATCGTTCCTGCATCTGCAACCGATGGTGGCGATTGTCATCAATATCGAAGCCGACCACATGGa                                                                         <  1:1629278/100‑1 (MQ=255)
                       cGGTACCTGATTGCCGAAGCAGATGAGAGTGATGCATCGTTCCTGCATCTGCAACCGATGGTGGCGATTGTCATCAATATCGAAGCCGACCACATGGAta                                                                       >  1:2048450/1‑100 (MQ=255)
                         ttACCGGATTGCCGAAGCAGATGAGATGGATGCATCGTTCCTGCATCTGCAACCGATGGTGGCGATTGTCATCAATATCGAAGCCGACCACATGGAtacc                                                                     <  1:930042/99‑1 (MQ=255)
                         gTACCTGATTGCCGAAGCAGATGAGAGTGATGCATCGTTCCTGCATCTGCAACCGATGGTGGCGATTGTCATCAATATCGAAGCCGACCACATGGAtacc                                                                     >  2:1278033/1‑100 (MQ=255)
                          tACCTGATTGCCGAAGCAGATGAGAGTGATGCATCGTTCCTGCATCTGCAACCGATGGTGGCGATTGTCATCAATATCGAAGCCGACCACATGGAtacct                                                                    >  1:2168851/1‑100 (MQ=255)
                          tACCTGATTGCCGAAGCAGATGAGAGTGATGCATCGTTCCTGCATCTGCAACCGATGGTGGCGATTGTCATCAATATCGAAGCCGACCACATGGAtacct                                                                    <  1:2597307/100‑1 (MQ=255)
                             cTGATTGCCGAAGCAGATGAGAGTGATGCATCGTTCCTGCATCTGCAACCGATGGTGGCGATTGTCATCAATATCGAAGCCGACCACATGGAtacctacc                                                                 <  2:706747/100‑1 (MQ=255)
                             cTGATTGCCGAAGCAGATGAGAGTGATGCATCGTTCCTGCATCTGCAACCGATGGTGGCGATTGTCATCAATATCGAAGCCGACCACATGGAtacctacc                                                                 <  2:46178/100‑1 (MQ=255)
                             cTGATTGCCGAAGCAGATGAGAGTGATGCATCGTTCCTGCATCTGCAACCGATGGTGGCGATTGTCATCAATATCGAAGCCGACCACATGGAtacctacc                                                                 <  1:1199331/100‑1 (MQ=255)
                              tGATTGCCGAAGCAGATGAGAGTGATGCATCGTTCCTGCATCTGCAACCGATGGTGGCGATTGTCATCAATATCGAAGCCGACCACATGGATACCTACCa                                                                >  1:2017104/1‑100 (MQ=255)
                               gATTGCCGAAGCAGATGAGAGTGATGCATCGTTCCTGCATCTGCAACCGATGGTGGCGATTGTCATCAATATCGAAGCCGACCACATGGATACCTACCAg                                                               >  1:293204/1‑100 (MQ=255)
                                aTTGCCGAAGCAGATGAGAGTGATGCATCGTTCCTGCATCTGCAACCGATGGTGGCGATTGTCATCAATATCGAAGCCGACCACATGGATACCTACCAgg                                                              >  1:1269503/1‑100 (MQ=255)
                                aTTGCCGAAGCAGATGAGAGTGATGCATCGTTCCTGCATCTGCAACCGATGGTGGCGATTGTCATCAATATCGAAGCCGACCACATGGATACCTACCAgg                                                              <  2:2584426/100‑1 (MQ=255)
                                aTTGCCGAAGCAGATGAGAGTGATGCATCGTTCCTGCATCTGCAACCGATGGTGGCGATTGTCATCAATATCGAAGCCGACCACATGGATACCTACCAgg                                                              <  2:2134786/100‑1 (MQ=255)
                                   gCCGAAGTAGATGAGAGTGATGCATGGTTCATGCATCTGCAACCGATGGTGGCGATTGTCATCAATATCGAAGCCGACCACATGGATACCTACCAGGGCg                                                           <  2:842232/100‑1 (MQ=255)
                                   gCCGAAGCAGATGAGAGTGATGCATCGTTCCTGCATCTGCAACCGATGGTGGCGATTTTCATCAATATCGAAGCCGACCACATGGATCCCTACCAGGGCg                                                           >  2:1265147/1‑100 (MQ=255)
                                   gCCGAAGCAGATGAGAGTGATGCATCGTTCCTGCATCTGCAACCGATGGTGGCGATTGTCATCAATATCGAAGCCGACCACATGGATACCTACCAGGGCg                                                           <  1:1688967/100‑1 (MQ=255)
                                   gCCGAAGCAGATGAGAGTGATGCATCGTTCCTGCATCTGCAACCGATGGTGGCGATTGTCATCAATATCGAAGCCGACCACATGGAAACCTACCAGGGCg                                                           >  1:2568504/1‑100 (MQ=255)
                                    ccGAAGCAGATGAGAGTGATGCATCGTTCCTGCATCTGCAACCGATGGTGGCGATTGTCATCAATATCGAAGCCGACCACATGGATACCTACCAGGGCGa                                                          >  1:598206/1‑100 (MQ=255)
                                     cGAAGCAGATGAGAGTGATGCATCGTTCCTGCATCTGCAACCGATGGTGGCGATTGTCATCAAAATCGAAGCCGGCCAGATGGATACCTACCAGGGCGAc                                                         >  2:2454304/1‑100 (MQ=255)
                                      gAAGCAGATGAGAGTGATGCATCGTTCCTGCATCTGCAACCGATGGTGGCGATTGTCATCAATATCGAAGCCGACCACATGGATACCTACCAGGGCGACt                                                        >  2:1324831/1‑100 (MQ=255)
                                           aGATGAGAGTGATGCATCGTTCCTGCATCTGCAACCGATGGTGGCGATTGTCATCAATATCGAAGCCGACCACATGGATACCTACCAGGGCGACTTTgag                                                   <  2:307670/100‑1 (MQ=255)
                                           aGATGAGAGTGATGCATCGTTCCTGCATCTGCAACCGATGGTGGCGATTGTCATCAATATCGAAGCCGACCACATGGATACCTACCAGGGCGACTTTgag                                                   <  1:1086330/100‑1 (MQ=255)
                                           aGATGAGAGTGATGCATCGTTCCTGCATCTGCAACCGATGGTGGCGATTGTCATCAATATCGAAGCCGACCACATGGATACCTACCAGGGCGACTTTgag                                                   <  2:2032189/100‑1 (MQ=255)
                                           aGATGAGAGTGATGCATCGTTCCTGCATCTGCAACCGATGGTGGCGATTGTCATCAATATCGAAGCCGACCACATGGATACCTACCAGGGCGACTTTgag                                                   <  1:1374992/100‑1 (MQ=255)
                                            gATGAGAGTGATGCATCGTTCCTGCATCTGCAACCGATGGTGGCGATTGTCATCAATATCGAAGCCGACCACATGGATACCTACCAGGGCGACTTTgaga                                                  >  1:23969/1‑100 (MQ=255)
                                                agagTGATGCATCGTTCCTGCATCTGCAACCGATGGTGGCGATTGTCATCAATATCGAAGCCGACCACATGGATACCTACCAGGGCGACTTTGAGAAttt                                              >  2:699966/1‑100 (MQ=255)
                                                 gagTGATGCATCGTTCCTGCATCTGCAACCGATGGTGGCGATTGTCATCAATATCGAAGCCGACCACATGGATACCTACCAGGGCGACTTTGAGAATTTa                                             >  1:1294053/1‑100 (MQ=255)
                                                 gagTGATGCAGCGTGCCTGCATCTGCAACCGATGGTGGCGATTGTCATCAATATCGAAGCCGACACATGGATACCTACCAGGGCGACTTTGAGAATTTaa                                            <  1:2461571/100‑1 (MQ=255)
                                                   gTGATGCATCGTTCCTGCATCTGCAACCGATGGTGGCGATTGTCATCAATATCGAAGCCGACCACATGGATACCTACCAGGGCGACTTTGAGAATTTaaa                                           <  2:1710960/100‑1 (MQ=255)
                                                   gTGATGCATCGTTCCTGCATCTGCAACCGATGGTGGCGATTGTCATCAATATCGAAGCCGACCACATGGATACCTACCAGGGCGACTTTGAGAATTTaaa                                           <  2:383418/100‑1 (MQ=255)
                                                   gTGATGCATCGTTCCTGCATCTGCAACCGATGGTGGCGATTGTCATCAATATCGAAGCCGACCACATGGATACCTACCAGGGCGACTTTGAGAATTTaaa                                           <  2:2287750/100‑1 (MQ=255)
                                                    tGATGCATGGTTCCGGCATCTGCAACCGATGGTGGCGATTGTCATCAATATCGAAGCCGACCACATGGATACCTACCAGGGCGACTTTGAGAATTTaaaa                                          <  1:698302/100‑1 (MQ=255)
                                                    tGATGCATCGTTCCTGCATCTGCAACCGATGGTGGCGATTGTCATCAATATCGAAGCCGACCACATGGATACCTACCAGGGCGACTTTGAGAATTTaaaa                                          >  2:2123355/1‑100 (MQ=255)
                                                     gATGCATCGTTCCTGCATCTGCAACCGATGGTGGCGATTGTCACCAATATCGAAGCCGACCACATGGATACCTACCAGGGCGACTTTGAGAATTTAAAAc                                         <  2:662990/100‑1 (MQ=255)
                                                     gATGCAGCGGGCCGGCATCTGCAACCGATGGTGGCGATTGTCATCAATATCGAAGCCGACCACATGGATACCTACCAGGGCGACTTTGAGAATTTAAAAc                                         <  1:2156641/100‑1 (MQ=255)
                                                      aTGCATCGTTCCTGCATCTGCAACCGATGGTGGCGATTGTCATCAATATCGAAGCCGACCACATGGATACCTACCAGGGCGACTTTGAGAATTTAAAACa                                        >  1:115754/1‑100 (MQ=255)
                                                       tGCATCGTTCCTGCATCTGCAACCGATGGTGGCGATTGTCATCAATATCGAAGCCGACCACATGGATACCTACCAGGGCGACTTTGAGAATTTAAAACAg                                       >  1:432015/1‑100 (MQ=255)
                                                       tGCATCGTTCCTGCATCTGCAACCGATGGTGGCGATTGTCATCAATATCGAAGCCGACCACATGGATACCTACCAGGGCGACTTTGAGAATTTAAAACAg                                       >  2:1469918/1‑100 (MQ=255)
                                                       tGCATCGTTCCTGCATCTGCAACCGATGGTGGCGATTGTCATCAATATCGAAGCCGACCACATGGATACCTACCAGGGCGACTTTGAGAATTTAAAACAg                                       >  1:553508/1‑100 (MQ=255)
                                                       tGCATCGTTCCTGCATCTGCAACCGAGGGTGGCGATTGTCATCAATATCGAAGCCGACCACATGGCTACCTACCAGGGCGACTTTGAGAATTTAAAACAg                                       >  1:1478291/1‑100 (MQ=255)
                                                        gCATCGTTCCTGCATCTGCAACCGATGGTGGCGATTGTCATCAATATCGAAGCCGACCACATGGATACCTACCAGGGCGACTTTGAGAATTTAAAACAGa                                      >  2:1127092/1‑100 (MQ=255)
                                                          aTCGTTCCTGCATCTGCAACCGATGGTGGCGATTGTCATCAATATCGAAGCCGACCACATGGATACCTACCAGGGCGACTTTGAGAATTTAAAACAGACt                                    >  1:666518/1‑100 (MQ=255)
                                                           tCGTTCCTGCATCTGCAACCGATGGTGGCGATTGTCATCAATATCGAAGCCGACCACATGGATACCTACCAGGGCGACTTTGAGAATTTAAAACAGACtt                                   >  2:1903599/1‑100 (MQ=255)
                                                             gTTCCGGCGTCTGCGACCGATGGTGGCGATTGTCATCAATATCGAAGCCGACCACATGGATACCTACCAGGGCGACTTTGAGAATTTAAAACAGACtttt                                 <  1:33677/100‑1 (MQ=255)
                                                                ccTGCATCTGCAACCGATGGTGGCGATTGTCATCAATATCGAAGCCGACCACATGGATACCTACCAGGGCGACTTTGAGAATTTAAAACAGACTTttatt                              >  2:1178619/1‑100 (MQ=255)
                                                                ccTGCATCTGCAACCGATGGTGGCGATTGTCATCAATATCGAAGCCGACCACATGGATACCTACCAGGGCGACTTTGAGAATTTAAAACAGACTTttatt                              <  1:1871509/100‑1 (MQ=255)
                                                                ccTGCATCTGCAACCGATGGTGGCGATTGTCATCAATATCGAAGCCGACCACATGGATACCTACCAGGGCGACTTTGAGAATTTAAAACAGACTTttatt                              <  1:1871591/100‑1 (MQ=255)
                                                                 cTGCATCTGCAACCGATGGTGGCGATTGTCATCAATATCGAAGCCGACCACATGGCTAACTACCAGGGCGACTTTGAGAATTTAAAACAGACTTttctta                             <  2:641303/100‑1 (MQ=255)
                                                                       cTGCAACCGATGGTGGCGATTGTCATCAATATCGAAGCCGACCACATGGATACCTACCAGGGCGACTTTGAGAATTTAAAACAGACTTTTATTAATTTTc                       <  2:76703/100‑1 (MQ=255)
                                                                       cTGCAACCGATGGTGGCGATTGTCATCAATATCGAAGCCGACCACATGGATACCTACCAGGGCGACTTTGAGAATTTAAAACAGACTTTTATTAATTTTc                       <  2:76487/100‑1 (MQ=255)
                                                                       cTGCAACCGATGGTGGCGATTGTCATCAATATCGAAGCCGACCACATGGATACCTACCAGGGCGACTTTGAGAATTTAAAACAGACTTTTATTAATTTTc                       >  1:358569/1‑100 (MQ=255)
                                                                        tGCAACCGATGGTGGCGATTGTCATCAATATCGAAGCCGACCACATGGATACCTACCAGGGCGACTTTGAGAATTTAAAACAGACTTTTATTAATTTTCt                      >  1:74977/1‑100 (MQ=255)
                                                                        tGCAACCGATGGTGGCGATCGTCATCAATATCGAAGCCGACCACATGGATACCTACCAGGGCGACTTTGCGCATTTACAACAGACATTTATTAAtttttt                      >  2:2371238/1‑98 (MQ=255)
                                                                          cAACCGATGGTGGCGATTGTCATCAATATCGAAGCCGACCACATGGATACCTACCAGGGCGACTTTGAGAATTTAAAACAGACTTTTATTAATTTTCTGc                    >  2:2344254/1‑100 (MQ=255)
                                                                          cAACCGATGGTGGCGATTGTCATCAATATCGAAGCCGACCACATGGATACCTACCAGGGCGACTTTGAGAATTTAAAACAGACTTTTATTAATTTTCTGc                    <  2:908690/100‑1 (MQ=255)
                                                                           aaCCGATGGTGGCGATTGTCATCAATATCGAAGCCGACCACATGGATACCTACCAGGGCGACTTTGAGAATTTAAAACAGACTTTTATTAATTTTCTGca                   <  1:1062239/100‑1 (MQ=255)
                                                                            aCCGATGGGGGCGATTGTCATCAATATCGAAGCCGACCACATGGATACCTACCAGGGCGACTTTGAGAATTTAAAACAGACTTTTATTAATTTTCTGcac                  <  1:100922/100‑1 (MQ=255)
                                                                             ccGTTGGTGTTGATTGTCTTTAATATCTAATTTGATCACATGGATACCTACCTGGGCGACTTTGAGAATTTAAAACAGACTTTTATTAATTTTCTGcaca                 <  1:1202705/100‑1 (MQ=255)
                                                                                   gtggGGATTGTCATCAATATCGAAGCCGACCACATGGATACCTACCAGGGCGACTTTGAGAATTTAAAACAGCCTTTTATTAATTTTCTGCACAACCTGc           >  2:2198549/1‑100 (MQ=255)
                                                                                   gtggCGATTGTCATCAATATCGAAGCCGACCACATGGATACCTACCAGGGCGACTTTGAGAATTTAAAACAGACTTTTGTTAATTTTCTGCACAACCTGc           >  1:611873/1‑100 (MQ=255)
                                                                                    tggCGATTGTCATCAATATCGAAGCCGACCACNTGGATACCTACCAGGGCGACTTTGAGAATTTAAATCAGACTTTTATTAATTTTCTGCACAACCTGcn          <  2:1855990/100‑2 (MQ=255)
                                                                                    tggCGATTGTCATCAATATCGAAGCCGACCACATGGATACCTACCAGGGCGACTTTGAGCCTTTAAAACAGACTTTTATTACTTTTCTGCACAACCTgac          >  2:1354969/1‑98 (MQ=255)
                                                                                     ggCGATTGTCATCAATATCGAAGCCGACCACATGGATAACTACTAGGGCGACTTTGAGAATTTAAAACAGACCTTTATTAAGTTTCTGCACAACCTGCCg         >  2:2107304/1‑100 (MQ=255)
                                                                                     ggCGACTGTCATCAATATCGAAGCCGACCACATGGATACCTACCAGGACGACTTTGAGAATTTAAAACAGACTTTTATTAATTTTCTGCACAACCTGCCg         >  2:702641/1‑100 (MQ=255)
                                                                                         aTTGTCATCAATATCGAAGCCGACCACATGGATACCTACCAGGGCGACTTTGAGAATTTAAAACAGACTTTTATTAATTTTCTGCACAACCTGCCGtttt     <  2:1631979/100‑1 (MQ=255)
                                                                                           tGTCATCAATATCGAAGCCGACCACATGGATACCTACCAGGGCGACTTTGAGAATTTAAAACAGACTTTTATTAATTTTCTGCACAACCTGCCGTTTTAc   <  1:1196384/100‑1 (MQ=255)
                                                                                           tGTCATCAATATCGAAGCCGACCACATGGATACCTACCAGGGCGACTTTGAGAATTTAAAACAGACTTTTATTAATTTTCTGCACAACCTGCCGTTTTAc   >  1:1987949/1‑100 (MQ=255)
                                                                                            gTCATCAATATCGAAGCCGACCACATGGATACCTACCAGGGCGACTTTGAGAATTTAAAACAGACTTTTATTAATTTTCTGCAGAACCTGCCGTTTTACg  >  1:2437228/1‑100 (MQ=255)
                                                                                            gTCATCAATATCGAAGCCGACCACATGGATACCTACCAGGGCGACTTTGAGAATTTAAAACAGACTTTTATTAATTTTCTGCACAACCTGCCGTTTTACg  <  2:397951/100‑1 (MQ=255)
                                                                                            gTCATCAATATCGAAGCCGACCACATGGATACCTACCAGGGCGACTTTGAGAATTTAAAACAGACTTTTATTAATTTTCTGCACAACCTGCCGTTTTACg  <  2:1763149/100‑1 (MQ=255)
                                                                                                |                                                                                               
TTCATGCGCGTTTGGGGCATGGTCGGTACCTGATTGCCGAAGCAGATGAGAGTGATGCATCGTTCCTGCATCTGCAACCGATGGTGGCGATTGTCACCAATATCGAAGCCGACCACATGGATACCTACCAGGGCGACTTTGAGAATTTAAAACAGACTTTTATTAATTTTCTGCACAACCTGCCGTTTTACG  >  NC_000913/101246‑101437

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: