Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,289,043 C→A 100% W8L (TGG→TTG)  ccmD ← cytochrome c biogenesis protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,289,0430CA96.8% 570.1 / 9.4 189W8L (TGG→TTG) ccmDcytochrome c biogenesis protein
Reads supporting (aligned to +/- strand):  ref base C (4/2);  new base A (82/101);  total (86/103)
Fisher's exact test for biased strand distribution p-value = 4.14e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.25e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 194 total aligned reads displayed.

GCCCGCTGTTGCGCCACGCCACGCAGAATTGCGCGATGTTGCATCACCGAGTGCACGACCAAAACCACCAGCGGAATAACGGTCATCACCACCGCCAGCCAGACAAAAAAGGCGTAACCGCCCATTGCGAAAAATTCATTCCAGGAAGCAAATGCAGGGGTCATTTACGGCCTCTTTTCAGTATCAGTTCACTCACCCACGGACGGCGTTTTTCCATCAGCAAAATCAAATTACGCATCCGCATCAGCGTCAGCGTGGCAGACAGGAGCAGGAAGCCAAAAATCG  >  CP009273/2288902‑2289186
                                                                                                                                             |                                                                                                                                               
gcCCGCTGTTGCGCCACGCCACGCAGAATTGCGCGATGTTGCATCACCGAGTGCACGACCAAAACCACCAGCGGAATAACGGTCATCACCACCGCCAGCCAGACAAAAAAGGCGTAACCGCCCATTGCGAAAAATTCATTCAAGGAAg                                                                                                                                           <  2:240898/148‑1 (MQ=255)
gcCCGCTGTTGCGCCACGCCACGCAGAATTGCGCGATGTTGCATCACCGAGTGCACGACCAAAACCACCAGCGGAATAACGGTCATCACCACCGCCAGCCAGACAAAAAAGGCGTAACCGCCCATTGCGAAAAATTCATTCAAGGAAg                                                                                                                                           <  1:2933734/148‑1 (MQ=255)
   cGCTGTTGCGCCACGCCACGCAGAATTGCGCGATGTTGCATCACCGAGTGCACGACCAAAACCACCAGCGGAATAACGGTCATCACCACCGCCAGCCAGACAAAAAAGGCGTAACCGCCCATTGCGAAAAATTCATTCAAGGAAGCaa                                                                                                                                        >  1:3119882/1‑148 (MQ=255)
    gCTGTTGCGCCACGCCACGCAGAATTGCGCGATGTTGCATCACCGAGTGCACGACCAAAACCACCAGCGGAATAACGGTCATCACCACCGCCAGCCAGACAAAAAAGGCGTAACCGCCCATTGCGAAAAATTCATTCAAGGAAGCaaa                                                                                                                                       <  2:1569057/148‑1 (MQ=255)
      tGTTGCGCCACGCCACGCAGAATTGCGCGATGTTGCATCACCGAGTGCACGACCAAAACCACCAGCGGAATAACGGTCATCACCACCGCCAGCCAGACAAAAAAGGCGTAACCGCCCATTGCGAAAAATTCATTCAAGGAAGCAAATg                                                                                                                                     <  1:91263/148‑1 (MQ=255)
        ttGCGCCACGCCACGCAGAATTGCGCGATGTTGCATCACCGAGTGCACGACCAAAACCACCAGCGGAATAACGGTCATCACCACCGCCAGCCAGACAAAAAAGGCGTAACCGCCCATTGCGTAAAATTCATTCAAGGAAGCAAATGCa                                                                                                                                   <  2:1720762/148‑1 (MQ=255)
        ttGCGCCACGCCACGCAGAATTGCGCGATGTTGCATCACCGAGTGCACGACCAAAACCACCAGCGGAATAACGGTCATCACCACCGCCAGCCAGACAAAAAAGGCGTAACCGCCCATTGCGAAAAATTCATTCAAGGAAGCAAATGCa                                                                                                                                   <  1:2698831/148‑1 (MQ=255)
        ttGCGCCACGCCACGCAGAATTGCGCGATGTTGCATCACCGAGTGCACGACCAAAACCACCAGCGGAATAACGGTCATCACCACCGCCAGCCAGACAAAAAAGGCGTAACCGCCCATTGCGAAAAATTCATTCAAGGAAGCAAATGCa                                                                                                                                   <  2:1992783/148‑1 (MQ=255)
         tGCGCCACGCCACGCAGAATTGCGCGATGTTGCATCACCGAGTGCACGACCAAAACCACCAGCGGAATAACGGTCATCACCACCGCCAGCCAGACAAAAAAGGCGTAACCGCCCATTGCGAAAAATTCATTCAAGGAAGCAAATGCAg                                                                                                                                  <  1:2595429/148‑1 (MQ=255)
           cgccacgccaCGCAGAATTGCGCGATGTTGCATCACCGAGTGCACGACCAAAACCACCAGCGGAATAACGGTCATCACCACCGCCAGCCAGACAAAAAAGGCGTAACCGCCCATTGCGAAAAATTCATTCAAGGAAGCAAATGCAggg                                                                                                                                <  1:118832/148‑1 (MQ=255)
             ccacgccacgCAGAATTGCGCGATGTTGCATCACCGAGTGCACGACCAAAACCACCAGCGGAATAACGGTCATCACCACCGCCAGCCAGACAAAAAAGGCGTAACCGCCCATTGCGAAAAATTCATTCAAGGAAGCAAATGCAGGGGt                                                                                                                              <  2:2145958/148‑1 (MQ=255)
             ccacgccacgCAGAATTGCGCGATGTTGCATCACCGAGTGCACGACCAAAACCACCAGCGGAATAACGGTCATCACCACCGCCAGCCAGACAAAAAAGGCGTAACCGCCCATTGCGAAAAATTCATTCAAGGAAGCAAATGCAGGGGt                                                                                                                              <  1:2059192/148‑1 (MQ=255)
             ccacgccacgCAGAATTGCGCGATGTTGCATCACCGAGTGCACGACCAAAACCACCAGCGGAATAACGGTCATCACCACCGCCAGCCAGACAAAAAAGGCGTAACCGCCCATTGCGAAAAATTCATTCAAGGAAGCAAATGCAGGGGt                                                                                                                              <  2:2146009/148‑1 (MQ=255)
               acgccacgcAGAATTGCGCGATGTTGCATCACCGAGTGCACGACCAAAACCACCAGCGGAATAACGGTCATCACCACCGCCAGCCAGACAAAAAAGGCGTAACCGCCCATTGCGAAAAATTCATTCAAGGAAGCAAATGCAGGGGTCa                                                                                                                            <  2:1026379/148‑1 (MQ=255)
                cgccacgcAGAATTGCGCGATGTTGCATCACCGAGTGCACGACCAAAACCACCAGCGGAATAACGGTCATCACCACCGCCAGCCAGACAAAAAAGGCGTAACCGCCCATTGCGAAAAATTCATTCAAGGAAGCAAATGCAGGGGTCAt                                                                                                                           <  1:746988/148‑1 (MQ=255)
                 gccacgcAGAATTGCGCGATGTTGCATCACCGAGTGCACGACCAAAACCACCAGCGGAATAACGGTCATCACCACCGCCAGCCAGACAAAAAAGGCGTAACCGCCCATTGCGAAAAATTCATTCAAGGAAGCAAATGCAGGGGTCAtt                                                                                                                          <  2:1073245/148‑1 (MQ=255)
                 gccacgcAGAATTGCGCGATGTTGCATCACCGAGTGCACGACCAAAACCACCAGCGGAATAACGGTCATCACCACCGCCAGCCAGACAAAAAAGGCGTAACCGCCCATTGCGAAAAATTCATTCAAGGAAGCAAATGCAGGGGTCAtt                                                                                                                          <  1:145856/148‑1 (MQ=255)
                 gccacgcAGAATTGCGCGATGTTGCATCACCGAGTGCACGACCAAAACCACCAGCGGAATAACGGTCATCACCACCGCCAGCCAGACAAAAAAGGCGTAACCGCCCATTGCGAAAAATTCATTCAAGGAAGCAAATGCAGGGGTCAtt                                                                                                                          <  2:495888/148‑1 (MQ=255)
                       cAGAATTGCGCGATGTTGCATCACCGAGTGCACGACCAACACCACCAGCGGAATAACGGTCATCACCACCGCCAGCCAGACAAAAAAGGCGTAACCGCCCATTGCGAAAAATTCATTCAAGGAAGCAAATGCAGGGGTCATTTACGGc                                                                                                                    <  1:588552/148‑1 (MQ=255)
                         gAATTGCGCGATGTTGCATCACCGAGTGCACGACCAAAACCACCAGCGGAATAACGGTCATCACCACCGCCAGCCAGACAAAAAAGGCGTAACCGCCCATTGCGAAAAATTCATTCAAGGAAGCAAATGCAGGGGTCATTTACGGCct                                                                                                                  <  2:2609969/148‑1 (MQ=255)
                          aaTTGCGCGATGTTGCATCACCGAGTGCACGACCAAAACCACCAGCGGAATAACGGTCATCACCACCGCCAGCCAGACAAAAAAGGCGTAACCGCCCATTGCGAAAAATTCATTCAAGGAAGCAAATGCAGGGGTCATTTACGGCctc                                                                                                                 >  2:3255702/1‑148 (MQ=255)
                           aTTGCGCGATGTTGCATCACCGAGTGCACGACCAAAACCACCAGCGGAATAACGGTCATCACCACCGCCAGCCAGACAAAAAAGGCGTAACCGCCCATTGCGAAAAATTCATTCAAGGAAGCAAATGCAGGGGTCATTTACGGCctct                                                                                                                <  1:2418073/148‑1 (MQ=255)
                           aTTGCGCGATGTTGCATCACCGAGTGCACGACCAAAACCACCAGCGGAATAACGGTCATCACCACCGCCAGCCAGACAAAAAAGGCGTAACCGCCCATTGCGAAAAATTCATTCAAGGAAGCAAATGCAGGGGTCATTTACGGCctct                                                                                                                <  1:759775/148‑1 (MQ=255)
                            ttGCGCGATGTTGCATCACCGAGTGCACGACCAAAACCACCAGCGGAATAACGGTCATCACCACCGCCAGCCAGACAAAAAAGGCGTAACCGCCCATTGCGAAAAATTCATTCAAGGAAGCAAATGCAGGGGTCATTTACGGCCTCtt                                                                                                               >  2:638956/1‑148 (MQ=255)
                            ttGCGCGATGTTGCATCACCGAGTGCACGACCAAAACCACCAGCGGAATAACGGTCATCACCACCGCCAGCCAGACAAAAAAGGCGTAACCGCCCATTGCGAAAAATTCATTCAAGGAAGCAAATGCAGGGGTCATTTACGGCCTCtt                                                                                                               <  2:838029/148‑1 (MQ=255)
                             tGCGCGATGTTGCATCACCGAGTGCACGACCAAAACCACCAGCGGAATAACGGTCATCACCACCGCCAGCCAGACAAAAAAGGCGTAACCGCCCAGTGCGAAAAATTCATTCAAGGAAGCAAATGCAGGGGTCATTTACGGCCTCttt                                                                                                              >  2:2426747/1‑148 (MQ=255)
                                 cgATGTTGCATCACCGAGTGCACGACCAAAACCACCAGCGGAATAACGGTCATCACCACCGCCAGCCAGACAAAAAAGGCGTAACCGCCCATTGCGAAAAATTCATTCAAGGAAGCAAATGCAGGGGTCATTTACGGCCTCTTTTCAg                                                                                                          <  1:2027049/148‑1 (MQ=255)
                                 cgATGTTGCATCACCGAGTGCACGACCAAAACCACCAGCGGAATAACGGTCATCACCACCGCCAGCCAGACAAAAAAGGCGTAACCGCCCATTGCGAAAAATTCATTCAAGGAAGCAAATGCAGGGGTCATTTACGGCCTCTTTTCAg                                                                                                          >  2:2769207/1‑148 (MQ=255)
                                    tGTTGCATCACCGAGTGCACGACCAAAACCACCAGCGGAATAACGGTCATCACCACCGCCAGCCAGACAAAAAAGGCGTAACCGCCCATTGCGAAAAATTCATTCAAGGAAGCAAATGCAGGGGTCATTTACGGCCTCTTTTCAGtct                                                                                                       >  2:2829326/1‑146 (MQ=255)
                                     gTTGCATCACCGAGTGCACGACCAAAACCACCAGCGGAATAACGGTCATCACCACCGCCAGCCAGACAAAAAAGGCGTAACCGCCCATTGCGAAAAATTCATTCAAGGAAGCAAATGCAGGGGTCATTTACGGCCTCTTTTCAGTATc                                                                                                      <  2:178353/148‑1 (MQ=255)
                                     gTTGCATCACCGAGTGCACGACCAAAACCACCAGCGGAATAACGGTCATCACCACCGCCAGCCAGACAAAAAAGGCGTAACCGCCCATTGCGAAAAATTCATTCAAGGAAGCAAATGCAGGGGTCATTTACGGCCTCTTTTCAGTATc                                                                                                      <  1:461702/148‑1 (MQ=255)
                                      ttGCATCACCGAGTGCACGACCAAAACCACCAGCGGAATAACGGTCATCACCACCGCCAGCCAGACAAAAAAGGCGTAACCGCCCATTGCGAAAAATTCATTCAAGGAAGCAAATGCAGGGGTCATTTACGGCCTCTTTTCAGTATCa                                                                                                     >  2:385998/1‑148 (MQ=255)
                                           tCACCGAGTGCACGACCAAAACCACCAGCGGAATAACGGTCATCACCACCGCCAGCCAGACAAAAAAGGCGTAACCGCCCATTGCGAAAAATTCATTCAAGGAAGCAAATGCAGGGGTCATTTACGGCCTCTTTTCAGTATCAGTtca                                                                                                <  2:3251633/148‑1 (MQ=255)
                                           tCACCGAGTGCACGACCAAAACCACCAGCGGAATAACGGTCATCACCACCGCCAGCCAGACAAAAAAGGCGTAACCGCCCATTGCGAAAAATTCATTCAAGGAAGCAAATGCAGGGGTCATTTACGGCCTCTTTTCAGTATCAGTtca                                                                                                <  2:223916/148‑1 (MQ=255)
                                             aCCGAGTGCACGACCAAAACCACCAGCGGAATAACGGTCATCACCACCGCCAGCCAGACAAAAAAGGCGTAACCGCCCATTGCGAAAAATTCATTCAAGGAAGCAAATGCAGGGGTCATTTACGGCCTCTTTGCAGTATCAGTtcact                                                                                              >  1:1654915/1‑148 (MQ=255)
                                                gAGTGCACGACCAAAACCACCAGCGGAATAACGGTCATCACCACCGCCAGCCAGACAAAAAAGGCGTAACCGCCCATTGCGAAAAATTCATTCAAGGAAGCAAATGCAGGGGTCATTTACGGCCTCTTTTCAGTATCAGTtcactcac                                                                                           <  2:1449812/148‑1 (MQ=255)
                                                   tGCACGACCAAAACCACCAGCGGAATAACGGTCATCACCACCGCCAGCCAGACAAAAAAGGCGTAACCGCCCATTGCGAAAAATTCATTCAAGGAAGCAAATGCAGGGGTCATTTACGGCCTCTTTTCAGTATCAGTTCACTCACCCa                                                                                        >  1:1252103/1‑148 (MQ=255)
                                                   tGCACGACCAAAACCACCAGCGGAATAACGGTCATCACCACCGCCAGCCAGACAAAAAAGGCGTAACCGCCCATTGCGAAAAATTCATTCAAGGAAGCAAATGCAGGGGTCATTTACGGCCTCTTTTCAGTATCAGTTCACTCACCCa                                                                                        >  2:1557906/1‑148 (MQ=255)
                                                   tGCACGACCAAAACCACCAGCGGAATAACGGTCATCACCACCGCCAGCCAGACAAAAAAGGCGTAACCGCCCATTGCGAAAAATTCATTCAAGGAAGCAAATGCAGGGGTCATTTACGGCCTCTTTTCAGTATCAGTTCACTCACCCa                                                                                        >  1:130079/1‑148 (MQ=255)
                                                   tGCACGACCAAAACCACCAGCGGAATAACGGTCATCACCACCGCCAGCCAGACAAAAAAGGCGTAACCGCCCATTGCGAAAAATTCATTCAAGGAAGCAAATGCAGGGGTCATTTACGGCCTCTTTTCAGTATCAGTTCACTCACCCa                                                                                        >  2:3198285/1‑148 (MQ=255)
                                                      aCGACCAAAACCACCAGCGGAATAACGGTCATCACCACCGCCAGCCAGACAAAAAAGGCGTAACCGCCCATTGCGAAAAATTCATTCAAGGAAGCAAATGCAGGGGTCATTTACGGCCTCTTTTCAGTATCAGTTCACTCACCCacgg                                                                                     <  1:3238083/148‑1 (MQ=255)
                                                          ccAAAACCACCAGCGGAATAACGGTCATCACCACCGCCAGCCAGACAAAAAAGGCGTAACCGCCCATTGCGAAAAATTCATTCAAGGAAGCAAATGCAGGGGTCATTTACGGCCTCTTTTCAGTATCAGTTCACTCACCCacggacgg                                                                                 >  2:2622633/1‑148 (MQ=255)
                                                            aaaaCCACCAGCGGAATAACGGTCATCACCACCGCCAGCCAGACAAAAAAGGCGTAACCGCCCATTGCGAAAAATTCATTCAAGGAAGCAAATGCAGGGGTCATTTACGGCCTCTTTTCAGTATCAGTTCACTCACCCACGGACGGCg                                                                               <  1:1385898/148‑1 (MQ=255)
                                                             aaaCCACCAGCGGAATAACGGTCATCACCACCGCCAGCCAGACAAAAAAGGCGTAACCGCCCATTGCGAAAAATTCATTCAAGGAAGCAAATGCAGGGGTCATTTACGGCCTCTTTTCAGTATCAGTTCACTCACCCACGGACGGCGt                                                                              <  2:1380460/148‑1 (MQ=255)
                                                             aaaCCACCAGCGGAATAACGGTCATCACCACCGCCAGCCAGACAAAAAAGGCGTAACCGCCCATTGCGAAAAATTCATTCAAGGAAGCAAATGCAGGGGTCATTTACGGCCTCTTTTCAGTATCAGTTCACTCACCCACGGACGGCGt                                                                              >  2:1004249/1‑148 (MQ=255)
                                                              aaCCACCAGCGGAATAACGGTCATCACCACCGCCTGCCAGACAAAAAAGGCGTAACCGCACATTGCGAAAAATTCATTCAAGGAAGCAATTGCAGGGGTCATTTACGGCCTCTTTTCAGTATCAGATCACTCACCCACGGACGGCGtt                                                                             >  2:1544738/1‑148 (MQ=255)
                                                              aaCCACCAGCGGAATAACGGTCATCACCACCGCCAGCCAGACAAAAAAGGCGTAACCGCCCATTGCGAAAAATTCATTCAAGGAAGCAAATGCAGGGGTCATTTACGGCCTCTTTTCAGTATCAGTTCACTCACCCACGGACGGCGtt                                                                             >  2:2515334/1‑148 (MQ=255)
                                                                 caccaGCGGAATAACGGTCATCACCACCGCCAGCCAGACAAAAAAGGCGTAACCGCCCATTGCGAAAAATTCATTCAAGGAAGCAAATGCAGGGGTCATTTACGGCCTCTTTTCAGTATCAGTTCACTCACCCACGGACGGCGttttt                                                                          <  1:963581/148‑1 (MQ=255)
                                                                 caccaGCGGAATAACGGTCATCACCACCGCCAGCCAGACAAAAAAGGCGTAACCGCCCATTGCGAAAAATTCATTCAAGGAAGCAAATGCAGGGGTCATTTACGGCCTCTTTTCAGTATCAGTTCACTCACCCACGGACGGCGttttt                                                                          <  1:2618400/148‑1 (MQ=255)
                                                                   ccaGCGGAATAACGGTCATCACCACCGCCAGCCAGACAAAAAAGGCGTAACCGCCCATTGCGAAAAATTCATTCAAGGAAGCAAATGCAGGGGTCATTTACGGCCTCTTTTCAGTATCAGTTCACTCACCCACGGACGGCGTTTTTcc                                                                        >  2:1474940/1‑148 (MQ=255)
                                                                    caGCGGAATAACTGTCATCACCAACGCCAGCCAGACAAAAAAGGCGTAACCGCCCATTGCGAAAAATTCATTCAAGGAAGCAAATGCAGGGGTCATTTACGGCCTCTTTTCAGTATCAGTTCACTCACCCACGGACGGCGTTTTTCCa                                                                       <  2:2213374/148‑1 (MQ=255)
                                                                      gCGGAATAACGGTCATCACCACCGCCAGCCAGACAAAAAAGGCGTAACCGCCCATTGCGAAAAATTCATTCAAGGAAGCAAATGCAGGGGTCATTTACGGCCTCTTTTCAGTATCAGTTCACTCACCCACGGACGGCGTTTTTCCATc                                                                     >  1:1889254/1‑148 (MQ=255)
                                                                       cGGAATAACGGTCATCACCACCGCCAGCCAGACAAAAAAGGCGTAACCGCCCATTGCGAAAAATTCATTCAAGGAAGCAAATGCAGGGGTCATTTACGGCCTCTTTTCAGTATCAGTTCACTCACCCACGGACGGCGTTTTTCCATCa                                                                    >  2:1927430/1‑148 (MQ=255)
                                                                        ggAATAACGGTCATCACCACCGCCAGCCAGACAAAAAAGGCGTAACCGCCCATTGCGAAAAATTCATTCAAGGAAGCAAATGCAGGGGTCATTTACGGCCTCTTTTCAGTATCAGTTCACTCACCCACGGACGGCGTTTTTCCATCAg                                                                   <  1:655499/148‑1 (MQ=255)
                                                                        ggAATAACGGTCATCACCACCGCCAGCCAGACAAAAAAGGCGTAACCGCCCATTGCGAAAAATTCATTCAAGGAAGCAAATGCAGGGGTCATTTACGGCCTCTTTTCAGTATCAGTTCACTCACCCACGGACGGCGTTTTTCCATCAg                                                                   <  1:654583/148‑1 (MQ=255)
                                                                         gAATAACGGTCATCACCACCGCCAGCCAGACAAAAAAGGCGTAACCGCCCATTGCGAAAAATTCATTCAAGGAAGCAAATGCAGGGGTCATTTACGGCCTCTTTTCAGTATCAGTTCACTCACCCACGGACGGCGTTTTTCCATCAGc                                                                  >  2:317385/1‑148 (MQ=255)
                                                                          aaTAACGGTCATCACCACCGCCAGCCAGACAAAAAAGGCGTAACCGCCCATTGCGAAAAATTCATTCCAGGAAGCAAATGCAGGGGTCATTTACGGCCTCTTTTCAGTATCAGTTCACTCACCCACGGACGGCGTTTTTCCATCAGCa                                                                 >  1:2426060/1‑148 (MQ=255)
                                                                           aTAACGGTCATCACCACCGCCAGCCAGACAAAAAAGGCGTAACCGCCCATTGCGAAAAATTCATTCAAGGAAGCAAATGCAGGGGTCATTTACGGCCTCTTTTCAGTATCAGTTCACTCACCCACGGACGGCGTTTTTCCATCAGCaa                                                                >  2:1230150/1‑148 (MQ=255)
                                                                            tAACGGTCATCACCACCGCCAGCCAGACAAAAAAGGCGTAACCGCCCATTGCGAAAAATTCATTCAAGGAAGCAAATGCAGGGGTCATTTACGGCCTCTTTTCAGTATCAGTTCACTCACCCACGGACGGCGTTTTTCCATCAGCaaa                                                               <  1:56878/148‑1 (MQ=255)
                                                                             aaCGGTCATCACCACCGCCAGCCAGACAAAAAAGGCGTAACCGCCCATTGCGAAAAATTCATTCAAGGAAGCAAATGCAGGGGTCATTTACGGCCTCTTTTCAGTATCAGTTCACTCACCCACGGACGGCGTTTTTCCATCAGCaaaa                                                              <  1:3092861/148‑1 (MQ=255)
                                                                                 gTCATCACCACCGCCAGCCAGACAAAAAAGGCGTAACCGCCCATTGCGAAAAATTCATTCAAGGAAGCAAATGCAGGGGTCATTTACGGCCTCTTTTCAGTATCAGTTCACTCACCCACGGACGGCGTTTTTCCATCAGCAAAATCaa                                                          <  1:2921227/148‑1 (MQ=255)
                                                                                  tcatcaCCACCGCCAGCCAGACAAAAAAGGCGTAACCGCCCATTGCGAAAAATTCATTCAAGGAAGCAAATGCAGGGGTCATTTACGGCCTCTTTTCAGTATCAGTTCACTCACCCACGGACGGCGTTTTTCCATCAGCAAAATCaaa                                                         <  1:990483/148‑1 (MQ=255)
                                                                                   catcaCCACCGCCAGCCAGACAAAAAAGGCGTAACCGCCCATTGCGAAAAATTCATTCAAGGAAGCAAATGCAGGGGTCATTTACGGCCTCTTTTCAGTATCAGTTCACTCACCCACGGACGGCGTTTTTCCATCAGCAAAATCAAAt                                                        <  2:442599/148‑1 (MQ=255)
                                                                                     tcaCCACCGCCAGCCAGACAAAAAAGGCGTAACCGCCCATTGCGAAAAATTCATTCAAGGAAGCAAATGCAGGGGTCATTTACGGCCTCTTTTCAGTATCAGTTCACTCACCCACGGACGGCGTTTTTCCATCAGCAAAATCAAATTa                                                      >  1:463345/1‑148 (MQ=255)
                                                                                     tcaCCACCGCCAGCCAGACAAAAAAGGCGTAACCGCCCATTGCGAAAAATTCATTCAAGGAAGCAAATGCAGGGGTCATTTACGGCCTCTTTTCAGTATCAGTTCACTCACCCACGGACGGCGTTTTTCCATCAGCAAAATCAAATTa                                                      <  1:1980381/148‑1 (MQ=255)
                                                                                           ccGCCAGCCAGACAAAAAAGGCGTAACCGCCCATTGCGAAAAATTCATTCAAGGAAGCAAATGCAGGGGTCATTTACGGCCTCTTTTCAGTATCAGTTCACTCACCCACGGACGGCGTTTTTCCATCAGCAAAATCAAATTAcgcatc                                                >  1:1429656/1‑148 (MQ=255)
                                                                                                  ccagACAAAAAAGGCGTAACCGCCCATTGCGAAAAATTCATTCAAGGAAGCAAATGCAGGGGTCATTTACGGCCTCTTTTCAGTATCAGTTCACTCACCCACGGACGGCGTTTTTCCATCAGCAAAATCAAATTACGCATCCGCAtca                                         >  2:1649160/1‑148 (MQ=255)
                                                                                                    agACAAAAAAGGCGTAACCGCCCATTGCGAAAAATTCATTCAAGGAAGCAAATGCAGGGGTCATTTACGGCCTCTTTTCAGTATCAGTTCACTCACCCACGGACGGCGTTTTTCCATCAGCAAAATCAAATTACGCATCCGCAtcagc                                       >  2:3070055/1‑148 (MQ=255)
                                                                                                     gACAAAAAAGGCGTAACCGCCCATTGCGAAAAATTCATTCAAGGAAGCAAATGCAGGGGTCATTTACGGCCTCTTTTCAGTATCAGTTCACTCACCCACGGACGGCGTTTTTCCATCAGCAAAATCAAATTACGCATCCGCAtcagcg                                      <  1:2482693/148‑1 (MQ=255)
                                                                                                     gACAAAAAAGGCGTAACCGCCCATTGCGAAAAATTCATTCAAGGAAGCAAATGCAGGGGTCATTTACGGCCTCTTTTCAGTATCAGTTCACTCACCCACGGACGGCGTTTTTCCATCAGCAAAATCAAATTACGCATCCGCAtcagcg                                      <  1:1116385/148‑1 (MQ=255)
                                                                                                         aaaaaGGCGTAACCGCCCATTGCGAAAAATTCATTCAAGGAAGCAAATGCAGGGGTCATTTACGGCCTCTTTTCAGTATCAGTTCACTCACCCACGGACGGCGTTTTTCCATCAGCAAAATCAAATTACGCATCCGCAtcagcgtcag                                  <  1:40740/148‑1 (MQ=255)
                                                                                                          aaaaGGCGTAACCGCCCATTGCGAAAAATTCATTCAAGGAAGCAAATGCAGGGGTCATTTACGGCCTCTTTTCAGTATCAGTTCACTCACCCACGGACGGCGTTTTTCCATCAGCAAAATCAAATTACGCATCCGCAtcagcgtcagc                                 >  2:2641618/1‑148 (MQ=255)
                                                                                                           aaaGGCGTAACCGCCCATTGCGAAAAATTCATTCAAGGAAGCAAATGCAGGGGTCATTTACGGCCTCTTTTCAGTATCAGTTCACTCACCCACGGACGGCGTTTTTCCATCAGCAAAATCAAATTACGCATCCGCAtcagcgtcagcg                                >  1:2391300/1‑148 (MQ=255)
                                                                                                           aaaGGCGTAACCGCCCATTGCGAAAAATTCATTCAAGGAAGCAAATGCAGGGGTCATTTACGGCCTCTTTTCAGTATCAGTTCACTCACCCACGGACGGCGTTTTTCCATCAGCAAAATCAAATTACGCATCCGCAtcagcgtcagcg                                <  2:3007896/148‑1 (MQ=255)
                                                                                                            aaGGCGTAACCGCCCATTGCGAAAAATTCATTCAAGGAAGCAAATGCAGGGGTCATTTACGGCCTCTTTTCAGTATCAGTTCACTCACCCACGGACGGCGTTTTTCCATCAGCAAAATCAAATTACGCATCCGCATcagcgtcagcgt                               >  2:1191491/1‑148 (MQ=255)
                                                                                                              ggCGTAACCGCCCATTGCGAAAAATTCAATCAAGGAAGCAAATGCAGGGGTCATGTACGGCCTCTTGTCAGTATCAGTGCACTCACCCACGGACGGCGTTTTTCCATCAGCAAAATCACATTACGCATCCGCATCAGCGTCAGCGTgg                             >  2:279528/1‑148 (MQ=255)
                                                                                                                  tAACCGCCCATTGCGAAAAATTCATTCCAGGAAGCAAATGCAGGGGTCATTTACGGCCTCTTTTCAGTATCAGTTCACTCACCCACGGACGGCGTTTTTCCATCAGCAAAATCAAATTACGCATCCGCATCAGCGTCAGCGTGGCAGa                         >  1:1729840/1‑148 (MQ=255)
                                                                                                                    aCCGCCCATTGCGAAAAATTCATTCAAGGAAGCAAATGCAGGGGTCATTTACGGCCTCTTTTCAGTATCAGTTCACTCACCCACGGACGGCGTTTTTCCATCAGCAAAATCAAATTACGCATCCGCATCAGCGTCAGCGTGGCAGACa                       >  2:2474369/1‑148 (MQ=255)
                                                                                                                     ccGCCCATTGCGAAAAATTCATTCAAGGAAGCAAATGCAGGGGTCATTTACGGCCTCTTTTCAGTATCAGTTCACTCACCCACGGACGGCGTTTTTCCATCAGCAAAATCAAATTACGCATCCGCATCAGCGTCAGCGTGGCAGACAg                      >  2:2324534/1‑148 (MQ=255)
                                                                                                                      cGCCCATTGCGAAAAATTCATTCAAGGAAGCACATGCAGGGGTCATTTACGGCCTCTTTTCAGTATCAGTTCACTCACCCACGGACGGCGTTTTTCCATCAGCAAAATCAAATTACGCATCCGCATCAGCGTCAGCGTGGCAGACAgg                     <  2:245849/148‑1 (MQ=255)
                                                                                                                      cGCCCATTGCGAAAAATTCATTCAAGGAAGCAAATGCAGGGGTCATTTACGGCCTCTTTTCAGTATCAGTTCACTCACCCACGGACGGCGTTTTTCCATCAGCAAAATCAAATTACGCATCCGCATCAGCGTCAGCGTGGCAGACAgg                     >  1:1561306/1‑148 (MQ=255)
                                                                                                                      cGCCCATTGCGAAAAATTCATTCAAGGAAGCAAATGCAGGGGTCATTTACGGCCTCTTTTCAGTATCAGTTCACTCACCCACGGACGGCGTTTTTCCATCAGCAAAATCAAATTACGCATCCGCATCAGCGTCAGCGTGGCAGACAgg                     <  2:245864/148‑1 (MQ=255)
                                                                                                                      cGCCCATTGCGAAAAATTCATTCAAGGAAGCAAATGCAGGGGTCATTTACGGCCTCTTTTCAGTATCAGTTCACTCACCCACGGACGGCGTTTTTCCATCAGCAAAATCAAATTACGCATCCGCATCAGCGTCAGCGTGGCAGACAgg                     >  1:696075/1‑148 (MQ=255)
                                                                                                                       gCCCATTGCGAAAAATTCATTCAAGGAAGCAAATGCAGGGGTCATTTACGGCCTCTTTTCAGTATCAGTTCACTCACCCACGGACGGCGTTTTTCCATCAGCAAAATCAAATTACGCATCCGCATCAGCGTCAGCGTGGCAGACAGGa                    <  1:188494/148‑1 (MQ=255)
                                                                                                                       gCCCATTGCGAAAAATTCATTCAAGGAAGCAAATGCAGGGGTCATTTACGGCCTCTTTTCAGTATCAGTTCACTCACCCACGGACGGCGTTTTTCCATCAGCAAAATCAAATTACGCATCCGCATCAGCGTCAGCGTGGCAGACAGGa                    <  1:2979483/148‑1 (MQ=255)
                                                                                                                         ccATTGCGAAAAATTCATTCAAGGAAGCAAATGCAGGGGTCATTTACGGCCTCTTTTCAGTATCAGTTCACTCACCCACGGACGGCGTTTTTCCATCAGCAAAATCAAATTACGCATCCGCATCAGCGTCAGCGTGGCAGACAGGAGc                  >  1:2224664/1‑148 (MQ=255)
                                                                                                                          cATTGCGAAAAATTCATTCAAGGAAGCAAATGCAGGGGTCATTTACGGCCTCTTTTCAGTATCAGTTCACTCACCCACGGACGGCGTTTTTCCATCAGCAAAATCAAATTACGCATCCGCATCAGCGTCAGCGTGGCAGACAGGAGCa                 >  2:386468/1‑148 (MQ=255)
                                                                                                                           aTTGCGAAAAATTCATTCAAGGAAGCAAATGCAGGGGTCATTTACGGCCTCTTTTCAGTATCAGTTCACTCACCCACGGACGGCGTTTTTCCATCAGCAAAATCAAATTACGCATCCGCATCAGCGTCAGCGTGGCAGACAGGAGCAg                >  2:676600/1‑148 (MQ=255)
                                                                                                                             tGCGAAAAATTCATTCAAGGAAGCAAATGCAGGGGTCATTTACGGCCTCTTTTCAGTATCAGTTCACTCACCCACGGACGGCGTTTTTCCATCAGCAAAATCAAATTACGCATCCGCATCAGCGTCAGCGTGGCAGACAGGAGCAGGa              >  1:274683/1‑148 (MQ=255)
                                                                                                                                gAAAAATTCATTCAAGGAAGCAAATGCAGGGGTCATTTACGGCCTCTTTTCAGTATCAGTTCACTCACCCACGGACGGCGTTTTTCCATCAGCAAAATCAAATTACGCATCCGCATCAGCGTCAGCGTGGCAGACAGGAGCAGGAAGc           <  2:162290/148‑1 (MQ=255)
                                                                                                                                 aaaaaTTCATTCAAGGAAGCAAATGCAGGGGTCATTTACGGCCTCTTTTCAGTATCAGTTCACTCACCCACGGACGGCGTTTTTCCATCAGCAAAATCAAATTACGCATCCGCATCAGCGTCAGCGTGGCAGACAGGAGCAGGAAGcc          >  2:1168890/1‑148 (MQ=255)
                                                                                                                                  aaaaTTCATTCAAGGAAGCAAATGCAGGGGTCATTTACGGCCTCTTTTCAGTATCAGTTCACTCACCCACGGACGGCGTTTTTCCATCAGCAAAATCAAATTACGCATCCGCATCAGCGTCAGCGTGGCAGACAGGAGCAGGAAGCCa         >  2:2905448/1‑148 (MQ=255)
                                                                                                                                   aaaTTCATTCAAGGAAGCAAATGCAGGGGTCATTTACGGCCTCTTTTCAGTATCAGTTCACTCACCCACGGACGGCGTTTTTCCATCAGCAAAATCAAATTACGCATCCGCATCAGCGTCAGCGTGGCAGACAGGAGCAGGAAGCCaa        <  2:2078048/148‑1 (MQ=255)
                                                                                                                                    aaTTCATTCAAGGAAGCAAATGCAGGGGTCATTTACGGCCTCTTTTCAGTATCAGTTCACTCACCCACGGACGGCGTTTTTCCATCAGCAAAATCAAATTACGCATCCGCATCAGCGTCAGCGTGGCAGACAGGAGCAGGAAGCCaaa       >  2:3091292/1‑148 (MQ=255)
                                                                                                                                     attcattcAAGGAAGCAAATGCAGGGGTCATTTACGGCCTCTTTTCAGTATCAGTTCACTCACCCACGGACGGCGTTTTTCCATCAGCAAAATCAAATTACGCATCCGCATCAGCGTCAGCGTGGCAGACAGGAGCAGGAAGCCaaaa      >  2:550725/1‑148 (MQ=255)
                                                                                                                                     attcattcAAGGAAGCAAATGCAGGGGTCATTTACGGCCTCTTTTCAGTATCAGTTCACTCACCCACGGACGGCGTTTTTCCATCAGCAAAATCAAATTACGCATCCGCATCAGCGTCAGCGTGGCAGACAGGAGCAGGAAGCCaaaa      <  1:168198/148‑1 (MQ=255)
                                                                                                                                     attcattcAAGGAAGCAAATGCAGGGGTCATTTACGGCCTCTTTTCAGTATCAGTTCACTCACCCACGGACGGCGTTTTTCCATCAGCAAAATCAAATTACGCATCCGCATCAGCGTCAGCGTGGCAGACAGGAGCAGGAAGCCaaaa      >  2:897964/1‑148 (MQ=255)
                                                                                                                                      ttcattcaAGGAAGCAAATGCAGGGGTCATTTACGGCCTCTTTTCAGTATCAGTTCACTCACCCACGGACGGCGTTTTTCCATCAGCAAAATCAAATTACGCATCCGCATCAGCGTCAGCGTGGCAGACAGGAGCAGGAAGCCaaaaa     >  1:3048255/1‑148 (MQ=255)
                                                                                                                                      ttcattcaAGGAAGCAAATGCAGGGGTCATTTACGGCCTCTTTTCAGTATCAGTTCACTCACCCACGGACGGCGTTTTTCCATCAGCAAAATCAAATTACGCATCCGCATCAGCGTCAGCGTGGCAGACAGGAGCAGGAAGCCaaaaa     >  1:1508644/1‑148 (MQ=255)
                                                                                                                                      ttcattcaAGGAAGCAAATGCAGGGGTCATTTACGGCCTCTTTTCAGTATCAGTTCACTCACCCACGGACGGCGTTTTTCCATCAGCAAAATCAAATTACGCATCCGCATCAGCGTCAGCGTGGCAGACAGGAGCAGGAAGCCaaaaa     >  1:1501323/1‑148 (MQ=255)
                                                                                                                                             |                                                                                                                                               
GCCCGCTGTTGCGCCACGCCACGCAGAATTGCGCGATGTTGCATCACCGAGTGCACGACCAAAACCACCAGCGGAATAACGGTCATCACCACCGCCAGCCAGACAAAAAAGGCGTAACCGCCCATTGCGAAAAATTCATTCCAGGAAGCAAATGCAGGGGTCATTTACGGCCTCTTTTCAGTATCAGTTCACTCACCCACGGACGGCGTTTTTCCATCAGCAAAATCAAATTACGCATCCGCATCAGCGTCAGCGTGGCAGACAGGAGCAGGAAGCCAAAAATCG  >  CP009273/2288902‑2289186

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