Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,830,747 T→C L99P (CTG→CCG)  xthA → exonuclease III

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,830,7470TC98.9% 579.7 / NA 188L99P (CTG→CCG) xthAexonuclease III
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (107/79);  minor base A (0/2);  total (107/81)
Fisher's exact test for biased strand distribution p-value = 1.84e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 210 total aligned reads displayed.

CCAAAGAGACGCCGATTGCCGTGCGTCGCGGCTTTCCCGGTGACGACGAAGAGGCGCAGCGGCGGATTATTATGGCGGAAATCCCCTCACTGCTGGGTAATGTCACCGTGATCAACGGTTACTTCCCGCAGGGTGAAAGCCGCGACCATCCGATAAAATTCCCGGCAAAAGCGCAGTTTT  >  NC_000913/1830657‑1830836
                                                                                          |                                                                                         
ccAAAGAGACGCCGATTGCCGTGCGTCGCGGCTTTCCCGGTGACGACGAAGAGGCGCAGCGGCGGATTATTATGGCGGAAATCCCCTCAccgc                                                                                         >  2:3042012/1‑93 (MQ=255)
ccAAAGAGACGCCGATTGCCGTGCGTCGCGGCTTTCCCGGTGACGACGAAGAGGCGCAGCGGCGGATTATTATGGCGGAAATCCCCTCAccgc                                                                                         <  1:928045/93‑1 (MQ=255)
ccAAAGAGACGCCGATTGCCGTGCGTCGCGGATTTCCCGGTGACGACGAAGAGGCGCAGCGGCGGATTATTATGGCGGAAAGCCCCCCAccgc                                                                                         >  1:2020344/1‑93 (MQ=255)
    agagACGCCGATTGCCGTGCGTCGCGGCTTTCCCGGTGACGACGAAGAGGCGCAGCGGCGGATTATTATGGCGGAAATCCCCTCACCGCTggg                                                                                     >  2:4957299/1‑93 (MQ=255)
     gagaCGCCGATTGCCGTGCGTCGCGGCTTTCCCGGTGACGACGAAGAGGCGCAGCGGCGGATTATTATGGCGGAAATCCCCTCACCGCTGGGt                                                                                    >  2:7236365/1‑93 (MQ=255)
     gagaCGCCGATTGCCGTGCGTCGCGGCTTTCCCGGTGACGACGAAGAGGCGCAGCGGCGGATTATTATGGCGGAAATCCCCTCACCGCTGGGt                                                                                    >  1:7213850/1‑93 (MQ=255)
     gagaCGCCGATTGCCGTGCGTCGCGGCTTTCCCGGTGACGACGAAGAGGCGCAGCGGCGGATTATTATGGCGGAAATCCCCTCACCCCTGGGt                                                                                    >  1:4585848/1‑93 (MQ=255)
     gagaCGCCGATTGCCGTGCGTCGCGGCTTTCCCGGTGACGACGAAGAGGCGCAGCGGCGGATTATTATGGCGGAAATACCCTCACCGCTgggg                                                                                    >  2:3723722/1‑92 (MQ=255)
       gaCGCCGATTGCCGTGCGTCGCGGCTTTCCCGGTGACGACGAAGAGGCGCAGCGGCGGATTATTATGGCGGAAATCCCCTCACCGCTGGGTaa                                                                                  <  1:7723664/93‑1 (MQ=255)
       gaCGCCGATTGCCGTGCGTCGCGGCTTTCCCGGTGACGACGAAGAGGCGCAGCGGCGGATTATTATGGCGGAAATCCCCTCACCGCTGGGTaa                                                                                  >  1:3703836/1‑93 (MQ=255)
           ccGATTGCCGTGCGTCGCGGCTTTCCCGGTGACGACGAAGAGGCGCAGCGGCGGATTATTATGGCGGAAATCCCCTCACCGCTGGGTAATGta                                                                              >  1:3918697/1‑92 (MQ=255)
            cGATTGCCGTGCGTCGCGGCTTTCCCGGTGACGACGAAGAGGCGCAGCGGCGGATTATTATGGCGGAAATCCCCTCACCGCTGGGTAATGTCa                                                                             >  1:1068251/1‑93 (MQ=255)
            cGATTGCCGTGCGTCGCGGCTTTCCCGGTGACGACGAAGAGGCGCAGCGGCGGATTATTATGGCGGAAATCCCCTCACCGCTGGGGAATGTCa                                                                             >  2:6074573/1‑93 (MQ=255)
              aTTGCCGTGCGTCGCGGCTTTCCCGGTGACGACGAAGAGGCGCAGCGGCGGATTATTATGGCGGAAATCCCCTCACCGCTGGGTAATGTCAcc                                                                           >  1:6273213/1‑93 (MQ=255)
              aTTGCCGTGCGTCGCGGCTTTCCCGGTGACGACGAAGAGGCGCAGCGGCGGATTATTATGGCGGAAATCCCCTCACCGCTGGGTAATGTCAcc                                                                           >  2:4255339/1‑93 (MQ=255)
              aTTGCCGTGCGTCGCGGCTTTCCCGGTGACGACGAAGAGGCGCAGCGGCGGATTATTATGGCGGAAATCCCCTCACCGCTGGGTAATGTCAcc                                                                           >  2:4163107/1‑93 (MQ=255)
               ttGCCGTGCGTCGCGGCTTTCCCGGTGACGACGAAGAGGCGCAGCGGCGGATTATTATGGCGGAAATCCACTCAACGCTGGGGAATGTCGCCg                                                                          >  1:2707818/1‑93 (MQ=255)
                tGCCGTGCGTCGCGGCTTTCCCGGTGACGACGAAGAGGCGCAGCGGCGGATTATTATGGCGGAAATCCCCTCACCGCTGGGTAATGTCACCGt                                                                         <  1:7899013/93‑1 (MQ=255)
                 gCCGTGGGTCGCGGCTTTCCCGGTGACGACGAAGAGGCGCAGCGGCGGATTATTATGGCGGAAATCCCCTCACCGCTGGGTAATGTCACCGTg                                                                        >  2:7678598/1‑93 (MQ=255)
                   cGTGCGTCGCGGCTTTCCCGGTGACGACGAAGAGGCGCAGCGGCGGATTATTATGGCGGAAATCCCCTCACCGCTGGGTAATGTCACCGTGAt                                                                      <  1:2814653/93‑1 (MQ=255)
                   cGTGCGTCGCGGCTTTCCCGGTGACGACGAAGAGGCGCAGCGGCGGATTATTATGGCGGAAATCCCCTCACCGCTGGGTAATGTCACCGTGAt                                                                      >  2:176221/1‑93 (MQ=255)
                     tGCGTCGCGGCTTTCCCGGTGACGACGAAGAGGCGCAGCGGCGGATTATTATGGCGGAAATCCCCTCACCGCTGGGTAATGTCACCGTGATCa                                                                    >  1:2293913/1‑93 (MQ=255)
                     tGCGTCGCGGCTTTCCCGGTGACGACGAAGAGGCGCAGCGGCGGATTATTATGGCGGAAATCCCCTCACCGCTGGGTAATGTCACCGTGATCa                                                                    >  2:122485/1‑93 (MQ=255)
                      gCGTCGCGGCTTTCCCGGTGACGACGAAGAGGCGCAGCGGCGGATTATTATGGCGGAAATCCCCTCACCGCTGGGTAATGTCACCGTGATCaa                                                                   <  1:3367284/93‑1 (MQ=255)
                       cGTCGCGGCTTTCCCGGTGACGACGAAGAGGCGCAGCGGCGGATTATTATGGCGGAAATCCCCTCACCGCTGGGTAATGTCACCGTGATCaat                                                                  <  2:4793513/93‑2 (MQ=255)
                       cGTCGCGGCTTTCCCGGTGACGACGAAGAGGCGCAGCGGCGGATTATTATGGCGGAAATCCCCTCACCGCTGGGTAATGTCACCGTGATCAAc                                                                  <  2:2728890/93‑1 (MQ=255)
                       cGTCGCGGCTTTCCCGGTGACGACGAAGAGGCGCAGCGGCGGATTATTATGGCGGAAATCCCCTCACCGCTGGGTAATGTCACCGTGATCAAc                                                                  <  2:5465907/93‑1 (MQ=255)
                        gTCGCGGCTTTCCCGGTGACGACGAAGAGGCGCAGCGGCGGATTATTATGGCGGAAATCCCCTCACCGCTGGGTAATGTCACCGTGATCAACg                                                                 >  1:1196852/1‑93 (MQ=255)
                         tCGCGGCTTTCCCGGTGACGACGAAGAGGCGCAGCGGCGGATTATTATGGCGGAAATCCCCTCACCGCTGGGTAATGTCACCGTGATCAACgg                                                                >  2:7796292/1‑93 (MQ=255)
                          cgcgGCTTTCCCGGTGACGACGAAGAGGCGCAGCGGCGGATTATTATGGCGGAAATCCCCTCACCGCTGGGTAATGTCACCGTGATCAACGGt                                                               <  2:1889724/93‑1 (MQ=255)
                           gcgGCTTCCCCGGTGACGACGAAGAGGCGCCGCGGCGGATTATTATGGCGGAAATCCCCTCGCCGCTGGGTAATGTCACCGTGATCAACGGtt                                                              <  2:7716497/93‑1 (MQ=255)
                            cgGCTTTCCCGGTGACGACGAAGAGGCGCAGCGGCGGATTATTATGGCGGAAATCCCCTCACCGCTGGGTAATGTCACCGTGATCAACGGTTa                                                             >  1:778877/1‑93 (MQ=255)
                            cgGCTTTCCCGGTGACGACGAAGAGGCGCAGCGGCGGATTATTATGGCGGAAATCCCCTCACCGCTGGGTAATGTCACCGTGATCAACGGTTa                                                             >  1:3879085/1‑93 (MQ=255)
                            cgGCTTTCCCGGTGACGACGAAGAGGCGCAGCGGCGGATTATTATGGCGGAAATCCCCTCACCGCTGGGTAATGTCACCGTGATCAACGGTTa                                                             <  2:6660853/93‑1 (MQ=255)
                             ggCTTTCCCGGTGACGACGAAGAGGCGCAGCGGCGGATTATTATGGCGGAAATCCCCTCACCGCTGGGTAATGTCACCGTGATCAACGGTTAc                                                            >  2:4828393/1‑93 (MQ=255)
                              gCTTTCCCGGTGACGACGAAGAGGCGCAGCGGCGGATTATTATGGCGGAAATCCCCTCACCGCTGGGTAATGTCACCGTGATCAACTGCTACt                                                           >  2:2448708/1‑93 (MQ=255)
                              gCTTTCCCGGTGACGACGAAGAGGCGCAGCGGCGGATTATTATGGCGGAAATCCCCTCACCGCTGGGTAATGTCACCGTGATCAACGGTTACt                                                           >  1:1166318/1‑93 (MQ=255)
                                 ttCCCGGTGACGACGAAGAGGCGCAGCGGCGGATTATTATGGCGGAAATCCCCTCACCGCTGGGTAATGTCACCGTGATCAACGGTTACTTcc                                                        <  2:1610085/93‑1 (MQ=255)
                                  tCCCGGTGACGACGAAGAGGCGCAGCGGCGGATTATTATGGCGGAAATCCCCTCACCGCTGGGTAATGTCACCGTGATCAACGGTTACTTccc                                                       <  1:4496510/93‑1 (MQ=255)
                                    ccGGTGACGACGAAGAGGCGCAGCGGCGGATTATTATGGCGGAAATCCCCTCACCGCTGGGGAATGTCACCGTGATCAACGGTTACTTCCCGc                                                     >  2:3517491/1‑93 (MQ=255)
                                     cGGTGACGACGAAGAGGCGCAGCGGCGGATTATTATGGCGGAAATCCCCTCACCGCTGGGTAATGTCACCGTGATCAACGGTTACTTCCCGCa                                                    <  2:7619900/93‑1 (MQ=255)
                                      ggTGACGACGAAGAGGCGCAGCGGCGGATTATTATGGCGGAAATCCCCTCACCGCTGGGTAATGTCACCGTGATCAACGGTTACTTCCCGCAg                                                   >  2:6156963/1‑93 (MQ=255)
                                      ggTGACGACGAAGAGGCGCAGCGGCGGATTATTATGGCGGAAATCCCCTCACCGCTGGGTAATGTCACCGTGATCAACGGTTACTTCCCGCAg                                                   <  2:1273025/93‑1 (MQ=255)
                                      ggTGACGACGAAGAGGCGCAGCGGCGGATTATTATGGCGGAAATCCCCTCACCGCTGGGTAATGTCACCGTGATCAACGGTTACTTCCCGCAg                                                   <  1:6285815/93‑1 (MQ=255)
                                      ggTGACGACGAAGAGGCGCAGCGGCGGATTATTATGGCGGAAATCCCCTCACCGCTGGGTAATGTCACCGTGATCAACGGTTACTTCCCGCAg                                                   <  2:574833/93‑1 (MQ=255)
                                       gTGACGACGAAGAGGCGCAGCGGCGGATTATTATGGCGGAAATCCCCTCACCGCTGGGTAATGTCACCGTGATCAACGGTTACTTCCCGCAgg                                                  >  2:4932592/1‑93 (MQ=255)
                                        tGACGACGAAGAGGCGCAGCGGCGGATTATTATGGCGGAAATCCCCTCACCGCTGGGTAATGTCACCGTGATCAACGGTTACTTCCCGCAggg                                                 >  1:4859908/1‑93 (MQ=255)
                                            gacgaAGAGGCGCAGCGGCGGATTATTATGGCGGAAATCCCCTCACCGCTGGGTAATGTCACCGTGATCAACGGTTACTTCCCGCAGGGTGaa                                             <  1:5904734/93‑1 (MQ=255)
                                              cgaAGAGGCGCAGCGGCGGATTATTATGGCGGAAATCCCCTCACCGCTGGGTAATGTCACCGTGATCAACGGTTACTTCCCGCAGGGTGAAAg                                           >  2:406871/1‑93 (MQ=255)
                                                aaGAGGCGCAGCGGCGGATTATTATGGCGGAAATCCCCTCACCGCTGGGTAATGTCACCGTGATCAACGGTTACTTCCCGCAGGGTGAAAGcc                                         >  1:3784504/1‑93 (MQ=255)
                                                aaGAGGCGCAGCGGCGGATTATTATGGCGGAAATCCCCTCACCGCTGGGTAATGTCACCGTGATCAACGGTTACTTCCCGCAGGGTGAAAGcc                                         >  2:4494761/1‑93 (MQ=255)
                                                aaGAGGCGCAGCGGCGGATTATTATGGCGGAAATCCCCTCACCGCTGGGTAATGTCACCGTGATCAACGGTTACTTCCCGCAGGGTGAAAGcc                                         >  2:4418210/1‑93 (MQ=255)
                                                 agagGCGCAGCGGCGGATTATTATGGCGGAAATCCCCTCACCGCTGGGTAATGTCACCGTGATCAACGGTTACTTCCCGCAGGGTGAAAGCcg                                        <  2:1680460/93‑1 (MQ=255)
                                                 agagGCGCAGCGGCGGATTATTATGGCGGAAATCCCCTCACCGCTGGGTAATGTCACCGTGATCAACGGTTACTTCCCGCAGGGTGAAAGCcg                                        >  1:6560187/1‑93 (MQ=255)
                                                         aGCGGCGGATTATTATGGCGGAAATCCCCTCACCGCTGGGTAATGTCACCGTGATCAACGGTTACTTCCCGCAGGGTGAAAGCCGCGACCATc                                <  2:3728983/93‑1 (MQ=255)
                                                          gcggcgGATTATTATGGCGGAAATCCCCTCACCGCTGGGTAATGTCACCGTGATCAACGGTTACTTCCCGCAGGGTGAAAGCCGCGACCATcc                               <  2:4408379/93‑1 (MQ=255)
                                                          gcggcgGATTATTATGGCGGAAATCCCCTCACCGCTGGGTAATGTCACCGTGATCAACGGTTACTTCCCGCAGGGTGAAAGCCGCGACCATcc                               >  2:1401599/1‑93 (MQ=255)
                                                          gcggcgGATTATTATGGCGGAAATCCCCTCACCGCTGGGTAATGTCACCGTGATCAACGGTTACTTCCCGCAGGGTGAAAGCCGCGACCATcc                               <  1:3411454/93‑1 (MQ=255)
                                                            ggcggATTATTATGGCGGAAATCCCCTCACCGCTGGGTAATGTCACCGTGATCAACGGTTACTTCCCGCAGGGTGAAAGCCGCGACCATCCGa                             <  2:6526346/93‑1 (MQ=255)
                                                            ggcggATTATTATGGCGGAAATCCCCTCACCGCTGGGTAATGTCACCGTGATCAACGGTTACTTCCCGCAGGGTGAAAGCCGCGACCATCCGa                             <  2:3077776/93‑1 (MQ=255)
                                                            ggcggATTATTATGGCGGAAATCCCCTCACCGCTGGGTAATGTCACCGTGATCAACGGTTACTTCCCGCAGGGTGAAAGCCGCGACCATCCGa                             <  2:6777841/93‑1 (MQ=255)
                                                            ggcggATTATTATGGCGGAAATCCCCTCACCGCTGGGTAATGTCACCGTGATCAACGGTTACTTCCCGCAGGGTGAAAGCCGCGACCATCCGa                             <  1:4219423/93‑1 (MQ=255)
                                                              cggATTATTATGGCGGAAATCCCCTCACCGCTGGGTAATGTCACCGTGATCAACGGTTACTTCCCGCAGGGTGAAAGCCGCGACCATCCGATa                           <  2:2606644/93‑1 (MQ=255)
                                                               ggATTATTATGGCGGAAATCCCCTCACCGCTGGGTAATGTCACCGTGATCAACGGTTACTTCCCGCAGGGTGAAAGCCGCGACCATCCGATaa                          >  2:1147995/1‑93 (MQ=255)
                                                                gATTATTATGGCGGAAATCCCCTCACCGCTGGGTAATGTCACCGTGATCAACGGTTACTTCCCGCAGGGTGAAAGCCGCGACCATCCGATaaa                         >  2:2767016/1‑93 (MQ=255)
                                                                gATTATTATGGCGGAAATCCCCTCACCGCTGGGTAATGTCACCGTGATCAACGGTTACTTCCCGCAGGGTGAAAGCCGCGACCATCCGATaaa                         >  2:5262191/1‑93 (MQ=255)
                                                                gATTATTATGGCGGAAATCCCCTCACCGCTGGGTAATGTCACCGTGATCAACGGTTACTTCCCGCAGGGTGAAAGCCGCGACCATCCGATaaa                         >  1:7569248/1‑93 (MQ=255)
                                                                gATTATTATGGCGGAAATCCCCTCACCGCTGGGTAATGTCACCGTGATCAACGGTTACTTCCCGCAGGGTGAAAGCCGCGACCATCCGATaaa                         >  1:5891389/1‑93 (MQ=255)
                                                                 attattATGGCGGAAATCCCCTCACCGCTGGGTAATGTCACCGTGATCAACGGTTACTTCCCGCAGGGTGAAAGCCGCGACCATCCGATaaaa                        <  1:2014029/93‑1 (MQ=255)
                                                                   tattatGGCGGAAATCCCCTCACCGCTGGGTAATGTCACCGTGATCAACGGTTACTTCCCGCAGGGTGAAAGCCGCGACCATCCGATAAAAtt                      <  1:1679011/93‑1 (MQ=255)
                                                                    attatGGCGGAAATCCCCTCACCGCTGGGTAATGTCACCGTGATCAACGGTTACTTCCCGCAGGGTGAAAGCCGCGACCATCCGATAAAATTc                     >  1:5241822/1‑93 (MQ=255)
                                                                     ttatGGCGGAAATCCCCTCACCGCTGGGTAATGTCACCGTGATCAACGGTTACTTCCCGCAGGGTGAAAGCCGCGACCATCCGATAAAATTcc                    >  1:1530559/1‑93 (MQ=255)
                                                                       atGGCGGAAATCCCCTCCCCGCTGGTTAATGTCACCGTGATCAACGGTTACTTCCCGCAGGGTGAAAGCCGCGACCATCCGATAAAATTCCCg                  <  1:1604064/93‑1 (MQ=255)
                                                                       atGGCGGAAATCCCCTCACCGCTGGGTAATGTCACCGTGATCAACGGTTACTTCCCGCAGGGGGAAAGCCGCGACCATCCGATAAAATTCCCg                  >  1:2517707/1‑93 (MQ=255)
                                                                       atGGCGGAAATCCCCTCACAGCTGGGTAATGTCACCGTGATCAACGGTTACTTCCCGCAGGGTGAAAGCCGCGCCCATCCGATAAAATTCCCg                  <  1:4050459/93‑1 (MQ=255)
                                                                         ggCGGAAATCCCCTCACCGCTGGGTAATGTCACCGTGATCAACGGTTACTTCCCGCAGGGTGAAGGCCGCGACCATCCGATAAAATTCCCGGc                <  1:1382450/93‑1 (MQ=255)
                                                                         ggCGGAAATCCCCTCACCGCTGGGTAATGTCACCGTGATCAACGGTTACTTCCCGCAGGGTGAAAGCCGCGACCATCCGATAAAATTCCCGGc                <  1:5386723/93‑1 (MQ=255)
                                                                          gCGGAAATCCCCTCACCGCTGGGTAATGTCACCGTGATCAACGGTTACTTCCCGCAGGGTGAAAGCCGCGACCATCCGATAAAATTCCCGGCa               <  2:1289985/93‑1 (MQ=255)
                                                                          gCGGAAATCCCCTCACAGCTGGGTAATGTAACAGTGATCAACGGTTACTTCCCGCAGGGTGAAAGCCGCGACCATCCGATAAAATTCCCGGCa               <  1:7036312/93‑1 (MQ=255)
                                                                            ggAAATCCCCTCACCGCTGGGTAATGTCACCGTGATCAACGGTTACTTCCCGCAGGGTGAAAGCCGCGACCATCCGATAAAATTCCCGGCaaa             <  2:3435343/93‑1 (MQ=255)
                                                                             gAAATCCCCTCACCGCTGGGTAATGTCACCGTGATCAACGGTTACTTCCCGCAGGGTGAAAGCCGCGACCATCCGATAAAATTCCCGGCaaaa            <  1:5732655/93‑1 (MQ=255)
                                                                             gAAATCCCCTCACCGCTGGGTAATGTCACCGTGATCAACGGTTACTTCCCGCAGGGTGAAAGCCGCGACCATCCGATAAAATTCCCGGCaaaa            <  1:7261383/93‑1 (MQ=255)
                                                                             gAAATCCCCTCACCGCTGGGTAATGTCACCGTGATCAACGGTTACTTCCCGCAGGGTGAAAGCCGCGACCATCCGATAAAATTCCCGGCaaaa            <  1:4439788/93‑1 (MQ=255)
                                                                              aaaTCCCCTCACCGCTGGGTAATGTCACCGTGATCAACGGTTACTTCCCGCAGGGTGAAAGCCGCGACCATCCGATAAAATTCCCGGCaaaac           >  1:653395/1‑92 (MQ=255)
                                                                               aaTCCCCTCACCGCTGGGTAATGTCACCGTGATCAACGGTTACTTCCCGCAGGGTGAAAGCCGCGACCATCCGATAAAATTCCCGGCAAAAgc          <  1:6890686/93‑1 (MQ=255)
                                                                               aaTCCCCTCACCGCTGGGTAATGTCACCGTGATCAACGGTTACTTCCCGCAGGGTGAAAGCCGCGACCATCCGATAAAATTCCCGGCAAAAgc          <  2:3771832/93‑1 (MQ=255)
                                                                                aTCCCCTCACCGCTGGGTAATGTCACCGTGATCAACGGTTACTTCCCGCAGGGTGAAAGCCGCGACCATCCGATAAAATTCCCGGCAAAAgcg         >  2:71127/1‑93 (MQ=255)
                                                                                 tCCCCTCACCGCTGGGTAATGTCTCCGTGATCAACGGTTACTTCCCGCAGTGTGAAAGCCGCGACCATCCGATAAAATTCCCGGCAAAAgcgc        <  1:3495374/93‑1 (MQ=255)
                                                                                 tCCCCTCACCGCTGGGTAATGTCACCGTGATCAACGGTTACTTCCCGCAGGGTGAAAGCCGCGACCATCCGATAAAATTCCCGGCAAAAgcgc        >  1:4094699/1‑93 (MQ=255)
                                                                                 tCCCCTCACCGCTGGGTAATGTCACCGTGATCAACGGTTACTTCCCGCAGGGTGAAAGCCGCGACCATCCGATAAAATTCCCGGCAAAAgcgc        >  1:2228590/1‑93 (MQ=255)
                                                                                  ccccTCACCGCTGGGTAATGTCACCGTGATCAACGGTTACTTCCCGCAGGGTGAAAGCCGCGACCATCCGATAAAATTCCCGGCAAAAGCGCa       <  2:2994098/93‑1 (MQ=255)
                                                                                  ccccTCACCGCTGGGTAATGTCACCGTGATCAACGGTTACTTCCCGCAGGGTGAAAGCCGCGACCATCCGATAAAATTCCCGGCAAAAGCGCa       <  1:2938617/93‑1 (MQ=255)
                                                                                   cccTCACCGCTGGGTAATGTCACCGTGATCAACGGTTACTTCCCGCAGGGTGAAAGCCGCGACCATCCGATAAAATTCCCGGCAAAAGCGCAg      >  1:5600512/1‑93 (MQ=255)
                                                                                    ccTCACCGCTGGGTAATGTCACCGTGATCAACGGTTACTTCCCGCAGGGTGAAAGCCGCGACCATCCGATAAAATTCCCGGCAAAAGCGCAGt     >  2:4137928/1‑93 (MQ=255)
                                                                                    ccTCACCGCTGGGTAATGTCACCGTGATCAACGGTTACTTCCCGCAGGGTGAAAGCCGCGACCATCCGATAAAATTCCCGGCAAAAGCGCAGt     >  2:2449481/1‑93 (MQ=255)
                                                                                     cTCACCGCTGGGTAATGTCACCGTGATCAACGGTTACTTCCCGCAGGGTGAAAGCCGCGACCATCCGATAAAATTCCCGGCAAAAGCGCAGtt    <  1:1934273/93‑1 (MQ=255)
                                                                                     cTCACCGCTGGGTAATGTCACCGTGATCAACGGTTACTTCCCGCAGGGTGAAAGCCGCGACCATCCGATAAAATTCCCGGCAAAAGCGCAGtt    <  2:2175040/93‑1 (MQ=255)
                                                                                      tCACCGCTGGGTAATGTCACCGTGATCAACGGTTACTTCCCGCAGGGTGAAAGCCGCGACCATCCGATAAAATTCCCGGCAAAAGCGCAGttt   >  1:7219490/1‑93 (MQ=255)
                                                                                      tCACCGCTGGGTAATGTCACCGTGATCAACGGTTACTTCCCGCAGGGTGAAAGCCGCGACCATCCGATAAAATTCCCGGCAAAAGCGCAGttt   >  2:803598/1‑93 (MQ=255)
                                                                                       cACCGCTGGGTAATGTCACCGTGATCAACGGTTACTTCCCGCAGGGTGAAAGCCGCGACCATCCGATAAAATTCCCGGCAAAAGCGCAGtttt  >  2:6365747/1‑93 (MQ=255)
                                                                                          |                                                                                         
CCAAAGAGACGCCGATTGCCGTGCGTCGCGGCTTTCCCGGTGACGACGAAGAGGCGCAGCGGCGGATTATTATGGCGGAAATCCCCTCACTGCTGGGTAATGTCACCGTGATCAACGGTTACTTCCCGCAGGGTGAAAGCCGCGACCATCCGATAAAATTCCCGGCAAAAGCGCAGTTTT  >  NC_000913/1830657‑1830836

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