Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 356,492 A→G 46.8% T108A (ACG→GCG)  codA → cytosine/isoguanine deaminase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913356,4920AG46.8% 18.7 / 248.7 186T108A (ACG→GCG) codAcytosine/isoguanine deaminase
Reads supporting (aligned to +/- strand):  ref base A (53/46);  new base G (47/40);  total (100/86)
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.34e-01

Only 100 of 190 total aligned reads displayed.

CTGGAATCAGTCCGGCACGCTGTTTGAAGGCATTGAACGCTGGGCCGAGCGCAAAGCGTTATTAACCCATGACGATGTGAAACAACGCGCATGGCAAACGCTGAAATGGCAGATTGCCAACGGCATTCAGCATGTGCGTACCCATGTCGATGTTTCGGATGCAACGCTAACTGCGCTGAAAGCAATGCTGGAAGTGA  >  NC_000913/356395‑356591
                                                                                                 |                                                                                                   
cTGGAATCAGTCCGGCACGCTGTTTGAAGGCATTGAACGCTGGGCCGAGCGCAAAGCGTTATTAACCCATGACGATGTGAAACAACGCGCATGGCAAGCg                                                                                                   >  1:3370067/1‑100 (MQ=255)
cTGGAATCAGTCCGGCACGCTGTTTGAAGGCATTGAACGCTGGGCCGAGCGCAAAGCGTTATTAACCCATGACGATGTGAAACAACGCGCATGGCAAACg                                                                                                   >  2:3260031/1‑100 (MQ=255)
  ggAATCAGTCCGGCACGCTGTTTGAAGGCATTGAACGCTGGGCCGAGCGCAAAGCGTTATTAACCCATGACGATGTGAAACAACGCGCATGGCAAGCGCt                                                                                                 >  1:2937160/1‑100 (MQ=255)
    aaTCAGTCCGGCACGCTGTTTGAAGGCATTGAACGCTGGGCCGAGCGCAAAGCGTTATTAACCCATGACGATGTGAAACAACGCGCATGGCAAGCGCTGa                                                                                               <  2:3719209/100‑1 (MQ=255)
    aaTCAGTCCGGCACGCTGTTTGAAGGCATTGAACGCTGGGCCGAGCGCAAAGCGTTATTAACCCATGACGATGTGAAACAACGCGCATGGCAAACGCTGa                                                                                               >  2:2911144/1‑100 (MQ=255)
       cAGTCCGGCACGCTGTTTGAAGGCATTGAACGCTGGGCCGAGAGCAAAGCGTTATTAACCCATGACGATGTGAAACAACGCGCATGGCAAACGCTGAAAt                                                                                            <  1:2993073/100‑1 (MQ=255)
        aGTCCGGCACGCTGTTTGAAGGCATTGAACGCTGGGCCGAGCGCAAAGCGTTATTAACCCATGACGATGTGAAACAACGCGCATGGCAAACGCTGAAATg                                                                                           <  1:260963/100‑1 (MQ=255)
         gTCCGGCACGCTGTTTGAAGGCATTGAACGCTGGGCCGAGCGCAAAGCGTTATTAACCCATGACGATGTGAAACAACGCGCATGGCAAGCGCTGAAATgg                                                                                          <  1:1685894/100‑1 (MQ=255)
          tCCGGCACGCTGTTTGAAGGCATTGAACGCTGGGCCGAGCGCAAAGCGTTATTAACCCATGACGATGTGAAACAACGCGCATGGCAAACGCTGAAATGGc                                                                                         <  2:1129664/100‑1 (MQ=255)
           ccGTCACGCGGTTTGAAGGCATTGAACGCTGGGCCGAGCGCAACGCGTTATTAACCCATGACGATGTGAAACAACGCGCATGGCAAACGCTGAAATGGCa                                                                                        >  1:1587692/1‑100 (MQ=255)
           ccGGCACGCTGTTTGAAGGCATTGGACGCTGGGCCGAGCGCAAAGCGTTATTAACCCATGACGATGTGAAACAACGCGCATGGCAAACGCTGAAATg                                                                                           <  1:4003622/97‑1 (MQ=255)
           ccGGCACGCTGTTTGAAGGCATTGAACGCTGGGCCGAGCGCAAAGCGTTATTAACCCATGACGATGTGAAACAACGCGCATGGCAAGCGCTGAAATGGCa                                                                                        <  1:538043/100‑1 (MQ=255)
           ccGGCACGCTGTTTGAAGGCATTGAACGCTGGGCCGAGCGCAAAGCGTTATTAACCCATGACGATGTGAAACAACGCGCATGGCAAACGCTGAAATGGCa                                                                                        >  1:1625110/1‑100 (MQ=255)
           ccGGCACGCTGTTTGAAGGCATTGAACGCTGGGCCGAGCGCAAAGCGTTATTAACCCATGACGATGTGAAACAACGCGCATGGCAAACGCTGAAATGGCa                                                                                        >  1:3178128/1‑100 (MQ=255)
           ccGGCACGCTGTTTGAAGGCATTGAACGCTGGGCCGAGCGCAAAGCGTTATTAACCCATGACGATGTGAAACAACGCGCATGGCAAACGCTGAAATGGCa                                                                                        >  2:1825688/1‑100 (MQ=255)
            cGGCACGCTGTTTGAAGGCATTGAGCGCTGGGCCGAGCGCAAAGCGGTATTAACCCATGACGATGTGAAACAACGCGCATGGCAAGCGCTGAAATGGCAg                                                                                       <  1:3438903/100‑1 (MQ=255)
            cGGCACGCTGTTTGAAGGCATTGAACGCTGGGCCGAGCGCAAAGCGTTATTAACCCATGACGATGTGAAACAACGCGCATGGCAAACGCTGAAATGGCAg                                                                                       >  1:2076429/1‑100 (MQ=255)
             ggCACGCTGTTTGAAGGCATTGAACGCTGGGCCGAGCGCAAAGCGTTATTAACCCATGACGATGTGAAACAACGCGCATGGCAAGCGCTGAAATGGCAGa                                                                                      >  1:2323573/1‑100 (MQ=255)
             ggCACGCTGTTTGAAGGCATTGAACGCTGGGCCGAGCGCAAAGCGTTATTAACCCATGACGATGTGAAACAACGCGCATGGCAAACGCTGAAATGGCAGa                                                                                      <  1:3279832/100‑1 (MQ=255)
              gCACGCTGTTTGAAGGCATTGAACGCTGGGCCGAGCGCAAAGCGTTATTAACCCATGACGATGTGAAACAACGCGCATGGCAAGCGCTGAAATGGCAGAt                                                                                     >  2:1516156/1‑100 (MQ=255)
              gCACGCTGTTTGAAGGCATTGAACGCTGGGCCGAGCGCAAAGCGTTATTAACCCATGACGATGTGAAACAACGCGCATGGCAAGCGCTGAAATGGCAGAt                                                                                     >  2:1769884/1‑100 (MQ=255)
              gCACGCTGTTTGAAGGCATTGAACGCTGGGCCGAGCGCAAAGCGTTATTAACCCATGACGATGTGAAACAACGCGCATGGCAAACGCTGAAATGGCAGAt                                                                                     >  1:863353/1‑100 (MQ=255)
                tcGCTGTTTGAAGGCATTGAACGCTGGGCCGAGCGCAAAGCGTTATTAACCCATGACGATGTGAAACAACGCGCATGGCAAGCGCTGAAATGGCAGATTg                                                                                   <  1:2519170/99‑1 (MQ=255)
                aCGCTGTTTGAAGGCATTGAACGCTGGGCCGAGCGCAAAGCGTTATTAACCCATGACGATGTGAAACAACGCGCATGGCAAACGCTGAAATGGCAGATTg                                                                                   <  1:811006/100‑1 (MQ=255)
                 cGCTGTTTGAAGGCATTGAACGCTGGGCCGAGCGCAAAGCGTTATTAACCCATGACGATGTGAAACAACGCGCATGGCAAACGCTGAAATGGCAGATTGc                                                                                  <  2:832129/100‑1 (MQ=255)
                  gCTGTTTGAAGGCATTGAACGCTGGGCCGAGCGCAAAGCGTTATTAACCCATGACGATGTGAAACAACGCGCATGGCAAGCGCTGAAATGGCAGATTGcc                                                                                 >  2:3481951/1‑100 (MQ=255)
                   cTGTTTGAAGGCATTAAACGCTGGGCCGAGCGCAAAGAGTTATTAACCCATGACGATGTGAAACAACGCGCATGGCAAGCTCTGAAATGGCAGATTGcc                                                                                 >  1:2018639/1‑99 (MQ=255)
                     gtatGAAGGCATTGAACGCTGGGCCTAGCGCAAAGCGTTATTAACACATGACGATGTGAAACAACGGGCATGGCAAGCGCTGAAATGGCAGATTGCTaa                                                                               >  2:257976/4‑99 (MQ=255)
                     gTTTGAAGGCAGTGAACGCTGGGCCGAGCGCAAAGCGTTATTAACCCATGACGATGTGAAACAACGCGCATGGCAAACGCTGAAATGGCAGATTGCCaa                                                                               <  1:3478003/99‑1 (MQ=255)
                      tttGAAGGCATTGAACGCTGGGCCGAGCGCAAAGCGTTATTAACCCATGACGATGTGAAACAACGCGCATGGCAAACGCTGAAATGGCAGATTGCCAACg                                                                             >  1:3608472/1‑100 (MQ=255)
                        tGAAGGCATTGAACGCTGGGCCGAGCGCAAAGCGTTATTAACCCATGACGATGTGAAACAACGCGCATGGCAAGCGCTGAAATGGCAGATTGCCAACGGc                                                                           <  1:482700/100‑1 (MQ=255)
                           aGGCATTGAACGCTGGGCCGAGCGCAAAGCGTTATTAACCCATGACGATGTGAAACAACGCGCATGGCAAACGCTGAAATGGCAGATTGCCAACGGCAtt                                                                        >  2:420769/1‑100 (MQ=255)
                               aTTGAACGCTGGGCCGAGCGCAAAGCGTTATTAACCCATGACGATGTGAAACAACGCGCATGGCAAACGCTGAAATGGCAGATTGCCAACGGCATTCAGc                                                                    >  2:3439451/1‑100 (MQ=255)
                                ttGAACGCTGGGCCGAGCGCAAAGCGTTATTAACCCATGACGATGTGAAACAACGCGCATGGCAAGCGCTGAAATGGCAGATTGCCAACGGCATTCAGCa                                                                   <  1:1139216/100‑1 (MQ=255)
                                ttGAACGCTGGGCCGAGCGCAAAGCGTTATTAACCCATGACGATGTGAAACAACGCGCATGGCAAACGCTGAAATGGCAGATTGCCAACGGCATTCAGCa                                                                   <  2:4141000/100‑1 (MQ=255)
                                 tGAACGCTGGGCCGAGCGCAAAGCGTTATTAACCCATGACGATGTGAAACAACGCGCATGGCAAACGCTGAAATGGCAGATTGCCAACGGCATTCAGCAt                                                                  >  2:2221119/1‑100 (MQ=255)
                                   aaCGCTGGGCCGAGCGCAAAGCGTTATTAACCCATGACGATGTGAAACAACGCGCATGGCAAACGCTGAAATGGCAGATTGCCAACGGCATTCAGCAtgt                                                                <  2:1095321/100‑1 (MQ=255)
                                    aCGCTGGGCCGAGCGCAAAGCGTTATTAACCCATTACGATGTGAAACAACGCGCATGGCAAGCGCTGAAATGGCAGATTGCCAACGGCATTCAGCAtgtg                                                               <  2:4008857/100‑1 (MQ=255)
                                    aCGCTGGGCCGAGCGCAAAGCGTTATTAACCCATGACGATGTGAAACAACGCGCATGGCAAGCGCTGAAATGGCAGATTGCCAACGGCATTCAGCAtgtg                                                               <  1:560639/100‑1 (MQ=255)
                                    aCGCTGGGCCGAGCGCAAAGCGTTATTAACCCATGACGATGTGAAACAACGCGCATGGCAAGCGCTGAAATGGCAGATTGCCAACGGCATTCAGCAtgtg                                                               <  1:317194/100‑1 (MQ=255)
                                      gCTGGGCCGAGCGCAAAGCGTTATTAACCCATGACGATGTGAAACAACGCGCATGGCAAGCGCTGAAATGGCAGATTGCCAACGGCATTCAGCATGTGCg                                                             <  1:734253/100‑1 (MQ=255)
                                      gCTGGGCCGAGCGCAAAGCGTTATTAACCCATGACGATGTGAAACAACGCGCATGGCAAGCGCTGAAATGGCAGATTGCCAACGGCATTCAGCATGTGCg                                                             <  1:27029/100‑1 (MQ=255)
                                        tGGGCCGAGCGCAAAGCGTTATTAACCCATGACGATGTGAAACAACGCGCATGGCAAACGCTGAAATGGCAGATTGCCAACGGCATTCAGCATGTGCGTa                                                           <  2:3237613/100‑1 (MQ=255)
                                         gggCCGAGCGCAAAGCGTTATTAACCCATGACGATGTGAAACAACGCGCATGGCAAGCGCTGAAATGGCAGATTGCCAACGGCATTCAGCATGTGCGTAc                                                          >  2:3930418/1‑100 (MQ=255)
                                         gggCCGAGCGCAAAGCGTTATTAACCCATGACGATGTGAAACAACGCGCATGGCAAACGCTGAAATGGCAGATTGCCAACGGCATTCAGCATGTGCGTAc                                                          >  2:160705/1‑100 (MQ=255)
                                          ggCCGAGCGCAAAGCGTTATTAACCCATGACGATGTGAAACAACGCGCATGGCAAACGCTGAAATGGCAGATTGCCAACGGCATTCAGCATGTGCGTAcc                                                         <  2:2076433/100‑1 (MQ=255)
                                          ggCCGAGCGCAAAGCGTTATTAACCCATGACGATGTGAAACAACGCGCATGGCAAACGCTGAAATGGCAGATTGCCAACGGCATTCAGCATGTGCGTAcc                                                         >  2:1409070/1‑100 (MQ=255)
                                           gCCGAGCGCAAAGCGTTATTAACCCATGACGATGTGAAACAACGCGCATGGCAAACGCTGAAATGGCAGATTGCCAACGGCATTCAGCATGTGCGTAccc                                                        >  2:195780/1‑100 (MQ=255)
                                            ccGAGCGCAAAGCGTTATTAACCCATGACGATGTGACACAACGCGCATGGCAAGCGCTGAAATGGCAGATTGCCAACGGCATTCAGCATGTGCGTACCCa                                                       >  2:3020493/1‑100 (MQ=255)
                                            ccGAGCGCAAAGCGTTATTAACCCATGACGATGTGAAACAACGCGCATGGCAAGCGCTGAAATGGCAGATTGCCAACGGCATTCAGCATGTGCGTACCCa                                                       <  1:2218890/100‑1 (MQ=255)
                                            ccGAGCGCAAAGCGTTATTAACCCATGACGATGTGAAACAACGCGCATGGCAAACGCTGAAATGGCAGATTGCCAACGGCATTCAGCATGTGCGTACCCa                                                       >  1:3824599/1‑100 (MQ=255)
                                            ccGAGCGCAAAGCGTTATTAACCCATGACGATGTGAAACAACGCGCATGGCAAACGCTGAAATGGCAGATTGCCAACGGCATTCAGCATGTGCGTACCCa                                                       <  1:2776632/100‑1 (MQ=255)
                                             cGAGCGCAAAGCGTTATTAACCCATGACGATGTGAAACAACGCGCATGGCAAGCGCTGAAATGGCAGATTGCCAACGGCATTCAGCATGTGCGTACCCAt                                                      >  2:4271007/1‑100 (MQ=255)
                                              gAGCGCAAAGCGTTATTAACCCATGACGATGTGAAACAACGCGCATGGCAAACGCTGAAATGGCAGATTGCCAACGGCATTCAGCATGTGCGTACCCATg                                                     >  2:1296527/1‑100 (MQ=255)
                                                 cgcAAAGCGTTATTAACCCATGACGATGTGAAACAACGCGCATGGCAAACGCTGAAATGGCAGATTGCCAACGGCATTCAGCATGTGCGTACCCATGTCg                                                  >  1:262773/1‑100 (MQ=255)
                                                   cAAAGCGTTATTAACCCATGACGATGTGAAACAACGCGCATGGCAAGCGCTGCAATGGCAGATTGCCAACGGCATTCAGCATGTGCGTACCCATGTCGAt                                                >  2:65009/1‑100 (MQ=255)
                                                   cAAAGCGTTATTAACCCATGACGATGTGAAACAACGCGCATGGCAAACGCTGAAATGGCAGATTGCCAACGGCATTCAGCATGTGCGTACCCATGTCGAt                                                >  2:2631912/1‑100 (MQ=255)
                                                    aaaGCGTTATTAACCCATGACGATGTGAAACAACGCGCATGGCAAGCGCTGAAATGGCAGATTGCCAACGGCATTCAGCATGTGCGTACCCATGTCGATg                                               <  2:1858775/100‑1 (MQ=255)
                                                      aGCGTTATTAACCCATGACGATGTGAAACAACGCGCATGGCAAGCGCTGAAATGGCAGATTGCCAACGGCATTCAGCATGTGCGTACCCATGTCGATGt                                              <  1:2148923/99‑1 (MQ=255)
                                                         gTTATTAACCCATGACGATGTGAAACAACGCGCATGGCAAACGCTGAAATGGCAGATTGCCAACGGCATTCAGCATGTGCGTACCCATGTCGATGTTTCg                                          <  2:1082145/100‑1 (MQ=255)
                                                            attaACCCATGACGATGTGAAACAACGCGCATGGCAAACGCTGAAATGGCAGATTGCCAACGGCATTCAGCATGTGCGTACCCATGTCGATGTTTCGGAt                                       <  1:2424964/100‑1 (MQ=255)
                                                            attaACCCATGACGATGTGAAACAACGCGCATGGCAAACGCTGAAATGGCAGATTGCCAACGGCATTCAGCATGTGCGTACCCATGTCGATGTTTCGGAt                                       <  1:1495168/100‑1 (MQ=255)
                                                             ttaACCCATGACGATGTGAAACAACGCGCATGGCAAGCGCTGAAATGGCAGATTGCCAACGGCATTCAGCATGTGCGTACTCATGTCGATGTTTCGGATg                                      <  1:1043996/100‑1 (MQ=255)
                                                             ttaACCCATGACGATGTGAAACAACGCGCATGGCAAGCGCTGAAATGGCAGATTGCCAACGGCATTCAGCATGTGCGTACCCATGTCGATGTTTCGGATg                                      <  2:1735780/100‑1 (MQ=255)
                                                              taACCCATGACGATGTGAAACAACGCGCATGGCAAGCGCTGAAATGGCAGATTGCCAACGGCATTCAGCATGTGCGTACCCATGTCGATGTTTCGGATGc                                     <  2:3265111/100‑1 (MQ=255)
                                                              taACCCATGACGATGTGAAACAACGCGCATGGCAAACGCTGAAATGGCAGATTGCCAACGGCATTCAGCATGTGCGTACCCATGTCGATGTTTCGGATGc                                     >  2:2689675/1‑100 (MQ=255)
                                                                aCCCATGACGATGTGAAACAACGCGCATGGCAAACGCTGAAATGGCAGATTGCCAACGGCATTCAGCATGTGCGTACCCATGTCGATGTTTCGGATGCaa                                   <  1:1029542/100‑1 (MQ=255)
                                                                 cccATGACGATGTGAAACAACGCGCATGGCAAACGCTGAAATGGCAGATTGCCAACGGCATTCAGCATGTGCGTACCCATGTCGATGTTTCGGATGCAAc                                  <  2:3912565/100‑1 (MQ=255)
                                                                  ccATGACGATGTGAAACAACGCGCATGGCAAGCGCTGAAATGGCAGATTGCCAACGGCATTCAGCATGTGCGTACCCATGTCGATGTTTCGGATGCAACg                                 >  1:1389931/1‑100 (MQ=255)
                                                                   cATGACGATGTGAAACAACGCGCATGGCAAGCGCTGAAATGGCAGATTGCCAACGGCATTCAGCATGTGCGTACCCATGTCGATGTTTCGGATGCAACGc                                >  1:4179236/1‑100 (MQ=255)
                                                                       aCGATGTGAAACAACGCGCATGGCAAACGCTGAAATGGCAGATTGCCAACGGCATTCAGCATGTGCGTACCCATGTCGATGTTTCGGATGCAACGCTAAc                            <  2:456837/100‑1 (MQ=255)
                                                                            gtgAAACAACGCGCATGGCAAGCGCTGAAATGGCAGATTGCCAACGGCATTCAGCATGTGCGTACCCATGTCGATGTTTCGGATGCAACGCTAACTgcgc                       <  2:61818/100‑1 (MQ=255)
                                                                            gtgAAACAACGCGCATGGCAAGCGCTGAAATGGCAGATTGCCAACGGCATTCAGCATGTGCGTACCCATGTCGATGTTTCGGATGCAACGCTAACTgcgc                       <  2:4130185/100‑1 (MQ=255)
                                                                            gtgAAACAACGCGCATGGCAAACGCTGAAATGGCAGATTGCCAACGGCATTCAGCATGTGCGTACCCATGTCGATGTTTCGGATGCAACGCTAACTgcgc                       <  2:2511066/100‑1 (MQ=255)
                                                                             tgAAACAACGCGCATGGCAAGCGCTGAAATGGCAGATTGCCAACGGCATTCAGCATGTGCGTACCCATGTCGATGTTTCGGATGCAACGCTAACTGCGCt                      >  2:2739807/1‑100 (MQ=255)
                                                                              gAAACAGCGCGCATGGCAAGCGCTGAAATGGCAGATTGCCAACGGCATTCAGCATGTGCGTACCCATGTCGATGTTTCGGATGCAACGCTAACTGCGCTg                     <  2:3807483/100‑1 (MQ=255)
                                                                              gAAACAACGCGCATGGCAAACGCTGAAATGGCAGATTGCCAACGGCATTCAGCATGTGCGTACCCATGTCGATGTTTCGGATGCAACGCTAACTgcgc                       >  1:1057657/1‑98 (MQ=255)
                                                                                aacaacGCGCATGGCAAGCGCTGAAATGGCAGATTGCCAACGGCATTCAGCATGTGCGTACCCATGTCGATGTTTCGGATGCAACGCTAACTGCGCTGaa                   >  2:2618916/1‑100 (MQ=255)
                                                                                aacaacGCGCATGGCAAGCGCTGAAATGGCAGATTGCCAACGGCATTCAGCATGTGCGTACCCATGTCGATGTTTCGGATGCAACGCTAACTGCGCTGaa                   >  1:2889962/1‑100 (MQ=255)
                                                                                aacaacGCGCATGGCAAACGCTGAAATGGCAGATTGCCAACGGCATTCAGCATGTGCGTACCCATGTCGATGTTTCGGATGCAACGCTAACTGCGCTGaa                   <  1:204055/100‑1 (MQ=255)
                                                                                 acaacGCGCATGGCAAGCGCTGAAATGGCAGATTGCCAACGGCATTCAGCATGTGCGTACCCATGTCGATGTTTCGGATGCAACGCTAACTGCGCTGaaa                  >  2:1277846/1‑100 (MQ=255)
                                                                                  caacGCGCATGGCAAGCGCTGAAATGGCAGATTGCCAACGGCATTCAGCATGTGCGTACCCATGTCGATGTTTCGGATGCAACGCTAACTGCGCTGAAAg                 >  2:768550/1‑100 (MQ=255)
                                                                                  caacGCGCATGGCAAACGCTGAAATGGCAGATTGCCAACGGCATTCAGCATGTGCGTACCCATGTCGATGTTTCGGATGCAACGCTAACTGCGCTGAAAg                 <  2:2828028/100‑1 (MQ=255)
                                                                                   aacGCGCATGGCAAACGCTGAAATGGCAGATTGCCAACGGCATTCAGCATGTGCGTACCCATGTCGATGTTTCGGATGCAACGCTAACTGCGCTGAAAGc                >  1:395631/1‑100 (MQ=255)
                                                                                   aacGCGCATGGCAAACGCTGAAATGGCAGATTGCCAACGGCATTCAGCATGTGCGTACCCATGTCGATGTTTCGGATGCAACGCTAACTGCGCTGAAAGc                <  1:1116613/100‑1 (MQ=255)
                                                                                    acGCGCATGGCAAGCGCTGAAATGGCAGATTGCCAACGGCATTCAGCATGTGCGTACCCATGTCGATGTTTCGGATGCAACGCTAACTGCGCTGAAAGCa               >  1:2675580/1‑100 (MQ=255)
                                                                                       cgcATGGCAAGCGCTGAAATGGCAGATTGCCAACGGCATTCAGCATGTGCGTACCCATGTCGATGTTTCGGATGCAACGCTAACTGCGCTGAAAGCAATg            >  1:246662/1‑100 (MQ=255)
                                                                                       cgcATGGCAAACGCTGAAATGGCAGATTGAACACGGCATTAAGCATGTGCGTCCCCATGTCGATATTTCGCCTCCAACGCTACCTACGCTGAAAGCCAt             >  1:3461152/1‑99 (MQ=255)
                                                                                        gcATGGCTAACACTGAAATGGCAGATTGCCAACGGCATTCAGCATGTGCGTACCCATGTCGATGTCTCGGATGCAACGCTAACTGCGCTGAAAGCAATGc           >  1:2237374/1‑100 (MQ=255)
                                                                                            ggCAAGCGCTGAAATGGCAGATTGCCAACGGCATTCAGCATGTGCGTACCCATGTCGATGTTTCGGATGCAACGCTAACTGCGCTGAAAGCAATGCTGGa       <  1:4268535/100‑1 (MQ=255)
                                                                                             gCAAGCGCTGAAATGGCAGATTGCCAACGGCATTCAGCATGTGCGTACCCATGTCGATGTTTCGGATGCAACGCTAACTGCGCTGAAAGCAATGCTGGaa      <  2:4187956/100‑1 (MQ=255)
                                                                                             gCAAACGCTGAAATGGCAGATTGCCAACGGCATTCAGCATGTGCGTACCCATGTCGATGTTTCGGATGCAACGCTAACTGCGCTGAAAGCAATGCTGGaa      >  2:2292138/1‑100 (MQ=255)
                                                                                              cAAGCGCTGAAATGGCAGATTGCCAACGGCATTCAGTATGTGCGTACCCATGTCGATGTTTCGGATGCAACGCTAACTGCGCTGAAAGCAATGCTGGAAg     >  2:3007966/1‑100 (MQ=255)
                                                                                              cAAACGCTGAAATGGCAGATTGCCAACGGCATTCAGCATGTGCGTACCCATGTCGATGTTTCGGATGCAACGCTAACTGCGCTGAAAGCAATGCTGGAAg     <  1:2503601/100‑1 (MQ=255)
                                                                                              cAAACGCTGAAATGGCAGATTGCCAACGGCATTCAGCATGTGCGTACCCATGTCGATGTTTCGGATGCAACGCTAACTGCGCTGAAAGCAATGCTGGAAg     <  2:2412261/100‑1 (MQ=255)
                                                                                              cAAACGCTGAAATGGCAGATTGCCAACGGCATTCAGCATGTGCGTACCCATGTCGATGTTTCGGATGCAACGCTAACTGCGCTGAAAGCAATGCTGGAAg     <  2:4063323/100‑1 (MQ=255)
                                                                                               aaaCGCTGAAATGGCAGATTGCCAACGGCATTCAGCATGTGCGTACCCATGTCGATGTTTCGGATGCAACGCTAACTGCGCTGAAAGCAATGCTGGAAGt    >  1:3021518/1‑100 (MQ=255)
                                                                                                aaCGCTGAAATGGCAGATTGCCAACGGCATTCAGCATGTGCGTACCCATGTCGATGTTTCGGATGCAACGCTAACTGCGCTGAAAGCAATGCTGGAAGTg   <  2:2989783/100‑1 (MQ=255)
                                                                                                aaCGCTGAAATGGCAGATTGCCAACGGCATTCAGCATGTGCGTACCCATGTCGATGTTTCGGATGCAACGCTAACTGCGCTGAAAGCAATGCTGGAAGTg   >  2:1016106/1‑100 (MQ=255)
                                                                                                 aCGCTGAAATGGCAGATTGCCAACGGCATTCAGCATGTGCGTACCCATGTCGATGTTTCGGATGCAACGCTAACTGCGCTGAAAGCAATGCTGGAAGTGa  <  2:2279523/100‑1 (MQ=255)
                                                                                                 |                                                                                                   
CTGGAATCAGTCCGGCACGCTGTTTGAAGGCATTGAACGCTGGGCCGAGCGCAAAGCGTTATTAACCCATGACGATGTGAAACAACGCGCATGGCAAACGCTGAAATGGCAGATTGCCAACGGCATTCAGCATGTGCGTACCCATGTCGATGTTTCGGATGCAACGCTAACTGCGCTGAAAGCAATGCTGGAAGTGA  >  NC_000913/356395‑356591

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