Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,950,955 Δ1 bp 16.8% coding (124/567 nt) yecD → isochorismatase family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,950,9540G.16.8% 320.0 / 70.7 119coding (123/567 nt)yecDisochorismatase family protein
Reads supporting (aligned to +/- strand):  ref base G (44/52);  new base . (9/11);  total (56/63)
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.50e-02

Only 100 of 127 total aligned reads displayed.

CACTGCGCTGGTGGTGATCGATTTACAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGATGCGCCC  >  NC_000913/1950855‑1951050
                                                                                                   |                                                                                                 
cACTGCGCTGGTGGTGATCGATTTACAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTggcg                                                                                                   <  2:2102505/100‑1 (MQ=255)
 aCTGCGCTGGTGGTGATCGATTTACAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTggcgg                                                                                                  <  2:1668210/100‑1 (MQ=255)
  cTGCGCTGGTGGTGATCGATTTACAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTggcggc                                                                                                 <  1:1032237/100‑1 (MQ=255)
  cTGCGCTGGTGGTGATCGATTTACAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTggcggc                                                                                                 <  2:312503/100‑1 (MQ=255)
  cTGCGCTGGTGGTGATCGATTTACAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGAAGAGGTAGTTAATCGCGCCGGGAAGCTggcggc                                                                                                 >  1:993102/1‑100 (MQ=255)
  cTGCGCTGGTGGTGATCGATTTACAAGAAGGCATCTTACCGTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTggcggc                                                                                                 <  2:2613424/100‑1 (MQ=255)
       ntggtggTGATCGATTTACAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAAt                                                                                            >  1:1318999/2‑100 (MQ=255)
       cTGGTGGTGATCGATTTACAAGAAGGCATCTTACCTTTTGTCGGTGGTCTTAATACTGGAAATGCTGTAGTTAATCGCGCCGGTACGCTGTCGGCGCAAt                                                                                            >  1:1015476/1‑100 (MQ=255)
         ggtggtGATCGATTTACAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAAttt                                                                                          <  2:1128149/100‑1 (MQ=255)
            ggtgATCGATTTACAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTcgc                                                                                      <  1:2542324/100‑1 (MQ=255)
            ggtgATCGATTTACAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTcgc                                                                                      <  1:2680209/100‑1 (MQ=255)
            ggtgATCGATTTACAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTc                                                                                         >  1:141557/1‑98 (MQ=255)
             gtgATCGATTTACAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTcgcg                                                                                     <  2:830779/100‑1 (MQ=255)
                 tCGATTTACAAGAAGGCATCTTACCTTTTGGCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGCGAAATTTcgccaccag                                                                                <  1:659618/100‑7 (MQ=255)
                 tCGATTTACAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAATCTGGCGCGAAATTTCGCCGCCAg                                                                                 <  2:2262453/100‑1 (MQ=255)
                  cGATTTACAAGAAGGCATCTTTCCTTTTGCCGGTGGTCCACATACTGCCGATGATGTAGTTAATCGCGCCGGGAAGCTGGCGCGAAATTTCGCCGCCAGc                                                                                <  1:2441335/100‑1 (MQ=255)
                  cGATTTACAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGCGAAATTTCGCCGCCAGc                                                                                <  2:2202324/100‑1 (MQ=255)
                  cGATTTACAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGCGAAATTTCGCCGCCAGc                                                                                <  2:1923038/100‑1 (MQ=255)
                  cGATTTACAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGCGAAATTTCGCCGCCAGc                                                                                <  2:1227393/100‑1 (MQ=255)
                  cGATTTACAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGCGAAATTTCGCCGCCAGc                                                                                <  2:2245190/100‑1 (MQ=255)
                   gagttACAAGAAGGCGTCTTAGCTTTTGCCGGTGGTCCACATACGGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGCGAAATTTCGCCGCCAGc                                                                                <  1:1884472/96‑1 (MQ=37)
                   gATTTACAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAAACGCCCCGGGAAGCTGGCGCGAAATTTCGCCGCCAGc                                                                                <  2:1283063/99‑1 (MQ=38)
                    aTTTACAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCTGCGAAATTTCGAGCCAGCgg                                                                              >  1:1315/1‑100 (MQ=255)
                    aTTTACAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCgg                                                                              <  2:2547760/100‑1 (MQ=255)
                    aTTTACAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCgg                                                                              >  2:411353/1‑100 (MQ=255)
                    aTTTACAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCgg                                                                              <  1:77240/100‑1 (MQ=255)
                     tttACAAGAAGGCATCTTAACTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCgg                                                                              <  2:1362465/99‑1 (MQ=255)
                       tACAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCCGGCGGCGAAATTTCGCGCCAGCGGTCa                                                                           >  1:490888/1‑100 (MQ=255)
                        aCAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAg                                                                          <  2:1840321/100‑1 (MQ=255)
                         cAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGACCGCGAAATTTCGCGCCAGCGCTCa                                                                           >  1:314220/1‑98 (MQ=255)
                          aagaagGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGcc                                                                        <  2:617957/100‑1 (MQ=255)
                           agaagGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGcc                                                                        >  1:145674/1‑99 (MQ=255)
                            gaagGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTACTCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCg                                                                      >  1:2339810/1‑100 (MQ=255)
                            gaagGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCg                                                                      >  2:607611/1‑100 (MQ=255)
                             aagGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGGTGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCgt                                                                     >  2:2896354/1‑100 (MQ=255)
                             aagGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGTAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCgt                                                                     <  1:1563211/100‑1 (MQ=255)
                                gCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGtt                                                                  <  2:2875098/100‑1 (MQ=255)
                                  aTCTTGCCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGTTTc                                                                <  1:951024/100‑1 (MQ=255)
                                  aTCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGtt                                                                  <  1:2271340/98‑1 (MQ=255)
                                  aTCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGTTTc                                                                >  2:1882122/1‑100 (MQ=255)
                                   tCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCt                                                               <  2:2590196/100‑1 (MQ=255)
                                     ttACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTgg                                                             <  2:2274960/100‑1 (MQ=255)
                                        ccTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGTgc                                                          <  2:495634/100‑1 (MQ=255)
                                          ttttGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGTgcgc                                                        >  1:908792/1‑100 (MQ=255)
                                          ttttGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGTgcgc                                                        >  1:1984573/1‑100 (MQ=255)
                                             tGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGtt                                                     <  2:629419/100‑1 (MQ=255)
                                              gCCGGTGGTCCACATACTGCCGATGAGGTAGTTACTCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTg                                                    >  1:557675/1‑100 (MQ=255)
                                               ccGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTgg                                                   <  1:293730/100‑1 (MQ=255)
                                                   tggtCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTgg                                               >  2:2367042/1‑100 (MQ=255)
                                                    ggtCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTgn                                               <  1:2539355/99‑2 (MQ=255)
                                                     gtCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTc                                             >  2:2335679/1‑100 (MQ=255)
                                                     ggccACATACTGCCGAGGAGGTAGTTAACCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTc                                             <  2:342902/98‑1 (MQ=255)
                                                      tCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCt                                            >  2:2387978/1‑100 (MQ=255)
                                                      tCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCt                                            <  2:961272/100‑1 (MQ=255)
                                                              cTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTa                                    <  2:540563/100‑1 (MQ=255)
                                                              cTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTa                                    >  1:1322803/1‑100 (MQ=255)
                                                                gCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACg                                  >  1:2364101/1‑100 (MQ=255)
                                                                gCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACg                                  >  1:2892614/1‑100 (MQ=255)
                                                                gCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGTGTCTGGTGCGCGTTGGCTGGTCTGCCGATTACg                                  >  1:250815/1‑100 (MQ=255)
                                                                gCCGATGAGGTAGTTAATCGCGCCGGCAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACg                                  >  2:1776416/1‑100 (MQ=255)
                                                                 ccGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGc                                 <  2:469183/100‑1 (MQ=255)
                                                                  cGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGc                                 >  1:202382/1‑99 (MQ=255)
                                                                      gAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCTTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGagg                            <  2:1971091/100‑3 (MQ=255)
                                                                      gAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAg                            <  1:1712877/100‑1 (MQ=255)
                                                                      gAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAg                            <  1:416487/100‑1 (MQ=255)
                                                                         gTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCAt                         >  1:2245993/1‑100 (MQ=255)
                                                                          tAGTTAATCGCGACGGGAAGCTGGCCGCGAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGTCTGGTCTGCCGATTACGCTGAAGCAtt                        >  2:142868/1‑100 (MQ=255)
                                                                           aGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTa                       >  1:2702104/1‑100 (MQ=255)
                                                                           aGTTAATAGCGCCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTa                       >  2:1532671/1‑100 (MQ=255)
                                                                            gTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTaa                      >  2:640719/1‑100 (MQ=255)
                                                                            gTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTaa                      <  1:10622/100‑1 (MQ=255)
                                                                             ttAATCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTaaa                     <  2:2212254/100‑1 (MQ=255)
                                                                             ttAATCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTaaa                     <  1:2168342/100‑1 (MQ=255)
                                                                              tAATCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTTGTGCGCTTTGGCTGGTCTGCCGATTACGCTAAAGCATTaaaa                    >  1:1555615/1‑100 (MQ=255)
                                                                              tAATCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTaaaa                    <  2:2513288/100‑1 (MQ=255)
                                                                               aaTCGCTCCGGGAATCTGGCGGCTAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTaaa                     <  2:1994995/98‑1 (MQ=255)
                                                                               aaTCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAAc                   >  1:1744330/1‑100 (MQ=255)
                                                                                aTCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGTCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACa                  <  1:657099/100‑1 (MQ=255)
                                                                                  cgcgTCGGGAAGCTGGCGGCGAAATTTGGCGCCAGCGGTCGGCCCGTGTTTCTGGTGCGCTTTGGCTGGTTTGCCGATTTCGCTGAAGCATTAAAACAGc                >  1:1154632/1‑100 (MQ=255)
                                                                                  cgcgCCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGc                >  2:2204545/1‑100 (MQ=255)
                                                                                  cgcgCCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGc                >  1:2663533/1‑100 (MQ=255)
                                                                                   gcgcCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGcc               <  1:896404/100‑1 (MQ=255)
                                                                                    cgcCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCg              <  2:1327102/100‑1 (MQ=255)
                                                                                       ctgaccgCTGGCGCGAAATTTCGCCGCCAGAGGGCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTc           >  1:1912898/7‑100 (MQ=37)
                                                                                       cGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTc           <  1:2355161/100‑1 (MQ=255)
                                                                                       cGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTc           <  1:484241/100‑1 (MQ=255)
                                                                                        gggAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCg          <  2:1311798/100‑1 (MQ=255)
                                                                                           ccgCTGGCGCGAAATTTCGCCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGATg       >  2:1070349/3‑100 (MQ=39)
                                                                                           ccgCTGGCGCGAAATTTCGCCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGATg       >  2:2287581/3‑100 (MQ=39)
                                                                                           ccgCTGGCGCGAAATTTCGCCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGATg       >  2:2799747/3‑100 (MQ=39)
                                                                                           ccgCTGGCGCGAAATTTCGCCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGATg       >  2:2584300/3‑100 (MQ=39)
                                                                                           ccgCTGGCGCGAAATTTCGCCGCCAGCGGTCAGCACGTGTTTCTGGTGCGCGTTGGCTGCTCTGCCGATTACGCTGAAGAATTAAAACAGCCGGTCGATg       >  1:276824/3‑100 (MQ=37)
                                                                                           aaGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGATg       <  2:803137/100‑1 (MQ=255)
                                                                                           aaGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGATg       <  2:2045277/100‑1 (MQ=255)
                                                                                            aGCTGGCGGCGGAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTTGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGATgc      <  1:784937/100‑1 (MQ=255)
                                                                                             gCTGGCGCGAAATTTCGCCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGATgcg     >  2:2854242/1‑100 (MQ=255)
                                                                                              cTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGATgcgc    <  1:1606966/100‑1 (MQ=255)
                                                                                              cTGGCGCGAAATTTCGCCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGATgcgc    >  1:2153976/1‑100 (MQ=255)
                                                                                                ggcggcGAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGATGCGccc  >  2:1152175/1‑100 (MQ=255)
                                                                                                ggcggcGAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGATGCGccc  <  1:1093601/100‑1 (MQ=255)
                                                                                                   |                                                                                                 
CACTGCGCTGGTGGTGATCGATTTACAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGATGCGCCC  >  NC_000913/1950855‑1951050

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: