Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,396,687 T→G I208S (ATC→AGC)  amiB → N‑acetylmuramoyl‑l‑alanine amidase II

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,396,6870TG100.0% 215.4 / NA 62I208S (ATC→AGC) amiBN‑acetylmuramoyl‑l‑alanine amidase II
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base G (32/30);  total (32/30)

AATCCGTTTAAAACGGAAAGTAACCGCACTACGGGTGTTATCAGCAGTAATACGGTAACGCGTCCGGCAGCGCGCGCGACGGCTAACACTGGCGATAAAATTATCATCGCTATTGATGCCGGACACGGCGGTCAGGACCCTGGCGCTATCGGCCCCGGTGGTACGCGGGAGAAAAATGTCACCATCGCCATCGCGCGTAAATTGCGTACTTTGCTCAATGACGATCCGATGTTTAAAGGCGTTTTAACCCGTGACGGGGATTACTTTATCTCGGTGATGGGGCGCAGTGATG  >  NC_000913/4396539‑4396830
                                                                                                                                                    |                                                                                                                                               
aaTCCGTTTAAAACGGAAAGTAACCGCACTACGGGTGTTATCAGCAGTAATACGGTAACGCGTCCGGCAGCGCGCGCGACGGCTAACACTGGCGATAAAATTATCATCGCTATTGATGCCGGACACGGCGGTCAGGACCCTGGCGCTAGCg                                                                                                                                               <  2:636523/151‑1 (MQ=255)
aaTACGTTTAAAACGGAAAGTAACCGCACTACGGGTGTTATCAGCAGTAATACGGTAACGCGTCCGGCAGCGCGCGCGACGGCTAACACTGGCGATAAAATTATCATCGCTATTGATGCCGGACACGGCGGTCAGGACCCTGGCGCTAGCg                                                                                                                                               >  2:609803/1‑151 (MQ=255)
   ccGTTTAAAACGGAAAGTAACCGCACTACGGGTGTTATCAGCAGTAATACGGTAACGCGTCCGGCAGCGCGCGCGACGGCTAACACTGGCGATAAAATTATCATCGCTATTGATGCCGGACACGGCGGTCAGGACCCTGGCGCCAGCgggc                                                                                                                                            >  2:300418/1‑149 (MQ=255)
    cGTTTAAAACGGAAAGTAACCGCACTACGGGTGTTATCAGCAGTAATACGGTAACGCGTCCGGCAGCGCGCGCGACGGCTAACACTGGCGATAAAATTATCATCGCTATTGATGCCGGACACGGCGGTCAGGACCCTGGCGCTAGCGGccc                                                                                                                                           >  2:772675/1‑151 (MQ=255)
      tttAAAACGGAAAGTAACCGCACTACGGGTGTTATCAGCAGTAATACGGTAACGCGTCCGGCAGCGCGCGCGACGGCTAACACTGGCGATAAAATTATCATCGCTATTGATGCCGGACACGGCGGTCAGGACCCTGGCGCTAGCGGcccc                                                                                                                                          >  2:908067/1‑150 (MQ=255)
      tttAAAACGGAAAGTAACCGCACTACGGGTGTTATCAGCAGTAATACGGTAACGCGTCCGGCAGCGCGCGCGACGGCTAACACTGGCGATAAAATTATCATCGCTATTGATGCCGGACACGGCGGTCAGGACCCTGGCGCTAGCGGCCCCg                                                                                                                                         >  1:877886/1‑151 (MQ=255)
       ttAAAACGGAAAGTAACCGCACTACGGGTGTTATCAGCAGTAATACGGTAACGCGTCCGGCAGCGCGCGCGACGGCTAACACTGGCGATAAAATTATCATCGCTATTGATGCCGGACACGGCGGTCAGGACCCTGGCGCTAGCGGCCCCg                                                                                                                                         <  2:708849/150‑1 (MQ=255)
        tAAAACGGAAAGTAACCGCACTACGGGTGTTATCAGCAGTAATACGGTAACGCGTCCGGCAGCGCGCGCGACGGCTAACACTGGCGATAAAATTATCATCGCTATTGATGCCGGACACGGCGGTCAGGACCCTGGCGCTAGCGGCCCCgg                                                                                                                                        >  1:551934/1‑150 (MQ=255)
           aaCGGAAAGTAACCGCACTACGGGTGTTATCAGCAGTAATACGGTAACGCGTCCGGCAGCGCGCGCGACGGCTAACACTGGCGATAAAATTATCATCGCTATTGATGCCGGACACGGCGGTCAGGACCCTGGCGCTAGCGGCCCCggtgg                                                                                                                                     >  2:173900/1‑150 (MQ=255)
                 aaGTAACCGCACTACGGGTGTTATCAGCAGTAATACGGTAACGCGTCCGGCAGCGCGCGCGACGGCTAACACTGGCGATAAAATTATCATCGCTATTGATGCCGGACACGGCGGTCAGGACCCTGGCGCTAGCGGCCCCGGTGGTACGCgg                                                                                                                              <  1:61462/151‑1 (MQ=255)
                  aGTAACCGCACTACGGGTGTTATCAGCAGTAATACGGTAACGCGTCCGGCAGCGCGCGCGACGGCTAACACTGGCGATAAAATTATCATCGCTATTGATGCCGGACACGGCGGTCAGGACCCTGGCGCTAGCGGCCCCGGTGGTACGCggg                                                                                                                             >  1:713180/1‑151 (MQ=255)
                   gTAACCGCACTACGGGTGTTATCAGCAGTAATACGGTAACGCGTCCGGCAGCGCGCGCGACGGCTAACACTGGCGATAAAATTATCATCGCTATTGATGCCGGACACGGCGGTCAGGACCCTGGCGCTAGCGGCCCCGGTGGTACGCGGga                                                                                                                            <  2:773334/151‑1 (MQ=255)
                    tAACCGCACTACGGGTGTTATCAGCAGTAATACGGTAACGCGTCCGGCAGCGCGCGCGACGGCTAACACTGGCGATAAAATTATCATCGCTATTGATGCCGGACACGGCGGTCAGGACCCTGGCGCTAGCGGCCCCGGTGGTACGCggg                                                                                                                             >  1:303219/1‑149 (MQ=255)
                         gCACTACGGGTGTTATCAGCAGTAATACGGTAACGCGTCCGGCAGCGCGCGCGACGGCTAACACTGGCGATAAAATTATCATCGCTATTGATGCCGGACACGGCGGTCAGGACCCTGGCGCTAGCGGCCCCGGTGGTACGCGGGAGaaaaa                                                                                                                      <  2:723108/151‑1 (MQ=255)
                         gCACTACGGGTGTTATCAGCAGTAATACGGTAACGCGTCCGGCAGCGCGCGCGACGGCTAACACTGGCGATAAAATTATCATCGCTATTGATGCCGGACACGGCGGTCAGGACCCTGGCGCTAGCGGCCCCGGTGGTACGCGGGAGaaaa                                                                                                                       >  2:782033/1‑150 (MQ=255)
                          cACTACGGGTGTTATCAGCAGTAATACGGTAACGCGTCCGGCAGCGCGCGCGACGGCTAACACTGGCGATAAAATTATCATCGCTATTGATGCCGGACACGGCGGTCAGGACCCTGGCGCTAGCGGCCCCGGTGGTACGCGGGAGaaaaa                                                                                                                      >  1:95543/1‑150 (MQ=255)
                           aCTACGGGTGTTATCAGCAGTAATACGGTAACGCGTCCGGCAGCGCGCGCGACGGCTAACACTGGCGATAAAATTATCATCGCTATTGATGCCGGACACGGCGGTCAGGACCCTGGCGCTAGCGGCCCCGGTGGTACGCGGGAGAAAAAt                                                                                                                     >  1:663152/1‑150 (MQ=255)
                             tACGGGTGTTATCAGCAGTAATACGGTAACGCGTCCGGCAGCGCGCGCGACGGCTAACACTGGCGATAAAATTATCATCGCTATTGATGCCGGACACGGCGGTCAGGACCCTGGCGCTAGCGGCCCCGGTGGTACGCGGGAGAAAAATGTc                                                                                                                  >  2:675466/1‑151 (MQ=255)
                              aCGGGTGTTATCAGCAGTAATACGGTAACGCGTCCGGCAGCGCGCGCGACGGCTAACACTGGCGATAAAATTATCATCGCTATTGATGCCGGACACGGCGGTCAGGACCCTGGCGCTAGCGGCCCCGGTGGTACGCGGGAGAAAAATGTc                                                                                                                  <  1:707042/150‑1 (MQ=255)
                               cGGGTGTTATCAGCAGTAATACGGTAACGCGTCCGGCAGCGCGCGCGACGGCTAACACTGGCGATAAAATTATCATCGCTATTGATGCCGGACACGGCGGTCAGGACCCTGGCGCTAGCGGCCCCGGTGGTACGCGGGAGAAAAATGTCAc                                                                                                                >  2:751418/1‑151 (MQ=255)
                                         cagcagTAATACGGTAACGCGTCCGGCAGCGCGCGCGACGGCTAACACTGGCGATAAAATTATCATCGCTATTGATGCCGGACACGGCGGTCAGGACCCTGGCGCTAGCGGCCCCGGTGGTACGCGGGAGAAAAATGTCAccatcgccatc                                                                                                      >  1:189826/1‑151 (MQ=255)
                                             agTAATACGGTAACGCGTCCGGCAGCGCGCGCGACGGCTAACACTGGCGATAAAATTATCATCGCTATTGATGCCGGACACGGCGGTCAGGACCCTGGCGCTAGCGGCCCCGGTGGTACGCGGGAGAAAAATGTCACCATCGCCATcgcg                                                                                                   >  2:614225/1‑150 (MQ=255)
                                              gTAATACGGTAACGCGTCCGGCAGCGCGCGCGACGGCTAACACTGGCGATAAAATTATCATCGCTATTGATGCCGGACACGGCGGTCAGGACCCTGGCGCTAGCGGCCCCGGTGGTACGCGGGAGAAAAATGTCACCATCGCCATcgcgcg                                                                                                 >  2:885793/1‑151 (MQ=255)
                                               tAATACGGTAACGCGTCCGGCAGCGCGCGCGACGGCTAACACTGGCGATAAAATTATCATCGCTATTGATGCCGGACACGGCGGTCAGGACCCTGGCGCTAGCGGCCCCGGGGGTACGCGGGAGAAAAATGTCACCATCGCCATCGCGCGt                                                                                                >  2:205823/1‑151 (MQ=255)
                                               tAATACGGTAACGCGTCCGGCAGCGCGCGCGACGGCTAACACTGGCGATAAAATTATCATCGCTATTGATGCCGGACACGGCGGGCAGGACCCTTGCGCTAGCGGCCCCGGGGGTAAGCGGGAGAAAAATGTCACCATCGCCCTcgcgc                                                                                                  >  2:435951/1‑149 (MQ=255)
                                                   aCGGTAACGCGTCCGGCAGCGCGCGCGACGGCTAACACTGGCGATAAAATTATCATCGCTATTGATGCCGGACACGGCGGTCAGGACCCTGGCGCTAGCGGCCCCGGTGGTACGCGGGAGAAAAATGTCACCATCGCCATCGCGCGTaaa                                                                                             <  1:761104/150‑1 (MQ=255)
                                                    cGGTAACGCGTCCGGCAGCGCGCGCGACGGCTAACACTGGCGATAAAATTATCATCGCTATTGATGCCGGACACGGCGGTCAGGACCCTGGCGCTAGCGGCCCCGGTGGTACGCGGGAGAAAAATGTCACCATCGCCATCGCGCGTAAAtt                                                                                           <  2:772341/151‑1 (MQ=255)
                                                    cGGTAACGCGTCCGGCAGCGCGCGCGACGGCTAACACTGGCGATAAAATTACCATCGCTATTGATGCCGGACACGGCGGTCAGGACCCTGGCGCTAGCGGCCCCGGTGGTACGCGGGAGAAAAATGTCACCATCGCCATCGCGCGTAAAtt                                                                                           <  2:228586/151‑1 (MQ=255)
                                                      gTAACGCGTCCGGCAGCGCGCGCGACGGCTAACACTGGCGATAAAATTATCATCGCTATTGATGCCGGACACGGCGGTCAGGACCCTGGCGCTAGCGGCCCCGGTGGTACGCGGGAGAAAAATGTCACCATCGCCATCGCGCGTAAAtt                                                                                           <  2:651352/149‑1 (MQ=255)
                                                         aCGCGTCCGGCAGCGCGCGCGACGGCTAACACTGGCGATAAAATTATCATCGCTATTGATGCCGGACACGGCGGTCAGGACCCTGGCGCTAGCGGCCCCGGTGGTACGCGGGAGAAAAATGTCACCATCGCCATCGCGCGTAAATTGCg                                                                                        >  2:21412/1‑149 (MQ=255)
                                                             gTCCGGCAGCGCGCGCGACGGCTAACACTGGCGATAAAATTATCATCGCTATTGATGCCGGACACGGCGGTCAGGACCCTGGCGCTAGCGGCCCCGGTGGTACGCGGGAGAAAAATGTCACCATCGCCATCGCGCGTAAATTGCGTACttt                                                                                  >  2:451324/1‑151 (MQ=255)
                                                               ccGGCAGCGCGCGCGACGGCTAACACTGGCGATAAAATTATCATCGCTATTGATGCCGGACACGGCGGTCAGGACCCTGGCGCTAGCGGCCCCGGTGGTACGCGGGAGAAAAAAATCACCATCGCCATCGCGCGGAAATTGCGATCttt                                                                                  >  2:489959/1‑149 (MQ=255)
                                                                cGGCAGCGCGCGCGACGGCTAACACTGGCGATAAAATTATCATCGCTATTGATGCCGGACACGGCGGTCAGGACCCTGGCGCTAGCGGCCCCGGTGGTACGCGGGAGAAAAATGTCACCATCGCCATCGCGCGTAAATTGCGTACTTTGCt                                                                               >  2:141677/1‑151 (MQ=255)
                                                                    aGCGCGCGCGACGGCTAACACTGGCGATAAAATTATCATCGCTATTGATGCCGGACACGGCGGTCAGGACCCTGGCGCTAGCGGCCCCGGTGGTACGCGGGAGAAAAATGTCACCATCGCCATCGCGCGTAAATTGCGTACTTTGCTCAAt                                                                           >  1:442301/1‑151 (MQ=255)
                                                                    aGCGCGCGCGACGGCTAACACTGGCGATAAAATTATCATCGCTATTGATGCCGGACACGGCGGTCAGGACCCTGGCGCTAGCGGCCCCGGTGGTACGCGGGAGAAAAATGTCACCATCGCCATCGCGCGTAAATTGCGTACTTTGCTCAAt                                                                           <  1:206462/151‑1 (MQ=255)
                                                                      cgcgcgcgACGGCTAACACTGGCGATAAAATTATCATCGCTATTGATGCCGGACACGGCGGTCAGGACCCTGGCGCTAGCGGCCCCGGTGGTACGCGGGAGAAAAATGTCACCATCGCCATCGCGCGTAAATTGCGTACTTTGCTCAATGa                                                                         <  2:619889/151‑1 (MQ=255)
                                                                       gcgcgcgACGGCTAACACTGGCGATAAAATTATCATCGCTATTGATGCCGGACACGGCGGTCAGGACCCTGGCGCTAGCGGCCCCGGTGGTACGCGGGAGAAAAATGTCACCATCGCCATCGCGCGTAAATTGCGTACTTTGCTCAATGAc                                                                        >  1:825450/1‑151 (MQ=255)
                                                                         gcgcgACGGCTAACACTGGCGATAAAATTATCATCGCTATTGATGCCGGACACGGCGGTCAGGACCCTGGCGCTAGCGGCCCCGGTGGTACGCGGGAGAAAAATGTCACCATCGCCATCGCGCGTAAATTGCGTACTTTGCTCAATGACGa                                                                      <  2:831289/151‑1 (MQ=255)
                                                                            cgACGGCTAACACTGGCGATAAAATTATCATCGCTATTGATGCCGGACACGGCGGTCAGGACCCTGGCGCTAGCGGCCCCGGTGGTACGCGGGAGAAAAATGTCACCATCGCCATCGCGCGTAAATTGCGTACTTTGCTCAATGACGATcc                                                                   <  1:725820/151‑1 (MQ=255)
                                                                                 gCTAACACTGGCGATAAAATTATCATCGCTATTGATGCCGGACACGGCGGTCAGGACCCTGGCGCTAGCGGCCCCGGTGGTACGCGGGAGAAAAATGTCACCATCGCCATCGCGCGTAAATTGCGTACTTTGCTCAATGACGATCCGATGt                                                              <  1:93505/151‑1 (MQ=255)
                                                                                    aaCACTGGCGATAAAATTATCATCGCTATTGATGCCGGACACGGCGGTCAGGACCCTGGCGCTAGCGGCCCCGGTGGTACGCGGGAGAAAAATGTCACCATCGCCATCGCGCGTAAATTGCGTACTTTGCTCAATGACGATCCGATGtt                                                             <  1:413854/149‑1 (MQ=255)
                                                                                              aTAAAATTATCATCGCTATTGATGCCGGACACGGCGGTCAGGACCCTGGCGCTAGCGGCCCCGGTGGTACGCGGGAGAAAAATGTCACCATCGCCATCGCGCGTAAATTGCGTACTTTGCTCAATGACGATCCGATGTTTAAAGGCGtttt                                                 >  1:431206/1‑151 (MQ=255)
                                                                                                aaaaTTATCATCGCTATTGATGCCGGACACGGCGGTCAGGACCCTGGCGCTAGCGGCCCCGGTGGTACGCGGGAGAAAAATGTCACCATCGCCATCGCGCGTAAATTGCGTACTTTGCTCAATGACGATCCGATGTTTAAAGGCGTTTTaa                                               <  1:60569/151‑1 (MQ=255)
                                                                                                  aaTTATCATCGCTATTGATGCCGGACACGGCGGTCAGGACCCTGGCGCTAGCGGCCCCGGTGGTACGCGGGAGAAAAATGTCACCATCGCCATCGCGCGTAAATTGCGTACTTTGCTCAATGACGATCCGATGTTTAAAGGCGTTTTAAcc                                             >  1:10393/1‑151 (MQ=255)
                                                                                                  aaTTATCATCGCTATTGATGCCGGACACGGCGGTAAGGACCCTGGCGCTAGCGGCCCCGGTGGTACGCGGGAGAAAAATGTCACCATCGCCATCGCGCGTAAATTGCGTACTTTGCTCAATGACGATCCGATGTTTAAAGGCGTTTTAAcc                                             <  1:627363/151‑1 (MQ=255)
                                                                                                        catcGCTATTGATGCCGGACACGGCGGTCAGGACCCTGGCGCTAGCGGCCCCGGTGGTACGCGGGAGAAAAATGTCACCATCGCCATCGCGCGTAAATTGCGTACTTTGCTCAATGACGATCCGATGTTTAAAGGCGTTTTAACCCGTGAc                                       >  1:802980/1‑151 (MQ=255)
                                                                                                           cGCTATTGATGCCGGACACGGCGGTCAGGACCCTGGCGCTAGCGGCCCCGGTGGTACGCGGGAGAAAAATGTCACCATCGCCATCGCGCGTAAATTGCGTACTTTGCTCAATGACGATCCGATGTTTAAAGGCGTTTTAACCCGTGACggg                                    <  1:536891/151‑1 (MQ=255)
                                                                                                               aTTGATGCCGGACACGGCGGTCAGGACCCTGGCGCTAGCGGCCCCGGTGGTACGCGGGAGAAAAATGTCACCATCGCCATCGCGCGTAAATTGCGTACTTTGCTCAATGACGATCCGATGTTTAAAGGCGTTTTAACCCGTGACGGGGAtt                                <  1:909635/151‑1 (MQ=255)
                                                                                                                   aTGCCGGACACGGCGGTCAGGACCCTGGCGCTAGCGGCCCCGGTGGTACGCGGGAGAAAAATGTCACCATCGCCATCGCGCGTAAATTGCGTACTTTGCTCAATGACGATCCGATGTTTAAAGGCGTTTTAACCCGTGACGGGGATTACtt                            >  2:616959/1‑151 (MQ=255)
                                                                                                                        ggACACGGCGGTCAGGACCCTGGCGCTAGCGGCCCCGGTGGTACGCGGGAGAAAAATGTCACCATCGCCATCGCGCGTAAATTGCGTACTTTGCTCAATGACGATCCGATGTTTAAAGGCGTTTTAACCCGTGACGGGGATTACTTTAtct                       <  2:280396/151‑1 (MQ=255)
                                                                                                                        ggACACGGCGGTCAGGACCCTGGCGCTAGCGGCCCCGGTGGTACGCGGGAGAAAAATGTCACCATCGCCATCGCGCGTAAATTGCGTACTTTGCTCAATGACGATCCGATGTTTAAAGGCGTTTTAACCCGTGACGGGGATTACTTTAtct                       <  1:130175/151‑1 (MQ=255)
                                                                                                                         gACACGGCGGTCAGGACCCTGGCGCTAGCGGCCCCGGTGGTACGCGGGAGAAAAATGTCACCATCGCCATCGCGCGTAAATTGCGTACTTTGCTCAATGACGATCCGATGTTTAAAGGCGTTTTAACCCGTGACGGGGATTACTTTAtctc                      >  1:807693/1‑151 (MQ=255)
                                                                                                                          acacGGCGGTCAGGACCCTGGCGCTAGCGGCCCCGGTGGTACGCGGGAGAAAAATGTCACCATCGCCATCGCGCGTAAATTGCGTACTTTGCTCAATGACGATCCGATGTTTAAAGGCGTTTTAACCCGTGACGGGGATTACTTTATCTCg                     <  2:410418/151‑1 (MQ=255)
                                                                                                                             cggcggTCAGGACCCTGGCGCTAGCGGCCCCGGTGGTACGCGGGAGAAAAATGTCACCATCGCCATCGCGCGTAAATTGCGTACTTTGCTCAATGACGATCCGATGTTTAAAGGCGTTTTAACCCGTGACGGGGATTACTTTATCTCGGTg                  <  2:783631/151‑1 (MQ=255)
                                                                                                                             cggcggTCAGGACCCTGGCGCTAGCGGCCCCGGTGGTACGCGGGAGAAAAATGTCACCATCGCCATCGCGCGTAAATTGCGTACTTTGCTCAATGACGATCCGATGTTTAAAGGCGTTTTAACCCGTGACGGGGATTACTTTATCTCGGTg                  <  2:890081/151‑1 (MQ=255)
                                                                                                                                cggTCAGGACCCTGGCGCTAGCGGCCCCGGTGGTACGCGGGAGAAAAATGTCACCATCGCCATCGCGCGTAAATTGCGTACTTTGCTCAATGACGATCCGATGTTTAAAGGCGTTTTAACCCGTGACGGGGATTACTTTATCTCGGTGATg               <  1:376872/151‑1 (MQ=255)
                                                                                                                                    cAGGACCCTGGCGCTAGCGGCCCCGGTGGTACGCGGGAGAAAAATGTCACCATCGCCATCGCGCGTAAATTGCGTACTTTGCTCAATGACGATCCGATGTTTAAAGGCGTTTTAACCCGTGACGGGGATTACTTTATCTCGGTGATGGGgc           <  1:326040/151‑1 (MQ=255)
                                                                                                                                       gACCCTGGCGCTAGCGGCCCCGGTGGTACGCGGGAGAAAAATGTCACCATCGCCATCGCGCGTAAATTGCGTACTTTGCTCAATGACGATCCGATGTTTAAAGGCGTTTTAACCCGTGACGGGGATTACTTTATCTCGGTGATGGGgcgc         >  1:252233/1‑150 (MQ=255)
                                                                                                                                        aCCCTGGCGCTAGCGGCCCCGGTGGTACGCGGGAGAAAAATGTCACCATCGCCATCGCGCGTAAATTGCGTACTTTGCTCAATGACGATCCGATGTTTAAAGGCGTTTTAACCCGTGACGGGGATTACTTTATCTCGGTGATGGGgcgc         >  2:231052/1‑149 (MQ=255)
                                                                                                                                             ggCGCTAGCGGCCCCGGTGGTACGCGGGAGAAAAATGTCACCATCGCCATCGCGCGTAAATTGCGTACTTTGCTCAATGACGATCCGATGTTTAAAGGCGTTTTAACCCGTGACGGGGATTACTTTATCTCGGTGATGGGGCGCAGTGa    <  2:51643/149‑1 (MQ=255)
                                                                                                                                             ggCGCTAGCGGCCCCGGTGGTACGCGGGAGAAAAATGTCACCATCGCCATCGCGCGTAAATTGCGTACTTTGCTCAATGACGATCCGATGTTTAAAGGCGTTTTAACCCGTGACGGGGATTACTTTATCTCGGTGATGGGGCGCAGTGa    <  2:879819/149‑1 (MQ=255)
                                                                                                                                             ggCGCTAGCGGCCCCGGTGGTACGCGGGAGAAAAATGTCACCATCGCCATCGCGCGTAAATTGCGTACTTTGCTCAATGACGATCCGATGTTTAAAGGCGTTTTAACCCGTGACGGGGATTACTTTATCTCGGTGATGGGGCGCAGTGAtg  <  1:614240/151‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                               
AATCCGTTTAAAACGGAAAGTAACCGCACTACGGGTGTTATCAGCAGTAATACGGTAACGCGTCCGGCAGCGCGCGCGACGGCTAACACTGGCGATAAAATTATCATCGCTATTGATGCCGGACACGGCGGTCAGGACCCTGGCGCTATCGGCCCCGGTGGTACGCGGGAGAAAAATGTCACCATCGCCATCGCGCGTAAATTGCGTACTTTGCTCAATGACGATCCGATGTTTAAAGGCGTTTTAACCCGTGACGGGGATTACTTTATCTCGGTGATGGGGCGCAGTGATG  >  NC_000913/4396539‑4396830

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 18 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: