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.

CACTGCGCTGGTGGTGATCGATTTACAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑CGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGATGCGCCC  >  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)
            ggtgATCGATTTACAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTcg‑c                                                                                      <  1:2542324/100‑1 (MQ=255)
            ggtgATCGATTTACAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTcg‑c                                                                                      <  1:2680209/100‑1 (MQ=255)
            ggtgATCGATTTACAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTc                                                                                         >  1:141557/1‑98 (MQ=255)
             gtgATCGATTTACAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTcg‑cg                                                                                     <  2:830779/100‑1 (MQ=255)
                 tCGATTTACAAGAAGGCATCTTACCTTTTGGCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGC‑GCGAAATTTcg‑ccaccag                                                                                <  1:659618/100‑7 (MQ=255)
                 tCGATTTACAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAATCTGGC‑GCGAAATTTCGCCGCCAg                                                                                 <  2:2262453/100‑1 (MQ=255)
                  cGATTTACAAGAAGGCATCTTTCCTTTTGCCGGTGGTCCACATACTGCCGATGATGTAGTTAATCGCGCCGGGAAGCTGGC‑GCGAAATTTCGCCGCCAGc                                                                                <  1:2441335/100‑1 (MQ=255)
                  cGATTTACAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGC‑GCGAAATTTCGCCGCCAGc                                                                                <  2:2202324/100‑1 (MQ=255)
                  cGATTTACAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGC‑GCGAAATTTCGCCGCCAGc                                                                                <  2:1923038/100‑1 (MQ=255)
                  cGATTTACAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGC‑GCGAAATTTCGCCGCCAGc                                                                                <  2:1227393/100‑1 (MQ=255)
                  cGATTTACAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGC‑GCGAAATTTCGCCGCCAGc                                                                                <  2:2245190/100‑1 (MQ=255)
                   gagttACAAGAAGGCGTCTTAGCTTTTGCCGGTGGTCCACATACGGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGC‑GCGAAATTTCGCCGCCAGc                                                                                <  1:1884472/96‑1 (MQ=37)
                   gATTTACAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAAACGCCCCGGGAAGCTGGC‑GCGAAATTTCGCCGCCAGc                                                                                <  2:1283063/99‑1 (MQ=38)
                    aTTTACAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCTGCGAAATTTCG‑AGCCAGCgg                                                                              >  1:1315/1‑100 (MQ=255)
                    aTTTACAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑CGCCAGCgg                                                                              <  2:2547760/100‑1 (MQ=255)
                    aTTTACAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑CGCCAGCgg                                                                              >  2:411353/1‑100 (MQ=255)
                    aTTTACAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑CGCCAGCgg                                                                              <  1:77240/100‑1 (MQ=255)
                     tttACAAGAAGGCATCTTAACTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑CGCCAGCgg                                                                              <  2:1362465/99‑1 (MQ=255)
                       tACAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCCGGCGGCGAAATTTCG‑CGCCAGCGGTCa                                                                           >  1:490888/1‑100 (MQ=255)
                        aCAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑CGCCAGCGGTCAg                                                                          <  2:1840321/100‑1 (MQ=255)
                         cAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGACCGCGAAATTTCG‑CGCCAGCGCTCa                                                                           >  1:314220/1‑98 (MQ=255)
                          aagaagGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑CGCCAGCGGTCAGcc                                                                        <  2:617957/100‑1 (MQ=255)
                           agaagGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑CGCCAGCGGTCAGcc                                                                        >  1:145674/1‑99 (MQ=255)
                            gaagGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTACTCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑CGCCAGCGGTCAGCCCg                                                                      >  1:2339810/1‑100 (MQ=255)
                            gaagGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑CGCCAGCGGTCAGCCCg                                                                      >  2:607611/1‑100 (MQ=255)
                             aagGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGGTGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑CGCCAGCGGTCAGCCCgt                                                                     >  2:2896354/1‑100 (MQ=255)
                             aagGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGTAGCTGGCGGCGAAATTTCG‑CGCCAGCGGTCAGCCCgt                                                                     <  1:1563211/100‑1 (MQ=255)
                                gCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑CGCCAGCGGTCAGCCCGTGtt                                                                  <  2:2875098/100‑1 (MQ=255)
                                  aTCTTGCCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑CGCCAGCGGTCAGCCCGTGTTTc                                                                <  1:951024/100‑1 (MQ=255)
                                  aTCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑CGCCAGCGGTCAGCCCGTGtt                                                                  <  1:2271340/98‑1 (MQ=255)
                                  aTCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑CGCCAGCGGTCAGCCCGTGTTTc                                                                >  2:1882122/1‑100 (MQ=255)
                                   tCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑CGCCAGCGGTCAGCCCGTGTTTCt                                                               <  2:2590196/100‑1 (MQ=255)
                                     ttACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑CGCCAGCGGTCAGCCCGTGTTTCTgg                                                             <  2:2274960/100‑1 (MQ=255)
                                        ccTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑CGCCAGCGGTCAGCCCGTGTTTCTGGTgc                                                          <  2:495634/100‑1 (MQ=255)
                                          ttttGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑CGCCAGCGGTCAGCCCGTGTTTCTGGTgcgc                                                        >  1:908792/1‑100 (MQ=255)
                                          ttttGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑CGCCAGCGGTCAGCCCGTGTTTCTGGTgcgc                                                        >  1:1984573/1‑100 (MQ=255)
                                             tGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑CGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGtt                                                     <  2:629419/100‑1 (MQ=255)
                                              gCCGGTGGTCCACATACTGCCGATGAGGTAGTTACTCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑CGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTg                                                    >  1:557675/1‑100 (MQ=255)
                                               ccGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑CGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTgg                                                   <  1:293730/100‑1 (MQ=255)
                                                   tggtCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑CGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTgg                                               >  2:2367042/1‑100 (MQ=255)
                                                    ggtCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑CGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTgn                                               <  1:2539355/99‑2 (MQ=255)
                                                     gtCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑CGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTc                                             >  2:2335679/1‑100 (MQ=255)
                                                     ggccACATACTGCCGAGGAGGTAGTTAACCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑CGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTc                                             <  2:342902/98‑1 (MQ=255)
                                                      tCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑CGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCt                                            >  2:2387978/1‑100 (MQ=255)
                                                      tCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑CGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCt                                            <  2:961272/100‑1 (MQ=255)
                                                              cTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑CGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTa                                    <  2:540563/100‑1 (MQ=255)
                                                              cTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑CGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTa                                    >  1:1322803/1‑100 (MQ=255)
                                                                gCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑CGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACg                                  >  1:2364101/1‑100 (MQ=255)
                                                                gCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑CGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACg                                  >  1:2892614/1‑100 (MQ=255)
                                                                gCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑CGCCAGCGGTCAGCCCGTGTGTCTGGTGCGCGTTGGCTGGTCTGCCGATTACg                                  >  1:250815/1‑100 (MQ=255)
                                                                gCCGATGAGGTAGTTAATCGCGCCGGCAAGCTGGCGGCGAAATTTCG‑CGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACg                                  >  2:1776416/1‑100 (MQ=255)
                                                                 ccGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑CGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGc                                 <  2:469183/100‑1 (MQ=255)
                                                                  cGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑CGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGc                                 >  1:202382/1‑99 (MQ=255)
                                                                      gAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑CGCCAGCGGTCAGCCCTTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGagg                            <  2:1971091/100‑3 (MQ=255)
                                                                      gAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑CGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAg                            <  1:1712877/100‑1 (MQ=255)
                                                                      gAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑CGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAg                            <  1:416487/100‑1 (MQ=255)
                                                                         gTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑CGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCAt                         >  1:2245993/1‑100 (MQ=255)
                                                                          tAGTTAATCGCGACGGGAAGCTGGCCGCGAAATTTCG‑CGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGTCTGGTCTGCCGATTACGCTGAAGCAtt                        >  2:142868/1‑100 (MQ=255)
                                                                           aGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑CGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTa                       >  1:2702104/1‑100 (MQ=255)
                                                                           aGTTAATAGCGCCGGGAAGCTGGCGGCGAAATTTCG‑CGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTa                       >  2:1532671/1‑100 (MQ=255)
                                                                            gTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑CGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTaa                      >  2:640719/1‑100 (MQ=255)
                                                                            gTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑CGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTaa                      <  1:10622/100‑1 (MQ=255)
                                                                             ttAATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑CGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTaaa                     <  2:2212254/100‑1 (MQ=255)
                                                                             ttAATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑CGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTaaa                     <  1:2168342/100‑1 (MQ=255)
                                                                              tAATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑CGCCAGCGGTCAGCCCGTGTTTCTTGTGCGCTTTGGCTGGTCTGCCGATTACGCTAAAGCATTaaaa                    >  1:1555615/1‑100 (MQ=255)
                                                                              tAATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑CGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTaaaa                    <  2:2513288/100‑1 (MQ=255)
                                                                               aaTCGCTCCGGGAATCTGGCGGCTAAATTTCG‑CGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTaaa                     <  2:1994995/98‑1 (MQ=255)
                                                                               aaTCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑CGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAAc                   >  1:1744330/1‑100 (MQ=255)
                                                                                aTCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑CGCCAGCGGTCAGTCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACa                  <  1:657099/100‑1 (MQ=255)
                                                                                  cgcgTCGGGAAGCTGGCGGCGAAATTTGG‑CGCCAGCGGTCGGCCCGTGTTTCTGGTGCGCTTTGGCTGGTTTGCCGATTTCGCTGAAGCATTAAAACAGc                >  1:1154632/1‑100 (MQ=255)
                                                                                  cgcgCCGGGAAGCTGGCGGCGAAATTTCG‑CGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGc                >  2:2204545/1‑100 (MQ=255)
                                                                                  cgcgCCGGGAAGCTGGCGGCGAAATTTCG‑CGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGc                >  1:2663533/1‑100 (MQ=255)
                                                                                   gcgcCGGGAAGCTGGCGGCGAAATTTCG‑CGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGcc               <  1:896404/100‑1 (MQ=255)
                                                                                    cgcCGGGAAGCTGGCGGCGAAATTTCG‑CGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCg              <  2:1327102/100‑1 (MQ=255)
                                                                                       ctgaccgCTGGC‑GCGAAATTTCGCCGCCAGAGGGCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTc           >  1:1912898/7‑100 (MQ=37)
                                                                                       cGGGAAGCTGGCGGCGAAATTTCG‑CGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTc           <  1:2355161/100‑1 (MQ=255)
                                                                                       cGGGAAGCTGGCGGCGAAATTTCG‑CGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTc           <  1:484241/100‑1 (MQ=255)
                                                                                        gggAAGCTGGCGGCGAAATTTCG‑CGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCg          <  2:1311798/100‑1 (MQ=255)
                                                                                           ccgCTGGC‑GCGAAATTTCGCCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGATg       >  2:1070349/3‑100 (MQ=39)
                                                                                           ccgCTGGC‑GCGAAATTTCGCCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGATg       >  2:2287581/3‑100 (MQ=39)
                                                                                           ccgCTGGC‑GCGAAATTTCGCCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGATg       >  2:2799747/3‑100 (MQ=39)
                                                                                           ccgCTGGC‑GCGAAATTTCGCCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGATg       >  2:2584300/3‑100 (MQ=39)
                                                                                           ccgCTGGC‑GCGAAATTTCGCCGCCAGCGGTCAGCACGTGTTTCTGGTGCGCGTTGGCTGCTCTGCCGATTACGCTGAAGAATTAAAACAGCCGGTCGATg       >  1:276824/3‑100 (MQ=37)
                                                                                           aaGCTGGCGGCGAAATTTCG‑CGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGATg       <  2:803137/100‑1 (MQ=255)
                                                                                           aaGCTGGCGGCGAAATTTCG‑CGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGATg       <  2:2045277/100‑1 (MQ=255)
                                                                                            aGCTGGCGGCGGAATTTCG‑CGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTTGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGATgc      <  1:784937/100‑1 (MQ=255)
                                                                                             gCTGGC‑GCGAAATTTCGCCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGATgcg     >  2:2854242/1‑100 (MQ=255)
                                                                                              cTGGCGGCGAAATTTCG‑CGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGATgcgc    <  1:1606966/100‑1 (MQ=255)
                                                                                              cTGGC‑GCGAAATTTCGCCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGATgcgc    >  1:2153976/1‑100 (MQ=255)
                                                                                                ggcggcGAAATTTCG‑CGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGATGCGccc  >  2:1152175/1‑100 (MQ=255)
                                                                                                ggcggcGAAATTTCG‑CGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGATGCGccc  <  1:1093601/100‑1 (MQ=255)
                                                                                                   |                                                                                                 
CACTGCGCTGGTGGTGATCGATTTACAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑CGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGATGCGCCC  >  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: ‑