Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,950,966:1 +C 15.8% coding (135/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,9651.C15.8% 397.0 / 78.7 120coding (134/567 nt)yecDisochorismatase family protein
Reads supporting (aligned to +/- strand):  ref base . (48/53);  new base C (9/10);  total (57/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 = 2.83e-01

Only 100 of 127 total aligned reads displayed.

GGTGATCGATTTACAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑CGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGATGCGCCCTCACCGGCTAAAGT  >  NC_000913/1950867‑1951064
                                                                                                   |                                                                                                   
ggtgATCGATTTACAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTcg‑c                                                                                                    <  1:2542324/100‑1 (MQ=255)
ggtgATCGATTTACAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTcg‑c                                                                                                    <  1:2680209/100‑1 (MQ=255)
 gtgATCGATTTACAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTcg‑cg                                                                                                   <  2:830779/100‑1 (MQ=255)
  tgATCGATTTACAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATAGCGCCGGGAAGCTGGCGGAGAAATTTcc‑cg                                                                                                   >  1:1831731/1‑99 (MQ=255)
   gATCGATTTACAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑CGcc                                                                                                 >  2:230546/1‑100 (MQ=255)
   gATCGATTTACAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCAGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑CGcc                                                                                                 >  1:2811781/1‑100 (MQ=255)
     tCGATTTACAAGAAGGCATCTTACCTTTTGGCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGC‑GCGAAATTTcg‑ccaccag                                                                                              <  1:659618/100‑7 (MQ=255)
     tCGATTTACAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAATCTGGC‑GCGAAATTTCGCCGCCAg                                                                                               <  2:2262453/100‑1 (MQ=255)
      cGATTTACAAGAAGGCATCTTGCCTTTTGCCGGTGGTCCACATACGGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGC‑GCGAAATTTCGCCGCCAg                                                                                               <  1:2334000/99‑1 (MQ=255)
      cGATTTACAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGC‑GCGAAATTTCGCCGCCAGc                                                                                              <  2:2202324/100‑1 (MQ=255)
      cGATTTACAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGC‑GCGAAATTTCGCCGCCAGc                                                                                              <  2:1923038/100‑1 (MQ=255)
       gagttACAAGAAGGCGTCTTAGCTTTTGCCGGTGGTCCACATACGGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGC‑GCGAAATTTCGCCGCCAGc                                                                                              <  1:1884472/96‑1 (MQ=37)
       gATTTACAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAAACGCCCCGGGAAGCTGGC‑GCGAAATTTCGCCGCCAGc                                                                                              <  2:1283063/99‑1 (MQ=38)
        aTTTACAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑CGCCAGCgg                                                                                            <  2:2547760/100‑1 (MQ=255)
        aTTTACAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑CGCCAGCgg                                                                                            <  1:77240/100‑1 (MQ=255)
           tACAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCCGGCGGCGAAATTTCG‑CGCCAGCGGTCa                                                                                         >  1:490888/1‑100 (MQ=255)
            aCAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑CGCCAGCGGTCAg                                                                                        <  2:1840321/100‑1 (MQ=255)
            aCAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑CGCCAGCGGTCAg                                                                                        <  2:1193859/100‑1 (MQ=255)
             cAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑CGCCAGCGGTCAGc                                                                                       >  1:1596033/1‑100 (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)
                 aagGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGGTGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑CGCCAGCGGTCAGCCCgt                                                                                   >  2:2896354/1‑100 (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)
                         ttACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCTGGAAGCTGGCGGCGAAATTTCG‑CGCCAGCGGTCAGCCCGTGTTTCTgg                                                                           >  1:1370329/1‑100 (MQ=255)
                         ttACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑CGCCAGCGGTCAGCCCGTGTTTCTgg                                                                           <  2:2274960/100‑1 (MQ=255)
                          tACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑CGCCAGCGGTCAGCCCGTGTTTCTGGt                                                                          >  1:2368966/1‑100 (MQ=255)
                            ccTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑CGCCAGCGGTCAGCCCGTGTTTCTGGTgc                                                                        <  2:495634/100‑1 (MQ=255)
                              ttttGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑CGCCAGCGGTCAGCCCGTGTTTCTGGTgcgc                                                                      >  1:908792/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)
                                         ggccACATACTGCCGAGGAGGTAGTTAACCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑CGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTc                                                           <  2:342902/98‑1 (MQ=255)
                                          tCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑CGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCt                                                          <  2:961272/100‑1 (MQ=255)
                                          tCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑CGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCt                                                          >  2:2387978/1‑100 (MQ=255)
                                                  cTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑CGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTa                                                  <  2:540563/100‑1 (MQ=255)
                                                  cTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑CGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTa                                                  >  1:1322803/1‑100 (MQ=255)
                                                    gCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑CGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACg                                                >  1:2892614/1‑100 (MQ=255)
                                                    gCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑CGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACg                                                >  1:2364101/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)
                                                     ccGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑CGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGc                                               <  2:37188/100‑1 (MQ=255)
                                                      cGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑CGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGc                                               >  1:202382/1‑99 (MQ=255)
                                                      cGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑CGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCt                                              <  2:2230461/100‑1 (MQ=255)
                                                          gAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑CGCCAGCGGTCAGCCCTTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGagg                                          <  2:1971091/100‑3 (MQ=255)
                                                          gAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑CGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAg                                          <  1:1712877/100‑1 (MQ=255)
                                                          gAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑CGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAg                                          <  2:59771/100‑1 (MQ=255)
                                                          gAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑CGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAg                                          <  1:416487/100‑1 (MQ=255)
                                                            ggTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑CGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCa                                        <  2:1067374/100‑1 (MQ=255)
                                                             gTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑CGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCAt                                       >  1:2245993/1‑100 (MQ=255)
                                                              tAGTTAATCGCGCCGGGATGCTGGCGGTGAATGTTTG‑TGCTAGCGGTCAGCCCGTGTTTCTGATTCGCGTGGGCTGGTCTGCCGATTATGCTGAAGCAtt                                      <  1:1122092/100‑1 (MQ=255)
                                                              tAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑CGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCAtt                                      >  1:2545082/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                                    <  1:10622/100‑1 (MQ=255)
                                                                gTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑CGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTaa                                    >  2:640719/1‑100 (MQ=255)
                                                                 ttAATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑CGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTaaa                                   >  2:135755/1‑100 (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)
                                                                       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:484241/100‑1 (MQ=255)
                                                                           cGGGAAGCTGGCGGCGAAATTTCG‑CGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTc                         <  1:2355161/100‑1 (MQ=255)
                                                                            gggAAGCTGGCGGCGAAATTTCG‑CGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCg                        <  2:1311798/100‑1 (MQ=255)
                                                                             gaccgCTGGC‑GCGAAATTTCGCCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGa                       >  2:651148/5‑100 (MQ=38)
                                                                               ccgCTGGC‑GCGAAATTTCGCCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGATg                     >  2:2287581/3‑100 (MQ=39)
                                                                               ccgCTGGC‑GCGAAATTTCGCCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGATg                     >  2:2584300/3‑100 (MQ=39)
                                                                               ccgCTGGC‑GCGAAATTTCGCCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGATg                     >  2:2799747/3‑100 (MQ=39)
                                                                               ccgCTGGC‑GCGAAATTTCGCCGCCAGCGGTCAGCACGTGTTTCTGGTGCGCGTTGGCTGCTCTGCCGATTACGCTGAAGAATTAAAACAGCCGGTCGATg                     >  1:276824/3‑100 (MQ=37)
                                                                               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)
                                                                                   tGGCGGCGAAATTTCG‑CGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGATGCGcc                 <  1:2029826/100‑1 (MQ=255)
                                                                                    ggcggcGAAATTTCG‑CGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGATGCGccc                <  1:1093601/100‑1 (MQ=255)
                                                                                    ggcggcGAAATTTCG‑CGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGATGCGccc                >  2:1152175/1‑100 (MQ=255)
                                                                                        gcgAAATTTCG‑CGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGATGCGCCCTCAc            >  2:1173354/1‑100 (MQ=255)
                                                                                         ggAAATTTCG‑CGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGATGGTCTGCCGGGTACGCTGAAGCATTAAAACAGCCGGTCGATGCGCCCTCAcc           >  2:2605828/2‑100 (MQ=255)
                                                                                         cgAAATTTCG‑CGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGATGCGCCCTCAcc           <  2:2364116/100‑1 (MQ=255)
                                                                                           aaaTTTCG‑CGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGATGCGCCCTCACCgg         <  2:1410554/100‑1 (MQ=255)
                                                                                              tttCG‑CGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGATGCGCCCTCACCGGCTa      <  2:1514875/100‑1 (MQ=255)
                                                                                               ttCG‑CGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGATGCGCCCTCACCGGCTaa     >  2:832495/1‑100 (MQ=255)
                                                                                                 cg‑cgCCAGCGGTCAGCCCGTTTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGATGCGCCCTCACCGGCTAAAg   <  1:2061018/100‑1 (MQ=255)
                                                                                                   |                                                                                                   
GGTGATCGATTTACAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑CGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGATGCGCCCTCACCGGCTAAAGT  >  NC_000913/1950867‑1951064

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: ‑