Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,427,336 T→A 12.6% D272V (GAC→GTC)  rhlE‑I ← ATP‑dependent RNA helicase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,427,3360TA12.6% 214.7 / 13.9 76D272V (GAC→GTC) rhlE‑IATP‑dependent RNA helicase
Reads supporting (aligned to +/- strand):  ref base T (29/37);  new base A (5/5);  total (34/42)
Fisher's exact test for biased strand distribution p-value = 7.45e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCAGGTCGAGGTTGACGACCAGCGGCAAGTCATCGATATCCAGGCCACGGGCCGCCACATCGGTCGCCACCAGCACTTGCACCTCACGGGCCTTGAAGCTGTCCAGCGCACGCTGGCGGGTTGCCTGTGGGCGGTCACCGTGGATGCCGTCGGCATTCACGCCTTCGGCCAGCAAGCGCTCCACCAGCTGATCGACGCCATTGCGGGTCTTGGCGAACACCAGCACCTGCTTCCAGCGCTGCTTGCGCAGCAGGTGGCAGAACAGGTCAGCCTTGCGCTTCTTGTCCACAGGCACC  >  NC_002947/5427187‑5427482
                                                                                                                                                     |                                                                                                                                                  
gCAGGTCGAGGTTGACGACCAGCGGCAAGTCATCGATATCCAGGCCACGGGCCGCCACATCGGTCGCCACCAGCACTTGCACCTCACGGGCCTTGAAGCTGTCCAGCGCACGCTGGCGGGTTGCCTGTGGGCGGTCACCGTGGATGCCGt                                                                                                                                                    >  1:231617/1‑150 (MQ=255)
      cGAGGTTGACGACCAGCGGCAAGTCATCGATATCCAGGCCACGGGCCGCCACATCGGTCGCCACCAGCACTTGCACCTCACGGGCCTTGAAGCTGTCCAGCGCACGCTGGCGGGTTGCCTGTGGGCGGTCACCGTGGATGCCGTCGGCAt                                                                                                                                              <  2:688989/150‑1 (MQ=255)
      cGAGGTTGACGACCAGCGGCAAGTCATCGATATCCAGGCCACGGGCCGCCACATCGGTCGCCACCAGCACTTGCACCTCACGGGCCTTGAAGCTGTCCAGCGCACGCTGGCGGGTTGCCTGTGGGCGGTCACCGTGGATGCCGTCGGCAt                                                                                                                                              <  1:19196/150‑1 (MQ=255)
       gAGGTTGACGACCAGCGGCAAGTCATCGATATCCAGGCCACGGGCCGCCACATCGGTCGCCACCAGCACTTGCACCTCACGGGCCTTGAAGCTATCCAGCGCACGCTGGCGGGTTGCCTGTGGGCGGTCACCGTGGATGCCGTCGGCAtt                                                                                                                                             >  1:1222286/1‑150 (MQ=255)
           ttGACGACCAGCGGCAAGTCATCGATATCCAGGCCACGGGCCGCCACATCGGTCGCCACCAGCACTTGCACCTCACGGGCCTTGAAGCTGTCCAGCGCACGCTGGCGGGTTGCCTGTGGGCGGTCACCGTGAATGCCGACGGCATCCACg                                                                                                                                         <  1:186779/150‑1 (MQ=255)
           ttGACGACCAGCGGCAAGTCATCGATATCCAGGCCACGGGCCGCCACATCGGTCGCCACCAGCACTTGCACCGCACGGGCCTTGAAGCTGTCCAGCGCACGCTGGCGGGTTGCCTGTGGGCGGTCACCGTGAATGCCGACGGCATCCACg                                                                                                                                         <  1:93474/150‑1 (MQ=39)
             gacgacCAGCGGCAAGTCATCGATATCCAGGCCACGGGCCGCCACATCGGTCGCCACCAGCACTTGCACCTCACGGGCCTTGAAGCTGTCCAGCGCACGCTGGCGGGTTGCCTGTGGGCGGTCACCGTGGATGCCGTCGGCATTCACGcc                                                                                                                                       >  2:1068580/1‑150 (MQ=255)
             gacgacCAGCGGCAAGTCATCGATATCCAGGCCACGGGCCGCCACATCGGTCGCCACCAGCACTTGCACCTCACGGGCCTTGAAGCTGTCCAGCGCACGCTGGCGGGTTGCCTGTGGGCGGTCACCGTGGATGCCGTCGGCATTCACGcc                                                                                                                                       >  1:463967/1‑150 (MQ=255)
             gacgacCAGCGGCAAGTCATCGATATCCAGGCCACGGGCCGCCACATCGGTCGCCACCAGCACTTGCACCTCACGGGCCTTGAAGCTGTCCAGCGCACGCTGGCGGGTTGCCTGTGGGCGGTCACCGTGAATGCCGACGGCATCCACggt                                                                                                                                       <  1:605550/150‑3 (MQ=255)
                  ccAGCGGCAAGTCATCGATATCCAGGCCACGGGCCGCCACATCGGTCGCCACCAGCACTTGCACCTCACGGGCCTTGAAGCTGTCCAGCGCACGCTGGCGGGTTGCCTGTGGGCGGTCACCGTGGATGCCGTCGGCATTCACGCCTTCgg                                                                                                                                  >  2:616000/1‑150 (MQ=255)
                     gCGGCAAGTCATCGATATCCAGGCCACGGGCCGCCACATCGGTCGCCACCAGCACTTGCACCTCACGGGCCTTGAAGCTGTCCAGCGCACGCTGGCGGGTTGCCTGTGGGCGGTCACCGTGGATGCCGTCGGCATTCACGCCTTCGGCCa                                                                                                                               <  2:231617/150‑1 (MQ=255)
                        gaaaGTGATCGATAGGCAGGCCACGGGCCGCCACATCGGTCGCCACCAGCACTTGCACCTCACGGGCCTTGAAGCTGTCCAGCGCACGCTGGCGGGTTGCCTGTGGGCGGTCACCGTGGATGCCGTCGGCATTCACGCCTTCGGCCAGCa                                                                                                                            >  2:707449/3‑150 (MQ=255)
                             tCATCGATATCCAGGCCACGGGCCGCCACATCGGTCGCCACCAGCACTTGCACCTCACGGGCCTTGAAGCTGTCCAGCGCACGCTGGCGGGTTGCCTGTGGGCGGTCACCGTGGATGCCGTCGGCATTCACGCCTTCGGCCAGCAAgcgc                                                                                                                       <  2:1276514/150‑1 (MQ=255)
                             tCATCGATATCCAGGCCACGGGCCGCCACATCGGTCGCCACCAGCACTTGCACCTCACGGGCCTTGAAGCTGTCCAGCGCACGCTGGCGGGTTGCCTGTGGGCGGTCACCGTGGATGCCGTCGGCATTCACGCCTTCGGCCAGCAAgcgc                                                                                                                       <  1:257378/150‑1 (MQ=255)
                                tCGATATCCAGGCCACGGGCCGCCACATCGGTCGCCACCAGCACTTGCACCTCACGGGCCTTGAAGCTGTCCAGCGCACGCTGGCGGGTTGCCTGTGGGCGGTCACCGTGGATGCCGTCGGCATTCACGCCTTCGGCCAGCAAGCGCTcc                                                                                                                    >  2:826613/1‑150 (MQ=255)
                                      tCCAGGCCACGGGCCGCCACATCGGTCGCCACCAGCACTTGCACCTCACGGGCCTTGAAGCTGTCCAGCGCACGCTGGCGGGTTGCCTGTGGGCGGTCACCGTGGATGCCGTCGGCATTCACGCCTTCGGCCAGCAAGCGCTCCACCAGc                                                                                                              <  2:818009/150‑1 (MQ=255)
                                              aCGGGCCGCCACATCGGTCGCCACCAGCACTTGCACCTCACGGGCCTTGAAGCTGTCCAGCGCACGCTGGCGGGTTGCCTGTGGGCGGTCACCGTGGATGCCGTCGGCATTCACGCCTTCGGCCAGCAAGCGCTCCACCAGCTGATCGAc                                                                                                      <  2:33935/150‑1 (MQ=255)
                                              aCGGGCCGCCACATCGGTCGCCACCAGCACTTGCACCTCACGGGCCTTGAAGCTGTCCAGCGCACGCTGGCGGGTTGCCTGTGGGCGGTCACCGTGGATGCCGTCGGCATTCACGCCTTCGGCCAGCAAGCGCTCCACCAGCTGATCGAc                                                                                                      <  1:616000/150‑1 (MQ=255)
                                                      ccACATCGGTCGCCACCAGCACTTGCACCTCACGGGCCTTGAAGCTGTCCAGCGCACGCTGGCGGGTTGCCTGTGGGCGGTCACCGTGGATGCCGTCGGCATTCACGCCTTCGGCCAGCAAGCGCTCCACCAGCTGATCGACGCCATTGc                                                                                              >  1:223308/1‑150 (MQ=255)
                                                      ccACATCGGTCGCCACCAGCACTTGCACCTCACGGGCCTTGAAGCTGTCCAGCGCACGCTGGCGGGTTGCCTGTGGGCGGTCACCGTGGATGCCGTCGGCATTCACGCCTTCGGCCAGCAAGCGCTCCACCAGCTGATCGACGCCATTGc                                                                                              >  2:229146/1‑150 (MQ=255)
                                                       cacaTCGGTCGCCACCAGCACTTGCACCTCACGGGCCTTGAAGCTGTCCAGCGCACGCTGGCGGGTTGCCTGTGGGCGGTCACCGTGGATGCCGTCGGCATTCACGCCTTCGGCCAGCAAGCGCTCCACCAGCTGATCGACGCCATTGCg                                                                                             <  2:288542/150‑1 (MQ=255)
                                                       cacaTCGGTCGCCACCAGCACTTGCACCTCACGGGCCTTGAAGCTGACCAGCGCACGCTGGCGGGTTGCCTGTGGGCGGTCACCGTGGCTGCCGTCGGCATTCACGCCTTCGGCCAGCAAGCGCTCCACCAGCTGATCGACGCCATTGCg                                                                                             >  2:252690/1‑150 (MQ=255)
                                                       cacaTCGGTCACCACCAGCACTTGCACCTCACGGGCCTTGAAGCTGTCCAGCGCACGCTGGCGGGTTGCCTGTGGGCGGTCACCGTGGATGCCGTCGGCATTCACGCCTTCGGCCAGCAAGCGCTCCACCAGCTGATCGACGCCATTGCg                                                                                             <  1:668796/150‑1 (MQ=255)
                                                        acaTCGGTCGCCACCAGCACTTGCACCTCACGGGCCTTGAAGCTCTCCAGCGCACGCTGGCGGGTTGCCTGTGGGCGGTCACCGTGGATGCCGTCGGCATTCACGCCTTCGGCCAGCAAGCGCTCCACCAGCTGATCGACGCCATTGCgg                                                                                            <  2:1028659/150‑1 (MQ=255)
                                                         caTCGGTCGCCACCAGCACTTGCACCTCACGGGCCTTGAAGCTGTCCAGCGCACGCTGGCGGGTTGCCTGTGGGCGGTCACCGTGGATGCAATCGTCATTCCCGACTTCGGCCAGCAAGCGCTCCCCCAGCTGACCGACGCCATTGCcgg                                                                                           <  2:760477/150‑1 (MQ=255)
                                                          aTCGGTCGCCACCAGCACTTGCACCTCACGGGCCTTGAAGCTGTCCAGCGCACGCTGGCGGGTTGCCTGTGGGCGGTCACCGTGGATGCCGTCGGCATTCACGCCTTCGGCCAGCAAGCGCTCCACCAGCTGATCGACGCCATTGCGGGt                                                                                          <  1:229146/150‑1 (MQ=255)
                                                              gTCGCCACCAGCACTTGCACCTCACGGGCCTTGAAGCTGTCCAGCGCACGCTGGCGGGTTGCCTGTGGGCGGTCACCGTGGATGCCGTCGGCATTCACGCCTTCGGCCAGCAAGCGCTCCACCAGCTGATCGACGCCATTGCGGGTCTTg                                                                                      <  1:1068580/150‑1 (MQ=255)
                                                               tCGCCACCAGCACTTGCACCTCACGGGCCTTGAAGCTGTCCAGCGCACGCTGGCGGGTTGCCTGTGGGCGGTCACCGTGGATGCCGTCGGCATTCACGCCTTCGGCCAGCAAGCGCTCCACCAGCTGATCGACGCCATTGCGGGTCTTgg                                                                                     >  2:284150/1‑150 (MQ=255)
                                                                   caccaGCACTTGCACCTCACGGGCCTTGAAGCTATCCAGCGCACGCTGGCGGGTTGCCTGTGGGCGGTCACCGTGGATGCCGTCGGCATTCACGCCTTCGGCCAGCAAGCGCTCCACCAGCTGATCGACGCCATTGCGGGTCTTGGCGaa                                                                                 <  2:1222286/150‑1 (MQ=255)
                                                                    accaGCACTTGCACCTCACGGGCCTTGAAGCTGTCCAGCGCACGCTGGCGGGTTGCCTGTGGGCGGTCACCGTGAATGCCGACGGCATCCACggt                                                                                                                                       <  2:46549/95‑3 (MQ=38)
                                                                    accaGCACTTGCACCTCACGGGCCTTGAAGCTGTCCAGCGCACGCTGGCGGGTTGCCTGTGGGCGGTCACCGTGAATGCCGACGGCATCCACggt                                                                                                                                       >  1:46549/1‑93 (MQ=39)
                                                                         cACTTGCACCTCACGGGCCTTGAAGCTGTCCAGCGCACGCTGGCGGGTTGCCTGTGGGCGGTCACCGTGGATGCCGTCGGCATTCACGCCTTCGGCCAGCAAGCGCTCCACCAGCTGATCGACGCCATTGCGGGTCTTGGCGAACACCAg                                                                           <  2:722958/150‑1 (MQ=255)
                                                                          aCTTGCACCTCACGGGCCTTGAAGCTGTCCAGCGCACGCTGGCGGGTTGCCTGTGGGCGGTCACCGTGGATGCCGTCGGCATTCACGCCTTCGGCCAGCAAGCGCTCCACCAGCTGATCGACGCCATTGCGGGTCTTGGCGAACACCAGc                                                                          >  2:844953/1‑150 (MQ=255)
                                                                           cTTGCACCTCACGGGCCTTGAAGCTGTCCAGCGCACGCTGGCGGGTTGCCTGTGGGCGGTCACCGTGGATGCCGTCGGCATTCACGCCTTCGGCCAGCAAGCGCTCCACCAGCTGATCGACGCCATTGCGGGTCTTGGCGAACACCAGc                                                                          <  2:993612/149‑1 (MQ=255)
                                                                           cTTGCACCTCACGGGCCTTGAAGCTGTCCAGCGCACGCTGGCGGGTTGCCTGTGGGCGGTCACCGTGGATGCCGTCGGCATTCACGCCTTCGGCCAGCAAGCGCTCCACCAGCTGATCGACGCCATTGCGGGTCTTGGCGAACACCAGc                                                                          >  1:993612/1‑149 (MQ=255)
                                                                            ttGCACCTCACGGGCCTTGAAGCTGTCCAGCGCACGCTGGCGGGTTGCCTGTGGGCGGTCACCGTGGATGCCGTCGGCATTCACGCCTTCGGCCAGCAAGCGCTCCACCAGCTGATCGACGCCATTGCGGGTCTTGGCGAACACCAGCAc                                                                        <  1:252690/150‑1 (MQ=255)
                                                                                  cTCACGGGCCTTGAAGCTGTCCAGCGCACGCTGGCGGGTTGCCTGTGGGCGGTCACCGTGGATGCCGTCGGCATTCACGCCTTCGGCCAGCAAGCGCTCCACCAGCTGATCGACGCCATTGCGGGTCTTGGCGAACACCAGCACCTGCtt                                                                  >  2:1295171/1‑150 (MQ=255)
                                                                                     aCGGGCCTTGAAGCTGTCCAGCGCACGCTGGCGGGTTGCCTGTGGGCGGTCACCGTGGATGCCGTCGGCATTCACGCCTTCGGCCa                                                                                                                               <  1:1275279/86‑1 (MQ=255)
                                                                                     aCGGGCCTTGAAGCTGTCCAGCGCACGCTGGCGGGTTGCCTGTGGGCGGTCACCGTGGATGCCGTCGGCATTCACGCCTTCGGCCa                                                                                                                               >  2:1275279/1‑86 (MQ=255)
                                                                                     aCGGGCCTTGAAGCTGTCCAGCGCACGCTGGCGGGTTGCCTGTGGGCGGTCACCGTGGATGCCGTCGGCATTCACGCCTTCGGCCAGCAAGCGCTCCACCAGCTGATACAGCCCATTGCGGGTCTTGGCGAACACCAGCACCTGCTTCCa                                                               >  1:727650/1‑150 (MQ=255)
                                                                                      cGGGCCTTGAAGCTGTCCAGCGCACGCTGGCGGGTTGCCTGTGGGCGGTCACCGTGGATGCCGTCGGCATTCACGCCTTCGGCCAGCAAGCGCTCCACCAGCTGATCGACGCCATTGCGGGTCTTGGCGAACACCAGCACCTGCTTCCAg                                                              >  2:1363976/1‑150 (MQ=255)
                                                                                           cTTGAAGCTGTCCAGCGCACGCTGGCGGGTTGCCTGTGGGCGGTCACCGTGGATGCCGTCGGCATTCACGCCTTCGGCCAGCAAGCGCTCCACCAGCTGATCGACGCCATTGCGGGTCTTGGCGAACACCAGCACCTGc                                                                    >  2:829217/1‑139 (MQ=255)
                                                                                           cTTGAAGCTGTCCAGCGCACGCTGGCGGGTTGCCTGTGGGCGGTCACCGTGGATGCCGTCGGCATTCACGCCTTCGGCCAGCAAGCGCTCCACCAGCTGATCGACGCCATTGCGGGTCTTGGCGAACACCAGCACCTGc                                                                    <  1:829217/139‑1 (MQ=255)
                                                                                           cTTGAAGCTGTCCAGCGCACGCTGGCGGGTTGCCTGTGGGCGGTCACCGTGGATGCCGTCGGCATTCACGCCTTCGGCCAGCAAGCGCTCCACCAGCTGATCGACGCCATTGCGGGTCTTGGCGAACACCAGCACCTGCTTCCAGCgctg                                                         >  1:686311/1‑150 (MQ=255)
                                                                                            ttGAAGCTGTCCAGCGCACGCTGGCGGGTTGCCTGTGGGCGGTCACCGTGGATGCCGTCGGCATTCACGCCTTCGGCCAGCAAGCGCTCCACCAGCTGATCGAc                                                                                                      >  2:388596/1‑104 (MQ=255)
                                                                                            ttGAAGCTGTCCAGCGCACGCTGGCGGGTTGCCTGTGGGCGGTAACCGTGGATGCCGTCGGCATTCACGCCTTCGGCCAGCAAGCGCTCCACCAGCTGATCGAc                                                                                                      <  1:388596/104‑1 (MQ=255)
                                                                                                     tCCAGCGCACGCTGGCGGGTTGCCTGTGGGCGGTCACCGTGGATGCCGTCGGCATTCACGCCTTCGGCCAGCAAGCGCTCCACCAGCTGATCGACGCCATTGCGGGTCTTGGCGAACACCAGCACCTGCTTCCAGCGCTGCTTGCgcagc                                               >  2:860129/1‑150 (MQ=255)
                                                                                                      ccAGCGCACGCTGGCGGGTTGCCTGTGGGCGGTCACCGTGGATGCCGTCGGCATTCACGCCTTCGGCCAGCAAGCGCTCCACCTGCTGATCGACGCCATTGCGGGTCTTGCCGAACACCAGCGCCTGCTTCCAGCGCTGCTTGCgcagca                                              <  1:763759/150‑1 (MQ=255)
                                                                                                       cAGCGCACGCTGGCGGGTTGCCTGTGGGCGGTCACCGTGGATGCCGTCGGCATTCACGCCTTCGGCCAGCAAGCGCTCCACCAGc                                                                                                              >  2:772347/1‑85 (MQ=255)
                                                                                                       cAGCGCACGCTGGCGGGTTGCCTGTGGGCGGTCACCGTGGATGCCGTCGGCATTCACGCCTTCGGCCAGCAAGCGCTCCACCAGc                                                                                                              <  1:772347/85‑1 (MQ=255)
                                                                                                       cAGCGCACGCTGGCGGGTTGCCTGTGGGCGGTCACCGTGGATGCCGTCGGCATTCACGCCTTCGGCCAGCAAGCGCTCCACCAGCTGATCGACGCCATTGCGGGTCTTGGCGaa                                                                                 <  2:1120075/114‑1 (MQ=255)
                                                                                                       cAGCGCACGCTGGCGGGTTGCCTGTGGGCGGTCACCGTGGATGCCGTCGGCATTCACGCCTTCGGCCAGCAAGCGCTCCACCAGCTGATCGACGCCATTGCGGGTCTTGGCGaa                                                                                 >  1:1120075/1‑114 (MQ=255)
                                                                                                           gcACGCTGGCGGGTTGCCTGTGGGCGGTCACCGTGGATGCCGTCGGCATTCACGCCTTCGGCCAGCAAGCGCTCCACCAGCTGATCGACGCCATTGCGGGTCTTGGCGAACACCAGCACCTGCTTCCAGCGCTGCTTGCGCAGCAGGTgg                                         <  2:391950/150‑1 (MQ=255)
                                                                                                                cTGGCGGGTTGCCTGTGGGCGGTCACCGTGGATGCCGTCGGCATTCACGCCTTCGGCCAGCAAGCGCTCCACCAGCTGATCGACGCCATTGCGGGTCTTGGCGAACACCAGCACCTGCTTCCAGCGCTGCTTGCGCAGCAGGTGGCAGaa                                    <  1:826613/150‑1 (MQ=255)
                                                                                                                   gCGGGTTGCCTGTGGGCGGTCACCGTGGATGCCGTCGGCATTCACGCCTTCGGCCAGCAAGCGCTCCACCAGCTGATCGACGCCATTGCGGGTCTTGGCGAACACCAGCACCTGCTTCCAGCGCTGCTTGCGCAGCAGGTGGCAGAACAg                                 <  1:707449/150‑1 (MQ=255)
                                                                                                                    cGGGTTGCCTGTGGGCGGTCACCGTGGATGCCGTCGGCATTCACGCCTTCGGCCAGCAAGCGCTCCACCAGCTGATCGACGCCATTGCGGGTCTTGGCGAACACCAGCACCTGCTTCCAGCGCTGCTTGCGCAGCAGGTGGCAGAACAgg                                <  2:223308/150‑1 (MQ=255)
                                                                                                                    cGGGTTGCCTGTGGGCGGTCACCGTGGATGCCGTCGGCATTCACGCCTTCGGCCAGCAAGCGCTCCACCAGCTGATCGACGCCATTGCGGGTCTTGGCGAACACCAGCACCTGCTTCCAGCGCTGCTTGCGCAGCAGGTGGCAGAACAgg                                <  1:1295171/150‑1 (MQ=255)
                                                                                                                       gTTGCCTGTGGGCGGTCACCGTGGATGCCGTCGGCATTCACGCCTTCGGCCAGCAAGCGCTCCACCAGCTGATCGACGCCATTGCgg                                                                                            >  1:945016/1‑87 (MQ=255)
                                                                                                                       gTTGCCTGTGGGCGGTCACCGTGGATGCCGTCGGCATTCACGCCTTCGGCCAGCAAGCGCTCCACCAGCTGATCGACGCCATTGCGGGTCTTGGCGAACACCAGCACCTGCTTCCAGCGCTGCTTGCGCAGCAGGTGGCAGAACAGGTCa                             <  1:284150/150‑1 (MQ=255)
                                                                                                                       gTTGCCTGTGGGCGGTCACCGTGGATGCCGTCGGCACTCACGCCTTCGGCCAGCAAGCGCTCCACCAGCTGATCGACGCCATTGCgg                                                                                            <  2:945016/87‑1 (MQ=255)
                                                                                                                            cTGTGGGCGGTCACCGTGGATGCCGTCGGCATTCACGCCTTCGGCCAGCAAGCGCTCCACCAGCTGATCGACGCCATTGCGGGTCTTGGCGAACACCAGCACCTGCTTCCAGCGCTGCTTGCGCAGCAGGTGGCAGAACAGGTCAGCCtt                        >  2:962078/1‑150 (MQ=255)
                                                                                                                                    cgaaggcGTGAATGCCGACGGCATCCACGCCTTCGGCCAGCAAGCGCTCCACCAGCTGATCGACGCCATTGCGGGTCTTGGCGAACACCAGCACCTGCTTCCAGCGCTGCTTGCGCAGCAGGTGGCAGAACAGGTCAGCCTTGCGcttc                 >  1:1029471/7‑149 (MQ=255)
                                                                                                                                    cgaaggcGTGAATGCCGACGGCATCCACGCCTTCGGCCAGCAAGCGCTCCACCAGCTGATCGACGCCATTGCGGGTCTTGGCGAACACCAGCACCTGCTTCCAGCGCTGCTTGCGCAGCAGGTGGCAGAACAGGTCAGCCTTGCGcttc                 <  2:1029471/143‑1 (MQ=255)
                                                                                                                                     gTCACCGTGGATGCCGTCGGCATTCACGCCTTCGGCCAGCAAGCGCTCCACCAGCTGATCGACGCCATTGCGGGTCTTGGCGAACACCAGCACCTGCTTCCAGCGCTGCTTGCGCAGCAGGTGGCAGAACAGGTCAGCCTTGCGcttctt               <  2:842292/150‑1 (MQ=255)
                                                                                                                                     gTCACCGTGGATGCCGTCGGCATTCACGCCTTCGGCCAGCAAGCGCTCCACCAGCTGATCGACGCCATTGCGGGTCTTGGCGAACACCAGCACCTGCTTCCAGCGCTGCTTGCGCAGCAGGTGGCAGAACAGGTCAGCCTTGCGcttctt               <  1:844953/150‑1 (MQ=255)
                                                                                                                                      tCACCGTGGATGCCGTCGGCATTCACGCCTTCGGCCAGCAAGCGCTCCACCAGCTGATCGACGCCATTGCGGGTCTTGGCGAACACCAGCACCTGCTTCCAGCGCTGCTTGCGCAGCAGGTGGCAGAACAGGTCAGCCTTGCGCTTCTTg              >  1:1393838/1‑150 (MQ=255)
                                                                                                                                         ccGTGGATGCCGTCGGCATTCACGCCTTCGGCCAGCAAGCGCTCCACCAGCTGATCGACGCCATTGCGGGTCTTGGCGAACACCAGCACCTGCTTCCAGCGCTGCTTGCGCAGCAGGTGGc                                        >  1:1178573/1‑121 (MQ=255)
                                                                                                                                         ccGTGGATGCCGTCGGCATTCACGCCTTCGGCCAGCAAGCGCTCCACCAGCTGATCGACGCCATTGCGGGTCTTGGCGAACACCAGCACCTGCTTCCAGCGCTGCTTGCGCAGCAGGTGGc                                        <  2:1178573/121‑1 (MQ=255)
                                                                                                                                           gTGGATGCCGTCGGCATTCACGCCTTCGGCCAGCAAGCGCTCCACCAGCTGATCGACGCCAtt                                                                                                >  1:176363/1‑63 (MQ=255)
                                                                                                                                           gTGGATGCCGTCGGCATTCACGCCTTCGGCCAGCAAGCGCTCCACCAGCTGATCGACGCCAtt                                                                                                <  2:176363/63‑1 (MQ=255)
                                                                                                                                           gTGAATGCCGACGGCATCCACGCCTTCGGCCAGCAAGCGCTCCACCAGCTGATCGACGCCCTTGCGGGTTTTTGCGAACACCAGCACCTGCTTCCATCGCTGCTTGCGCAGCAGGTGGGAGAACAGGTCAGCCTTGCGCTTCTTGTCcac         >  2:979751/1‑150 (MQ=255)
                                                                                                                                           gTGAATGCCGACGGCATCCACGCCTTCGGCCAGCAAGCGCTCCACCAGCTGATCGACGCCATTGCGGGTCTTGGCGAACACCAGCACCTGCTTCCAGCGCTGCTTGCGCAGCAGGTGGCAGAACAGGTCAGCCTTGCGCTTCTTGTCcac         >  1:112499/1‑150 (MQ=255)
                                                                                                                                           gTGAATGCCGACGGCATCCACGCCTTCGGCCAGCAAGCGCTCCACCAGCTGATCGACGCCATTGCGGGTCTTGGCGAACACCAGCACCTGCTTCCAGCGCTGCTTGCGCAGCAGGTGGCAGAACAGGTCAGCCTTGCGCTTCTTGTCcac         >  1:915315/1‑150 (MQ=255)
                                                                                                                                               aTGCCGTCGGCATTCACGCCTTCGGCCAGCAAGCGCTCCACCAGCTGATCGACGCCATTGCGGGTCTTGGCGAACACCAGCACCTGCTTCCAGCGCTTCTTGCGCAGCAGGAGCCAGAACAGGTCAGCCTTGCGCTTCTTGTCCACAGGc     >  2:691266/1‑150 (MQ=255)
                                                                                                                                                tGCCGTCGGCATTCACGCCTTCGGCCAGCAAGCGCTCCACCAGCTGATCGACGCCATTGCGGGTCTTGGCGAACACCAGCACCTGCTTc                                                                 <  1:143902/89‑1 (MQ=255)
                                                                                                                                                tGCCGTCGGCATTCACGCCTTCGGCCAGCAAGCGCTCCACCAGCTGATCGACGCCATTGCGGGTCTTGGCGAACACCAGCACCTGCTTc                                                                 >  2:143902/1‑89 (MQ=255)
                                                                                                                                                  ccGTCGGCATTCACGCCTTCGGCCAGCAAGCGCTCCCCCAGCTGATCGACGCCATTGCGGGTCTTGGCGAACACCAGCACCTGCTTCCAGCGCTGCTTGCGCAGCAGGTGGCCGAACAGGTCAGGCTTGCGCTTCTTGTCCACAGGCAcc  <  1:691266/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                  
GCAGGTCGAGGTTGACGACCAGCGGCAAGTCATCGATATCCAGGCCACGGGCCGCCACATCGGTCGCCACCAGCACTTGCACCTCACGGGCCTTGAAGCTGTCCAGCGCACGCTGGCGGGTTGCCTGTGGGCGGTCACCGTGGATGCCGTCGGCATTCACGCCTTCGGCCAGCAAGCGCTCCACCAGCTGATCGACGCCATTGCGGGTCTTGGCGAACACCAGCACCTGCTTCCAGCGCTGCTTGCGCAGCAGGTGGCAGAACAGGTCAGCCTTGCGCTTCTTGTCCACAGGCACC  >  NC_002947/5427187‑5427482

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