Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,827,811 T→G 18.0% K150T (AAA→ACA)  astB ← succinylarginine dihydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,827,8110TG18.0% 314.8 / 64.1 137K150T (AAA→ACA) astBsuccinylarginine dihydrolase
Reads supporting (aligned to +/- strand):  ref base T (62/50);  new base G (16/9);  total (78/59)
Fisher's exact test for biased strand distribution p-value = 5.06e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.93e-01

Only 100 of 142 total aligned reads displayed.

CGAGACGATTGTGGTTTGCCGCCCCCTCATCACCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACTTTGCCATCCAGCGTATCGGCAG  >  NC_000913/1827713‑1827910
                                                                                                  |                                                                                                   
cgAGACGATTGTGGTTTGCCGCCCCCTCATCACCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCttt                                                                                                    >  2:3290476/1‑100 (MQ=255)
 gagaCGATTGTGGTTTGCCGCCCCCTCATCACCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCtttt                                                                                                   <  2:2697948/100‑1 (MQ=255)
  agaCGATTGTGGTTTGCCGCCCCCTCATCACCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTa                                                                                                  >  1:2798908/1‑100 (MQ=255)
  agaCGATTGTGGTTTGCCGCCCCCTCATCACCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTa                                                                                                  <  2:341539/100‑1 (MQ=255)
    aCGATTGTGGTTTGCCGCCCCCTCATCACCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAAc                                                                                                <  2:1071552/100‑1 (MQ=255)
     cGATTGTGGTTTGCCGCCCCCTCATCACCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACa                                                                                               >  2:1268235/1‑100 (MQ=255)
     cGATTGTGGTTTGCCGCCCCCTCATCACCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACa                                                                                               >  2:3620178/1‑100 (MQ=255)
     cGATTGTGGTTTGCCGCCCCCTCATCACCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACa                                                                                               <  1:3176256/100‑1 (MQ=255)
           tgGTTTGCCGCCCCCTCATCACCTAGCAACGCTACCTGTTTCAACGCCTAATGGACGCTAAAGTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGAt                                                                                          <  1:825728/99‑1 (MQ=255)
           tgGTTTGCCGCCCCCTCATCACCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTCAAAAATCGCTTTTAACAGCGAtt                                                                                         >  1:3965590/1‑100 (MQ=255)
            ggTTTGCCGCCCCCTCATCACCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTc                                                                                        <  2:1953730/100‑1 (MQ=255)
            ggTTTGCCGCCCCCTCATCACCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTc                                                                                        >  2:3962771/1‑100 (MQ=255)
             gTTTGCCGCCCCCTCATCACCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTGTTAAAAGCGAtttt                                                                                       <  2:3574411/100‑3 (MQ=255)
                tGCCGCCCCCTCATCACCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTGTTAAAAGCGAtttttaa                                                                                    <  2:3872743/100‑6 (MQ=255)
                 gccgccCCCTCATCACCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTGTTAAAAGCGAtttttaac                                                                                   <  1:2066496/100‑7 (MQ=39)
                 gccgccCCCTCATCACCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTGTTAAAAGCGAtttttaac                                                                                   <  1:1754472/100‑7 (MQ=39)
                  ccgccCCCTCATCACCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGAc                                                                                  >  1:2291622/1‑100 (MQ=255)
                  ccgccCCCTCATCACCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTGTTAAAAGCGAtttttaacg                                                                                  <  2:2245694/100‑8 (MQ=38)
                    gccCCCTCATCACCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACgg                                                                                >  2:951524/1‑100 (MQ=255)
                      ccccTCATCACCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGgc                                                                              >  1:1964391/1‑100 (MQ=255)
                        ccTCATCACCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGgcgc                                                                            >  2:2324496/1‑100 (MQ=255)
                        ccTCATCACCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGgcg                                                                             <  1:4065092/99‑1 (MQ=255)
                         cTCATCACCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCt                                                                           <  2:132995/100‑1 (MQ=255)
                         cTCATCACCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCt                                                                           <  2:960232/100‑1 (MQ=255)
                          tcatcaCCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCtt                                                                          >  2:1231445/1‑100 (MQ=255)
                           catcaCCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTc                                                                         >  2:3332598/1‑100 (MQ=255)
                            accaCCGAGCAACACTACCTGTGGCAACGCCGAATGGACGCTGAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTcc                                                                        >  2:2969983/3‑100 (MQ=255)
                             tcaCCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCa                                                                       >  2:177896/1‑100 (MQ=255)
                             tcaCCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCa                                                                       >  2:2817507/1‑100 (MQ=255)
                              caCCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAg                                                                      >  2:1615238/1‑100 (MQ=255)
                              caCCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAg                                                                      >  1:1257541/1‑100 (MQ=255)
                              caCCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAg                                                                      >  2:4113686/1‑100 (MQ=255)
                               aCCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGc                                                                     <  1:3543779/100‑1 (MQ=255)
                               aCCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGc                                                                     >  2:3333158/1‑100 (MQ=255)
                                      aaCGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGt                                                              >  1:766787/1‑100 (MQ=255)
                                       aCGCTACCTGTGGCAACGCCGGATGGACGCTAAATTTCTCTTCGGCGTTAAAAATCGCTTTTAACAGCGCTTCAGTGACGGGCGCTTCCGGCGAACgg                                                               >  1:3118939/1‑98 (MQ=255)
                                       aCGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTg                                                             <  1:2528572/100‑1 (MQ=255)
                                        cGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGa                                                            <  1:3872765/100‑1 (MQ=255)
                                         gCTACCTGTCGCAACGCCGAATGGACGCTAAATTTCTCGTCGTCGTTAAAAATCGCTTTTAACAGCGAGTCAGTGACGGGCGCTTCCAGCGAACGGTGaa                                                           <  2:4255217/100‑1 (MQ=255)
                                               tgtgGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGGAATTTAt                                                     <  1:3197614/100‑1 (MQ=255)
                                               tgtgGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTAt                                                     >  2:2158471/1‑100 (MQ=255)
                                                gtgGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTAtt                                                    >  1:1162923/1‑100 (MQ=255)
                                                 tgGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTg                                                   >  1:709448/1‑100 (MQ=255)
                                                 nggCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTGCCAGCGAACGGTGAAATTTATTg                                                   >  2:1371811/2‑100 (MQ=255)
                                                  ggCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGCTTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGt                                                  >  1:1377019/1‑100 (MQ=255)
                                                  ggCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGt                                                  >  2:4255523/1‑100 (MQ=255)
                                                   gCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGtt                                                 >  2:368334/1‑100 (MQ=255)
                                                       cGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTGTAACAGCGATTGAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCa                                               >  1:2820645/1‑98 (MQ=255)
                                                        gCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGt                                            <  1:4113657/100‑1 (MQ=255)
                                                        gCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGt                                            >  1:1397776/1‑100 (MQ=255)
                                                        gCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGt                                            <  2:2360985/100‑1 (MQ=255)
                                                          cGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTg                                          >  1:1816979/1‑100 (MQ=255)
                                                          cGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTg                                          <  1:3364032/100‑1 (MQ=255)
                                                           gAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTgg                                         <  2:599961/100‑1 (MQ=255)
                                                           gAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTgg                                         <  2:3327732/100‑1 (MQ=255)
                                                           gAATGGACGCTAAATTTCTCTTCGTAGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTgg                                         <  1:1313036/100‑1 (MQ=255)
                                                            aaTGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGc                                        <  1:2728814/100‑1 (MQ=255)
                                                            aaTGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTCCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGc                                        <  2:1130283/100‑1 (MQ=255)
                                                                 aCGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACgg                                   >  1:2680870/1‑100 (MQ=255)
                                                                 aCGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACgg                                   >  1:916214/1‑100 (MQ=255)
                                                                 aCGCTAAATTTCTCTTCGTCGTCAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACgg                                   <  2:2147530/100‑1 (MQ=255)
                                                                   gCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTg                                 >  2:803308/1‑100 (MQ=255)
                                                                    cTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTga                                >  1:2986096/1‑100 (MQ=255)
                                                                    cTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTga                                >  2:665190/1‑100 (MQ=255)
                                                                    cTAAATTTCTCTTCGTCGTTAAAAATCGCTGTTGACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTAGCAACGGTga                                >  1:2891518/1‑100 (MQ=255)
                                                                      aaaTTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTgaga                              >  1:2480186/1‑100 (MQ=255)
                                                                         ttgctctTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATGTATTGTTCAGGTTGGCAACGGTGAGATGc                           <  1:1639560/97‑1 (MQ=255)
                                                                          ttCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCa                          <  2:545509/100‑1 (MQ=255)
                                                                           tctcTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCAc                         >  2:3405481/1‑100 (MQ=255)
                                                                             tctTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGTTTGGCAACGGTGAGATGCACtt                       <  1:2557282/100‑1 (MQ=255)
                                                                                tcgtcgTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACTTTGc                    <  2:50816/100‑1 (MQ=255)
                                                                                tcgtcgTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACTTTGc                    <  2:2235557/100‑1 (MQ=255)
                                                                                  gtcgtTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACTTTGCCa                  <  1:2540755/100‑1 (MQ=255)
                                                                                  gtcgtTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACTTTGCCa                  <  1:1615233/100‑1 (MQ=255)
                                                                                   tcgtTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACTTTGCCAt                 <  1:2491707/100‑1 (MQ=255)
                                                                                   tcgtTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACTTTGCCAt                 >  1:3333439/1‑100 (MQ=255)
                                                                                    gtcactgaaTCGCTGTTAAAAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACTTTGCCATc                >  2:1330240/8‑100 (MQ=38)
                                                                                    cgtTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACTTTGCCATc                >  1:3490796/1‑100 (MQ=255)
                                                                                     gtTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACTTTGCCATcc               >  2:3955956/1‑100 (MQ=255)
                                                                                      ttAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACTTTGCCATCCa              <  1:3816279/100‑1 (MQ=255)
                                                                                      ttAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACTTTGCCATCCa              >  1:551406/1‑100 (MQ=255)
                                                                                       actgaaTCGCTGTTAAAAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACTTTGCCATCCAg             >  2:698635/5‑100 (MQ=255)
                                                                                        ctgaaTCGCTGTTAAAAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACTTTGCCATCCAg             >  2:1792295/4‑99 (MQ=255)
                                                                                        ctgaaTCGCTGTTAAAAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACTTTGCCATCCAGc            >  1:3330670/4‑100 (MQ=255)
                                                                                        ctgaaTCGCTGTTAAAAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACTTTGCCATCCAGc            >  2:3119478/4‑100 (MQ=255)
                                                                                        ctgaaTCGCTGTTAAAAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACTTTGCCATCCAGc            >  2:930539/4‑100 (MQ=255)
                                                                                        ctgaaTCGCTGTTAAAAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACTTTGCCATCCAGc            >  1:2292076/4‑100 (MQ=255)
                                                                                        ctgaaTCGCTGTTAAAAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACTTTGCCATCCAGc            >  2:2591471/4‑100 (MQ=255)
                                                                                        ctgaaTCGCGGTTAAAAGCGATTCAGTGACGGGCGCTTCCAGCGAACAGCGAAATTTACTGTTCAGGTTGTCAACGGTGAGATGCTCTTTGCCATCCa              >  1:481302/4‑98 (MQ=255)
                                                                                        aaaaaTCGCTGTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACTTTGCCATCCAg             <  2:1638037/99‑1 (MQ=255)
                                                                                         aaaaTCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCATGGAACGGTGAAATTTATTGTGCAGGTTGGCAACGGTGAGATGCACTTTGCCATCCAGCg           >  1:2394416/1‑100 (MQ=255)
                                                                                         aaaaTCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACTTTGCCATCCAGCg           <  1:311601/100‑1 (MQ=255)
                                                                                         aaaaTCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACTTTGCCATCCAGCg           >  1:709332/1‑100 (MQ=255)
                                                                                              cGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACTTTGCCATCCAGCGTATCg      <  1:2990175/100‑1 (MQ=255)
                                                                                               gCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACTTTGCCATCCAGCGTATCgg     <  2:1909206/100‑1 (MQ=255)
                                                                                                cTTTTAACAGCGCTTCAGTGACGGGCGCTTCCTGCGGACGGTGAAAGTTATTGTTCAGGTTGGCAACGGTGAGATGCACTTTGCCATCCAGCGTATCGGc    >  1:3695867/1‑100 (MQ=255)
                                                                                                cTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACTTTGCCATCCAGCGTATCGGc    <  2:209173/100‑1 (MQ=255)
                                                                                                 ttttAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACTTTGCCATCCAGCGTATCGGCa   >  2:2316504/1‑100 (MQ=255)
                                                                                                 ttttAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACTTTGCCATCCAGCGTATCGGCa   >  2:365856/1‑100 (MQ=255)
                                                                                                  tttAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACTTTGCCATCCAGCGTATCGGCAg  <  1:1203569/100‑1 (MQ=255)
                                                                                                  |                                                                                                   
CGAGACGATTGTGGTTTGCCGCCCCCTCATCACCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACTTTGCCATCCAGCGTATCGGCAG  >  NC_000913/1827713‑1827910

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: