Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,427,343 T→C 13.7% N270D (AAT→GAT)  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,3430TC13.7% 184.8 / 16.4 73N270D (AAT→GAT) rhlE‑IATP‑dependent RNA helicase
Reads supporting (aligned to +/- strand):  ref base T (28/35);  new base C (5/5);  total (33/40)
Fisher's exact test for biased strand distribution p-value = 7.47e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

GAGGTTGACGACCAGCGGCAAGTCATCGATATCCAGGCCACGGGCCGCCACATCGGTCGCCACCAGCACTTGCACCTCACGGGCCTTGAAGCTGTCCAGCGCACGCTGGCGGGTTGCCTGTGGGCGGTCACCGTGGATGCCGTCGGCATTCACGCCTTCGGCCAGCAAGCGCTCCACCAGCTGATCGACGCCATTGCGGGTCTTGGCGAACACCAGCACCTGCTTCCAGCGCTGCTTGCGCAGCAGGTGGCAGAACAGGTCAGCCTTGCGCTTCTTGTCCACAGGCACCAGCCACTGC  >  NC_002947/5427194‑5427491
                                                                                                                                                     |                                                                                                                                                    
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                                                                                                                                <  1:257378/150‑1 (MQ=255)
                      tCATCGATATCCAGGCCACGGGCCGCCACATCGGTCGCCACCAGCACTTGCACCTCACGGGCCTTGAAGCTGTCCAGCGCACGCTGGCGGGTTGCCTGTGGGCGGTCACCGTGGATGCCGTCGGCATTCACGCCTTCGGCCAGCAAgcgc                                                                                                                                <  2:1276514/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                                                                             <  1:829217/139‑1 (MQ=255)
                                                                                    cTTGAAGCTGTCCAGCGCACGCTGGCGGGTTGCCTGTGGGCGGTCACCGTGGATGCCGTCGGCATTCACGCCTTCGGCCAGCAAGCGCTCCACCAGCTGATCGACGCCATTGCGGGTCTTGGCGAACACCAGCACCTGc                                                                             >  2:829217/1‑139 (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                                                                                                                       <  1:772347/85‑1 (MQ=255)
                                                                                                cAGCGCACGCTGGCGGGTTGCCTGTGGGCGGTCACCGTGGATGCCGTCGGCATTCACGCCTTCGGCCAGCAAGCGCTCCACCAGc                                                                                                                       >  2:772347/1‑85 (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                        <  1:844953/150‑1 (MQ=255)
                                                                                                                              gTCACCGTGGATGCCGTCGGCATTCACGCCTTCGGCCAGCAAGCGCTCCACCAGCTGATCGACGCCATTGCGGGTCTTGGCGAACACCAGCACCTGCTTCCAGCGCTGCTTGCGCAGCAGGTGGCAGAACAGGTCAGCCTTGCGcttctt                        <  2:842292/150‑1 (MQ=255)
                                                                                                                               tCACCGTGGATGCCGTCGGCATTCACGCCTTCGGCCAGCAAGCGCTCCACCAGCTGATCGACGCCATTGCGGGTCTTGGCGAACACCAGCACCTGCTTCCAGCGCTGCTTGCGCAGCAGGTGGCAGAACAGGTCAGCCTTGCGCTTCTTg                       >  1:1393838/1‑150 (MQ=255)
                                                                                                                                  ccGTGGATGCCGTCGGCATTCACGCCTTCGGCCAGCAAGCGCTCCACCAGCTGATCGACGCCATTGCGGGTCTTGGCGAACACCAGCACCTGCTTCCAGCGCTGCTTGCGCAGCAGGTGGc                                                 <  2:1178573/121‑1 (MQ=255)
                                                                                                                                  ccGTGGATGCCGTCGGCATTCACGCCTTCGGCCAGCAAGCGCTCCACCAGCTGATCGACGCCATTGCGGGTCTTGGCGAACACCAGCACCTGCTTCCAGCGCTGCTTGCGCAGCAGGTGGc                                                 >  1:1178573/1‑121 (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                                                                          >  2:143902/1‑89 (MQ=255)
                                                                                                                                         tGCCGTCGGCATTCACGCCTTCGGCCAGCAAGCGCTCCACCAGCTGATCGACGCCATTGCGGGTCTTGGCGAACACCAGCACCTGCTTc                                                                          <  1:143902/89‑1 (MQ=255)
                                                                                                                                           ccGTCGGCATTCACGCCTTCGGCCAGCAAGCGCTCCCCCAGCTGATCGACGCCATTGCGGGTCTTGGCGAACACCAGCACCTGCTTCCAGCGCTGCTTGCGCAGCAGGTGGCCGAACAGGTCAGGCTTGCGCTTCTTGTCCACAGGCAcc           <  1:691266/150‑1 (MQ=255)
                                                                                                                                                    ttCACGCCTTCGGCCAGCAAGCGCTCCACCAGCTGATCGACGCCATTGCGGGTCTTGGCGAACACCAGCACCTGCTTCCAGCGCTGCTTGCGCAGCAGGTGGCAGAACAGGTCAGCCTTGCGcttc                          >  2:1889/1‑126 (MQ=255)
                                                                                                                                                    ttCACGCCTTCGGCCAGCAAGCGCTCCACCAGCTGATCGACGCCATTGCGGGTCTTGGCGAACACCAGCACCTGCTTCCAGCGCTGCTTGCGCAGCAGGTGGCAGAACAGGTCAGCCTTGCGcttc                          <  1:1889/126‑1 (MQ=255)
                                                                                                                                                    ttCACGCCTTCGGCCAGCAAGCGCTCCACCAGCTGATCGACGCCATTGCGGGTCTTGGCGAACACCAGCACCTGCTTCCAGCGCTGCTTGCGCAGCAGGTGGCAGAACAGGTCAGCCTTGCGCTTCTTGTCCACAGGCACCAGCCACTGc  >  2:250960/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
GAGGTTGACGACCAGCGGCAAGTCATCGATATCCAGGCCACGGGCCGCCACATCGGTCGCCACCAGCACTTGCACCTCACGGGCCTTGAAGCTGTCCAGCGCACGCTGGCGGGTTGCCTGTGGGCGGTCACCGTGGATGCCGTCGGCATTCACGCCTTCGGCCAGCAAGCGCTCCACCAGCTGATCGACGCCATTGCGGGTCTTGGCGAACACCAGCACCTGCTTCCAGCGCTGCTTGCGCAGCAGGTGGCAGAACAGGTCAGCCTTGCGCTTCTTGTCCACAGGCACCAGCCACTGC  >  NC_002947/5427194‑5427491

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: