Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,975,807 A→G L372P (CTG→CCG)  lysA ← diaminopimelate decarboxylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,975,8070AG98.1% 488.4 / ‑2.5 156L372P (CTG→CCG) lysAdiaminopimelate decarboxylase
Reads supporting (aligned to +/- strand):  ref base A (1/2);  new base G (76/77);  total (77/79)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

Only 100 of 158 total aligned reads displayed.

CGCCTGACCATTATCAAACAGAACTTCTGGTAACAGCGGACGGCTATTGTAGTTGGATGACATTGATGCGCCATATGCCCCTGTATCATGCAGTACCAGATAATCACCTGCCTTCACTTCCGGCAAGGCGCGGGTTTCAACATTTCCCCCTTCCTGCTGGGTAAAGACATCGCCCGATTCA  >  NC_000913/2975716‑2975896
                                                                                           |                                                                                         
cgcCTGACCATTATCAAACAGAACTTCTGGTAACAGCGGACGGCTATTGTAGTTGGATGACATTGATGCGCCATATGCCCCTGTATCATGCAg                                                                                          <  2:2007624/93‑1 (MQ=255)
   cTGACCATTATCAAACAGAACTTCTGGTAACAGCGGACGGCTATTGTAGTTGGATGACATTGATGCGCCATATGCCCCTGTATCATGCGGTAc                                                                                       <  1:933523/93‑1 (MQ=255)
    tGACCATTATCAAACAGAACTTCTGGTAACAGCGGACGGCTATTGTAGTTGGATGACATTGATGCGCCATATGCCCCTGTATCATGCGGTAcc                                                                                      >  1:3032640/1‑93 (MQ=255)
    tGACCATTATCAAACAGAACTTCTGGTAACAGCGGACGGCTATTGTAGTTGGATGACATTGATGCGCCATATGCCCCTGTATCATGCGGTAcc                                                                                      <  2:1845448/93‑1 (MQ=255)
    tGACCATTATCAAACAGAACTTCTGGTAACAGCGGACGGCTATTGTAGTTGGATGACATTGATGCGCCATATGCCCCTGTATCATGCGGTAcc                                                                                      >  1:2782587/1‑93 (MQ=255)
     gACCATTATCAAACAGAACTTCTGGTAACAGCGGACGGCTATTGTAGTTGGATGACATTGATGCGCCATATGCCCCTGTATCATGCGGTACCa                                                                                     >  2:4340800/1‑93 (MQ=255)
      aCCATTATCAAACAGAACTTCTGGTAACAGCGGACGGCTATTGTAGTTGGATGACATTGATGCGCCATATGCCCCTGTATCATGCGGTACCAg                                                                                    <  1:4286942/93‑1 (MQ=255)
      aCCATTATCAAACAGAACTTCTGGTAACAGCGGACGGCTATTGTAGTTGGATGACATTGATGCGCCATATGCCCCTGTATCATGCGGTACCAg                                                                                    >  1:156685/1‑93 (MQ=255)
       ccATTATCAAACAGAACTTCTGGTAACAGCGGACGGCTATTGTAGTTGGATGACATTGATGCGCCATATGCCCCTGTATCATTCGGTACCAGa                                                                                   <  2:466475/93‑1 (MQ=255)
       ccATTATCAAACAGAACTTCTGGTAACAGCGGACGGCTATTGTAGTTGGATGACATTGATGCGCCATATGCCCCTGTATCATGCGGTACCAGa                                                                                   <  1:4691797/93‑1 (MQ=255)
        cATTATCAAACAGAACTTCTGGTAACAGCGGACGGCTATTGTAGTTGGATGACATTGATGCGCCATATGCCCCTGTATCATGCGGTACCAGAt                                                                                  >  2:2037160/1‑93 (MQ=255)
         aTTATCAAACAGAACTTCTGGTAACAGCGGACGGCTATTGTAGTTGGATGACATTGATGCGCCATATGCCCCTGTATCATGCGGTACCAGATa                                                                                 >  1:38856/1‑93 (MQ=255)
          ttATCAAACAGAACTTCTGGTAACAGCGGACGGCTATTGTAGTTGGATGACATTGATGCGCCATATGCCCCTGTATCATGCGGTACCAGATaa                                                                                <  2:4193110/93‑1 (MQ=255)
          ttATCAAACAGAACTTCTGGTAACAGCGGACGGCTATTGTAGTTGGATGACATTGATGCGCCATATGCCCCTGTATCATGCGGTACCAGATaa                                                                                >  2:168070/1‑93 (MQ=255)
          ttATCAAAAAGAACTTCTGGTAACAGCGGACGGCTATTGTAGTTGGATGACATTGATGCGCCATATGCCCCTGTATCATGCGGTACCAGATaa                                                                                <  1:4525390/93‑1 (MQ=255)
           tATCAAACAGAACTTCTGGTAACAGCGGACGGCTATTGTAGTTGGATGACATTGATGCGCCATATGCCCCTGTATCATGCGGTACCAGATAAt                                                                               <  2:291595/93‑1 (MQ=255)
           tATCAAACAGAACTTCTGGTAACAGCGGACGGCTATTGTAGTTGGATGACATTGATGCGCCATATGCCCCTGTATCATGCGGTACCAGATAAt                                                                               <  1:526968/93‑1 (MQ=255)
            aTCAAACAGAACTTCTGGTAACAGCGGACGGCTATTGTAGTTGGATGACATTGATGCGCCATATGCCCCTGTATCATGCGGTACCAGATAATc                                                                              <  2:3036913/93‑1 (MQ=255)
             tCAAACAGAACTTCTGGTAACAGCGGACGGTTATTGTAGTTGGATGACATTGATGCGCCATATGCCCCTGTATCAGGCGGTACCAGATAATCa                                                                             <  1:3200269/93‑1 (MQ=255)
             tCAAACAGAACTTCTGGTAACAGCGGACGGCTATTGTAGTTGGATGACATTGATGCGCCATATGCCCCTGTATCATGCGGTACCAGATAATCa                                                                             >  1:2792518/1‑93 (MQ=255)
             tCAAACAGAACTTCTGGTAACAGCGGACGGCTATTGTAGTTGGATGACATTGATGCGCCATATGCCCCTGTATCATGCGGTACCAGATAATCa                                                                             <  1:3200267/93‑1 (MQ=255)
               aaaCAGAACTTCTGGTAACAGCGGACGGCTATTGTAGTTGGATGACATTGATGCGCCATATGCCCCTGTATCATGCGGTACCAGATAATCAcc                                                                           <  1:3939491/93‑1 (MQ=255)
               aaaCAGAACTTCTGGCAACAGCGGACGGCTATTGTAGTTGGATGACATTGATGCGCCATATGCCCCTGTATCATGCAGTACCAGATAATCAcc                                                                           <  2:2371014/93‑1 (MQ=255)
                aaCAGAACTTCTGGTAACAGCGGACGGCTATTGTAGTTGGATGACATTGATGCGCCATATGCCCCTGTATCATGCGGTACCAGATAATCACCt                                                                          <  1:4124107/93‑1 (MQ=255)
                 aCAGAACTTCTGGTAACAGCGGACGGCTATTGTAGTTGGATGACATTGATGCGCCATATGCCCCTGTATCATGCGGTACCAGATAATCACCTg                                                                         <  1:173216/93‑1 (MQ=255)
                  cAGAACTTCTGGTAACAGCGGACGGCTATTGTAGTTGGATGACATTGATGCGCCATATGCCCCTGTATCATGCGGTACCAGATAATCACCTGc                                                                        >  1:3606401/1‑93 (MQ=255)
                   aGAACTTCTGGTAACAGCGGACGGCTATTGTAGTTGGATGACATTGATGCGCCATATGCCCCTGTATCATGCGGTACCAGATAATCACCTGcc                                                                       <  1:3674772/93‑1 (MQ=255)
                     aaCTTCTGGTAACAGCGGACGGCTATTGTAGTTGGATGACATTGATGCGCCATATGCCCCTGTATCATGCGGTACCAGATAATCACCTGCCtt                                                                     <  1:4143250/93‑1 (MQ=255)
                     aaCTTCTGGTAACAGCGGACGGCTATTGAAGTTGGATGACATTGATGCGCCATATGCCCCTGTATCATGCGGTACCAGATAATCACCTGCCtt                                                                     <  1:4608761/93‑1 (MQ=255)
                        ttCTGGTAACAGCGGACGGCTATTGTAGTTGGATGACATTGATGCGCCATATGCCCCTGTATCATGCAGTACCAGATAATCACCTGCCTTCAc                                                                  >  1:4708436/1‑93 (MQ=255)
                           tGGTAACAGCGGACGGCTATTGTAGTTGGATGACATTGATGCGCCATATGCCCCTGTATCATGCGGTACCAGATAATCACCTGCCTTCACTTc                                                               >  1:2975652/1‑93 (MQ=255)
                           tGGTAACAGCGGACGGCTATTGTAGTTGGATGACATTGATGCGCCATATGCCCCTGTATCATGCGGTACCAGATAATCACCTGCCTTCACTTc                                                               >  1:638838/1‑93 (MQ=255)
                           tGGTAACAGCGGACGGCTATTGTAGTTGGATGACATTGATGCGCCATATGCCCCTGTATCATGCGGTACCAGATAATCACCTGCCTTCACTTc                                                               <  2:438760/93‑1 (MQ=255)
                           tGGTAACAGCGGACGGCTATTGTAGTTGGATGACATTGATGCGCCATATGCCCCTGTATCATGCGGTACCAGATAATCACCTGCCTTCACTTc                                                               <  2:663410/93‑1 (MQ=255)
                           tGGTAACAGCGGACGGCTATTGTAGTTGGATGACATTGATGCGCCATATGCCCCTGTATCATGCGGTACCAGATAATCACCTGCCTTCACTTc                                                               <  2:2852177/93‑1 (MQ=255)
                           tGGTAACAGCGGACGGCTATTGTAGTTGGATGACATTGATGCGCCATATGCCCCTGTATCATGCGGTACCAGATAATCACCTGCCTTCACTTc                                                               <  2:1673878/93‑1 (MQ=255)
                           tGGTAACAGCGGACGGCTATTGTAGTTGGATGACATTGATGCGCCATATGCCCCTGTATCATGCGGTACCAGATAATCACCTGCCTTCACTTc                                                               <  1:2798480/93‑1 (MQ=255)
                            ggTAACAGCGGACGGCTATTGTAGTTGGCTGACATTGATGCGCCATATGCCCCTGTATCATGCGGTACCAGATAATCACCTGCCTTCACTTcc                                                              >  2:4399315/1‑93 (MQ=255)
                            ggTAACAGCGGACGGCTATTGTAGTTGGATGACATTGATGCGCCATATGCCCCTGTATCATGCGGTACCAGATAATCACCTGCCTTCACTTcc                                                              >  2:2758490/1‑93 (MQ=255)
                            ggTAACAGCGGACGGCTATTGTAGTTGGATGACATTGATGCGCCATATGCCCCTGTATCATGCGGTACCAGATAATCACCTGCCTTCACTTcc                                                              >  2:544772/1‑93 (MQ=255)
                            ggTAACAGCGGACGGCTATTGTAGTTGGATGACATTGATGCGCCATATGCCCCTGTATCATGCGGTACCAGATAATCACCTGCCTTCACTTcc                                                              <  2:2134446/93‑1 (MQ=255)
                            ggTAACAGCGGACGGCTATTGTAGTTGGATGACATTGATGCGCCATATGCCCCTGTATCATGCGGTACCAGATAATCACCTGCCTTCACTTcc                                                              >  2:1940462/1‑93 (MQ=255)
                             gTAACAGCGGACGGCTATTGTAGTTGGATGACATTGATGCGCCATATGCCCCTGTATCATGCGGTACCAGATAATCACCTGCCTTCACTTCCg                                                             >  1:1141597/1‑93 (MQ=255)
                             gTAACAGCGGACGGCTATTGTAGTTGGATGACATTGATGCGCCATATGCCCCTGTATCATGCGGTACCAGATAATCACCTGCCTTCACTTCCg                                                             <  1:590639/93‑1 (MQ=255)
                                aCAGCGGACGGCTATTGTAGTTGGATGACATTGATGCGCCATATGCCCCTGTATCATGCGGTACCAGATAATCACCTGCCTTCACTTCCGGCa                                                          >  2:3292838/1‑93 (MQ=255)
                                aCAGCGGACGGCTATTGTAGTTGGATGACATTGATGCGCCATATGCCCCTGTATCATGCGGTACCAGATAATCACCTGCCTTCACTTCCGGCa                                                          >  1:1960679/1‑93 (MQ=255)
                                 cAGCGGACGGCTATTGTAGTTGGATGACATTGATGCGCCATATGCCCCTGTATCATGCGGTACCAGATAATCACCTGCCTTCACTTCCGGCaa                                                         <  1:4491072/93‑1 (MQ=255)
                                 cAGCGGACGGCTATTGTAGTTGGATGACATTGATGCGCCATATGCCCCTGTATCATGCGGTACCAGATAATCACCTGCCTTCACTTCCGGCaa                                                         >  2:1505433/1‑93 (MQ=255)
                                  aGCGGACGGCTATTGTAGTTGGATGACATTGATGCGCCATATGCCCCTGTATCATGCGGTACCAGATAATCACCTGCCTTCACTTCCGGCAAg                                                        <  2:3992260/93‑1 (MQ=255)
                                   gCGGACGGCTATTGTAGTTGGATGACATTGATGCGCCATATGCCCCTGTATCATGCGGTACCAGATAATCACCTGCCTTCACTTCCGGCAAgg                                                       <  2:3319296/93‑1 (MQ=255)
                                   gCGGACGGCTATTGTAGTTGGATGACATTGATGCGCCATATGCCCCTGTATCATGCGGTACCAGATAATCACCTGCCTTCACTTCCGGCAAgg                                                       <  1:798621/93‑1 (MQ=255)
                                   gCGGACGGCTATTGTAGTTGGATGACATTGATGCGCCATATGCCCCTGTATCATGCGGTACCAGATAATCACCTGCCTTCACTTCCGGCAAgg                                                       >  2:3752568/1‑93 (MQ=255)
                                    cGGACGGCTGTTGTAGTTGGATGACATTGATGCGCCATATGCCCCTGTATCATGCGGTACCAGATAATCACCTGGCTTCACTTCCGGCAAGgc                                                      >  1:378207/1‑93 (MQ=255)
                                    cGGACGGCTATTGTAGTTGGATGACATTGATGCGCCATATGCCCCTGTATCATGCGGTACCAGATAATCACCTGCCTTCACTTCCGGCAAGgc                                                      <  1:2318188/93‑1 (MQ=255)
                                     ggACGGCTATTGTAGTTGGATGACATTGATGCGCCATATGCCCCTGTATCATGCGGTACCAGATAATCACCTGCCTTCACTTCCGGCAAGgcg                                                     <  1:2579790/93‑1 (MQ=255)
                                      gACGGCTATTGTAGTTGGATGACATTGATGCGCCATATGCCCCTGTATCATGCGGTACCAGATAATCACCTGCCTTCACTTCCGGCAAGgcgc                                                    >  2:1330295/1‑93 (MQ=255)
                                       aCGGCTATTGTAGTTGGATGACATTGATGCGCCATATGCCCCTGTATCATGCGGTACCAGATAATCACCTGCCTTCACTTCCGGCAAGGcgcg                                                   <  2:4531053/93‑1 (MQ=255)
                                          gCTATTGTAGTTGGATGACATTGATGCGCCATATGCCCCTGTATCATGCGGTACCAGATAATCACCTGCCTTCACTTCCGGCAAGGCGCGGGt                                                >  2:1480512/1‑93 (MQ=255)
                                           cTATTGTAGTTGGATGACATTGATGCGCCATATGCCCCTGTATCATGCGGTACCAGATAATCACCTGCCTTCACTTCCGGCAAGGCGCGGGtt                                               <  2:3079153/93‑1 (MQ=255)
                                            tATTGTAGTTGGATGACATTGATGCGCCATATGCCCCTGTATCATGCGGTACCAGATAATCACCTGCCTTCACTTCCGGCAAGGCGCGGGttt                                              >  2:2056270/1‑93 (MQ=255)
                                            tATTGTAGTTGGATGACATTGATGCGCCATATGCCCCTGTATCATGCGGTACCAGATAATCACCTGCCTTCACTTCCGGCAAGGCGCGGGttt                                              >  1:182047/1‑93 (MQ=255)
                                            tATTGTAGTTGGATGACATTGATGCGCCATATGCCCCTGTATCATGCGGTACCAGATAATCACCTGCCTTCACTTCCGGCAAGGCGCGGGttt                                              >  2:2685057/1‑93 (MQ=255)
                                               tGTAGTTGGATGACATTGATGCGCCATATGCCCCTGTATCATGCGGTACCAGATAATCACCTGCCTTCACTTCCGGCAAGGCGCGGGTTTCaa                                           >  2:1057086/1‑93 (MQ=255)
                                                gTAGTTGGATGACATTGATGCGCCATATGCCCCTGTATCATGCGGTACCAGATAATCACCTGCCTTCACTTCCGGCAAGGCGCGGGTTTCAAc                                          <  1:1113665/93‑1 (MQ=255)
                                                   gTTGGATGACATTGATGCGCCATATGCCCCTGTATCATGCGGTACCAGATAATCACCTGCCTTCACTTCCGGCAAGGCGCGGGTTTCAACAtt                                       >  2:4127549/1‑93 (MQ=255)
                                                    ttGGATGACATTGATGCGCCATATGCCCCTGTATCATGCGGTACCAGATAATCACCTGCCTTCACTTCCGGCAAGGCGCGGGTTTCAACAttt                                      <  1:3896494/93‑1 (MQ=255)
                                                      ggATGACATTGATGCGCCATATGCCCCTGTATCATGCGGTACCAGATAATCACCTGCCTTCACTTCCGGCAAGGCGCGGGTTTCAACATTTcc                                    >  2:4452399/1‑93 (MQ=255)
                                                      ggATGACATTGATGCGCCATATGCCCCTGTATCATGCGGTACCAGATAATCACCTGCCTTCACTTCCGGCAAGGCGCGGGTTTCAACATTTcc                                    >  2:3412910/1‑93 (MQ=255)
                                                        aTGACATTGATGCGCCATATGCCCCTGTATCATGCGGTACCAGATAATCACCTGCCTTCACTTCCGGCAAGGCGCGGGTTTCAACATTTcccc                                  <  2:2501933/93‑1 (MQ=255)
                                                        aTGACATTGATGCGCCATATGCCCCTGTATCATGCGGTACCAGATAATCACCTGCCTTCACTTCCGGCAAGGCGCGGGTTTCAACATTTcccc                                  <  1:3713776/93‑1 (MQ=255)
                                                        aTGACATTGATGCGCCATATGCCCCTGTATCATGCGGTACCAGATAATCACCTGCCTTCACTTCCGGCAAGGCGCGGGTTTCAACATTTcccc                                  <  2:4615373/93‑1 (MQ=255)
                                                         tGACATTGATGCGCCATATGCCCCTGTATCATGCGGTACCAGATAATCACCTGCCTTCACTTCCGGCAAGGCGCGGGTTTCAACATTTccccc                                 <  2:47429/93‑1 (MQ=255)
                                                          gACATTGATGCGCCATATGCCCCTGTATCATGCGGTACCAGATAATCACCTGCCTTCACTTCCGGCAAGGCGCGGGTTTCAACATTTCCCCCt                                <  2:1191307/93‑1 (MQ=255)
                                                           aCATTGATGCGCCATATGCCCCTGTATCATGCGGTACCAGATAATCACCTGCCTTCACTTCCGGCAAGGCGCGGGTTTCAACATTTCCCCCtt                               >  1:2945358/1‑93 (MQ=255)
                                                             aTTGATGCGCCATATGCCCCTGTATCATGCGGTACCAGATAATCACCTGCCTTCACTTCCGGCAAGGCGCGGGTTTCAACATTTCCCCCTTcc                             <  2:936113/93‑1 (MQ=255)
                                                             aTTGATGCGCCATATGCCCCTGTATCATGCGGTACCAGATAATCACCTGCCTTCACTTCCGGCAAGGCGCGGGTTTCAACATTTCCCCCTTcc                             <  2:3087326/93‑1 (MQ=255)
                                                             aTTGATGCGCCATATGCCCCTGTATCATGCGGTACCAGATAATCACCTGCCTTCACTTCCGGCAAGGCGCGGGTTTCAACATTTCCCCCTTcc                             <  1:742410/93‑1 (MQ=255)
                                                             aTTGATGCGCCATATGCCCCTGTATCATGCGGTACCAGATAATCACCTGCCTTCACTTCCGGCAAGGCGCGGGGTTCAACATTTCCCCCTTcc                             >  1:511743/1‑93 (MQ=255)
                                                               tGATGCGCCATATGCCCCTGTATCATGCGGTACCAGATAATCACCTGCCTTCACTTCCGGCAAGGCGCGGGTTTCAACATTTCCCCCTTCctg                           >  2:4284849/1‑93 (MQ=255)
                                                                gATGCGCCATATGCCCCTGTATCATGCGGTACCAGATAATCACCTGCCTTCACTTCCGGCAAGGCGCGGGTTTCAACATTTCCCCCTTCctgc                          >  1:3097936/1‑93 (MQ=255)
                                                                 aTGCGCCATATGCCCCTGTATCATGCGGTACCAGATAATCACCTGCCTTCACTTCCGGCAAGGCGCGGGTTTCAACATTTCCCCCTTCctgct                         <  1:1206468/93‑1 (MQ=255)
                                                                  tGCGCCATATGCCCCTGTATCATGCGGTACCAGATAATCACCTGCCTTCACTTCCGGCAAGGCGCGGGTTTCAACATTTCCCCCTTCctgctg                        >  2:1957195/1‑93 (MQ=255)
                                                                  tGCGCCATATGCCCCTGTATCATGCGGTACCAGATAATCACCTGCCTTCACTTCCGGCAAGGCGCGGGTTTCAACATTTCCCCCTTCctgctg                        <  2:4413595/93‑1 (MQ=255)
                                                                    cgcCATATGCCCCTGTATCATGCGGTACCAGATAATCACCTGCCTTCACTTCCGGCAAGGCGCGGGTTTCAACATTTCCCCCTTCCTGCTggg                      <  1:1891357/93‑1 (MQ=255)
                                                                     cccATATGCCCCTGTATCATGCGGTACCAGATAATCACCTGCCTTCACTTCCGGCAAGGCGCGGGTTTCAACATTTCCCCCTTCCTGCTGGGt                     <  1:1655216/92‑1 (MQ=255)
                                                                        atatGCCCCTGTATCATGCGGTACCAGATAATCACCTGCCTTCACTTCCGGCAAGGCGCGGGTTTCAACATTTCCCCCTTCCTGCTGGGTaaa                  >  2:1205353/1‑93 (MQ=255)
                                                                        atatGCCCCTGTATCATGCGGTACCAGATAATCACCTGCCTTCACTTCCGGCAAGGCGCGGGCTTCAACATTTCCCCCTTCCTGCTGGGTaaa                  >  2:2789896/1‑93 (MQ=255)
                                                                         tatGCCCCTGTATCATGCGGTACCAGATAATCACCTGCCTTCACTTCCGGCAAGGCGCGGGTTTCAACATTTCTCCCTTCCTGCTGGGTAAAg                 >  2:1845149/1‑93 (MQ=255)
                                                                            gCCCCTGTATCATGCGGTACCAGATAATCACCTGCCTTCACTTCCGGCAAGGCGCGGGTTTCAACATTTCCCCCTTCCTGCTGGGTAAAGACa              >  2:717999/1‑93 (MQ=255)
                                                                            gCCCCTGTATCATGCGGTACCAGATAATCACCTGCCTTCACTTCCGGCAAGGCGCGGGTTTCAACATTTCCCCCTTCCTGCTGGGTAAAGACa              <  1:1909468/93‑1 (MQ=255)
                                                                              cccTGTATCATGCGGTACCAGATAATCACCTGCCTTCACTTCCGGCAAGGCGCGGGTTTCAACATTTCCCCCTTCCTGCTGGGTAAAGACATc            <  1:3132102/93‑1 (MQ=255)
                                                                              cccTGTATCATGCGGTACCAGATAATCACCTGCCTTCACTTCCGGCAAGGCGCGGGTTTCAACATTTCCCCCTTCCTGCTGGGTAAAGACATc            <  2:540050/93‑1 (MQ=255)
                                                                                cTGTATCATGCGGTACCAGATAATCACCTGCCTTCACTTCCGGCAAGGCGCGGGTTTCAACATTTCCCCCTTCCTGCTGGGTAAAGACATCGc          <  1:1031582/93‑1 (MQ=255)
                                                                                cTGTATCATGCGGTACCAGATAATCACCTGCCTTCACTTCCGGCAAGGCGCGGGTTTCAACATTTCCCCCTTCCTGCTGGGTAAAGACATCGc          <  1:2321322/93‑1 (MQ=255)
                                                                                cTGTATCATGCGGTACCAGATAATCACCTGCCTTCACTTCCGGCAAGGCGCGGGTTTCAACATTTCCCCCTTCCTGCTGGGTAAAGACATCGc          <  1:807261/93‑1 (MQ=255)
                                                                                   tATCATGCGGTACCAGATAATCACCTGCCTTCACTTCCGGCAAGGCGCGGGTTTCAACATTTCCCCCTTCCCGCTGGGTAAAGACATCGCCCg       <  2:3009479/93‑1 (MQ=255)
                                                                                      cATGCGGTACCAGATAATCACCTGCCTTCACTTCCGGCAAGGCGCGGGTTTCAACATTTCCCCCTTCCTGCTGGGTAAAGACATCGCCCGAtt    <  2:4481525/93‑1 (MQ=255)
                                                                                       aTGCGGTACCAGATAATCACCTGCCTTCACTTCCGGCAAGGCGCGGGTTTCAACATTTCCCCCTTCCTGCTGGGTAAAGACATCGCCCGATTc   >  2:4275504/1‑93 (MQ=255)
                                                                                        tGCGGTACCAGATAATCACCTGCCTTCACTTCCGGCAAGGCGCGGGTTTCAACATTTCCCCCTTCCTGCTGGGTAAAGACATCGCCCGATTca  <  1:2553494/93‑1 (MQ=255)
                                                                                        tGCGGTACCAGATAATCACCTGCCTTCACTTCCGGCAAGGCGCGGGTTTCAACATTTCCCCCTTCCTGCTGGGTAAAGACATCGCCCGATTca  <  2:3764151/93‑1 (MQ=255)
                                                                                           |                                                                                         
CGCCTGACCATTATCAAACAGAACTTCTGGTAACAGCGGACGGCTATTGTAGTTGGATGACATTGATGCGCCATATGCCCCTGTATCATGCAGTACCAGATAATCACCTGCCTTCACTTCCGGCAAGGCGCGGGTTTCAACATTTCCCCCTTCCTGCTGGGTAAAGACATCGCCCGATTCA  >  NC_000913/2975716‑2975896

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