Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,427,343 T→C 6.6% 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,3430TC 6.6% 452.2 / 14.3 152N270D (AAT→GAT) rhlE‑IATP‑dependent RNA helicase
Reads supporting (aligned to +/- strand):  ref base T (62/79);  new base C (3/7);  total (65/87)
Fisher's exact test for biased strand distribution p-value = 5.16e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.53e-01

Only 100 of 156 total aligned reads displayed.

AGGTTGACGACCAGCGGCAAGTCATCGATATCCAGGCCACGGGCCGCCACATCGGTCGCCACCAGCACTTGCACCTCACGGGCCTTGAAGCTGTCCAGCGCACGCTGGCGGGTTGCCTGTGGGCGGTCACCGTGGATGCCGTCGGCATTCACGCCTTCGGCCAGCAAGCGCTCCACCAGCTGATCGACGCCATTGCGGGTCTTGGCGAACACCAGCACCTGCTTCCAGCGCTGCTTGCGCAGCAGGTGGCAGAACAGGTCAGCCTTGCGCTTCTTGTCCACAGGCACCAGCCACTGC  >  NC_002947/5427195‑5427491
                                                                                                                                                    |                                                                                                                                                    
  gTTGACGACCAGCGGCAAGTCATCGATATCCAGGCCACGGGCCGCCACATCGGTCGCCACCAGCACTTGCACCTCACGGGCCTTGAAGCTGTCCAGCGCACGCTGGCGGGTTGCCTGTGGGCGGTCACCGTGGATGCCGTCGGCATTCAc                                                                                                                                                   >  2:2347143/1‑150 (MQ=255)
   ttGACGACCAGCGGCAAGTCATCGATATCCAGGCCACGGGCCGCCACTGCGGTCGCCACCAGCACTTGCACCTCACGGGCCTTGAAGCTGTCCAGCGCACGCTGGCGGGTTGCCTGTGGGCGGTCACCGTGAATGCCGACGGCATCCACg                                                                                                                                                  <  1:872010/150‑1 (MQ=255)
   ttGACGACCAGCGGCAAGTCATCGATATCCAGGCCACGGGCCGCCACATCGGTCGCCACCAGCACTTGCACCTCACGGGCCTTGAAGCTGTCCAGCGCACGCTGGCGGGTTGCCTGTGGGCGGTCACCGTGAATGCCGACGGCATCCACg                                                                                                                                                  <  1:1830462/150‑1 (MQ=255)
         acCAGCGGCAAGTCATCGATATCCAGGCCACGGGCCGCCACATCGGTCGCCACCAGCACTTGCACCTCACGGGCCTTGAAGCTGTCCAGCGCACGCTGGCGGGTTGCCTGTGGGCGGTCACCGTGGATGCCGTCGGCATTCACGCCTTCg                                                                                                                                            <  1:2310776/150‑1 (MQ=255)
         acCAGCGGCAAGTCATCGATATCCAGGCCACGGGCCGCCACATCGGTCGCCACCAGCACTTGCACCTCACGGGCCTTGAAGCTGTCCAGCGCACGCTGGCGGGTTGCCTGTGGGCGGTCACCGTGAATGCCGACGGCATCCACggtgacc                                                                                                                                            <  1:1049868/150‑7 (MQ=39)
           cAGCGGCAAGTCATCGATATCCAGGCCACGGGCCGCCACATCGGTCGCCACCAGCACTTGCACCTCACGGGCCTTGAAGCTGTCCAGCGCACGCTGGCGGGTTGCCTGTGGGCGGTCACCGTGAATGCCGACGGCATCCACgg                                                                                                                                                 >  1:318917/1‑142 (MQ=255)
           cAGCGGCAAGTCATCGATATCCAGGCCACGGGCCGCCACATCGGTCGCCACCAGCACTTGCACCTCACGGGCCTTGAAGCTGTCCAGCGCACGCTGGCGGGTTGCCTGTGGGCGGTCACCGTGAATGCCGACGGCATCCACgg                                                                                                                                                 <  2:318917/143‑2 (MQ=255)
           cAGCGGCAAGTCATCGATATCCAGGCCACGGGCCGCCACATCGGTCGCCACCAGCACTTGCACCTCACGGGCCTTGAAGCTGTCCAGCGCACGCTGGCGGGTTGCCTGTGGGCGGTCAACGTGGATGCCGTCGGCATTCACGCCGTCGGc                                                                                                                                          >  1:38096/1‑150 (MQ=255)
             gCGGCAAGTCATCGATATCCAGGCCACGGGCCGCCACATCGGTCGCCACCAGCACTTGCACCTCACGGGCCTTGAAGCTGTCCAGCGCACGCTGGCGGGTTGCCTGTGGGCGGTCACCGTGGATGCCGTCGGCATTCACGCCTTCGGCCa                                                                                                                                        >  1:2187793/1‑150 (MQ=255)
             gCGGCAAGTCATCGATATCCAGGCCACGGGCCGCCACATCGGTCGCCACCAGCACTTGCACCTCACGGGCCTTGAAGCTGTCCAGCGCACGCTGGCGGGTTGCCTGTGGGCGGTCACCGTGGATGCCGACGGCATTCACGCCTTCGGCCa                                                                                                                                        <  1:2022700/150‑1 (MQ=255)
                       aTCGATATCCAGGCCACGGGCCGCCACATCGGTCGCCACCAGCACTTGCACCTCACGGGCCTTGAAGCTGTCCAGCGCACGCTGGCGGGTTGCCTGTGGGCGGTCACCGTGGATGCCGTCGGCATTCACGCCTTCGGCCAGCAAGCGCTc                                                                                                                              <  1:778578/150‑1 (MQ=255)
                        tCGATATCCAGGCCACGGGCCGCCACATCGGTCGCCACCAGCACTTGCACCTCACGGGCCTTGAAGCTGTCCAGCGCACGCTGGCGGGTTGCCTGTGGGCGGTCACCGTGGATGCCGTCGGCATTCACGCCTTCGGCCAGCAAGCGCTcc                                                                                                                             <  1:1080885/150‑1 (MQ=255)
                           atatCCAGGCCACGGGCCGCCACATCGGTCGCCACCAGCACTTGCACCTCACGGGCCTTGAAGCTGTCCAGCGCACGCTGGCGGGTTGCCTGTGGGCGGTCACCGTGGATGCCGTCGGCATTCACGCCTTCGGCCAg                                                                                                                                       >  1:1758954/1‑137 (MQ=255)
                           atatCCAGGCCACGGGCCGCCACATCGGTCGCCACCAGCACTTGCACCTCACGGGCCTTGAAGCTGTCCAGCGCACGCTGGCGGGTTGCCTGTGGGCGGTCACCGTGGATGCCGTCGGCATTCACGCCTTCGGCCAg                                                                                                                                       <  2:1758954/137‑1 (MQ=255)
                             atCCAGGCCACGGGCCGCCACATCGGTCGCCACCAGCACTTGCACCTCACGGGCCTTGAAGCTGTCCAGCGCACGCTGGCGGGTTGCCTGTGGGCGGTCACCGTGGATGCCGTCGGCATTCACGCCTTCGGCCAGCAAGCGCTCCACCAg                                                                                                                        <  2:131050/150‑1 (MQ=255)
                                cAGGCCACGGGCCGCCACATCGGTCGCCACCAGCACTTGCACCTCACGGGCCTTGAAGCTGTCCAGCGCACGCTGGCGGGTTGCCTGTGGGCGGTCACCGTGGATGCCGTCGGCATTCACGc                                                                                                                                                 <  1:1256460/122‑1 (MQ=255)
                                cAGGCCACGGGCCGCCACATCGGTCGCCACCAGCACTTGCACCTCACGGGCCTTGAAGCTGTCCAGCGCACGCTGGCGGGTTGCCTGTGGGCGGTCACCGTGGATGCCGTCGGCATTCACGc                                                                                                                                                 <  1:1255456/122‑1 (MQ=255)
                                cAGGCCACGGGCCGCCACATCGGTCGCCACCAGCACTTGCACCTCACGGGCCTTGAAGCTGTCCAGCGCACGCTGGCGGGTTGCCTGTGGGCGGTCACCGTGGATGCCGTCGGCATTCACGc                                                                                                                                                 >  2:1255456/1‑122 (MQ=255)
                                cAGGCCACGGGCCGCCACATCGGTCGCCACCAGCACTTGCACCTCACGGGCCTTGAAGCTGTCCAGCGCACGCTGGCGGGTTGCCTGTGGGCGGTCACCGTGGATGCCGTCGGCATTCACGCCTTCGGCCAGCAAGCGCTCCACCAGCTg                                                                                                                     >  2:2024095/1‑150 (MQ=255)
                                cAGGCCACGGGCCGCCACATCGGTCGCCACCAGCACTTGCACCTCACGGGCCTTGAAGCTGTCCAGCGCACGCTGGCGGGTTGCCTGTGGGCGGTCACCGTGGATGCCGTCGGCATTCACGCCTTCGGCCAGCAAGCGCTCCACCAGCTg                                                                                                                     >  2:612064/1‑150 (MQ=255)
                                 aGGCCACGGGCCGCCACATCGGTCGCCACCAGACCCCGCACCTCACGGGCCTTGAAGCTGTCCAGCGCACGCTAGCGGGTTGCCTGTGGGCGGTCACCGTGGATGCCGTCGGCATTCACGCCTTCGGCCAGCAAGCGCTCCACCAGCTGa                                                                                                                    >  1:1172216/1‑150 (MQ=255)
                                   gCCACGGGCCGCCACATCGGTCGCCACCAGCACTTGCACCTCACGGGCCTTGAAGCTGTCCAGCGCACGCTGGCGGGTTGCCTGTGGGCGGTCACCGTGGATGCCGTCGGCATTCACGCCTTCGGCCAGCAAGCGCTCCACCAGCTGATc                                                                                                                  <  2:458989/150‑1 (MQ=255)
                                   cccACTGGCCGCCACATCGGTCGCCACCAGCACTTGCACCTCACGGGCCTTCAAGCTGTCCAGCGCACGCTGGCGGGTTGCCTGCGGGCGGTCACCGTGGATGCCGTCGGCATTCACGCCTTCGGCCAGCAAGCGCTCCACCAGCTGATc                                                                                                                  <  1:255346/149‑1 (MQ=255)
                                     cACGGGCCGCCACATCGGTCGCCACCAGCACTTGCACCTCACGGGCCTTGAAGATGTTGAACGCACGCTTGCGTGGTGCCTTTGGGCGGTCACCGTGGATGCCGTCGGCATTCACGCCTTCGGCCAGCAAGCGCTCCACCAGCTGATCGa                                                                                                                >  1:2100580/1‑150 (MQ=255)
                                      aCGGGCCGCCACATCGGTCGCCACCAGCACTTGCACCTCACGGGCCTTGAAGCTGTCCAGCGCACGCTGGCGGGTTGCCTGTGGGCGGTCACCGTGGATGCCGTCGGCATTCACGCCTTCGGCCAGCAAGCGCTCCACCAGCTGATCGAc                                                                                                               >  2:1918793/1‑150 (MQ=255)
                                              ccACATCGGTCGCCACCAGCACTTGCACCTCACGGGCCTTGAAGCTGTCCAGCGCACGCTGGCGGGTTGCCTGTGGGCGGTCACCGTGGATGCCGTCGGCATTCACGCCTTCGGCCAGCAAGCGCTCCACCAGCt                                                                                                                      >  1:38342/1‑135 (MQ=255)
                                              ccACATCGGTCGCCACCAGCACTTGCACCTCACGGGCCTTGAAGCTGTCCAGCGCACGCTGGCGGGTTGCCTGTGGGCGGTCACCGTGGATGCCGTCGGCATTCACGCCTTCGGCCAGCAAGCGCTCCACCAGCt                                                                                                                      <  2:38342/135‑1 (MQ=255)
                                                acctCGGTCGCCACCAGCACTTGCACCTCACGGGCCTTGAAGCTGTCCAGCGCACGCTGGCGGGTTGCCTGTGGGCGGTCACCGTGGATGCCGTCGGCATTCACGCCTTCGGCCAGCAAGCGCTCCACCAGCTGATCGACGCCATTGCgg                                                                                                     <  2:1019723/147‑1 (MQ=255)
                                                acaTCGGTCGCCACCAGCACTTGCACCTCACGGGCCTTGAAGCTGTCCAGCGCACGCTGGCGGGTTGCCTGTGGGCGGTCACCGTGGATGCCGTCGGCATTCACGCCTTCGGCCAGCAAGCGCTCCACCAGCTGATCGACGCCATTGCgg                                                                                                     >  2:1864684/1‑150 (MQ=255)
                                                      gTCGCCACCAGCACTTGCACCTCACGGGCCTTGAAGCTGTCCAGCGCACGCTGGCGGGTTGCCTGTGGGCGGTCACCGTGGATGCCGTCGGCATTCACGCCTTCGGCCAGCAAGCGCTCACACGTCTGATCGACGCCATTGCGGGTCTTg                                                                                               >  1:1811836/1‑150 (MQ=255)
                                                       tCGCCACCAGCACTTGCACCTCACGGGCCTTGAAGCTGTCCAGCGCACGCTGGCGGGTTGCCTGTGGGCGGTCACCGTGGATGCCGTCGGCATTCACGCCTTCGGCCAGCAAGCGCTCCACCAGCTGATCGACGCCATTGCGGGTCTTgg                                                                                              <  2:2187793/150‑1 (MQ=255)
                                                            accaGCACTTGCACCTCACGGGCCTTGAAGCTGTCCAGCGCACGCTGGCGGGTTGCCTGTGGGCGGTCACCGTGGATGCCGTCGGCATTCACGCCTTCg                                                                                                                                            >  1:55680/1‑99 (MQ=255)
                                                             ccaGCACTTGCACCTCACGGGCCTTGAAGCTGTCCAGCGCACGCTGGCGGGTTGCCTGTGGGCGGTCACCGTGGATGCCGTCGGCATTCACGCCTTCg                                                                                                                                            >  1:774774/1‑98 (MQ=255)
                                                              caGCACTTGCACCTCACGGGCCTTGAAGCTGTCCAGCGCACGCTGGCGGGTTGCCTGTGGGCGGTCACCGTGGATGCCGTCGGCATTCACGCCTTCGGCCAGCAAGCGCTCCACCAGCTGATCGACGCCATTGCGGGTCTTGGCGAagac                                                                                       <  2:1200036/150‑1 (MQ=255)
                                                                gCACTTGCACCTCACGGGCCTTGAAGCTGTCCAGCGCACGCTGGCGGGTTGCCTGTGGGCGGTCACCGTGGATGCCGTCGGCATTCACGCCTTCGGCCAGCAAGCGCTCCACCAGCTGATCGACGCCATTGCGGGTCTTGGCgccaaaca                                                                                     <  2:1243244/150‑8 (MQ=255)
                                                                   cTTGCACCTCACGGGCCTTGAAGCTGTCCAGCGCACGCTGGCGGGTTGCCTGTGGGCGGTCACCGTGGATGCCGTCGGCATTCACGCCTTCGGCCAGCAAGCGCTCCACCAGCTGATCGACGCCATTGCGGGTCTTGGCGAACAACAGCa                                                                                  <  2:1202608/150‑1 (MQ=255)
                                                                    ttGCACCTCACGGGCCTTGAAGCTGTCCAGCGCACGCTGGCGGGTTGCCTGTGGGCGGTCACCGTGGATGCCGTCGGCATTCACGCCTTCGGCCAGCAAGCGCTCCACCAGCTGATCGACGCCATTGCgg                                                                                                     >  2:856844/1‑130 (MQ=255)
                                                                    ttGCACCTCACGGGCCTTGAAGCTGTCCAGCGCACGCTGGCGGGTTGCCTGTGGGCGGTCACCGTGGATGCCGTCGGCATTCACGCCTTCGGCCAGCAAGCGCTCCACCAGCTGATCGACGCCATTGCgg                                                                                                     <  1:168352/130‑1 (MQ=255)
                                                                    ttGCACCTCACGGGCCTTGAAGCTGTCCAGCGCACGCTGGCGGGTTGCCTGTGGGCGGTCACCGTGGATGCCGTCGGCATTCACGCCTTCGGCCAGCAAGCGCTCCACCAGCTGATCGACGCCATTGCgg                                                                                                     >  2:168352/1‑130 (MQ=255)
                                                                    ttGCACCTCACGGGCCTTGAAGCTGTCCAGCGCACGCTGGCGGGTTGCCTGTGGGCGGTCACCGTGGATGCCGTCGGCATTCACGCCTTCGGCCAGCAAGCGCTCCACCAGCTGATCGACGCCATTGCGGGTCTTGTCGAACACCAGCAc                                                                                 <  2:2394989/150‑1 (MQ=255)
                                                                    ttGCACCTCACGGGCCTTGAAGCTGTCCAGCGCACGCTGGCGGGTTGCCTGTGGGCGGTCACCGTGGATGCCGTCGGCATTCACGCCTTCGGCCAGCAAGCGCTCCACCAGCTGATCGACGCCATTGCGGGTCTTGGCGAACACCAGCAc                                                                                 >  1:940824/1‑150 (MQ=255)
                                                                      gCACCTCACGGGCCTTGAAGCTGTCCAGCGCACGCTGGCGGGTTGCCTGTGGGCGGTCACCGTGGATGCCGTCGGCATTCACGCCTTCGGCCAGCAAGCGCTCCACCAGCTGATCGACGCCATTTCGGGTCTTGGCGAACACCAGCACCt                                                                               >  2:446492/1‑150 (MQ=255)
                                                                         ccTCACGGGCCTTGAAGCTGTCCAGCGCACGCTGGCGGGTTGCCTGTGGGCGGTCACCGTGGATGCCGTCGGCATTCACGCCTTCGGCCAGCAAGCGCTCCACCAGCTGATCGACGCCATTGCGGGTCTTGGCGAACACCAGCACCTGCt                                                                            >  1:122534/1‑150 (MQ=255)
                                                                          cTCACGGGCCTTGAAGCTGTCCAGCGCACGCTGGCGGGTTGCCTGTGGGCGGTCACCGTGGATGCCGTCGGCATTCACGCCTTCGGCCAGCAAGCGCTCCACCAGCTGATCGACGCCATTGCGGGTCTTGGCGAACACCAGCACCTGCtt                                                                           >  1:2025885/1‑150 (MQ=255)
                                                                            cACGGGCCTTGAAGCTGTCCAGCGCACGCTGGCGGGTTGCCTGTGGGCGGTCACCGTGGATGCCGTCGGCATTCACGCCTTCGGCCAGCAAGCGCTCCACCAGCTGATCGACGCCATTGCGGGTCTTGGCGAACACCAGCACCTGCTTcc                                                                         <  1:446492/150‑1 (MQ=255)
                                                                             aCGGGCCTTGAAGCTGTCCACCGCACGCTGGCGGGTTGCCGGTGGGCGGTCACCGTGGATGCCGTCGGCATTCACGCCTTCGGCCAGCAAGCGCTCCACCAGCTGATCGACGCCATTGCGGGTCTTGGCGAACACCAGCACCTGCTGCCa                                                                        <  1:1350989/150‑1 (MQ=255)
                                                                              cGGGCCTTGAAGCTGTCCAGCGCACGCTGGCGGGTTGCCTGTGGGCGGTCACCGTGGATGCCGTCGGCATTCACGCCTTCGGCCAGCAAGCGCTCCACCAGCTGATCGACGCCATTGCGGGTCTTGGCGAACACCAGCACCTGCTTCCAg                                                                       <  1:1236804/150‑1 (MQ=255)
                                                                                 gCCTTGAAGCTGTCCAGCGCACGCTGGCGGGTTGCCTGTGGGCGGTCACCGTGGATGCCGTCGGCATTCACGCCTTCGGCCAGCAAGCGCTCCACCAGCTGATCGACGc                                                                                                             >  2:1056861/1‑109 (MQ=255)
                                                                                      gAAGCTGTCCAGCGCACGCTGGCGGGTTGCCTGTGGGCGGTCACCGTGGATGCCGTCGGCATTCACGCCTTCGGCCAGCAAGCGCTCCACCAGCTGATCGACGCCATTGCGGGTCTTGGCGAACACCAGCACCTGCTTCCAGCGCTGCtt                                                               >  2:1444170/1‑150 (MQ=255)
                                                                                          cccgccAGCGTGCGCTGGCGGGTTGCCTGTGGGCGGTCACCGTGGATGCCGTCGGCATTCACGCCTTCGGCCAGCAAGCGCTCCACCAGCTGATCGACGCCATTGCGGGTCTTGGCGAACACCAGCACCTGCTTCCAGCGCTGCTTgcgc                                                           >  2:40270/5‑150 (MQ=255)
                                                                                             tCCAGCGCACGCTGGCGGGTTGCCTGTGGGCGGTCACCGTGGATGCCGTCGGCATTCACGCCTTCGGCCAGCAAGCGCTccac                                                                                                                           <  1:1189967/83‑1 (MQ=255)
                                                                                             acttgcgcACGCTGGCGGGTTGCCTGTGGGCGGTCACCGTGGATGCCGTCGGCATTCACGCCTTCGGCCAGCAAGCGCTCCACCAGCTGATCGACGCCATTGCGGGTCTTGGCGAACACCAGCACCTGCTTCCAGCGCTGCTTGCgcagc                                                        >  2:1553081/5‑150 (MQ=255)
                                                                                              ccAGCGCACGCTGGCGGGTTGCCTGTGGGCGGTCACCGTGGATGCCGTCGGCATTCACGCCTTCGGCCAGCAAGCGCGCCACCAGCTGATCGACGCCATTGCGGGTCTTGGCGAACACCAGCACCTGCTTCCAGCGCTGCTTGCgcagca                                                       <  2:754883/150‑1 (MQ=255)
                                                                                                aGCGCACGCTGGCGGGTTGCCTGTGGGCGGTCACCGTGGATGCCGTCGGCATTCACGCCTTCGGCCAGCAAGCGCTCCACCAGCTGATCGACGCCATTGCGGGTCTTGGCGAACACCAGCACCTGCTTCCAGCGCTGCTTGCGCAGCAgg                                                     <  1:2335605/150‑1 (MQ=255)
                                                                                                 gcgcACGCTGGCGGGTTGCCTGTGGGCGGTCACCGTGGATGCCGTCGGCATTCACGCCTTCGGCCAGCAAGCGCTCCACCAGCTGATCGACGCCATTGCGGGTCTTGGCGAACACCAGCACCTGCTTCCAGCGCTGCTTGCGCAGCAGGt                                                    >  1:857053/1‑150 (MQ=255)
                                                                                                  cgcACGCTGGCGGGTTGCCTGTGGGCGGTCACCGTGGATGCCGTCGGCATTCACGCCTTCGGCCAGCAAGCGCTCCACCAGCTGATCGACGCCATTGCGGGTCTTGGCGAACACCAGCACCTGCTTCCAGCGCTGCTTGCGCAGCAGGTg                                                   <  1:2361233/150‑1 (MQ=255)
                                                                                                   gcACGCTGGCGGGTTGCCTGTGGGCGGTCACCGTGGATGCCGTCGGCATTCACGCCTTCGGCCAGACCGTGCTCCACCAGCTGATCGACGCCATTGCGGGTCTTGGCGAACACCAGCACCTGCTTCCAGCGCTGCTTGCGCAGCAGGTgg                                                  <  1:1956558/150‑1 (MQ=255)
                                                                                                     aCGCTGGCGGGTTGCCTGTGGGCGGTCACCGTGGATGCCGTCGGTATTCACCCCTTCGGCCAGCAAGCGCTCCACCAGCTGCTCGACGCCATTGCGGGTCTTNGCGAACACCAGCACCTGCTTCCAGCGCTGCTTGCGCAGCAGGTGGCa                                                <  1:677362/150‑1 (MQ=255)
                                                                                                       gCTGGCGGGTTGCCTGTGGGCGGTCACCGTGGATGCCGTCGGCATTCACGCCTTCGGCCAGCAAGCGCTCCACCAGCTGATCGACGCCATTGCGGGTCTTGGCGAACACCAGCACCTGCTTCCAGCGCTGCTTGCGCAGCAGGTGCCAGa                                              >  2:1376571/1‑150 (MQ=255)
                                                                                                         tGGCGGGTTGCCTGTGGGCGGTCACCGTGGATGCCGTCGGCATTCACGCCTTCGGCCAGCAAGCGCTCCACCAGCTGATCGACGCCATTGCGGGTCTTGGCGAACACCAGCACCTGCTTCCAGCGCTGCTTGCGCAGCAGGTGGCAGAAc                                            <  1:474228/150‑1 (MQ=255)
                                                                                                           gCGGGTTGCCTGTGGGCGGTCACCGTGGATGCCGTCGGCATTCACGCCTTCGGCCa                                                                                                                                        <  1:1146120/56‑1 (MQ=255)
                                                                                                           gCGGGTTGCCTGTGGGCGGTCACCGTGGATGCCGTCGGCATTCACGCCTTCGGCCa                                                                                                                                        >  2:1146120/1‑56 (MQ=255)
                                                                                                           gCGGGTTGCCTGTGGGCGGTCACCGTGGATGCCGTCGGCATTCACGCCTTCGGCCAGCAAGCGCTCCACCAGCTGATCGACGCCATTGCGGGTCTTGGCGAACACCAGCACCTGCTTCCAGCGCTGCTTGCGCAGCAGGTGGCAGAACAg                                          >  1:1568018/1‑150 (MQ=255)
                                                                                                           gCGGGTTGCCTGTGGGCGGTCACCGTGGATGCCGTCGGCATTCACGCCTTCGGCCAGCAAGCGCTCCACCAGCTGATCGACGCCATTGCGGGTCTTGGCGAACACCAGCACCTGCTTCCAGCGCTGCTTGCGCAGCAGGTGGCAGAACAg                                          <  1:699159/150‑1 (MQ=255)
                                                                                                               gTTGCCTGTGGGCGGTCACCGTGGATGCCGTCGGCATTCACGCCTTCGGCCAGCAAGCGCTCCACCAGCTGATCGACGCCATTGCGGGTCTTGGCGAACACCAGCACCTGCTTCCAGCGCTGCTTGCGCAGCAGGTGGCAGAACAg                                          <  1:69435/146‑1 (MQ=255)
                                                                                                               gTTGCCTGTGGGCGGTCACCGTGGATGCCGTCGGCATTCACGCCTTCGGCCAGCAAGCGCTCCACCAGCTGATCGACGCCATTGCGGGTCTTGGCGAACACCAGCACCTGCTTCCAGCGCTGCTTGCGCAGCAGGTGGCAGAACAg                                          >  2:69435/1‑146 (MQ=255)
                                                                                                               gTTGCCTGTGGGCGGTCACCGTGGATGCCGTCGGCATTCACGCCTTCGGCCAGCAAGCGCTCCACCAGCTGATCGACGCCATTGCGGGTCTTGGCGAACACCAGCACCTGCTTCCAGCGCTGCTTGCGCAGCAGGTGGCAGAACAGGTCa                                      <  2:1841770/150‑1 (MQ=255)
                                                                                                               gTTGCCTGTGGGCGGTCACCGTGGATGCCGTCGGCATTCACGCCTTCGGCCAGCAAGCGCTCCACCAGCTGATCGACGCCATTGCGGGTCTTGGCGAACACCAGCACCTGCTTCCAGCGCTGCTTGCGCAGCAGGTGGCAGAACAGGTCa                                      >  2:766608/1‑150 (MQ=255)
                                                                                                                    cTGTGGGCGGTCACCGTGGATGCCGTCGGCATTCACGCCTTCAGCCACCAAGCGCTCCACCAGCTGATCGACGCCATTGCGGGTCTTGGCGAACACCAGCACCTGCTTCCAGCGCTGCTTGCGCAGCAGGTGGCAGAACAGGTCAGCCtt                                 <  1:1467433/150‑1 (MQ=255)
                                                                                                                     tgtgGGCGGTTACCGTTGATGCCGTCGGCATTCACGCCTTCGGCCAGCAAGCGCTCCACCAGCTGATCGACGCCATTGCGGGTCTTGGCGAACACCAGCACCTGCTTCCAGCGCTGCTTGCGCAGCAGGTGGCAGAACAGTTCAGCCTTg                                <  2:1079711/150‑1 (MQ=255)
                                                                                                                      gtgGGCGGTCACCGTGGATGCCGTCGGCATTCACGCCTTCGGCCAGCAAGCGCTCCACCAGCTGATCGACGCCATTGCGGGTCTTGGCGAACACCAGCACCTGCTTCCAGCGCTGCTTGCGCAGCAGGTGGCAGAACAGGTCAGCCTTgc                               <  2:1568018/150‑1 (MQ=255)
                                                                                                                          gCGGTCACCGTGGATGCCGTCGGCATTCACGCCTTCGGCCAGCAAGCGCTCCACCAGCTGATCGACGCCATTGCGGGTCTTGGCGAACACCAGCACCTGCTTCCAGCGCTGCTTGCGCAGCAGGTGGCAGAACAGGTCAGCCTTGCGctt                           <  2:1809402/150‑1 (MQ=255)
                                                                                                                            ggTCACCGTGGATGCCGTCGGCATTCACGCCTTCGGCCAGCAAGCGCTCCACCAGCTGATCGACGCCATTGCGGGTCTTGGCGAACACCAGCAcc                                                                                >  2:586498/1‑95 (MQ=255)
                                                                                                                             gTCACCGTGGATGCCGTCGGCATTCACGCCTTCGGCCAGAAAGCGCTCCACCAGCTGATCGACGCCATTGCGGGTCTTGGCGAACACCa                                                                                     <  2:1751047/89‑1 (MQ=255)
                                                                                                                             gTCACCGTCGATGCCGTCGGCATTCACGCCTTCGGCCAGCAAGCGCTCCACCAGCTGATCGACGCCATTGCGGGTCTTGGCGAACACCAGCACCTGCTTCCAGCGCTGCTTGCGCAGCAGGTGGCAGAACAGGTCAGCCTTGCGcttctt                        <  1:1037922/150‑1 (MQ=255)
                                                                                                                             gTCACCGGGGATGCCGTCGGCATTCACGCCTTCGGCCAGAAAGCGCTCCACCAGCTGATCGACGCCATTGCGGGTCTTGGCGAACACCa                                                                                     >  1:1751047/1‑89 (MQ=255)
                                                                                                                               cACCGTGGATGCCGTCGGCATTCACGCCTTCGGCCAGCAAGCGCTCCACCAGCTGATCGACGCCATTGCGGGTCTTGGCGAACACCAGCACCTGCTTCCAGCGCTGCTTGCGCAGCAGGTGGCAGAACAGGTCAGCCTTGCGCTTCTTGt                      <  2:2100580/150‑1 (MQ=255)
                                                                                                                                aCCGTGGATGCCGTCGGCATTCACGCCTTCGGCCAGCAAGCGCTCCACCAGCTGATCGACGCCATTGCGGGTCTTGGCGAACACCAGCACCTGCTTCCAGCGCTGCTTGCGCAGCAGGTGGCAGAACAGGTCAGCCTTGCGCTTCTTGTc                     <  1:1553081/150‑1 (MQ=255)
                                                                                                                                aCCGTGGATGCCGTCGGCATTCACGCCTTCGGCCAGCAAGCGCTCCACCAGCTGATCGACGCCATTGCGGGTCTTGGCGAACACCAGCACCTGCTTCCAGCGCTGCTTGCGCAGCAGGTGGCAGAACAGGTCAGCCTTGCGCTTCTTGTc                     >  2:1317348/1‑150 (MQ=255)
                                                                                                                                aCCGTGGATGCCGTCGGCATTCACGCCTTCGGCCAGCAAGCGCTCCACCAGCTGATCGACGCCATTGCGGGTCTTGGCGAACACCAGCACCTGCTTCCAGCGCTGCTTGCGCAGCAGGTGGCAGAACAGGTCAGCCTTGCGCTTCTTGTc                     <  2:1969830/150‑1 (MQ=255)
                                                                                                                                  cGTGGATGCCGTCGGCATTCACGCCTTCGGCCAGCAAGCGCTCCACCAGCTGATCGACGCCATTGCGGGTCTTGGCGAACACCAGCACCTGCTTCCAGCGCTGCTTGCGCAGCAGGTGGCAGAACAGGTCAGCCTTGCGCTTCTTGTCca                   <  2:680483/150‑1 (MQ=255)
                                                                                                                                   gTGGATGCCGTCGGCATTCACGCCTTCGGCCAGCAAGCGCTCCACCAGCTGATCGACGCCATTGCGGGTCTTGGCGAACACCAGCACCTGCTTCCAGCGCTGCTTGCGCAGCAGGTGGCAGAACAGGTCAGCCTTGCGCTTCTTg                       >  1:326937/1‑145 (MQ=255)
                                                                                                                                   gTGGATGCCGTCGGCATTCACGCCTTCGGCCAGCAAGCGCTCCACCAGCTGATCGACGCCATTGCGGGTCTTGGCGAACACCAGCACCTGCTTCCAGCGCTGCTTGCGCAGCAGGTGGCAGAACAGGTCAGCCTTGCGCTTCTTg                       <  2:326937/145‑1 (MQ=255)
                                                                                                                                   gTGAATGCCGACGGCATCCACGCCTTCGGCCAGCAAGCGCTCCACCAGCTGATCGACGCCATTGCGGGTCTTGGCGAACACCAGCACCTGc                                                                             <  2:1174711/91‑1 (MQ=255)
                                                                                                                                   gTGAATGCCGACGGCATCCACGCCTTCGGCCAGCAAGCGCTCCACCAGCTGATCGACGCCATTGCGGGTCTTGGCGAACACCAGCACCTGc                                                                             >  1:1174711/1‑91 (MQ=255)
                                                                                                                                     ggATGCCGTCGGCATTCACGCCTTCGGCCAGCAAGCGCTCCACCAGCTGATCGACGCCATTGCGGGTCTTGGCGAACACCAGCACCTGCTTCCAGCGCTGCTTGCGCAGCAGGt                                                    <  1:2369684/114‑1 (MQ=255)
                                                                                                                                      gATGCCGTCGGCATTCACGCCTTCGGCCAGCAAGCGCTCCACCAGCTGATCGACGCCATTGCGGGTCTTGGCGAACACCa                                                                                     >  1:1329108/1‑80 (MQ=255)
                                                                                                                                       aTGCCGTCGGCATTCACGCCTTCGGCCAGCAAGCGCTCCACCAGCTGATCGACGCCATTGCGGGTCTTGGCGAACACCAGCACCTGCTTCCAGCGCTGCTTGCGCAGCAGGTGGCAGAACAGGTCAGCCTTGCGCTTCTTGTCCACAGGc              <  1:1444170/150‑1 (MQ=255)
                                                                                                                                           cGTCGGCATTCACGCCTTCGGCCAGCAAGCGCTCCACCAGCTGATCGACGCCATTGCGGGTCTTGGCGAACACCAGCACCTGCTTCCAGCGCTGCTTGCGCAGCAGGTGGCAGAACAGGTCAGCCTTGCGcttctt                        >  2:2111191/1‑136 (MQ=255)
                                                                                                                                           cGTCGGCATTCACGCCTTCGGCCAGCAAGCGCTCCACCAGCTGATCGACGCCATTGCGGGTCTTGGCGAACACCAGCACCTGCTTCCAGCGCTGCTTGCGCAGCAGGTGGCAGAACAGGTCAGCCTTGCGcttctt                        <  1:2111191/136‑1 (MQ=255)
                                                                                                                                           cGTCGGCATGCACGCCTTCGGCCAGCAAGCGCTCCACCAGCTGATCGACGCCATTGCGGGTCTTGGCGAACACCAGCACCTGCTTCCAGCGCTGCTTGCGCAGCAGGTGGCAGAACAGGTCAGCCTTGCGCTTCTTGTCCACAGGCACCa          <  2:1811836/150‑1 (MQ=255)
                                                                                                                                            gTCGGCATTCACGCCTTCGGCCAGCAAGCGCTCCACCAGCTGATCGACGCCATTGCGGGTCTTGGCGAACACCAGCACCTGCTTCCAGCGCTGCTTGCGCAGCAGGTGGCAGAACAGGTCAGCCTTGCGCTTCTTGTCCACAGGCACCAg         >  2:2209782/1‑150 (MQ=255)
                                                                                                                                               ggCATTCACGCCTTCGGCCAGCAAGCGCTCCACCAGCTGATCGACGCCATTGCGGGTCTTGGCGAACACCa                                                                                     >  2:1234906/1‑71 (MQ=255)
                                                                                                                                               ggCATTCACGCCTTCGGCCAGCAAGCGCTCCACCAGCTGATCGACGCCATTGCGGGTCTTGGCGAACACCAGCACCTGCTTCCAGCGCTGCTTGCGCAGCAGGTGGCAGAACAGGTCAGCCTTGCGCTTCTTGTCcaca                 >  1:2219605/1‑139 (MQ=255)
                                                                                                                                               ggCATTCACGCCTTCGGCCAGCAAGCGCTCCACCAGCTGATCGACGCCATTGCGGGTCTTGGCGAACACCAGCACCTGCTTCCAGCGCTGCTTGCGCAGCAGGTGGCAGAACAGGTCAGCCTTGCGCTTCTTGTCcaca                 <  2:2219605/139‑1 (MQ=255)
                                                                                                                                                  agtCACGCCTTCGGCCAGCAAGCGCTCCACCAGCTGATCGACGCCATTGCGGGTCTTGGCGAACACCAGCACCTGCTTCCAGCGCTGCTTGCGCAGCAGGTGGCAGAACAGGTCAGCCTTGCGCTTCTTGTCCACAGGCACCAGCCACTg   >  2:2212825/3‑150 (MQ=255)
                                                                                                                                                  aTTCACGCCTTCGGCCAGCAAGCGCAAAACAAGCTGATCAAAACAACAGAAGCACTTGGCAAACACCAGCAACTACTGCCAGGGCTGCTTGCGCAGCAGGTGGCAGAACAGGTCAGCCTTGCGCTTCTTGTCCACAGGCACCAGCCACTg   >  2:275903/1‑150 (MQ=255)
                                                                                                                                                   ttCACGCCTTCGGCCAGCAAGCGCTCCACCAGCTGATCGACGCCATTGCGGGTCTTGGCGAACACCAGCACCTGCTTCCAGCGCTGCTTGCGCAGCAGGTGGCAGAACAGGTCAGCCTTGCGCTTCTTGTCCACAGGCACCAGCCACTGc  <  2:1180443/150‑1 (MQ=255)
                                                                                                                                                    tCACGCCTTCGGCCAGCAAGCGCTCCACCAGCTGATCGACGCCATTGCGGGTCTTGGCGAACACCAg                                                                                    >  1:541671/1‑67 (MQ=255)
                                                                                                                                                    tCACGCCTTCGGCCAGCAAGCGCTCCACCAGCTGATCGACGCCATTGCGGGTCTTGGCGAACACCAg                                                                                    <  2:541671/67‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
AGGTTGACGACCAGCGGCAAGTCATCGATATCCAGGCCACGGGCCGCCACATCGGTCGCCACCAGCACTTGCACCTCACGGGCCTTGAAGCTGTCCAGCGCACGCTGGCGGGTTGCCTGTGGGCGGTCACCGTGGATGCCGTCGGCATTCACGCCTTCGGCCAGCAAGCGCTCCACCAGCTGATCGACGCCATTGCGGGTCTTGGCGAACACCAGCACCTGCTTCCAGCGCTGCTTGCGCAGCAGGTGGCAGAACAGGTCAGCCTTGCGCTTCTTGTCCACAGGCACCAGCCACTGC  >  NC_002947/5427195‑5427491

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