Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 626,572 C→A 6.4% L208L (CTG→CTT dedA ← inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947626,5720CA 6.4% 295.7 / 5.8 93L208L (CTG→CTTdedAinner membrane protein
Reads supporting (aligned to +/- strand):  ref base C (45/42);  new base A (3/3);  total (48/45)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.62e-01

GCCGTCCAGCATTGGCAGCCGCTCGCGTACGGCCTGCCATTGCTGTGGCGTGATCGGGTAACGCGCCAGAGTACGCTTGCGCCGCCAGGCGCTGAATGACCACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCAGGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACAGGATGATGATGCCCACCACCATCAGCGACAGGTGCTGCTTGATGAATGGCACGTTGCCGAAGAAGTAGCCCAAGGTGACCAGGCCGCCGACCCAAAGCAGCGAAC  >  NC_002947/626423‑626720
                                                                                                                                                     |                                                                                                                                                    
gCCGTCCAGCATTGGCAGCCGCTCGCGTACGGCCTGCCATTGCTGTGGCGTGATCGGGTAACGCGCCAGAGTACGCTTGCGCCGCCAGGCGCTGAATGACCACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTACGGccc                                                                                                                                                      >  2:1140584/1‑150 (MQ=255)
 ccGTCCAGCATTGGCAGCCGCTCGCGTACGGCCTGCCATTGCTGTGGCGTGATCGGGTAACGCGCCAGAGTACGCTTGCGCCGCCAGGCGCTGAATGACCACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCa                                                                                                                                                     <  1:822428/150‑1 (MQ=255)
  cGTCCAGCATTGGCAGCCGCTCGCGTACGGCCTGCCATTGCTGTGGCGTGATCGGGTAACGCGCCAGAGTACGCTTGCGCCGCCAGGCGCTGAATGACCACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCAg                                                                                                                                                    <  1:1140584/150‑1 (MQ=255)
   gTCCAGCATTGGCAGCCGCTCGCGTACGGCCTGCCATTGCTGTGGCGTGATCGGGTAACGCGCCAGAGTACGCTTGCGCCGCCAGGCGCTGAATGACCACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCAgg                                                                                                                                                   >  2:2337084/1‑150 (MQ=255)
   gTCCAGCATTGGCAGCCGCTCGCGTACGGCCTGCCATTGCTGTGGCGTGATCGGGTAACGCGCCAGAGTACGCTTGCGCCGCCAGGCGCTGAATGACCACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCAAgg                                                                                                                                                   <  1:458775/150‑1 (MQ=255)
   gTCCAGCATTGGCAGCCGCTCGCGTACGGCCTGCCATTGCTGTGGCGTGATCGGGTAACGCGCCAGAGTACGCTTGCGCCGCCAGGCGCTGAATGACCACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTGGGCCGTACGGCCAAgg                                                                                                                                                   <  1:1596867/150‑1 (MQ=255)
    tCCAGCATTGGCAGCCGCTCGCGTACGGCCTGCCACTGCTGTGGCGTGATCGGGTAACGCGCCAGAGTACGCTTGCGCCGCCAGGCGCTGAATGACCACATGCCGTTCGTGGTTGACCTGTTCAGTGAGCCTTGGCCGTACGGCCCAGgc                                                                                                                                                  >  1:739211/1‑150 (MQ=255)
    cacagcaatggcaggccgtacgcgagcggctgccaatgctgTGGCGTGATCGGGTAACGCGCCAGAGTACGCTTGCGCCGCCAGGCGCTGAATGACCACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCAGgc                                                                                                                                                  >  1:2371468‑M2/42‑150 (MQ=255)
        GCATTGGCAGCCGCTCGCGTACGGCCTGCCATTGCTGTGGCGTGATCGGGTAACGCGCCAGAGTACGCTTGCGCCGCCAGGCGCTGAATGACCACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTGGGCCGTACGGCcaaggctcac                                                                                                                                              >  2:400390‑M1/1‑146 (MQ=255)
                 gCCGCTCGCGTACGGCCTGCCATTGCTGTGGCGTGATCGGGTAACGCGCCAGAGTACGCTTGCGCCGCCAGGCGCTGAATGACCACATGCCGTTCGTTCTTTACCTGTTCCGTGAGCCTTGGCCGTACGGCCCAGGCGCCCACGCAGCaa                                                                                                                                     >  1:1231428/1‑150 (MQ=255)
                 gCCGCTCGCGTACGGCCTGCCATTGCTGTGGCGTGATCGGGTAACGCGCCAGAGTACGCTTGCGCCGCCAGGCGCTGAATGACCACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCAGGCGCCCACGCAGCaa                                                                                                                                     <  2:1151292/150‑1 (MQ=255)
                  ccGCTCGCGTACGGCCTGCCATTGCTGTGGCGTGATCGGGTAACGCGCCAGAGTACGCTTGCGCCGCCAGGCGCTGAATGACCACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCAGGCGCCCACGCAGCAAg                                                                                                                                    <  1:2198611/150‑1 (MQ=255)
                    gCTCGCGTACGGCCTGCCATTGCTGTGGCGTGATCGGGTAACGCGCCAGAGTACGCTTGCGCCGCCAGGCGCTGAATGACCACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCAGGCGCCCACGCAGCAAGcc                                                                                                                                  <  2:1790007/150‑1 (MQ=255)
                    gCTCGCGTACGGCCTGCCATTGCTGTGGCGTGATCGGGTAACGCGCCAGAGTACGCTTGCGCCGCCAGGCGCTGAATGACCACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCAGGCGCCCACGCAGCAAGcc                                                                                                                                  <  1:617628/150‑1 (MQ=255)
                              ggCCTGCCATTGCTGTGGCGTGCTAGGGTGCCCAGGCAGATGACTCTTGCGCCGCCATTCGCTGAATGACCACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCAGGCGCCCACGCCGCAAGCCGAGAatcatc                                                                                                                        <  2:1231428/150‑1 (MQ=255)
                              ggCCTGCCATTGCTGTGGCGTGATCGGGTAACGCGCCAGAGTACGCTTGCGCCGCCAGGCGCTGAATGACCACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCAGGCGCCCACGCAGCAAGCCGAGAatcatc                                                                                                                        >  1:2197910/1‑150 (MQ=255)
                                 cTGCCATTGCTGTGGCGTGATCGGGTAACGCGCCAGAGTACGCTTGCGCCGCCAGGCGCTGAATGACCACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCAGGCGCCCACGCAGCAAGCCGAGAATCATCGGc                                                                                                                     >  2:2049025/1‑150 (MQ=255)
                                   gCCATTGCTGTGGCGTGATCGGGTAACGCGCCAGAGTACGCTTGCGCCGCCAGGCGCTGAATGACCACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCAGGCGCCCACGCAGCAAGCCGAGAATCATCGGCAc                                                                                                                   >  2:2222912/1‑150 (MQ=255)
                                       ttGCTGTGGCGTGATCGGGTAACGCGCCAGAGTACGCTTGCGCCGCCAGGCGCTGAATGACCACATGCCTTTCGTTCTTTACCTGTTCGGTGAGCCTTGGCCGTACGGCCCAGGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGc                                                                                                               >  1:816159/1‑150 (MQ=255)
                                         gctgTGGCGTGATCGGGTAACGCGCCAGAGTACGCTTGCGCCGCCAGGCGCTGAATGACCACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCAGGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGa                                                                                                             <  1:1055312‑M2/146‑1 (MQ=255)
                                          ctgTGGCGTGATCGGGTAACGCGCCAGAGTACGCTTGCGCCGCCAGGCGCTGAATGACCACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCAGGCGCCCACgca                                                                                                                                         >  2:2399358‑M2/4‑121 (MQ=255)
                                          ctgTGGCGTGATCGGGTAACGCGCCAGAGTACGCTTGCGCCGCCAGGCGCTGAATGACCACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCAGGCGCCCACgca                                                                                                                                         <  1:2399358‑M2/118‑1 (MQ=255)
                                           tgTGGCGTGATCGGGTAACGCGCCTGAGTACGCTTGCGCCGCCAGGCGCTGAATGACCACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGCGGTACGGCCCAGGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACa                                                                                                           >  1:1055114‑M2/3‑150 (MQ=255)
                                           tgTGGCGTGATCGGGTAACGCGCCAGAGTACGCTTGCGCCGCCAGGCGCTGAATGACCACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCAGGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACa                                                                                                           <  2:2051002‑M2/148‑1 (MQ=255)
                                            gTGGCGTGATCGGGTAACGCGCCAGAGTACGCTTGCGCCGCCAGGCGCTGAATGACCACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCAGGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACAg                                                                                                          <  1:692635‑M2/149‑1 (MQ=255)
                                                 gTGATCGGGTAACGCGCCAGAGTACGCTTGCGCCGCCAGGCGCTGAATGACCACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCAGGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACAGgatg                                                                                                      >  1:812235/1‑149 (MQ=255)
                                                 gTGATCGGGTAACGCGCCAGAGTACGCTTGCGCCGCCAGGCGCTGAATGACCACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCAGGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACAGgatg                                                                                                      <  2:812235/149‑1 (MQ=255)
                                                    aTCGGGTAACGCGCCAGAGTACGCTTGCGCCGCCAGGCGCTGAATGACCACATGCCTTTCGTTCTTTACCTGTTCGGTGAGCCTTGGCCGTACGGCCCAGGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACAGgatgatga                                                                                                  <  2:816159/150‑1 (MQ=255)
                                                     tCGGGTAACGCGCCAGAGTACGCTTGCGCCGCCAGGCGCTGAATGACCACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCAGGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACAGgatgatgat                                                                                                 <  1:1332893/150‑1 (MQ=255)
                                                       gggTAACGCGCCAGAGTACGCTTGCGCCGCCAGGCGCTGAATGACCACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCAGGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACAGGATGATGATGc                                                                                               >  1:875774/1‑150 (MQ=255)
                                                           aaCGCGCCAGAGTACGCTTGCGCCGCCAGGCGCTGAATGACCACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCAGGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACAGGATGATGATGCccac                                                                                           >  1:1672169/1‑150 (MQ=255)
                                                               cgcCAGAGTACGCTTGCGCCGCCAGGCGCTGAATGACCACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCAGGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACAGGATGATGATGCccaccacc                                                                                       <  2:2197910/150‑1 (MQ=255)
                                                                gcCAGAGTACGCTTGCGCCGCCAGGCGCTGAATGACCACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCAGGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACAGGATGATGATGCccaccacca                                                                                      >  2:515279/1‑150 (MQ=255)
                                                                gcCAGAGTACGCTTGCGCCGCCAGGCGCTGAATGACCACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTGGGCCGTACGGCCAAgg                                                                                                                                                   >  1:30825/1‑89 (MQ=255)
                                                                gcCAGAGTACGCTTGCGCCGCCAGGCGCTGAATGACCACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTGGGCCGTACGGCCAAgg                                                                                                                                                   <  2:30825/89‑1 (MQ=255)
                                                                    gcgTACGCTTGCGCCGCCAGGCGCTGACTCAACACTCGTAGGTTGTTCTTTACCCGTTCCGTGAGCCTTGGCCGTCCGGCGCAGGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACAGGATGATGATGCCCACCACCATCAg                                                                                  <  2:803811/148‑1 (MQ=255)
                                                                            ttGCGCCGCCAGGCGCTGAATGACCACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCAGGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACAGGATGATGATGCCCACCACCATCAGCGACAGGt                                                                          >  1:1301904/1‑150 (MQ=255)
                                                                            ttGCGCCGCCAGGCGCTGAATGACCACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCAGGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACAGGATGATGATGCCCACCACCATCAGCGACAGGt                                                                          >  1:515965/1‑150 (MQ=255)
                                                                            ttGCGCCGCCAGGCGCTGAATGACCACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCAGGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACAGGATGATGATGCCCACCACCATCAGCGACAGGt                                                                          >  1:1491404/1‑150 (MQ=255)
                                                                              gcgcCGCCAGGCGCTGAATGACCCCATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCAGGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACAGGATGATGATGCCCACCACCATCAGCGACAGGtgc                                                                        <  2:739211/150‑1 (MQ=255)
                                                                                      aGGCGCTGAATGACCACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCAGGCGCCCACGCAGCAAGCCGAGAATCATCg                                                                                                                       >  1:849078/1‑95 (MQ=255)
                                                                                      aGGCGCTGAATGACCACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCAGGCGCCCACGCAGCAAGCCGAGAATCATCg                                                                                                                       <  2:849078/95‑1 (MQ=255)
                                                                                      aGGCGCTGAATGACCACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCAGGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACAGGATGATGATGCCCACCACCATCAGCGACAGGTGCTGCTtgat                                                                <  2:1301904/150‑1 (MQ=255)
                                                                                         cgcTGAATGACCACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCAGGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACAGGATGATGATGCCCACCACCATCAGCGACAGGTGCTGCTTGATGaa                                                             >  2:938827/1‑150 (MQ=255)
                                                                                          gcTGAATGACCACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCAGGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACAGGATGATGATGCCCACCACCATCAGCGACAGGTGCTGCTTGATGAAt                                                            <  2:1617875/150‑1 (MQ=255)
                                                                                          gcTGAATGACCACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCAGGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACAGGATGATGATGCCCACCACCATCAGCGACAGGTGCTGCTTGATGAAt                                                            <  1:2337084/150‑1 (MQ=255)
                                                                                          gcTGAATGACCACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCAGGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACAGGATGATGATGCCCACCACCATCAGCGACAGGTGCTGCTTGATGAAt                                                            <  1:515279/150‑1 (MQ=255)
                                                                                          gcTGAATGACCACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTGGGCCGTACGGCCCAGGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACAGGATGATGATGCCCACCACCATCAGCGACAGGTGCTGCTTGATGAAt                                                            <  1:1479112/150‑1 (MQ=255)
                                                                                            tGAATGACCACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCAGGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACAGGATGATGATGCCCACCACCATCAGCGACAGGTGCTGCTTGATGAATgg                                                          >  2:718509/1‑150 (MQ=255)
                                                                                                   ccACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCAGGCGCCCACGCAGCAAGCCGAGAATCATCGGCa                                                                                                                    >  2:1779204/1‑85 (MQ=255)
                                                                                                   ccACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCAGGCGCCCACGCAGCAAGCCGAGAATCATCGGCa                                                                                                                    <  1:1779204/85‑1 (MQ=255)
                                                                                                   ccACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCAGGCGCCCACGCAGCAAGCCGAGAAGCTTCGGCACCAGCGAAAGGTTGCTCATGCCATGGGCCAGCGGTTACAGGTGCTGCTTGATGAATGGCACGTTg                                                   >  1:2341037/1‑150 (MQ=255)
                                                                                                    cacaTGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCAGGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACAGGATGATGATGCCCACCACCATCAGCGACAGGTGCTGCTTGATGAATGGCACGTTGc                                                  <  1:718509/150‑1 (MQ=255)
                                                                                                      caTGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCAGGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACAGGATGATGATGCCCACCACCATCAGCGACAGGTGCTGCTTGATGAATGGCACGTTGCCg                                                >  1:581288/1‑150 (MQ=255)
                                                                                                       aTGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCAGGCGCCACCAGCGCAAGCCGAGAAGCAGCTGCACCAGCGACAGGATGATGATGCCCACCACCATCAGCGACAGGTGCTGCTTGATGAATGGCACGTTGCCga                                               >  1:2100527/1‑150 (MQ=255)
                                                                                                          ccGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTACGGTCCAGGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACAGGATGATGATGCCCACCACCATCAGCGACAg                                                                            <  1:250180/118‑1 (MQ=255)
                                                                                                          ccGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTACGGTCCAGGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACAGGATGATGATGCCCACCACCATCAGCGACAg                                                                            >  2:250180/1‑118 (MQ=255)
                                                                                                           cgttcgttCTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCAGGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACAGGATGATGATGCCCACCACCATCAGCGACAGGTGCTGCTTGATGAATGGCACGTTGCCgaagaa                                           <  2:1002370/150‑1 (MQ=255)
                                                                                                            gttcgttcTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCAGGCGCCCACGCAGCAAGCCGAGAATCATCg                                                                                                                       >  2:73729/1‑73 (MQ=255)
                                                                                                            gttcgttcTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCAGGCGCCCACGCAGCAAGCCGAGAATCATCg                                                                                                                       <  1:73729/73‑1 (MQ=255)
                                                                                                             ttcgttcTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCAGGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACAGGATGATGATGCCCACCACCATCAGCGACAg                                                                            >  2:2246058/1‑115 (MQ=255)
                                                                                                             ttcgttcTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCAGGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACAGGATGATGATGCCCACCACCATCAGCGACAg                                                                            <  1:2246058/115‑1 (MQ=255)
                                                                                                                gttcTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCAGGCGCCCACGCAGCAAGCCGAGAATCTTTGGCCTCGGCGACAGGATGATGATGCCCACCACCATCAGCGACAGGTGCTGCTTGATGAATGGCACGTTGCCGAAGAAGTAGc                                      <  1:2222912/150‑1 (MQ=255)
                                                                                                                gttcTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCAGGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACAGGATGATGATGGCGACTACCATCAGCGACAGGTGCTGCTTGATGAATGGCACGTTGCCGAAGAAGTAGc                                      >  1:24448/1‑150 (MQ=255)
                                                                                                                  tcTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCAGGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACAGGATGATGATGCCCACCACCATCAGCGACAGGTGCTGCTTGATGAATGGCACGTTGCCgaa                                              <  1:2231951/140‑1 (MQ=255)
                                                                                                                  tcTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCAGGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACAGGATGATGATGCCCACCACCATCAGCGACAGGTGCTGCTTGATGAATGGCACGTTGCCgaa                                              >  2:2231951/1‑140 (MQ=255)
                                                                                                                    tttACCTGTTCAGTGAGCCTTGGCCGTACGGCCCAGGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACAGGATGATGATGc                                                                                               <  2:1650243/89‑1 (MQ=255)
                                                                                                                    tttACCTGTTCAGTGAGCCTTGGCCGTACGGCCCAGGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACAGGATGATGATGc                                                                                               >  1:1650243/1‑89 (MQ=255)
                                                                                                                    tttACCTGTTCAGTGAGCCTTGGCCGTACGGCCCAGGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACAGGATGATGATGCCCACCACCATCAGCGACAGGTGCTGCTTGATGaa                                                             <  1:2235079/123‑1 (MQ=255)
                                                                                                                    tttACCTGTTCAGTGAGCCTTGGCCGTACGGCCCAGGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACAGGATGATGATGCCCACCACCATCAGCGACAGGTGCTGCTTGATGaa                                                             >  2:2235079/1‑123 (MQ=255)
                                                                                                                       aCCTGTTCAGTGAGCCTTGGCCGTACGGCCCAGGCGCCCACGCAGCAAGCCGAGAATCATCGGCAACAGCGACAGGATGATGCTGCCCACCACCATCAGCGACAGGTGCTGCTTGATGAATGGCACGTTGCCGAAGAAGTAGCCCAAGGt                               >  2:1279224/1‑150 (MQ=255)
                                                                                                                        ccTGTTCAGTGAGCCTTGGCCGTACGGCCCAAGAGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACAGGATGATGATGCCCACCACCATCAGCGACAGGTGCTGCTTGATGAATGGCACGTTGCCGAAGAAGTAGCCCAAGGTg                              >  1:315890/1‑150 (MQ=255)
                                                                                                                         cTGTTCAGTGAGCCTTGGCCGTACGGCCCAGGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACAGGATGATGATGCCCACCACCATCAGCGACAGGTGCTGCTTGATGAATGGCACGTTGCCGAAGAAGTAGCCCAAGGTGa                             >  1:1095200/1‑150 (MQ=255)
                                                                                                                         cTGTTCAGTGAGCCTTGGCCGTACGGCCCAGGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACAGGATGATGATGCCCACCACCATCAGCGACAGGTGCTGCTTGATGAATGGCACGTTGCCGAAGAAGTAGCCCAAGGTGa                             >  1:1050747/1‑150 (MQ=255)
                                                                                                                          tGTTCAGTGAGCCTTGGCCGTACGGCCCAGGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACAGGATGATGATGCCCACCACCATCAGCGACAGGTGCTGCTTGATGAATGGCACGTTGCCGAAGAAGTAGCCCAAGGTGAc                            >  2:700146/1‑150 (MQ=255)
                                                                                                                            ttCAGTGAGCCTTGGCCGTACGGCCCAGGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACAGGATGATGATGCCCACCACCATCAGCGACAGGTGCTGCTTGATGAATGGCACGTTGCCGAAGAAGTAGCCCAAGGTGACCa                          <  1:938827/150‑1 (MQ=255)
                                                                                                                               aGTGAGCCTTGGCCGTACGGCCCAGGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACAGGATGATGATGCCCACCACCATCAGCGACAGGTGCTGCTTGATGAAt                                                            >  2:1350508/1‑113 (MQ=255)
                                                                                                                               aGTGAGCCTTGGCCGTACGGCCCAGGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACAGGATGATGATGCCCACCACCATCAGCGACAGGTGCTGCTTGATGAAt                                                            <  1:1350508/113‑1 (MQ=255)
                                                                                                                               aGTGAGCCTTGGCCGTACGGCCCAGGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACAGGATGATGATGCCCACCACCATCAGCGACAGGTGCTGCTTGATGAATGGCACGTTGCCGAAGAAGTAGCCCAAGGAGACCAGgc                       >  1:576231/1‑150 (MQ=255)
                                                                                                                                 tGAGCCTTGGCCGTACGGCCCAGGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACAGTATGATGATGCCCACCACCATCAGCGACAGGTGCTGCTTGATGAATGGCACGTTGCCGAAGAAGAAGCCCAAGGTGACCAAgccg                     >  2:250674/1‑150 (MQ=255)
                                                                                                                                 tGAGCCTTGGCCGTACGGCCCAGGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACAGGATGATGATGCCCACCACCATCAGCGACAGGTGCTGCTTGATGAATGGCACGTTGCCGAAGAAGTAGCCCAAGGTGACCAGgccg                     >  1:1419255/1‑150 (MQ=255)
                                                                                                                                  gAGCCTTGGCCGTACGGCCCAGGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACAGGATGATGATGCCCACCACCATCAGCGACAGGTGCTGCTTGATGAATGGCACGTTGCCGAAGAAGTAGCCCAAGGTGACCAGgccgc                    <  2:581288/150‑1 (MQ=255)
                                                                                                                                   aGCCTTGGCCGTACGGCCCAGGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACAGGATGATGATGCCCACCACCATCAGc                                                                                 <  1:1810861/88‑1 (MQ=255)
                                                                                                                                   aGCCTTGGCCGTACGGCCCAGGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACAGGATGATGATGCCCACCACCATCAGc                                                                                 >  2:1810861/1‑88 (MQ=255)
                                                                                                                                     ccTTGGCCGTACGGCCCAGGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACAGGATGATGATGCCCACCACCATCAGCGACAGGTGCTGCTTGATGAATGGCACGTTGCCGAAGAAGTAGCCCAAGGTGACCAGGCCGCCGa                 <  2:1483171/150‑1 (MQ=255)
                                                                                                                                      ctgggCCGTACGGCCAAGGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACAGGATGATGATGCCCACCACCATCAGCGACAGGTGCTGCTTGATGAATGGCACGTTGCCGAAGAAGTAGCCCAAGGTGACCAGGCCGCCGAc                >  2:1897172/4‑150 (MQ=255)
                                                                                                                                       ttGGCCGTACGGCCCAGGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACAGGATGATGATGCCCACCATCATCAGCGACAGGTGCTCCTCGGCGACTGTCACGTGGCCTAAGAAGTAGCCCAAGGTGACCAGGCCGCCGAcc               >  1:1740753/1‑150 (MQ=255)
                                                                                                                                       ttGGCCGTACGGCCCAGGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACAGGATGATGATGCCCACCACCATCAGCGACAGGTGCTGCTTGATGAATGGCACGTTGCCGAAGAAGTAGCCCAAGGTGACCAGGCCGCCGAcc               >  1:1974954/1‑150 (MQ=255)
                                                                                                                                       tgggCCGTACGGCCAAGGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACAGGATGATGATGCCCACCACCATCAGCGACAGGTGCTGCTTGATGAATGGCACGTTGCCGAAGAAGTAGCCCAAGGTGACCAGGCCGCCGAcc               >  1:1895814/3‑150 (MQ=255)
                                                                                                                                           ccGTACGGCCCAGGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACAGGATGATGATGCCCACCACCATCAGCGACAGGTGCTGCTTGatggatgg                                                          >  2:1632779/1‑103 (MQ=255)
                                                                                                                                           ccGTACGGCCCAGGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACAGGATGATGATGCCCACCACCATCAGCGACAGGTGCTGCTTGatggatgg                                                          <  1:1632779/103‑1 (MQ=255)
                                                                                                                                           ccGTACGGCCCAGGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACAGGATGATGATGCCCACCACCATCAGCGACAGGTGCTGCTTGATGAATGGCACGTTGCCGAAGAAGTAGCCCAAGGTGACCAGGCCGCCGACCCaaa           >  1:1776648/1‑150 (MQ=255)
                                                                                                                                            cGTACGGCCCAGGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACAGGATGATGATGCCCACCACCATCAGCGACAGGTGCTGCTTGATGAATGGCACGTTGCCGAAGAAGTAGCCCAAGGTGACCAGGCCGCCGACCCAAag          <  2:1672169/150‑1 (MQ=255)
                                                                                                                                               aCGGCCCAGGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACAGGATGATGATGCCCACCACCATCAGCGACAGGTGCTGCTTGATGAATGGCACGTTGCCGAAGAAGTAGCCCAAGGTGACCAGGCCGCCGACCCAAagcag       <  2:2371468/150‑1 (MQ=255)
                                                                                                                                               aCGGCCCAGGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACAGGATGATGATGCCCACCACCATCAGCGACAGGTGCTGCTTGATGAATGGCACGTTGCCGAAGAAGTAGCCCAAGGTGACCAGGCCGCCGACCCAAagcag       <  2:2341037/150‑1 (MQ=255)
                                                                                                                                                    ccAGGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACAGGATGATGATGCCCACCACCATCAGCGACAGGTGCTGCTTGATGAATGGCACGTTGCCGAAGAAGTAGCCCAAGGTGACCAGGCCGCCGACCCAAAGCAGCGAAc  <  1:700146/150‑1 (MQ=255)
                                                                                                                                                    ccAGGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACAGGATGATGATGCCCACCACCATCAGCGACAGGTGCTGCTTGATGAATGGCACGTTGCCGAAGAAGTAGCCCAAGGTGACCAGGCCGCCGACCCAAAGCAGCGAAc  <  2:1776648/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
GCCGTCCAGCATTGGCAGCCGCTCGCGTACGGCCTGCCATTGCTGTGGCGTGATCGGGTAACGCGCCAGAGTACGCTTGCGCCGCCAGGCGCTGAATGACCACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCAGGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACAGGATGATGATGCCCACCACCATCAGCGACAGGTGCTGCTTGATGAATGGCACGTTGCCGAAGAAGTAGCCCAAGGTGACCAGGCCGCCGACCCAAAGCAGCGAAC  >  NC_002947/626423‑626720

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: