Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,182,712 G→A 32.1% A55T (GCC→ACC)  zupT → zinc transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,182,7120GA32.1% 114.3 / 67.2 90A55T (GCC→ACC) zupTzinc transporter
Reads supporting (aligned to +/- strand):  ref base G (29/32);  new base A (14/15);  total (43/47)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.47e-01

CGTTCTCGGGCAAAAACCCTCGAACCGCTTACTGGCGTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTGCCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTCTGCTTGGCTATTTTGGCCTGGACCGCATGTTGCCACATG  >  NC_000913/3182615‑3182811
                                                                                                 |                                                                                                   
cGTTCTCGGGCAAAAACCCTCGAACCGCTTACTGGCGTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTAcc                                                                                                   <  2:1289949/100‑1 (MQ=255)
 gTTCTCGGGCAAAAACCCTCGAACCGCTTACTGGCGTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTACCg                                                                                                  <  2:837509/100‑1 (MQ=255)
  ttCTCGGGCAAAAACCCTCGAACCGCTTACTGGCGTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTGCCGc                                                                                                 <  1:421647/100‑1 (MQ=255)
    ctcGGGCAAAAACCCTCGAACCGCTTACTGGCGTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTGCCGCAc                                                                                               >  1:449372/1‑100 (MQ=255)
     tcGGGCAAGAACCCTCGAACCGCTTACTGGCGTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTGCCGCACt                                                                                              <  1:408437/100‑1 (MQ=255)
     tcGGGCAAAAACCCTCGAACCGCTTACTGGCGTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTGCCGCACt                                                                                              <  1:1858154/100‑1 (MQ=255)
      cGGGCAAAAACCCTCGAACCGCTTACTGGCGTTTTCGCTATGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTGCCGCACTa                                                                                             >  1:2238210/1‑100 (MQ=255)
      cGGGCAAAAACCCTCGAACCGCTTACTGGCGTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTGCCGCACTa                                                                                             >  1:288193/1‑100 (MQ=255)
      cGGGCAAAAACCCTCGAACCGCTTACTGGCGTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTGCCGCACTa                                                                                             <  2:2182458/100‑1 (MQ=255)
         gCAAAAACCCTCGAACCGCTTACTGGCGTTTTCGCTAGGCTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTGCCGCACTagc                                                                                           >  2:428110/1‑99 (MQ=255)
          cAAAAACCCTCGAACCGCTTACTGGCGTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTGCCGCACTagcag                                                                                         >  2:1236244/1‑100 (MQ=255)
          cAAAAACCCTCGAACCGCTTACTGGCGTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTGCCGCACTagcag                                                                                         >  2:254724/1‑100 (MQ=255)
          cAAAAACCCTCGAACCGCTTACTGGCGTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTGCCGCACTagcag                                                                                         >  1:964566/1‑100 (MQ=255)
          cAAAAACCCTCGAACCGCTTACTGGCGTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTGCCGCACTagcag                                                                                         >  1:1708895/1‑100 (MQ=255)
           aaaaaCCCTCGAACCGCTTACTGGCGTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTGCCGCACTagcagc                                                                                        <  2:1794444/100‑1 (MQ=255)
            aaaaCCCTCGAACCGCTTACTGGCGTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTGCCGCACTAGCAGCt                                                                                       <  1:235954/100‑1 (MQ=255)
             aaaCCCTCGAACCGCTTACTGGCGTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTGCCGCACTAGCAGCTg                                                                                      >  1:1920225/1‑100 (MQ=255)
                cccTCGAACCGCTTACTGGCGTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTACCGCACTAGCAGCTGAAg                                                                                   >  2:870969/1‑100 (MQ=255)
                 ccTCGAACCGCTTACTGGCGTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTACCGCACTAGCAGCTGAAgg                                                                                  >  1:1657472/1‑100 (MQ=255)
                   tCGAACCGCTTACTGGCGTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTGCCGCACTAGCAGCTGAAGGaa                                                                                <  1:1001336/100‑1 (MQ=255)
                      aaCCGCTTACTGGCGTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTGCCGCACTAGCAGCTGAAGGAATg                                                                              <  2:655969/99‑1 (MQ=255)
                       aCCGCTTACTGGCGTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTGCCGCACTAGCAGCTGAAGGAATGTc                                                                            >  2:1260128/1‑100 (MQ=255)
                        ccGCTTACTGGCGTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTGCCGCACTAGCAGCTGAAGGAATGTCg                                                                           <  1:2161791/100‑1 (MQ=255)
                          gCTTACTGGCGTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTGCCGCACTAGCAGCTGAAGGAATGTCGcc                                                                         <  2:964564/100‑1 (MQ=255)
                          gCTTACTGGCGTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTACCGCACTAGCAGCTGAAGGAATGTCGcc                                                                         <  2:1657461/100‑1 (MQ=255)
                           cTTAGTGGCGTTTTCGCTCGGTTTTGCGGCGGTGATAATGCTGCTCTTCTCATTAATGGTGATGCTTCCTGCCGGACTAGCAGCTGAAGGAATg                                                                              <  2:636434/94‑1 (MQ=255)
                            ttACTGGCGTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTGCCGCACTAGCAGCTGAAGGAATATCGCCtg                                                                       <  1:1169820/100‑1 (MQ=255)
                            ttACTGGCGTTTTCGCTAGGTTGTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTACCGCACTAGCAGCTGAAGGAATGTCGCCtg                                                                       <  1:872901/100‑1 (MQ=255)
                               cTGGCGTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTGCCGCACTAGCAGCTGACGGAATGTCGCCTgtgt                                                                    >  2:768326/1‑100 (MQ=255)
                               cTGGCGTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTGCCGCACTAGCAGCTGAAGGAATGTCGCCTgtgt                                                                    >  2:1815545/1‑100 (MQ=255)
                                tGGCGTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTACCGCACTAGCAGCTGAAGGAATGTCGCCTGTGtt                                                                   >  1:1108416/1‑100 (MQ=255)
                                 ggCGTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTGCCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTg                                                                  >  2:348977/1‑100 (MQ=255)
                                 ggCGTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTACCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTg                                                                  >  1:2171572/1‑100 (MQ=255)
                                    gTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTGCCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGt                                                               >  1:574836/1‑100 (MQ=255)
                                       ttCGCTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTGCCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTAt                                                            <  2:1127244/100‑1 (MQ=255)
                                         cGCTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTGCCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATgg                                                          >  2:1476672/1‑100 (MQ=255)
                                         cGCTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTGCCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATg                                                           >  1:688329/1‑99 (MQ=255)
                                          gCTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTACCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGa                                                         >  1:903351/1‑100 (MQ=255)
                                           cTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTACCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGaa                                                        >  1:33948/1‑100 (MQ=255)
                                             aGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTGCCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATg                                                      <  1:24636/100‑1 (MQ=255)
                                               gTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTGCCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGtt                                                    <  1:56129/100‑1 (MQ=255)
                                               gTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTGCCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGtt                                                    <  1:67782/100‑1 (MQ=255)
                                               gTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTGCCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGtt                                                    <  1:251225/100‑1 (MQ=255)
                                               gTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTACCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGtt                                                    <  1:1412957/100‑1 (MQ=255)
                                                 tttGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTGCCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTa                                                  <  1:476060/100‑1 (MQ=255)
                                                 tttGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCCGCCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTa                                                  <  1:2021921/100‑1 (MQ=255)
                                                    gcggcgGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTGCCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCt                                               >  2:2167723/1‑100 (MQ=255)
                                                     cggcggGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTGCCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCt                                               >  1:1883505/1‑99 (MQ=255)
                                                      ggcggGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTACCGCACTTGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTc                                             <  2:373836/100‑1 (MQ=255)
                                                       gcggGGGTCATGTTGCTCATCTCATTAATGGAAATGCTTCCTGACGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTAGCTTCg                                            >  2:713588/1‑100 (MQ=255)
                                                       gcggGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTGCCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCg                                            <  2:3540/100‑1 (MQ=255)
                                                       gcggGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTACCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCg                                            <  1:798445/100‑1 (MQ=255)
                                                        cggGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTGCCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCgg                                           <  1:891715/100‑1 (MQ=255)
                                                           ggATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTGCCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTCt                                        <  2:1917424/100‑1 (MQ=255)
                                                           ggATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTACCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTCt                                        >  1:2156281/1‑100 (MQ=255)
                                                                 tGTTGCTCATCTCATTAATGGAAATGCTTCCTGCCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTCTGCTTgg                                  >  2:1850677/1‑100 (MQ=255)
                                                                 tGTTGCTCATCTCATTAATGGAAATGCTTCCTGCCGCACTAGCAGCTGAAGAAATGTCGCCTGTGTCGGGTTATGGAATGTTTATCTTCGGTCt                                        >  2:1850682/1‑94 (MQ=255)
                                                                 tGTTGCTCATCTCATTAATGGAAATGCTTCCTACCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTCTGCTTg                                   >  1:1045259/1‑99 (MQ=255)
                                                                  gTTGCTCATCTCATTAATGGAAATGCTTCCTGCCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTCTGCTTGGc                                 <  1:1627936/100‑1 (MQ=255)
                                                                  gTTGCTCATCTCATTAATGGAAATGCTTCCTACTGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGGATGGTTATCTTCGGTCTGCTTGGc                                 >  2:1449922/1‑100 (MQ=255)
                                                                   ttGCTCATCTCATTAATGGAAATGCTTCCTGCCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTCTGCTTGGCt                                >  1:1495598/1‑100 (MQ=255)
                                                                    tGCTCATCTCATTAATGGAAATGCTTCCTACCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTCTGCTTGGCTa                               >  2:787218/1‑100 (MQ=255)
                                                                     gCTCATCTCATTAATGGAAATGCTTCCTACCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTCTGCTTGGCTAt                              >  2:1041623/1‑100 (MQ=255)
                                                                      ctcatctcatTAATGGAAATGCTTCCTACCGCACTAGCAGCTTAAGGAATGTAGCCTGTGTTGGGTTATGGACTGTTTATCTTCGGTCTGTTTgg                                  >  2:2093231/1‑95 (MQ=255)
                                                                          tctcatTAATGGAAATGCTTCCTGCCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTCTGCTTGGCTATTTTgg                         <  1:624248/100‑1 (MQ=255)
                                                                          tctcatTAATGGAAATGCTTCCTGCCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTCTGCTTGGCTATTTTgg                         >  1:2248872/1‑100 (MQ=255)
                                                                            tcatTAATGGAAATGCTTCCTGCCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTCTGCTTGGCTATTTTGGcc                       <  1:878461/100‑1 (MQ=255)
                                                                            tcatTAATGGAAATGCTTCCTACCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTCTGCTTGGCTATTTTGGcc                       <  2:41154/100‑1 (MQ=255)
                                                                            tcatTAATGGAAATGCTTCCTACCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTCTGCTTGGCTATTTTGGcc                       <  1:1555095/100‑1 (MQ=255)
                                                                             catTAATGGAAATGCTTCCTGCCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTCTGCTTGGCTATTTTGGCCt                      >  1:1801821/1‑100 (MQ=255)
                                                                             catTAATGGAAATGCTTCCTACCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTCTGCTTGGCTATTTTGGCCt                      <  1:870970/100‑1 (MQ=255)
                                                                             catTAATGGAAATGCTTCCTACCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTCTGCTTGGCTATTTTGGCCt                      >  1:1010008/1‑100 (MQ=255)
                                                                              atTAATGGAAATGCTTCCTGCCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTCTGCTTGGCTATTTTGGCCTg                     <  1:1988018/100‑1 (MQ=255)
                                                                               ttAATGGAAATGCTTCCTACCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTCTGCTTGGCTATTTTGGCCTgg                    <  1:1507814/100‑1 (MQ=255)
                                                                               ttAATGGAAATGCTTCCTACCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTCTGCTTGGCTATTTTGGCCTgg                    >  1:514333/1‑100 (MQ=255)
                                                                                tAATGGAAATGCTTCCTGCCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTCTGCTTGGCTATTTTGGCCTGGa                   >  1:280135/1‑100 (MQ=255)
                                                                                    ggAAATGCTTCCTGCCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTCTGCTTGGCTATTTTGGCCTGGACCGc               >  1:874524/1‑100 (MQ=255)
                                                                                     gAAATGCTTCCTGCCGCACTAGCAGCTGAAGGAATGTTGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTCTGCTTGGCTATTTTGGCCTGGACCGCa              >  2:1194265/1‑100 (MQ=255)
                                                                                     gAAATGCTTCCTGCCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTCTGCTTGGCTATTTTGGCCTGGACCGCa              <  2:1034819/100‑1 (MQ=255)
                                                                                     gAAATGCTTCCTACCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTCTGCTTGGCTATTTTGGCCTGGACCGCa              <  2:590893/100‑1 (MQ=255)
                                                                                      aaaTGCTTCCTACCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTCTGCTTGGCTATTTTGGCCTGGACCGCAt             <  2:1617768/100‑1 (MQ=255)
                                                                                      aaaTGCTTCCTACCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTCTGCTTGGCTATTTTGGCCTGGACCACAt             <  1:1626438/100‑1 (MQ=255)
                                                                                        aTGCTTCCTGCCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTCTGCTTGGCTATTTTGGCCTGGACCGCATGt           <  1:792522/100‑1 (MQ=255)
                                                                                        aTGCTTCCTGCCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTCTGCTTGGCTATTTTGGCCTGGACCGCATGt           <  1:952327/100‑1 (MQ=255)
                                                                                          gCTTCCTGCCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTCTGCTTGGCTATTTTGGCCTGGACCGCATGTTg         <  1:225427/100‑1 (MQ=255)
                                                                                          gCTTCCTGCCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTCTGCTTGGCTATTTTGGCCTGGACCGCATGTTg         <  2:1708884/100‑1 (MQ=255)
                                                                                            ttCCTGCCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTCTGCTTGGCTATTTTGGCCTGGACCGCATGTTGcc       <  2:449373/100‑1 (MQ=255)
                                                                                            ttCCTACCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTCTGCTTGGCTATTTTGGCCTGGACCGCATGTTGcc       <  2:144670/100‑1 (MQ=255)
                                                                                               cTGCCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTCTGCTTGGCTATTTTGGCCTGGACCGCATGTTGCcaca    >  1:2160823/1‑100 (MQ=255)
                                                                                                tGCCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTCTGCTTGGCTATTTTGGCCTGGACCGCATGTTGCCACAt   >  1:1004420/1‑100 (MQ=255)
                                                                                                 gCCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTCTGCTTGGCTATTTTGGCCTGGACCGCATGTTGCCACATg  >  2:886080/1‑100 (MQ=255)
                                                                                                 |                                                                                                   
CGTTCTCGGGCAAAAACCCTCGAACCGCTTACTGGCGTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTGCCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTCTGCTTGGCTATTTTGGCCTGGACCGCATGTTGCCACATG  >  NC_000913/3182615‑3182811

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: