Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,950,955 Δ1 bp 19.9% 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.19.9% 387.3 / 102.0 141coding (123/567 nt)yecDisochorismatase family protein
Reads supporting (aligned to +/- strand):  ref base G (53/58);  new base . (18/10);  total (73/68)
Fisher's exact test for biased strand distribution p-value = 1.41e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.64e-01

Only 100 of 153 total aligned reads displayed.

CACTGCGCTGGTGGTGATCGATTTACAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGATGCGCCCTCA  >  NC_000913/1950855‑1951053
                                                                                                   |                                                                                                    
cACTGCGCTGGTGGTGATCGATTTACAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTggcg                                                                                                      <  1:987717/100‑1 (MQ=255)
 aCTGCGCTGGTGGTGATCGATTTACAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGCAGTTAATCGCGCCGGGAAGCTggcgg                                                                                                     <  1:2549488/100‑1 (MQ=255)
  cTGCGCTGGTGGTGATCGATTTACAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTggcggc                                                                                                    >  1:3822020/1‑100 (MQ=255)
  cTGCGCTGGTGGTGATCGATTTACAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTggcggc                                                                                                    <  1:2471848/100‑1 (MQ=255)
         ggtggtGATCGATTTACAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAAttt                                                                                             >  1:3703865/1‑100 (MQ=255)
              tgATCGATTTACAAGAAGGCATCTTACCTTTTGCCGGTGGCCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCCGCGAGATATCgcgc                                                                                       >  2:3667567/1‑100 (MQ=255)
               gATCGATTTACAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGcc                                                                                      >  2:3830998/1‑100 (MQ=255)
               gATCGATTTACAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGcc                                                                                      >  1:545161/1‑100 (MQ=255)
               gATCGATTTACAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGcc                                                                                      >  1:437000/1‑100 (MQ=255)
                aTCGATTTACAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGCGAAATTTcgccgcca                                                                                    <  1:50302/100‑6 (MQ=255)
                  cGATTTACAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGCGAAATTTCGCCGCCAGc                                                                                   <  2:283856/100‑1 (MQ=255)
                  cGATTTACAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGCGAAATTTCGCCGCCAGc                                                                                   <  1:446794/100‑1 (MQ=255)
                  cGATTTACAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGCGAAATTTCGCCGCCAGc                                                                                   <  1:3585111/100‑1 (MQ=255)
                  cGATTTACAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGCGAAATTTCGCCGCCAGc                                                                                   <  2:735125/100‑1 (MQ=255)
                  cGATTTACAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGCGAAATTTCGCCGCCAGc                                                                                   <  1:2090255/100‑1 (MQ=255)
                   gATTTACAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCg                                                                                  <  1:4032182/100‑1 (MQ=255)
                   gATTTACAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCg                                                                                  >  2:1233594/1‑100 (MQ=255)
                     tttACAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGt                                                                                >  1:2343984/1‑100 (MQ=255)
                     tttACAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGt                                                                                <  2:59487/100‑1 (MQ=255)
                      ttACAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGTc                                                                               <  2:3722358/100‑1 (MQ=255)
                       tACAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACCGCCGATGAGGTAGTTAATCGCGCCGGGACGCTGGCGGCGAAATTTCGCGCCAGCGGTc                                                                               >  2:857409/1‑99 (MQ=255)
                         cAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGc                                                                            >  1:1877481/1‑100 (MQ=255)
                         cAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGc                                                                            >  1:1019702/1‑100 (MQ=255)
                         cAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGCGAAATTTCGCCGCCAGcttcccgg                                                                            <  1:1590239/100‑8 (MQ=37)
                         cAAGAAGGCAGCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGc                                                                            <  1:178285/100‑1 (MQ=255)
                          aagaagGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGcc                                                                           <  1:1370575/100‑1 (MQ=255)
                          aagaagGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGcc                                                                           >  2:3406579/1‑100 (MQ=255)
                            taacgCATCGTACCTTTTGCCGGTTGTCCACATACTTCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGccct                                                                         <  2:4034852/96‑2 (MQ=255)
                              agGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCgtg                                                                       >  1:128840/1‑100 (MQ=255)
                                gCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGtt                                                                     <  1:318752/100‑1 (MQ=255)
                                  aTCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGTTTc                                                                   >  1:2176425/1‑100 (MQ=255)
                                  aTCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGTTTc                                                                   >  2:2252810/1‑100 (MQ=255)
                                  aTCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGTTTc                                                                   <  2:3196221/100‑1 (MQ=255)
                                      tACCTTTTGCCGGTGGTCCACATACCGCCATCGAGGTAGTTAATCGCCCCGCCTAGCTGGCCGCGAAATTTCGCGCCAGCGGTCAGCCCgtcgtcctg                                                                 >  2:3294103/1‑91 (MQ=255)
                                        ccTTTTGCCGGTGGTCCACGTACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGTgc                                                             <  2:1624834/100‑1 (MQ=255)
                                         cTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGTgcg                                                            >  1:694414/1‑100 (MQ=255)
                                          ttttGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGTgcgc                                                           >  1:674534/1‑100 (MQ=255)
                                              gCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTg                                                       <  2:399023/100‑1 (MQ=255)
                                                   tggtCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTgg                                                  <  1:1165798/100‑1 (MQ=255)
                                                    ggtCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGt                                                 >  2:1133779/1‑100 (MQ=255)
                                                    ggtCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGt                                                 <  2:1646105/100‑1 (MQ=255)
                                                     gtCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTc                                                <  2:547501/100‑1 (MQ=255)
                                                     gtCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTc                                                <  1:2252815/100‑1 (MQ=255)
                                                      tCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCt                                               >  1:1283673/1‑100 (MQ=255)
                                                      tCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCt                                               <  2:3029980/100‑1 (MQ=255)
                                                       ccACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTg                                              <  2:507412/100‑1 (MQ=255)
                                                       ccACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTg                                              >  1:3595588/1‑100 (MQ=255)
                                                          caTACTGGCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCg                                           <  2:200844/100‑1 (MQ=255)
                                                          caTACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCg                                           <  1:3408620/100‑1 (MQ=255)
                                                          caTACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCg                                           <  1:2268980/100‑1 (MQ=255)
                                                            tACTGCCGATGAGGTGGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGAt                                         >  2:3517402/1‑100 (MQ=255)
                                                              cTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTa                                       <  1:1090811/100‑1 (MQ=255)
                                                               tGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTAc                                      <  2:673906/100‑1 (MQ=255)
                                                                 ccGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGc                                    <  1:3512203/100‑1 (MQ=255)
                                                                   gATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTg                                  >  2:2448209/1‑100 (MQ=255)
                                                                       aGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGGCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGc                              <  2:823054/100‑1 (MQ=255)
                                                                       aGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGc                              <  2:2256344/100‑1 (MQ=255)
                                                                       aGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGc                              <  1:2448214/100‑1 (MQ=255)
                                                                          tAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCAtt                           <  1:1480327/100‑1 (MQ=255)
                                                                          tAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCAtt                           <  2:3940864/100‑1 (MQ=255)
                                                                           aGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTa                          <  2:128840/100‑1 (MQ=255)
                                                                           aGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTa                          >  1:786012/1‑100 (MQ=255)
                                                                           aGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTa                          <  1:3834494/100‑1 (MQ=255)
                                                                           aGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTa                          >  2:3008392/1‑100 (MQ=255)
                                                                            gTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTATGCTGAAGCATTaa                         <  1:2842052/100‑1 (MQ=255)
                                                                            gTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTaa                         >  1:1848476/1‑100 (MQ=255)
                                                                              tAATCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTaaaa                       >  2:739962/1‑100 (MQ=255)
                                                                               attCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAAc                      >  1:2850763/3‑100 (MQ=255)
                                                                               aaTCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAAc                      <  2:944827/100‑1 (MQ=255)
                                                                                aTCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACa                     <  1:102068/100‑1 (MQ=255)
                                                                                  cgcgCCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAGGCATTAAAACAGc                   >  1:2438315/1‑100 (MQ=255)
                                                                                  cgcgCCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGc                   >  1:383954/1‑100 (MQ=255)
                                                                                  cgcgCCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGc                   >  2:414164/1‑100 (MQ=255)
                                                                                  cgcgCCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGc                   >  2:1226691/1‑100 (MQ=255)
                                                                                    gggctgaccgCTGGCGCGAAATTTCGCCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCg                 >  1:3308136/10‑100 (MQ=37)
                                                                                     gcCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCgg                <  2:1058288/100‑1 (MQ=255)
                                                                                     gcCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCgg                <  1:1448664/100‑1 (MQ=255)
                                                                                        gggAATCTGGCGGCGAAATTTCGTGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCg             <  2:3831597/100‑1 (MQ=255)
                                                                                        gggAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCg             <  2:3435160/100‑1 (MQ=255)
                                                                                         ggAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCg             >  2:1013845/1‑99 (MQ=255)
                                                                                         gaccgCTGGCGCGAAATTTCGCCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGa            >  2:2433379/5‑100 (MQ=38)
                                                                                          gAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATGACGCTGAAGCATTAAAACAGCCGGTCGAt           >  2:2523413/1‑100 (MQ=255)
                                                                                          accgCTGGCGCGAAATTTCGCCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGAt           >  2:4148691/4‑100 (MQ=38)
                                                                                          accgCTGGCGCGAAATTTCGCCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGAt           >  2:273705/4‑100 (MQ=38)
                                                                                          accgCTGGCGCGAAATTTCGCCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGAt           <  2:1684699/97‑1 (MQ=38)
                                                                                           ccgCTGGCGCGAAATTTCGCCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGATg          >  2:3276916/3‑100 (MQ=39)
                                                                                           ccgCTGGCGCGAAATTTCGCCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGATg          >  1:3855001/3‑100 (MQ=39)
                                                                                           ccgCTGGCGCGAAATTTCGCCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGATg          >  1:2341857/3‑100 (MQ=39)
                                                                                           ccgCTGGCGCGAAATTTCGCCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGATg          >  1:1764160/3‑100 (MQ=39)
                                                                                           aaGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGATg          <  2:190163/100‑1 (MQ=255)
                                                                                            cgCTGGCGCGAAATTTCGCCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCAGACGATgc         >  2:771966/2‑100 (MQ=38)
                                                                                             gCTGGCGCGAAATTTCGCCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGCAGCATCAAAACAGCCGGTCGATgcg        >  2:2840210/1‑100 (MQ=255)
                                                                                             gCTGGCGCGAAATTTCGCCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGATgcg        >  1:408297/1‑100 (MQ=255)
                                                                                             gCTGGCGCGAAATTTCGCCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGATgcg        >  2:2686175/1‑100 (MQ=255)
                                                                                             gCTGGCGCGAAATTTCGCCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGATgcg        >  1:1978707/1‑100 (MQ=255)
                                                                                             gCTGGCGCGAAATTTCGCCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGATgcg        >  2:2729198/1‑100 (MQ=255)
                                                                                             gCTGGCGCGAAATTTCGCCGCCAGCGGTCAGCCCCTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGACGATgcg        >  2:2500569/1‑100 (MQ=255)
                                                                                                ggcggcGAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGATGCGccc     <  2:1338177/100‑1 (MQ=255)
                                                                                                 gcggcgAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTGTTGCGCGTGGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGATg          <  2:780285/94‑1 (MQ=255)
                                                                                                  cggcgAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGATGCGCCCTc   >  2:1714781/1‑100 (MQ=255)
                                                                                                   |                                                                                                    
CACTGCGCTGGTGGTGATCGATTTACAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGATGCGCCCTCA  >  NC_000913/1950855‑1951053

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