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,0430CA98.8% 1094.5 / 2.1 344W8L (TGG→TTG) ccmDcytochrome c biogenesis protein
Reads supporting (aligned to +/- strand):  ref base C (3/2);  new base A (142/196);  total (145/199)
Fisher's exact test for biased strand distribution p-value = 6.54e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 347 total aligned reads displayed.

TCACGCGCCCGCTGTTGCGCCACGCCACGCAGAATTGCGCGATGTTGCATCACCGAGTGCACGACCAAAACCACCAGCGGAATAACGGTCATCACCACCGCCAGCCAGACAAAAAAGGCGTAACCGCCCATTGCGAAAAATTCATTCCAGGAAGCAAATGCAGGGGTCATTTACGGCCTCTTTTCAGTATCAGTTCACTCACCCACGGACGGCGTTTTTCCATCAGCAAAATCAAATTACGCATCCGCATCAGCGTCAGCGTGGCAGACAGGAGCAGGAAGCCAAAAATCGACC  >  CP009273/2288896‑2289189
                                                                                                                                                   |                                                                                                                                                  
tCACGCGCCCGCTGTTGCGCCACGCCACGCAGAATTGCGCGATGTTGCATCACCGAGTGCACGACCAAAACCACCAGCGGAATAACGGTCATCACCACCGCCAGCCAGACAAAAAAGGCGTAACCGCCCATTGCGAAAAATTCATTCAAg                                                                                                                                                  <  2:908889/150‑1 (MQ=255)
 cACGCGCCCGCTGTTGCGCCACGCCACGCAGAATTGCGCGATGTTGCATCACCGAGTGCACGACCAAAACCACCAGCGGAATAACGGTCATCACCACCGCCAGCCAGACAAAAAAGGCGTAACCGCCCATTGCGAAAAATTCATTCAAgg                                                                                                                                                 <  2:2498417/150‑1 (MQ=255)
    gcgcCCGCTGTTGCGCCACGCCACGCAGAATTGCGCGATGTTGCATCACCGAGTGCACGACCAAAACCACCAGCGGAATAACGGTCATCACCACCGCCAGCCAGACAAAAAAGGCGTAACCGCCCATTGCGAAAAATTCATTCAAGGAAg                                                                                                                                              <  1:539824/150‑1 (MQ=255)
      gcCCGCTGTTGCGCCACGCCACGCAGAATTGCGCGATGTTGCATCACCGAGTGCACGACCAAAACCACCAGCGGAATAACGGTCATCACCACCGCCAGCCAGACAAAAAAGGCGTAACCGCCCATTGCGAAAAATTCATTCAAGGAAGCa                                                                                                                                            <  1:3404519/150‑1 (MQ=255)
      gcCCGCTGTTGCGCCACGCCACGCAGAATTGCGCGATGTTGCATCACCGAGTGCACGACCAAAACCACCAGCGGAATAACGGTCATCACCACCGCCAGCCAGACAAAAAAGGCGTAACCGCCCATTGCGAAAAATTCATTCAAGGAAGCa                                                                                                                                            <  2:4507604/150‑1 (MQ=255)
       cccGCTGTTGCGCCACGCCACGCAGAATTGCGCGATGTTGCATCACCGAGTGCACGACCAAAACCACCAGCGGAATAACGGTCATCACCACCGCCAGCCAGACAAAAAAGGCGTAACCGCCCATTGCGAAAAATTCATTCAAGGAAGCaa                                                                                                                                           >  2:3454298/1‑150 (MQ=255)
       cccGCTGTTGCGCCACGCCACGCAGAATTGCGCGATGTTGCATCACCGAGTGCACGACCAAAACCACCAGCGGAATAACGGTCATCACCACCGCCAGCCAGACAAAAAAGGCGTAACCGCCCATTGCGAAAAATTCATTCAAGGAAGCaa                                                                                                                                           <  1:900691/150‑1 (MQ=255)
        ccGCTGTTGCGCCACGCCACGCAGAATTGCGCGATGTTGCATCACCGAGTGCACGACCAAAACCACCAGCGGAATAACGGTCATCACCACCGCCAGCCAGACAAAAAAGGCGTAACCGCCCATTGCGAAAAATTCATTCAAGGAAGCaaa                                                                                                                                          <  1:1857848/150‑1 (MQ=255)
         cGCTGTTGCGCCACGCCACGCAGAATTGCGCGATGTTGCATCACCGAGTGCACGACCAAAACCACCAGCGGAATAACGGTCATCACCACCGCCAGCCAGACAAAAAAGGCGTAACCGCCCATTGCGAAAAATTCATTCAAGGAAGCAAAt                                                                                                                                         >  2:829231/1‑150 (MQ=255)
         cGCTGTTGCGCCACGCCACGCAGAATTGCGCGATGTTGCATCACCGAGTGCACGACCAAAACCACCAGCGGAATAACGGTCATCACCACCGCCAGCCAGACAAAAAAGGCGTAACCGCCCATTGCGAAAAATTCATTCAAGGAAGCAAAt                                                                                                                                         <  1:3574238/150‑1 (MQ=255)
          gCTGTTGCGCCACGCCACGCAGAATTGCGCGATGTTGCATCACCGAGTGCACGACCAAAACCACCAGCGGAATAACGGTCATCACCACCGCCAGCCAGACAAAAAAGGCGTAACCGCCCATTGCGAAAAATTCATTCAAGGAAGCAAATg                                                                                                                                        <  1:2313239/150‑1 (MQ=255)
           cTGTTGCGCCACGCCACGCAGAATTGCGCGATGTTGCATCACCGAGTGCACGACCAAAACCACCAGCGGAATAACGGTCATCACCACCGCCAGCCAGACAAAAAAGGCGTAACCGCCCATTGCGAAAAATTCATTCAAGGAAGCAAATGc                                                                                                                                       >  2:1344497/1‑150 (MQ=255)
              ttGCGCCACGCCACGCAGAATTGCGCGATGTTGCATCACCGAGTGCACGACCAAAACCACCAGCGGAATAACGGTCATCACCACCGCCAGCCAGACAAAAAAGGCGTAACCGCCCATTGCGAAAAATTCATTCAAGGAAGCAAATGCAgg                                                                                                                                    <  2:3865784/150‑1 (MQ=255)
               tGCGCCACGCCACGCAGAATTGCGCGATGTTGCATCACCGAGTGCACGACCAAAACCACCAGCGGAATAACGGTCATCACCACCGCCAGCCAGACAAAAAAGGCGTAACCGCCCATTGCGAAAAATTCATTCAAGGAAGCAAATGCAggg                                                                                                                                   <  2:4207835/150‑1 (MQ=255)
                 cgccacgccaCGCAGAATTGCGCGATGTTGCATCACCGAGTGCACGACCAAAACCACCAGCGGAATAACGGTCATCACCACCGCCAGCCAGACAAAAAAGGCGTAACCGCCCATTGCGAAAAATTCATTCAAGGAAGCAAATGCAGGGGt                                                                                                                                 >  1:4412434/1‑150 (MQ=255)
                   ccacgccacgCAGAATTGCGCGATGTTGCATCACCGAGTGCACGACCAAAACCACCAGCGGAATAACGGTCATCACCACCGCCAGCCAGACAAAAATGGCGTAACCGCCCATGGCGAAAATTTCATTCAAGGTAGAAAATGCAGGGGtaa                                                                                                                               >  2:2683779/1‑148 (MQ=255)
                   ccacgccacgCAGAATTGCGCGATGTTGCATCACCGAGTGCACGACCAAAACCACCAGCGGAATAACGGTCATCACCACCGCCAGCCAGACAAAAAAGGCGTAACCGCCCATTGCGAAAAATTCATTCAAGGAAGCAAATGCAGGGGTCa                                                                                                                               >  1:2419235/1‑150 (MQ=255)
                     acgccacgcAGAATTGCGCGATGTTGCATCACCGAGTGCACGACCAAAACCACCAGCGGAATAACGGTCATCACCACCGCCAGCCAGACAAAAAAGGCGTAACCGCCCATTGCGAAAAATTCATTCAAGGAAGCAAATGCAGGGGTCAtt                                                                                                                             <  1:3324660/150‑1 (MQ=255)
                     acgccacgcAGAATAGCGCGATGATGCATCACCGAGTGAACGACCAAAACCACCAGCGGAATAACGGTCATCACCACCGCCAGCCAGACAAAAAAGGCGTAACCGCCCAGTGCGAAAAATTCATTCAAGGAAGCAAAAGCATGGGTCAtt                                                                                                                             >  2:1285163/1‑150 (MQ=255)
                      cgccacgcAGAATTGCGCGATGTTGCATCACCGAGTGCACGACCAAAACCACCAGCGGAATAACGGTCATCACCACCGCCAGCCAGACAAAAAAGGCGTAACCGCCCATTGCGAAAAATTCATTCAAGGAAGCAAATGCAGGGGTCAttt                                                                                                                            <  1:267741/150‑1 (MQ=255)
                             cAGAATTGCGCGATGTTGCATCACCGAGTGCACGACCAAAACCACCAGCGGAATAACGGTCATCACCACCGCCAGCCAGACAAAAAAGGCGTAACCGCCCATTGCGAAAAATTCATTCAAGGAAGCAAATGCAGGGGTCATTTACGGCct                                                                                                                     <  1:833761/150‑1 (MQ=255)
                             cAGAATTGCGCGATGTTGCATCACCGAGTGCACGACCAAAACCACCAGCGGAATAACGGTCATCACCACCGCCAGCCAGACAAAAAAGGCGTAACCGCCCATTGCGAAAAATTCATTCAAGGAAGCAAATGCAGGGGTCATTTACGGCct                                                                                                                     <  2:1647181/150‑1 (MQ=255)
                                 aTTGCGCGATGTTGCATCACCGAGTGCACGACCAAAACCACCAGCGGAATAACGGTCATCACCACCGCCAGCCAGACAAAAAAGGCGTAACCGCCCATTGCGAAAAATTCATTCAAGGAAGCAAATGCAGGGGTCATTTACGGCCTCttt                                                                                                                 >  1:1173273/1‑150 (MQ=255)
                                    gcgcGATGTTGCATCACAGAGTGCACGACCAAAACCACCAGCGGAATAACGGTCATCACCACCGCCAGCCAGACAAAAAAGGCGTAACCGCCCATTGCGAAAAATTCATTCAAGGAAGCAAATGCAGGGGTCATTTACGGCCTCTTTTCa                                                                                                              <  1:5073839/150‑1 (MQ=255)
                                      gcgATGTTGCATCACCGAGTGCACGACCAAAACCACCAGCGGAATAACGGTCATCACCACCGCCAGCCAGACAAAAAAGGCGTAACCGCCCATTGCGAAAAATTCATTCAAGGAAGCAAATGCAGGGGTCATTTACGGCCTCTTTTCAGt                                                                                                            >  1:2947450/1‑150 (MQ=255)
                                      gcgATGTTGCATCACCGAGTGCACGACCAAAACCACCAGCGGAATAACGGTCATCACCACCGCCAGCCAGACAAAAAAGGCGTAACCGCCCATTGCGAAAAATTCATTCAAGGAAGCAAATGCAGGGGTCATTTACGGCCTCTTTTCAGt                                                                                                            <  1:3150604/150‑1 (MQ=255)
                                          tGTTGCATCACCGAGTGCACGACCAAAACCACCAGCGGAATAACGGTCATCACCACCGCCAGCCAGACAAAAAAGGCGTAACCGCCCATTGCGAAAAATTCATTCAAGGAAGCAAATGCAGGGGTCATTTACGGCCTCTTTTCAGTATCa                                                                                                        <  2:92446/150‑1 (MQ=255)
                                                 tCACCGAGTGCACGACCAAAACCACCAGCGGAATAACGGTCATCACCACCGCCAGCCAGACAAAAAAGGCGTAACCGCCCATTGCGAAAAATTCATTCAAGGAAGCAAATGCAGGGGTCATTTACGGCCTCTTTTCAGTATCAGTtcact                                                                                                 <  1:1333468/150‑1 (MQ=255)
                                                 tCACCGAGTGCACGACCAAAACCACCAGCGGAATAACGGTCATCACCACCGCCAGCCAGACAAAAAAGGCGTAACCGCCCATTGCGAAAAATTCATTCAAGGAAGCAAATGCAGGGGTCATTTACGGCCTCTTTTCAGTATCAGTtcact                                                                                                 >  2:3311364/1‑150 (MQ=255)
                                                     cGAGTGCACGACCAAAACCACCAGCGGAATAACGGTCATCACCACCGCCAGCCAGACAAAAAAGGCGTAACCGCCCATTGCGAAAAATTCATTCAAGGAAGCAAATGCAGGGGTCATTTACGGCCTCTTTTCAGTATCAGTTCACTCAcc                                                                                             <  2:354828/150‑1 (MQ=255)
                                                     cGAGTGCACGACCAAAACCACCAGCGGAATAACGGTCATCACCACCGCCAGCCAGACAAAAAAGGCGTAACCGCCCATTGCGAAAAATTCATTCAAGGAAGCAAATGCAGGGGTCATTTACGGCCTCTTTTCAGTATCAGTTCACTCAcc                                                                                             <  1:727652/150‑1 (MQ=255)
                                                      gAGTGCACGACCAAAACCACCAGCGGAATAACGGTCATCACCACCGCCAGCCAGACAAAAAAGGCGTAACCGCCCATTGCGAAAAATTCATTCAAGGAAGCAAATGCAGGGGTCATTTACGGCCTCTTTTCAGTATCAGTTCACTCAccc                                                                                            >  1:2060076/1‑150 (MQ=255)
                                                       aGTGCACGACCAAAACCACCAGCGGAATAACGGTCATCACCACCGCCAGCCAGACAAAAAAGGCGTAACCGCCCATTGCGAAAAATTCATTCAAGGAAGCAAATGCAGGGGTCATTTACGGCCTCTTTTCAGTATCAGTTCACTCACCCa                                                                                           >  2:3570980/1‑150 (MQ=255)
                                                        gTGCACGACCAAAACCACCAGCGGAATAACGGTCATCACCACCGCCAGCCAGACAAAAAAGGCGTAACCGCCCATTGCGAAAAATTCATTCAAGGAAGCAAATGCAGGGGTCATTTACGGCCTCTTTTCAGTATCAGTTCACTCACCCac                                                                                          <  1:959013/150‑1 (MQ=255)
                                                           cACGACCAAAACCACCAGCGGAATAACGGTCATCACCACCGCCAGCCAGACAAAAAAGGCGTAACCGCCCATTGCGAAAAATTCATTCAAGGAAGCAAATGCAGGGGTCATTTACGGCCTCTTTTCAGTATCAGTTCACTCACCCacgga                                                                                       >  2:4569713/1‑150 (MQ=255)
                                                           cACGACCAAAACCACCAGCGGAATAACGGTCATCACCACCGCCAGCCAGACAAAAAAGGCGTAACCGCCCATTGCGAAAAATTCATTCAAGGAAGCAAATGCAGGGGTCATTTACGGCCTCTTTTCAGTATCAGTTCACTCACCCacgga                                                                                       <  2:2133369/150‑1 (MQ=255)
                                                                    aaCCACCAGCGGAATAACGGTCATCACCACCGCCAGCCAGACAAAAAAGGCGTAACCGCCCATTGCGAAAAATTCATTCAAGGAAGCAAATGCAGGGGTCATTTACGGCCTCTTTTCAGTATCAGTTCACTCACCCACGGACGGCGtttt                                                                              <  2:4715418/150‑1 (MQ=255)
                                                                    aaCCACCAGCGGAATAACGGTCATCACCACCGCCAGCCAGACAAAAAAGGCGTAACCGCCCATTGCGAAAAATTCATTCAAGGAAGCAAATGCAGGGGTCATTTACGGCCTCTTTTCAGTATCAGTTCACTCACCCACGGACGGCGtttt                                                                              >  1:392446/1‑150 (MQ=255)
                                                                         ccaGCGGAATAACGGTCATCACCACCGCCAGCCAGACAAAAAAGGCGTAACCGCCCATTGCGAAAAATTCATTCAAGGAAGCAAATGCAGGGGTCATTTACGGCCTCTTTTCAGTATCAGTTCACTCACCCACGGACGGCGTTTTTCCAt                                                                         <  1:3092648/150‑1 (MQ=255)
                                                                          caGCGGAATAACGGTCATCACCACCGCCAGCCAGACAAAAAAGGCGTAACCGCCCATTGCGAAAAATTCATTCAAGGAAGCAAATGCAGGGGTCATTTACGGCCTCTTTTCAGTATCAGTTCACTCACCCACGGACGGCGTTTTTCCATc                                                                        <  1:4204799/150‑1 (MQ=255)
                                                                          caGCGGAATAACGGTCATCACCACCGCCAGCCAGACAAAAAAGGCGTAACCGCCCATTGCGAAAAATTCATTCAAGGAAGCAAATGCAGGGGTCATTTACGGCCTCTTTTCAGTATCAGTTCACTCACCCACGGACGGCGTTTTTCCATc                                                                        <  1:1790617/150‑1 (MQ=255)
                                                                                aaTAACGGTCATCACCACCGCCAGCCAGACAAAAAAGGCGTAACCGCCCATTGCGAAAAATTCATTCAAGGAAGCAAATGCAGGGGTCATTTACGGCCTCTTTTCAGTATCAGTTCACTCACCCACGGACGGCGTTTTTCCATCAGCaaa                                                                  >  2:2102156/1‑150 (MQ=255)
                                                                                   aaCGGTCATCACCACCGCCAGCCAGACAAAAAAGGCGTAACCGCCCATTGCGAAAAATTCATTCAAGGAAGCAAATGCAGGGGTCATTTACGGCCTCTTTTCAGTATCAGTTCACTCACCCACGGACGGCGTTTTTCCATCAGCAAAATc                                                               >  2:4092294/1‑150 (MQ=255)
                                                                                    aCGGTCATCACCACCGCCAGCCAGACAAAAAAGGAGTAACCGCCCTTTGCGAAAAATTCATTCAAGGAAGCAAATGCTGGGGTCATTTACGGCCTCTTTTCAGTATCAGTTCACTCACCCACGGACGGCGTTTTTCCATCAGCAAAATCa                                                              <  2:283340/150‑1 (MQ=255)
                                                                                      ggTCATCACCACCGCCAGCCAGACAAAAAAGGCGTAACCGCCCATTGCGAAAAATTCATTCAAGGAAGCAAATGCAGGGGTCATTTACGGCCTCTTTTCAGTATCAGTTCACTCACCCACGGACGGCGTTTTTCCATCAGCAAAATCaaa                                                            <  1:1549984/150‑1 (MQ=255)
                                                                                       gTCATCACCACCGCCAGCCAGACAAAAAAGGCGTAACCGCCCATTGCGAAAAATTCATTCAAGGAAGCAAATGCAGGGGTCATTTACGGCCTCTTTTCAGTATCAGTTCACTCACCCACGGACGGCGTTTTTCCATCAGCAAAATCAAAt                                                           <  2:2233144/150‑1 (MQ=255)
                                                                                       gTCATCACCACCGCCAGCCAGACAAAAAAGGCGTAACCGCCCATTGCGAAAAATTCATTCAAGGAAGCAAATGCAGGGGTCATTTACGGCCTCTTTTCAGTATCAGTTCACTCACCCACGGACGGCGTTTTTCCATCAGCAAAATCAAAt                                                           <  1:3657540/150‑1 (MQ=255)
                                                                                         catcaCCACCGCCAGCCAGACAAAAAAGGCGTAACCGCCCATTGCGAAAAATTCATTCAAGGAAGCAAATGCAGGGGTCATTTACGGCCTCTTTTCAGTATCAGTTCACTCACCCACGGACGGCGTTTTTCCATCAGCAAAATCAAATTa                                                         <  1:3545580/150‑1 (MQ=255)
                                                                                         catcaCCACCGCCAGCCAGACAAAAAAGGCGTAACCGCCCATTGCGAAAAATTCATTCAAGGAAGCAAATGCAGGGGTCATTTACGGCCTCTTTTCAGTATCAGTTCACTCACCCACGGACGGCGTTTTTCCATCAGCAAAATCAAATTa                                                         >  1:4838418/1‑150 (MQ=255)
                                                                                           tcaCCACCGCCAGCCAGACAAAAAAGGCGTAACCGCCCATTGCGAAAAATTCATTCAAGGAAGCAAATGCAGGGGTCATTTACGGCCTCTTTTCAGTATCAGTTCACTCACCCACGGACGGCGTTTTTCCATCAGCAAAATCAAATTAcg                                                       <  2:1901374/150‑1 (MQ=255)
                                                                                            caccacCGCCAGCCAGACAAAAAAGGCGTAACCGCCCATTGCGAAAAATTCATTCAAGGAAGCAAATGCAGGGGTCATTTACGGCCTCTTTTCAGTATCAGTTCACTCACCCACGGACGGCGTTTTTCCATCAGCAAAATCAAATTAcgc                                                      <  2:4026108/150‑1 (MQ=255)
                                                                                            caccacCGCCAGCCAGACAAAAAAGGCGTAACCGCCCATTGCGAAAAATTCATTCAAGGAAGCAAATGCAGGGGTCATTTACGGCCTCTTTTCAGTATCAGTTCACTCACCCACGGACGGCGTTTTTCCATCAGCAAAATCAAATTAcgc                                                      >  2:25653/1‑150 (MQ=255)
                                                                                            caccacCGCCAGCCAGACAAAAAAGGCGTAACCGCCCATTGCGAAAAATTCATTCAAGGAAGCAAATGCAGGGGTCATTTACGGCCTCTTTTCAGTATCAGTTCACTCACCCACGGACGGCGTTTTTCCATCAGCAAAATCAAATTAcgc                                                      <  2:621142/150‑1 (MQ=255)
                                                                                            caccacCGCCAGCCAGACAAAAAAGGCGTAACCGCCCATTGCGAAAAATTCATTCAAGGAAGCAAATGCAGGGGTCATTTACGGCCTCTTTTCAGTATCAGTTCACTCACCCACGGACGGCGTTTTTCCATCAGCAAAATCAAATTAcgc                                                      <  2:5124160/150‑1 (MQ=255)
                                                                                               caccGCCAGCCAGACAAAAAAGGCGTAACCGCCCATTGCGAAAAATTCATTCAAGGAAGCAAATGCAGGGGTCATTTACGGCCTCTTTTCAGTATCAGTTCACTCACCCACGGACGGCGTTTTTCCATCAGCAAAATCAAATTAcgcatc                                                   <  2:3669174/150‑1 (MQ=255)
                                                                                               caccGCCAGCCAGACAAAAAAGGCGTAACCGCCCATTGCGAAAAATTCATTCAAGGAAGCAAATGCAGGGGTCATTTACGGCCTCTTTTCAGTATCAGTTCACTCACCCACGGACGGCGTTTTTCCATCAGCAAAATCAAATTAcgcatc                                                   <  1:3481295/150‑1 (MQ=255)
                                                                                               caccGCCAGCCAGACAAAAAAGGCGTAACCGCCCATTGCGAAAAATTCATTCAAGGAAGCAAATGCAGGGGTCATTTACGGCCTCTTTTCAGTATCAGTTCACTCACCCACGGACGGCGTTTTTCCATCAGCAAAATCAAATTAcgcatc                                                   >  2:1057806/1‑150 (MQ=255)
                                                                                                accGCCAGCCAGACAAAAAAGGCGTAACCGCCCATTGCGAAAAATTCATTCAAGGAAGCAAATGCAGGGGTCATTTACGGCCTCTTTTCAGTATCAGTTCACTCACCCACGGACGGCGTTTTTCCATCAGCAAAATCAAATTAcgcatcc                                                  <  1:1619748/150‑1 (MQ=255)
                                                                                                accGCCAGCCAGACAAAAAAGGCGTAACCGCCCATTGCGAAAAATTCATTCAAGGAAGCAAATGCAGGGGTCATTTACGGCCTCTTTTCAGTATCAGTTCACTCACCCACGGACGGCGTTTTTCCATCAGCAAAATCAAATTAcgcatcc                                                  <  1:1481138/150‑1 (MQ=255)
                                                                                                 ccGCCAGCCAGACAAAAAAGGCGTAACCGCCCATTGCGAAAAATTCATTCAAGGAAGCAAATGCAGGGGTCATTTACGGCCTCTTTTCAGTATCAGTTCACTCACCCACGGACGGCGTTTTTCCATCAGCAAAATCAAATTAcgcatccg                                                 <  2:2956608/150‑1 (MQ=255)
                                                                                                 ccGCCAGCCAGACAAAAAAGGCGTAACCGCCCATTGCGAAAAATTCATTCAAGGAAGCAAATGCAGGGGTCATTTACGGCCTCTTTTCAGTATCAGTTCACTCACCCACGGACGGCGTTTTTCCATCAGCAAAATCAAATTAcgcatccg                                                 >  1:1353369/1‑150 (MQ=255)
                                                                                                   gccagccaGACAAAAAAGGCGTAACCGCCCATTGCGAAAAATTCATTCAAGGAAGCAAATGCAGGGGTCATTTACGGCCTCTTTTCAGTATCAGTTCACTCACCCACGGACGGCGTTTTTCCATCAGCAAAATCAAATTAcgcatccgca                                               >  1:3583759/1‑150 (MQ=255)
                                                                                                      agccagACAAAAAAGGCGTAACCGCCCATTGCGAAAAATTCATTCAAGGAAGCAAATGCAGGGGTCATTTACGGCCTCTTTTCAGTATCAGTTCACTCACCCACGGACGGCGTTTTTCCATCAGCAAAATCAAATTACGCATCCGCAtca                                            >  1:1702931/1‑150 (MQ=255)
                                                                                                      agccagACAAAAAAGGCGTAACCGCCCATTGCGAAAAATTCATTCAAGGAAGCAAATGCAGGGGTCATTTACGGCCTCTTTTCAGTATCAGTTCACTCACCCACGGACGGCGTTTTTCCATCAGCAAAATCAAATTACGCATCCGCAtca                                            >  1:1485994/1‑150 (MQ=255)
                                                                                                      agccagACAAAAAAGGCGTAACCGCCCATTGCGAAAAATTCATTCAAGGAAGCAAATGCAGGGGTCATTTACGGCCTCTTTTCAGTATCAGTTCACTCACCCACGGACGGCGTTTTTCCATCAGCAAAATCAAATTACGCATCCGCAtca                                            >  2:236614/1‑150 (MQ=255)
                                                                                                      agccagACAAAAAAGGCGTAACCGCCCATTGCGAAAAATTCATTCAAGGAAGCAAATGCAGGGGTCATTTACGGCCTCTTTTCAGTATCAGTTCACTCACCCACGGACGGCGTTTTTCCATCAGCAAAATCAAATTACGCATCCGCAtca                                            <  1:409953/150‑1 (MQ=255)
                                                                                                         cagACAAAAAAGGCGTAACCGCCCATTGCGAAAAATTCATTCAAGGAAGCAAATGCAGGGGTCATTTACGGCCTCTTTTCAGTATCAGTTCACTCACCCACGGACGGCGTTTTTCCATCAGCAAAATCAAATTACGCATCCGCAtcagcg                                         <  1:4992853/150‑1 (MQ=255)
                                                                                                           gACAAAAAAGGCGTAACCGCCCATTGCGAAAAATTCATTCAAGGAAGCAAATGCAGGGGTCATTTACGGCCTCTTTTCAGTATCAGTTCACTCACCCACGGACGGCGTTTTTCCATCAGCAAAATCAAATTACGCATCCGCAtcagcgtc                                       <  1:1652318/150‑1 (MQ=255)
                                                                                                           gACAAAAAAGGCGTAACCGCCCATTGCGAAAAATTCATTCAAGGAAGCAAATGCAGGGGTCATTTACGGCCTCTTTTCAGTATCAGTTCACTCACCCACGGACGGCGTTTTTCCATCAGCAAAATCAAATTACGCATCCGCAtcagcgtc                                       <  2:1923887/150‑1 (MQ=255)
                                                                                                            actaaaaaGTCGTAACCGCAGATTGCGAAAAATTCATTAAAGGAAGCAAATGCAGGGGTCATTTACGGCCTCTTTTCAGTATCAGTTCACTCACCCAAGGACGGCGTTTTTCAATCAGGAAAATCAAATTACGCATCAGCTtcagcgtct                                      <  2:2816281/147‑2 (MQ=255)
                                                                                                            aCAAAAAAGGCGTAACCGCCCATTGCGAAAAATTCATTCAAGGAAGCAAATGCAGGGGTCATTTACGGCCTCTTTTCAGTATCAGTTCACTCACCCACGGACGGCGTTTTTCCATCAGCAAAATCAAATTACGCATCCGCAtcagcgtca                                      <  1:5155753/150‑1 (MQ=255)
                                                                                                            aCAAAAAAGGCGTAACCGCCCATTGCGAAAAATTCATTCAAGGAAGCAAATGCAGGGGTCATTTACGGCCTCTTTTCAGTATCAGTTCACTCACCCACGGACGGCGTTTTTCCATCAGCAAAATCAAAATACGCATCCGCAtcagcgtca                                      >  2:2618833/1‑150 (MQ=255)
                                                                                                             cAAAAAAGGCGTAACCGCCCATTGCGAAAAATTCATTCAAGGAAGCAAATGCAGGGGTCATTTACGGCCTCTTTTCAGTATCAGTTCACTCACCCACGGACGGCGTTTTTCCATCAGCAAAATCAAATTACGCATCCGCAtcagcgtcag                                     <  1:3275080/150‑1 (MQ=255)
                                                                                                                aaaaGGCGCAACCGCCCATTGCGAAAAATTCATTCAAGGAAGCAAATGCAGGGGTCATTTACGGCCTCTTTTCAGTATCAGTTCACTCACCCACGGACGGCGTTTTTCCATCAGCAAAATCAAATTACGCATCCGCATcagcgtcagcgt                                  <  1:4295744/150‑1 (MQ=255)
                                                                                                                 aaaGGCGTAACCGCCCATTGCGAAAAATTCATTCCAGGAAGCAAATGCAGGGGTCATTTACGGCCTCTTTTCAGTATCAGTTCACTCACCCACGGACGGCGTTTTTCCATCAGCAAAATCAAATTACGCATCCGCATCAGCGTCAGCGTg                                 <  2:3109557/150‑1 (MQ=255)
                                                                                                                 aaaGGCGTAACCGCCCATTGCGAAAAATTCATTCAAGGAAGCAAATGCAGGGGTCATTTACGGCCTCTTTTCAGTATCAGTTCACTCACCCACGGACGGCGTTTTTCCATCAGCAAAATCAAATTACGCATCCGCATCAGCGTCAGCGTg                                 <  1:1769349/150‑1 (MQ=255)
                                                                                                                 aaaGGCGTAACCGCCCATTGCGAAAAATTCATTCAAGGAAGCAAATGCAGGGGTCATTTACGGCCTCTTTTCAGTATCAGTTCACTCACCCACGGACGGCGTTTTTCCATCAGCAAAATCAAATTACGCATCCGCATCAGCGTCAGCGTg                                 >  2:4354678/1‑150 (MQ=255)
                                                                                                                   aGGCGTAACCGCCCATTGCGAAAAATTCATTCAAGGAAGCAAATGCAGGGGTCATTTACGGCCTCTTTTCAGTATCAGTTCACTCACCCACGGACGGCGTTTTTCCATCAGCAAAATCAAATTACGCATCCGCATCAGCGTCAGCGTGGc                               >  2:2288156/1‑150 (MQ=255)
                                                                                                                        tAACCGCCCATTGCGAAAAATTCATTCAAGGAAGCAAATGCAGGGGTCATTTACGGCCTCTTTTCAGTATCAGTTCACTCACCCACGGACGGCGTTTTTCCATCAGCAAAATCAAATTACGCATCCGCATCAGCGTCAGCGTGGCAGACa                          <  2:3003809/150‑1 (MQ=255)
                                                                                                                        tAACCGCCCATTGCGAAAAATTCATTCAAGGAAGCAAATGCAGGGGTCATTTACGGCCTCTTTTCAGTATCAGTTCACTCACCCACGGACGGCGTTTTTCCATCAGCAAAATCAAATTACGCATCCGCATCAGCGTCAGCGTGGCAGACa                          >  2:3371412/1‑150 (MQ=255)
                                                                                                                        tAACCGCCCATTGCGAAAAATTCATTCAAGGAAGCAAATGCAGGGGTCATTTACGGCCTCTTTTCAGTATCAGTTCACTCACCCACGGACGGCGTTTTTCCATCAGCAAAATCAAATTACGCATCCGCATCAGCGTCAGCGTGGCAGACa                          <  1:982109/150‑1 (MQ=255)
                                                                                                                         aaCCGCCCATTGCGAAAAATTCATTCAAGGAAGCAAATGCAGGGGTCATTTACGGCCTCTTTTCAGTATCAGTTCACTCACCCACGGACGGCGTTTTTCCATCAGCAAAATCAAATTACGCATCCGCATCAGCGTCAGCGTGGCAGACAg                         >  1:4635879/1‑150 (MQ=255)
                                                                                                                          aCCGCCCATTGCGAAAAATTCATTCAAGGAAGCAAATGCAGGGGTCATTTACGGCCTCTTTTCAGTATCAGTTCACTCACCCACGGACGGCGTTTTTCCATCAGCAAAATCAAATTACGCATCCGCATCAGCGTCAGCGTGGCAGACAgg                        >  1:1060506/1‑150 (MQ=255)
                                                                                                                           ccGCCCATTGCGAAAAATTCATTCAAGGAAGCAAATGCAGGGGTCATTTACGGCCTCTTTTCAGTATCAGTTCACTCACCCACGGACGGCGTTTTTCCATCAGCAAAATCAAATTACGCATCCGCATCAGCGTCAGCGTGGCAGACAGGa                       <  1:3764011/150‑1 (MQ=255)
                                                                                                                           ccGCCCATTGCGAAAAATTCATTCAAGGAAGCAAATGCAGGGGTCATTTACGGCCTCTTTTCAGTATCAGTTCACTCACCCACGGACGGCGTTTTTCCATCAGCAAAATCAAATTACGCATCCGCATCAGCGTCAGCGTGGCAGACAGGa                       <  1:2818579/150‑1 (MQ=255)
                                                                                                                            cGCCCATTGCGAAAAATTCATTCAAGGAAGCAAATGCAGGGGTCATTTACGGCCTCTTTTCAGTATCAGTTCACTCACCCACGGACGGCGTTTTTCCATCAGCAAAATCAAATTACGCATCCGCATCAGCGTCAGCGTGGCAGACAGGAg                      >  2:3643423/1‑150 (MQ=255)
                                                                                                                             gCCCATTGCGAAAAATTCATTCAAGGAAGCAAATGCAGGGGTCATTTACGGCCTCTTTTCAGTATCAGTTCACTCACCCACGGACGGCGTTTTTCCATCAGCAAAATCAAATTACGCATCCGCATCAGCGTCAGCGTGGCAGACAGGAGc                     >  2:171126/1‑150 (MQ=255)
                                                                                                                              cccATTGCGAAAAATTCATTCAAGGAAGCAAATGCAGGGGTCATTTACGGCCTCTTTTCAGTATCAGTTCACTCACCCACGGACGGCGTTTTTCCATCAGCAAAATCAAATTACGCATCCGCATCAGCGTCAGCGTGGCAGACAGGAGCa                    <  1:4527079/150‑1 (MQ=255)
                                                                                                                               ccATTGCGAAAAATTCATTCAAGGAAGCAAATGCAGGGGTCATTTACGGCCTCTTTTCAGTATCAGTTCACTCACCCACGGACGGCGTTTTTCCATCAGCAAAATCAAATTACGCATCCGCATCAGCGTCAGCGTGGCAGACAGGAGCAg                   >  1:4383534/1‑150 (MQ=255)
                                                                                                                                cATTGCGAAAAATTCATTCAAGGAAGCAAATGCAGGGGTCATTTACGGCCTCTTTTCAGTATCAGTTCACTCACCCACGGACGGCGTTTTTCCATCAGCAAAATCAAATTACGCATCCGCATCAGCGTCAGCGTGGCAGACAGGAGCAgg                  >  1:1935807/1‑150 (MQ=255)
                                                                                                                                   tGCGAAAAATTCATTCAAGGAAGCAAATGCAGGGGTCATTTACGGCCTCTTTTCAGTATCAGTTCACTCACCCACGGACGGCGTTTTTCCATCAGCAAAATCAAATTACGCATCCGCATCAGCGTCAGCGTGGCAGACAGGAGCAGGAAg               <  2:1173273/150‑1 (MQ=255)
                                                                                                                                    gCGAAAAATTCATTCAAGGAAGCAAATGCAGGGGTCATTTACGGCCTCTTTTCAGTATCAGTTCACTCACCCACGGACGGCGTTTTTCCATCAGCAAAATCAAATTACGCATCCGCATCAGCGTCAGCGTGGCAGACAGGAGCAGGAAGc              >  1:1669937/1‑150 (MQ=255)
                                                                                                                                     cGAAAAATTCATTCAAGGAAGCAAATGCAGGGGTCATTTACGGCCTCTTTTCAGTATCAGTTCACTCACCCACGGACGGCGTTTTTCCATCAGCAAAATCAAATTACGCATCCGCATCAGCGTCAGCGTGGCAGACAGGAGCAGGAAGcc             >  1:3459406/1‑150 (MQ=255)
                                                                                                                                      gAAAAATTCATTCAAGGAAGCAAATGCAGGGGTCATTTACGGCCTCTTTTCAGTATCAGTTCACTCACCCACGGACGGCGTTTTTCCATCAGCAAAATCAAATTACGCATCCGCATCAGCGTCAGCGTGGCAGACAGGAGCAGGAAGCCa            >  2:4967585/1‑150 (MQ=255)
                                                                                                                                        aaaaTTCATTCAAGGAAGCAAATGCAGGGGTCATTTACGGCCTCTTTTCAGTATCAGTTCACTCACCCACGGACGGCGTTTTTCCATCAGCAAAATCAAATTACGCATCCGCATCAGCGTCAGCGTGGCAGACAGGAGCAGGAAGCCaaa          >  1:2027330/1‑150 (MQ=255)
                                                                                                                                         agattcattcAATGAATCAAATGCAGGTGTCATTTTCGGCCTATTTTCAGTATCAGTTCACTCACCCACGGACGTCGTTTTTCCATCAGCAAAATCAAATTACGCATCCGCATCAGCGTCAGCGTGGCAGACAGGAGCAGGAAGCCaaaa         <  2:4825935/148‑1 (MQ=255)
                                                                                                                                            ttcattcaAGGAAGCAAATGCAGGGGTCATTTACGGCCTCTTTTCAGTATCAGTTCACTCACCCACGGACGGCGTTTTTCCATCAGCAAAATCAAATTACGCATCCGCATCAGCGTCAGCGTGGCAGACAGGAGCAGGAAGCCAAAAATc      <  1:2618359/150‑1 (MQ=255)
                                                                                                                                                ttcAAGGAAGCAAATGCAGGGGTCATTTACGGCCTCTTTTCAGTATCAGTTCACTCACCCACGGACGGCGTTTTTCCATCAGCAAAATCAAATTACGCATCCGCATCAGCGTCAGCGTGGCAGACAGGAGCAGGAAGCCAAAAATCGAcc  >  1:4874926/1‑150 (MQ=255)
                                                                                                                                                ttcAAGGAAGCAAATGCAGGGGTCATTTACGGCCTCTTTTCAGTATCAGTTCACTCACCCACGGACGGCGTTTTTCCATCAGCAAAATCAAATTACGCATCCGCATCAGCGTCAGCGTGGCAGACAGGAGCAGGAAGCCAAAAATCGAcc  <  1:1363812/150‑1 (MQ=255)
                                                                                                                                                ttcAAGGAAGCAAATGCAGGGGTCATTTACGGCCTCTTTTCAGTATCAGTTCACTCACCCACGGACGGCGTTTTTCCATCAGCAAAATCAAATTACGCATCCGCATCAGCGTCAGCGTGGCAGACAGGAGCAGGAAGCCAAAAATCGAcc  >  1:4050313/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
TCACGCGCCCGCTGTTGCGCCACGCCACGCAGAATTGCGCGATGTTGCATCACCGAGTGCACGACCAAAACCACCAGCGGAATAACGGTCATCACCACCGCCAGCCAGACAAAAAAGGCGTAACCGCCCATTGCGAAAAATTCATTCCAGGAAGCAAATGCAGGGGTCATTTACGGCCTCTTTTCAGTATCAGTTCACTCACCCACGGACGGCGTTTTTCCATCAGCAAAATCAAATTACGCATCCGCATCAGCGTCAGCGTGGCAGACAGGAGCAGGAAGCCAAAAATCGACC  >  CP009273/2288896‑2289189

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: