Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,289,043 C→A 41.4% 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,0430CA41.4% 378.0 / inf 362W8L (TGG→TTG) ccmDcytochrome c biogenesis protein
Reads supporting (aligned to +/- strand):  ref base C (94/118);  new base A (73/77);  total (167/195)
Fisher's exact test for biased strand distribution p-value = 4.54e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 368 total aligned reads displayed.

CCTCACGCGCCCGCTGTTGCGCCACGCCACGCAGAATTGCGCGATGTTGCATCACCGAGTGCACGACCAAAACCACCAGCGGAATAACGGTCATCACCACCGCCAGCCAGACAAAAAAGGCGTAACCGCCCATTGCGAAAAATTCATTCCAGGAAGCAAATGCAGGGGTCATTTACGGCCTCTTTTCAGTATCAGTTCACTCACCCACGGACGGCGTTTTTCCATCAGCAAAATCAAATTACGCATCCGCATCAGCGTCAGCGTGGCAGACAGGAGCAGGAAGCCAAAAATCGACCAGC  >  CP009273/2288894‑2289192
                                                                                                                                                     |                                                                                                                                                     
   cACGCGCCCGCTGTTGCGCCACGCCACGCAGAATTGCGCGATGTTGCATCACCGAGTGCACGACCAAAACCACCAGCGGAATAACGGTCATCACCACCGCCAGCCAGACAAAAAAGGCGTAACCGCCCATTGCGAAAAATTCATTCCAgg                                                                                                                                                    <  2:473367/150‑1 (MQ=255)
   cACGCGCCCGCTGTTGCGCCACGCCACGCAGAATTGCGCGATGTTGCATCACCGAGTGCACGACCAAAACCACCAGCGGAATAACGGTCATCACCACCGCCAGCCAGACAAAAAAGGCGTAACCGCCCATTGCGAAAAATTCATTCCAgg                                                                                                                                                    <  2:1400169/150‑1 (MQ=255)
       cgcCCGCTGTTGCGCCACGCCACGCAGAATTGCGCGATGTTGCATCACCGAGTGCACGACCAAAACCACCAGCGGAATAACGGTCATCACCACCGCCAGCCAGACAAAAAAGGCGTAACCGCCCATTGCGAAAAATTCATTCAAGGAAGc                                                                                                                                                <  2:353856/150‑1 (MQ=255)
       cgcCCGCTGTTGCGCCACGCCACGCAGAATTGCGCGATGTTGCATCACCGAGTGCACGACCAAAACCACCAGCGGAATAACGGTCATCACCACCGCCAGCCAGACAAAAAAGGCGTAACCGCCCATTGCGAAAAATTCATTCAAGGAAGc                                                                                                                                                <  1:1538218/150‑1 (MQ=255)
        gcCCGCTGTTGCGCCACGCCACGCAGAATTGCGCGATGTTGCATCACCGAGTGCACGACCAAAACCACCAGCGGAATAACGGTCATCACCACCGCCAGCCAGACAAAAAAGGCGTAACCGCCCATTGCGAAAAATTCATTCAAGGAAGCa                                                                                                                                               >  2:2595635/1‑150 (MQ=255)
          ccGCTGTTGCGCCACGCCACGCAGAATTGCGCGATGTTGCATCACCGAGTGCACGACCAAAACCACCAGCGGAATAACGGTCATCACCACCGCCAGCCAGACAAAAAAGGCGTAACCGCCCATTGCGAAAAATTCATTCCAGGAAGCaaa                                                                                                                                             <  2:3406498/150‑1 (MQ=255)
           cGCTGTTGCGCCACGCCACGCAGAATTGCGCGATGTTGCATCACCGAGTGCACGACCAAAACCACCAGCGGAATAACGGTCATCACCACCGCCAGCCAGACAAAAAAGGCGTAACCGCCCATTGCGAAAAATTCATTCAAGGAAGCAAAt                                                                                                                                            <  2:3978747/150‑1 (MQ=255)
               gTTGCGCCACGCCACGCAGAATTGCGCGATGTTGCATCACCGAGTGCACGACCAAAACCACCAGCGGAATAACGGTCATCACCACCGCCAGCCAGACAAAAAAGGCGTAACCGCCCATTGCGAAAAATTCATTCCAGGAAGCAAATGCAg                                                                                                                                        <  1:3587113/150‑1 (MQ=255)
                ttGCGCCACGCCACGCAGAATTGCGCGATGTTGCATCACCGAGTGCACGACCAAAACCACCAGCGGAATAACGGTCATCACCACCGCCAGCCAGACAAAAAAGGCGTAACCGCCCATTGCGAAAAATTCATTCCAGGAAGCAAATGCAgg                                                                                                                                       >  1:1074002/1‑150 (MQ=255)
                ttGCGCCACGCCACGCAGAATTGCGCGATGTTGCATCACCGAGTGCACGACCAAAACCACCAGCGGAATAACGGTCATCACCACCGCCAGCCAGACAAAAAAGGCGTAACCGCCCATTGCGAAAAATTCATTCCAGGAAGCAAATGCAgg                                                                                                                                       >  1:4683959/1‑150 (MQ=255)
                ttGCGCCACGCCACGCAGAATTGCGCGATGTTGCATCACCGAGTGCACGACCAAAACCACCAGCGGAATAACGGTCATCACCACCGCCAGCCAGACAAAAAAGGCGTAACCGCCCATTGCGAAAAATTCATTCAAGGAAGCAAATGCAgg                                                                                                                                       <  2:3788625/150‑1 (MQ=255)
                   cgccacgccaCGCAGAATTGCGCGATGTTGCATCACCGAGTGCACGACCAAAACCACCAGCGGAATAACGGTCATCACCACCGCCAGCCAGACAAAAAAGGCGTAACCGCCCATTGCGAAAAATTCATTCAAGGAAGCAAATGCAGGGGt                                                                                                                                    >  1:1592337/1‑150 (MQ=255)
                    gccacgccacGCAGAATTGCGCGATGTTGCATCACCGAGTGCACGACCAAAACCACCAGCGGAATAACGGTCATCACCACCGCCAGCCAGACAAAAAAGGCGTAACCGCCCATTGCGAAAAATTCATTCCAGGAAGCAAATGCAGGGGTc                                                                                                                                   >  2:843119/1‑150 (MQ=255)
                       acgccacgcAGTATTGCGCGATGTTGCATCACCGAGTGCACGACCAAAACCACCAGCGGAATAACGGTCATCACCACCGCCAGCCAGACAAAAAAGGCGTAACCGCCCATTGCGAAAAATTCATTCAAGGAAGCAAATGCAGGGGTCAtt                                                                                                                                <  2:5020396/150‑1 (MQ=255)
                        cgccacgcAGAATTGCGCGATGTTGCATCACCGAGTGCACGACCAAAACCACCAGCGGAATAACGGTCATCACCACCGCCAGCCAGACAAAAAAGGCGTAACCGCCCATTGCGAAAAATTCATTCCAGGAAGCAAATGCAGGGGTCAttt                                                                                                                               <  1:3741832/150‑1 (MQ=255)
                        cgccacgcAGAATTGCGCGATGTTGCATCACCGAGTGCACGACCAAAACCACCAGCGGAATAACGGTCATCACCACCGCCAGCCAGACAAAAAAGGCGTAACCGCCCATTGCGAAAAATTCATTCCAGGAAGCAAATGCAGGGGTCAttt                                                                                                                               <  1:546130/150‑1 (MQ=255)
                                aGAATTGCGCGATGTTGCATCACCGAGTGCACGACCAAAACCACCAGCGGAATAACGGTCATCACCACCGCCAGCCAGACAAAAAAGGCGTAACCGCCCATTGCGAAAAATTCATTCCAGGAAGCACATGCAGGGGTCATTTACGGCctc                                                                                                                       >  2:3192062/1‑150 (MQ=255)
                                  aaTTGCGCGATGTTGCATCACCGAGTGCACGACCAAAACCACCAGCGGAATAACGGTCATCACCACCGCCAGCCAGACAAAAAAGGCGTAACCGCCCATTGCGAAAAATTCATTCCAGGAAGCAAATGCAGGGGTCATTTACGGCCTCtt                                                                                                                     <  1:123458/150‑1 (MQ=255)
                                      gcgcGATGTTGCATCACCGAGTGCACGACCAAAACCACCAGCGGAATAACGGTCATCACCACCGCCAGCCAGACAAAAAAGGCGTAACCGCCCATTGCGAAAAATTCATTCAAGGAAGCAAATGCAGGGGTCATTTACGGCCTCTTTTCa                                                                                                                 <  1:2989471/150‑1 (MQ=255)
                                        gcgATGTTGCATCACCGAGTGCACGACCAAAACCACCAGCGGAATAACGGTCATCACCACCGCCAGCCAGACAAAAAAGGCGTAACCGCCCATTGCGAAAAATTCATTCAAGGAAGCAAATGCAGGGGTCATTTACGGCCTCTTTTCAGt                                                                                                               <  1:2792073/150‑1 (MQ=255)
                                          gATGTTGCATCACCGAGTGCACGACCAAAACCACCAGCGGAATAACGGTCATCACCACCGCCAGCCAGACAAAAAAGGCGTAACCGCCCATTGCGAAAAATTCATTCCAGGAAGCAAATGCAGGGGTCATTTACGGCCTCTTTTCAGTAt                                                                                                             >  2:1633128/1‑150 (MQ=255)
                                           aTGTTGCATCACCGAGTGCACGACCAAAACCACCAGCGGAATAACGGTCATCACCACCGCCAGCCAGACAAAAAAGGCGTAACCGCCCATTGCGAAAAATTCATTCAAGGAAGCAAATGCAGGGGTCATTTACGGCCTCTTTTCAGTATc                                                                                                            <  1:381844/150‑1 (MQ=255)
                                              ttGCATCACCGAGTGCACGACCAAAACCACCAGCGGAATAACGGTCATCACCACCGCCAGCCAGACAAAAAAGGCGTAACCGCCCATTGCGAAAAATTCATTCCAGGAAGCAAATGCAGGGGTCATTTACGGCCTCTTTTCAGTATCAGt                                                                                                         >  1:2447589/1‑150 (MQ=255)
                                               tGCATCACCGAGTGCACGACCAAAACCACCAGCGGAATAACGGTCATCACCACCGCCAGCCAGACAAAAAAGGCGTAACCGCCCATTGCGAAAAATTCATTCAAGGAAGCAAATGCAGGGGTCATTTACGGCCTCTTTTCAGTATCAGtt                                                                                                        <  1:2134334/150‑1 (MQ=255)
                                                  aTCACCGAGTGCACGACCAAAACCACCAGCGGAATAACGGTCATCACCACCGCCAGCCAGACAAAAAAGGCGTAACCGCCCATTGCGAAAAATTCATTCCAGGAAGCAAATGCAGGGGTCATTTACGGCCTCTTTTCAGTATCAGTtcac                                                                                                     >  2:3702225/1‑150 (MQ=255)
                                                  aTCACCGAGTGCACGACCAAAACCACCAGCGGAATAACGGTCATCACCACCGCCAGCCAGACAAAAAAGGCGTAACCGCCCATTGCGAAAAATTCATTCCAGGAAGCAAATGCAGGGGTCATTTACGGCCTCTTTTCAGTATCAGTtcac                                                                                                     >  1:4852881/1‑150 (MQ=255)
                                                  aTCACCGAGTGCACGACCAAAACCACCAGCGGAATAACGGTCATCACCACCGCCAGCCAGACAAAAAAGGCGTAACCGCCCATTGCGAAAAATTCATTCCAGGAAGCAAATGCAGGGGTCATTTACGGCCTCTTTTCAGTATCAGTtaac                                                                                                     <  1:2000485/150‑1 (MQ=255)
                                                  aTCACCGAGTGCACGACCAAAACCACCAGCGGAATAACGGTCATCACCACCGCCAGCCAGACAAAAAAGGCGTAACCGCCCATTGCGAAAAATTCATTCAAGGAAGCACATGCAGGGGTCATTTACGGCCTCTTTTCAGTATCAGTtcac                                                                                                     >  2:4581883/1‑150 (MQ=255)
                                                    cACCGAGTGCACGACCAAAACCACCAGCGGAATAACGGTCATCACCACCGCCAGCCAGACAAAAAAGGCGTAACCGCCCATTGCGAAAAATTCATTCAAGGAAGCAAATGCAGGGGTCATTTACGGCCTCTTTTCAGTATCAGTtcactc                                                                                                   >  2:1746395/1‑150 (MQ=255)
                                                     aCCGAGTGCACGACCAAAACCACCAGCGGAATAACGGTCATCACCACCGCCAGCCAGACAAAAAAGGCGTAACCGCCCATTGCGAAAAATTCATTCAAGGAAGCAAATGCAGGGGTCATTTACGGCCTCTTTTCAGTATCAGTtcactca                                                                                                  <  2:1858390/150‑1 (MQ=255)
                                                       cGAGTGCACGACCAAAACCACCAGCGGAATAACGGTCATCACCACCGCCAGCCAGACAAAAAAGGCGTAACCGCCCATTGCGAAAAATTCATTCCAGGAAGCAAATGCAGGGGTCATTTACGGCCTCTTTTCAGTATCAGTTCACTCAcc                                                                                                >  1:2171697/1‑150 (MQ=255)
                                                        gAGTGCACGACCAAAACCACCAGCGGAATAACGGTCATCACCACCGCCAGCCAGACAAAAAAGGCGTAACCGCCCATTGCGAAAAATTCATTCCAGGAAGCAAATGCAGGGGTCATTTACGGCCTCTTTTCAGTATCAGTTCACTCAccc                                                                                               >  1:761104/1‑150 (MQ=255)
                                                        gAGTGCACGACCAAAACCACCAGCGGAATAACGGTCATCACCACCGCCAGCCAGACAAAAAAGGCGTAACCGCCCATTGCGAAAAATTCATTCCAGGAAGCAAATGCAGGGGTCATTTACGGCCTCTTTTCAGTATCAGTTAACTCAccc                                                                                               <  1:4220651/150‑1 (MQ=255)
                                                         aGTGCACGACCAAAACCACCAGCGGAATAACGGTCATCACCACCGCCAGCCAGACAAAAAAGGCGTAACCGCCCATTGCGAAAAATTCATTCAAGGAAGCAAATGCAGGGGTCATTTACGGCCTCTTTTCAGTATCAGTTCACTCACCCa                                                                                              <  1:2250747/150‑1 (MQ=255)
                                                           tGCACGACCAAAACCACCAGCGGAATAACGGTCATCACCACCGCCAGCCAGACAAAAAAGGCGTAACCGCCCATTGCGAAAAATTCATTCCAGGAAGCAAATGCAGGGGTCATTTACGGCCTCTTTTCAGTATCAGTTCACTCACCCacg                                                                                            >  1:4307346/1‑150 (MQ=255)
                                                              aCGACCAAAACCACCAGAGGAATAACGGTCATCACCACCGCCAGCCAGACAAAAAAGGCGTAACCGCCCATTGCGAAAAATTCATTCAAGGAAGCAAATGCAGGGGTCATTTACGGCCTCTTTTCAGTATCAGTTCACTCACCCacggac                                                                                         <  2:1706798/150‑1 (MQ=255)
                                                               cGACCAACACCACCAGCGGAATAACGGTCATCACCACCGCCAGCCAGACAAAAAAGGCGTAACCGCCCATTGCGAAAAATTCATTCCAGGAAGCAAATGCAGGGGTCATTTACGGCCTCTTTTCAGTATCAGTTAACTCACCCacggacg                                                                                        <  2:39044/150‑1 (MQ=255)
                                                                gACCAAAACCACCAGCGGAATAACGGTCATCACCACCGCCAGCCAGACAAAAAAGGCGTAACCGCCCATTGCGAAAAATTCATTCAAGGAAGCAAATGCAGGGGTCATTTACGGCCTCTTTTCAGTATCAGTTCACTCACCCacggacgg                                                                                       >  2:2501070/1‑150 (MQ=255)
                                                                 aCCAAAACCACCAGCGGAATAACGGTCATCACCACCGCCAGCCAGACAAAAAAGGCGTAACCGCCCATTGCGAAAAATTCATTCAAGGAAGCAAATGCAGGGGTCATTTACGGCCTCTTTTCAGTATCAGTTCACTCACCCACGGACGGc                                                                                      <  1:3616102/150‑1 (MQ=255)
                                                                    aaaaCCACCAGCGGAATAACGGTCATCACCACCGCCAGCCAGACAAAAAAGGCGTAACCGCCCATTGCGAAAAATTCATTCCAGGAAGCAAATGCAGGGGTCATTTACGGCCTCTTTTCAGTATCAGTTCACTCACCCACGGACGGCGtt                                                                                   <  1:650569/150‑1 (MQ=255)
                                                                      aaCCACCAGCGGAATAACGGTCATCACCACCGCCAGCCAGACAAAAAAGGCGTAACCGCCCATTGCGAAAAATTCATTCAAGGAAGCAAATGCAGGGGTCATTTACGGCCTCTTTTCAGTATCAGTTCACTCACCCACGGACGGCGtttt                                                                                 >  1:1581017/1‑150 (MQ=255)
                                                                         caccaGCGGAATAACGGTCATCACCACCGCCAGCCAGACAAAAAAGGCGTAACCGCCCATTGCGAAAAATTCATTCCAGGAAGCAAATGCAGGGGTCATTTACGGCCTCTTTTCAGTATCAGTTCACTCACCCACGGACGGCGTTTTTcc                                                                              <  2:1609035/150‑1 (MQ=255)
                                                                            caGCGGAATAACGGTCATCACCACCGCCAGCCAGACAACAAAGGCGTAACCGCCCATTGCGAAAAATTCATTCAAGGAAGCAAATGCAGGGGTCATTTACGGCCTCTTTTCAGTATCAGTTCACTCACCCACGGACGGCGTTTTTCCATc                                                                           <  2:355100/150‑1 (MQ=255)
                                                                              gCGGAATAACGGTCATCACCACCGCCAGCCAGACAAAAAAGGCGTAACCGCCCATTGCGAAAAATTCATTCCAGGAAGCAAATGCAGGGGTCATTTACGGCCTCTTTTCAGTATCAGTTCACTCACCCACGGACGGCGTTTTTCCATCAg                                                                         <  2:1217333/150‑1 (MQ=255)
                                                                               cGGAATAACGGTCATCACCACCGCCAGCCAGACAAAAAAGGCGTAACCGCCCATTGCGAAAAATTCATTCAAGGAAGCAAATGCAGGGGTCATTTACGGCCTCTTTTCAGTATCAGTTCACTCACCCACGGACGGCGTTTTTCCATCAGc                                                                        <  1:1684208/150‑1 (MQ=255)
                                                                                ggAATAACGGTCATCACCACCGCCAGCCAGACAAAAAAGGCGTAACCGCCCATTGCGAAAAATTCATTCCAGGAAGCAAATGCAGGGGTCATTTACGGCCTCTTTTCAGTATCAGTTCACTCACCCACGGACGGCGTTTTTCCATCAGCa                                                                       >  1:287490/1‑150 (MQ=255)
                                                                                 gAATAACGGTCATCACCACCGCCAGCCAGACAAAAAAGGCGTAACCGCCCATTGCGAAAAATTCATTCCAGGAAGCAAATGCAGGGGTCATTTACGGCCTCTTTTCAGTATCAGTTCACTCACCCACGGACGGCGTTTTTCCATCAGCaa                                                                      <  2:2239089/150‑1 (MQ=255)
                                                                                  aaTAACGGTCATCACCACCGCCAGCCAGACAAAAAAGGCGTAACCGCCCATTGCGAAAAATTCATTCCAGGAAGCAAATGCAGGGGTCATTTACGGCCTCTTTTCAGTATCAGTTAACTCACCCACGGACGGCGTTTTTCCATCAGCaaa                                                                     >  1:4930843/1‑150 (MQ=255)
                                                                                     aaCGGTCATCACCACCGCCAGCCAGACAAAAAAGGCGTAACCGCCCATTGCGAAAAATTCATTCAAGGAAGCAAATGCAGGGGTCATTTACGGCCTCTTTTCAGTATCAGTTCACTCACCCACGGACGGCGTTTTTCCATCAGCAAAATc                                                                  <  2:2710218/150‑1 (MQ=255)
                                                                                     aaCGGTCATCACCACCGCCAGCCAGACAAAAAAGGCGTAACCGCCCATTGCGAAAAATTCATTCAAGGAAGCAAATGCAGGGGTCATTTACGGCCTCTTTTCAGTATCAGTTCACTCACCCACGGACGGCGTTTTTCCATCAGCAAAATc                                                                  >  2:4538452/1‑150 (MQ=255)
                                                                                      aCGGTCATCACCACCGCCAGCCAGACAAAAAAGGCGTAACCGCCCATTGCGAAAAATTCATTCAAGGAAGCAAATGCAGGGGTCATTTACGGCCTCTTTTCAGTATCAGTTCACTCACCCACGGACGGCGTTTTTCCATCAGCAAAATCa                                                                 <  2:2696885/150‑1 (MQ=255)
                                                                                       cGGTCATCACCACCGCCAGCCAGACAAAAAAGGCGTAACCGCCCATTGCGAAAAATTCATTCCAGGAAGCAAATGCAGGGGTCATTTACGGCCTCTTTTCAGTATCAGTTCACTCACCCACGGACGGCGTTTTTCCATCAGCAAAATCaa                                                                <  1:1180005/150‑1 (MQ=255)
                                                                                             tcaCCACCGCCAGCCAGACAAAAAAGGCGTAACCGCCCATTGCGAAAAATTCATTCCAGGAAGCAAATGCAGGGGTCATTTACGGCCTCTTTTCAGTATCAGTTCACTCACCCACGGACGGCGTTTTTCCATCAGCAAAATCAAATTAcg                                                          <  2:1045866/150‑1 (MQ=255)
                                                                                             tcaCCACCGCCAGCCAGACAAAAAAGGCGTAACCGCCCATTGCGAAAAATTCATTCCAGGAAGCAAATGCAGGGGTCATTTACGGCCTCTTTTCAGTATCAGTTCACTCACCCACGGACGGCGTTTTTCCATCAGCAAAATCAAATTAcg                                                          <  1:1848089/150‑1 (MQ=255)
                                                                                              caccacCGCCAGCCAGACAAAAAAGGCGTAACCGCCCATTGCGAAAAATTCATTCCAGGAAGCAAATGCAGGGGTCATTTACGGCCTCTTTTCAGTATCAGTTCACTCACCCACGGACGGCGTTTTTCCATCAGCAAAATCAAATTAcgc                                                         <  2:2616327/150‑1 (MQ=255)
                                                                                               accaccGCCAGCCAGACAAAAAAGGCGTAACCGCCCATTGCGAAAAATTCATTCCAGGAAGCAAATGCAGGGGTCATTTACGGCCTCTTTTCAGTATCAGTTCACTCACCCACGGACGGCGTTTTTCCATCAGCAAAATCAAATTAcgca                                                        >  2:2270583/1‑150 (MQ=255)
                                                                                                ccaccGCCAGCCAGACAAAAAAGGCGTAACCGCCCATTGCGAAAAATTCATTCCAGGAAGCAAATGCAGGGGTCATTTACGGCCTCTTTTCAGTATCAGTTCACTCACCCACGGACGGCGTTTTTCCATCAGCAAAATCAAATTAcgcat                                                       >  2:1148719/1‑150 (MQ=255)
                                                                                                ccaccGCCAGCCAGACAAAAAAGGCGTAACCGCCCATTGCGAAAAATTCATTCCAGGAAGCAAATGCAGGGGTCATTTACGGCCTCTTTTCAGTATCAGTTAACTCACCCACGGACGGCGTTTTTCCATCAGCAAAATCAAATTAcgcat                                                       <  1:3956613/150‑1 (MQ=255)
                                                                                                 caccGCCAGCCAGACAAAAAAGGCGTAACCGCCCATTGCGAAAAATTCATTCCAGGAAGCAAATGCAGGGGTCATTTACGGCCTCTTTTCAGTATCAGTTCACTCACCCACGGACGGCGTTTTTCCATCAGCAAAATCAAATTAcgcatc                                                      >  2:1346054/1‑150 (MQ=255)
                                                                                                  accGCCAGCCAGACAAAAAAGGCGTAACCGCCCATTGCGAAAAATTCATTCCAGGAAGCAAATGCAGGGGTCATTTACGGCCTCTTTTCAGTATCAGTTCACTCACCCACGGACGGCGTTTTTCCATCAGCAAAATCAAATTAcgcatcc                                                     >  1:460041/1‑150 (MQ=255)
                                                                                                   ccGCCAGCCAGACAAACAAGGCGTAACCGCCCATTGCGAAAAATTCATTCCAGGAAGCAAATGCAGGGGTCATTTACGGCCTCTTTTCAGTATCAGTTCACTCACCCACGGACGGCGTTTTTCCATCAGCAAAATCAAATTAcgcatccg                                                    <  1:2895933/150‑1 (MQ=255)
                                                                                                   ccGCCAGCCAGACAAAAAAGGCGTAACCGCCCATTGCGAAAAATTCATTCCAGGAAGCAAATGCAGGGGTCATTTACGGCCTCTTTTCAGTATCAGTTCACTCACCCACGGACGGCGTTTTTCCATCAGCAAAATCAAATTAcgcatccg                                                    <  1:1124282/150‑1 (MQ=255)
                                                                                                   ccGCCAGCCAGACAAAAAAGGCGTAACCGCCCATTGCGAAAAATTCATTCCAGGAAGCAAATGCAGGGGTCATTTACGGCCTCTTTTCAGTATCAGTTCACTCACCCACGGACGGCGTTTTTCCATCAGCAAAATCAAATTAcgcatccg                                                    <  1:1302957/150‑1 (MQ=255)
                                                                                                    cGCCAGCCAGACAAAAAAGGCGTAACCGCCCATTGCGAAAAATTCATTCCAGGAAGCAAATGCAGGGGTCATTTACGGCCTCTTTTCAGTATCAGTTCACTCACCCACGGACGGCGTTTTTCCATCAGCAAAATCAAATTAcgcatccgc                                                   >  2:3267297/1‑150 (MQ=255)
                                                                                                      ccagccagACAAAAAAGGCGTAACCGCCCATTGCGAAAAATTCATTCCAGGAAGCAAATGCAGGGGTCATTTACGGCCTCTTTTCAGTATCAGTTCACTCACCCACGGACGGCGTTTTTCCATCAGCAAAATCAAATTAcgcatccgcat                                                 >  2:4425325/1‑150 (MQ=255)
                                                                                                        agccagACAAAAAAGGCGTAACCGCCCATTGCGAAAAATTCATTCCAGGAAGCAAATGCAGGGGTCATTTACGGCCTCTTTTCAGTATCAGTTCACTCACCCACGGACGGCGTTTTTCCATCAGCAAAATCAAATTACGCATCCGCAtca                                               <  1:229284/150‑1 (MQ=255)
                                                                                                           cagACAAAAAAGGCGTAACCGCCCATTGCGAAAAATTCATTCCAGGAAGCAAATGCAGGGGTCATTTACGGCCTCTTTTCAGTATCAGTTCACTCACCCACGGACGGCGTTTTTCCATCAGCAAAATCAAATTACGCATCCGCAtcagcg                                            <  1:3049909/150‑1 (MQ=255)
                                                                                                            agACAAAAAAGGCGTAACCGCCCATTGCGAAAAATTCATTCCAGGAAGCAAATGCAGGGGTCATTTACGGCCTCTTTTCAGTATCAGTTCACTCACCCACGGACGGCGTTTTTCCATCAGCAAAATCAAATTACGCATCCGCAtcagcgt                                           >  1:445894/1‑150 (MQ=255)
                                                                                                             gACAAAAAAGGCGTAACCGCCCATTGCGAAAAATTCATTCAAGGAAGCAAATGCAGGGGTCATTTACGGCCTCTTTTCAGTATCAGTTCACTCACCCACGGACGGCGTTTTTCCATCAGCAAAATCAAATTACGCATCCGCAtcagcgtc                                          <  2:4023724/150‑1 (MQ=255)
                                                                                                                aaaaaaGGCGTAACCGCCCATTGCGAAAAATTCATTCCAGGAAGCAAATGCAGGGGTCATTTACGGCCTCTTTTCAGTATCAGTTCACTCACCCACGGACGGCGTTTTTCCATCAGCAAAATCAAATTACGCATCCGCAtcagcgtcagc                                       >  1:1433901/1‑150 (MQ=255)
                                                                                                                aaaaaaGGCGTAACCGCCCATTGCGAAAAATTCATTCCAGGAAGCAAATGCAGGGGTCATTTACGGCCTCTTTTCAGTATCAGTTAACTCACCCACGGACGGCGTTTTTCCATCAGCAAAATCAAATTACGCATCCGCAtcagcgtcagc                                       >  2:2641508/1‑150 (MQ=255)
                                                                                                                   aaaGGCGTAACCGCCCATTGCGAAAAATTCATTCAAGGAAGCAAATGCAGGGGTCATTTACGGCCTCTTTTCAGTATCAGTTCACTCACCCACGGACGGCGTTTTTCCATCAGCAAAATCAAATTACGCATCCGCATCAGCGTCAGCGTg                                    <  1:258328/150‑1 (MQ=255)
                                                                                                                   aaaGGCGTAACCGCCCATGCGAAAAATTCATTCCAGGAAGCAAATGCAGGGGTCATTTACGGCCTCTTTTCAGTATCAGTTCACTCACCCACGGACGGCGTTTTTCCATCAGCAAAATCAAATTACGCATCCGCATCAGCGTCAGCGTgg                                   >  2:4885235/1‑150 (MQ=255)
                                                                                                                    aaGGCGTAACCGCCCATTGCGAAAAATTCATTCCAGGAAGCAAATGCAGGGGTCATTTACGGCCTCTTTTCAGTATCAGTTCACTCACCCACGGACGGCGTTTTTCCATCAGCAAAATCAAATTACGCATCCGCATCAGCGTCAGCGTgg                                   >  2:1517499/1‑150 (MQ=255)
                                                                                                                    aaGGCGTAACCGCCCATTGCGAAAAATTCATTCCAGGAAGCAAATGCAGGGGTCATTTACGGCCTCTTTTCAGTATCAGTTCACTCACCCACGGACGGCGTTTTTCCATCAGCAAAATCAAATTACGCATCCGCATCAGCGTCAGCGTgg                                   <  1:2650432/150‑1 (MQ=255)
                                                                                                                      ggCGTAACCGCCCATTGCGAAAAATTCATTCAAGGAAGCAAATGCAGGGGTCATTTACGGCCTCTTTTCAGTATCAGTTCACTCACCCACGGACGGCGTTTTTCCATCAGCAAAATCAAATTACGCATCCGCATCAGCGTCAGCGTGGCa                                 <  2:3336608/150‑1 (MQ=255)
                                                                                                                        cGTAACCGCCCATTGCGAAAAATTCATTCCAGGAAGCAAATGCAGGGGTCATTTACGGCCTCTTTTCAGTATCAGTTCACTCACCCACGGACGGCGTTTTTCCATCAGCAAAATCAAATTACGCATCCGCATCAGCGTCAGCGTGGCAGa                               <  2:761104/150‑1 (MQ=255)
                                                                                                                          tAACCGCCCATTGCGAAAAATTCATTCCAGGAAGCAAATGCAGGGGTCATTTACGGCCTCTTTTCAGTATCAGTTCACTCACCCACGGACGGCGTTTTTCCATCAGCAAAATCAAATTACGCATCCGCATCAGCGTCAGCGTGGCAGACa                             <  2:4853452/150‑1 (MQ=255)
                                                                                                                          tAACCGCCCATTGCGAAAAATTCATTCAAGGAAGCAAATGCAGGGGTCATTTACGGCCTCTTTTCAGTATCAGTTCACTCACCCACGGACGGCGTTTTTCCATCAGCAAAATCAAATTACGCATCCGCATCAGCGTCAGCGTGGCAGACa                             >  2:3122446/1‑150 (MQ=255)
                                                                                                                            aCCGCCCATTGCGAAAAATTCATTCCAGGAAGCAAATGCAGGGGTCATTTACGGCCTCTTTTCAGTATCAGTTAACTCACCCACGGACGGCGTTTTTCCATCAGCAAAATCAAATTACGCATCCGCATCAGCGTCAGCGTGGCAGACAgg                           >  1:3241009/1‑150 (MQ=255)
                                                                                                                            aCCGCCCATTGCGAAAAATTCATTCAAGGAAGCAAATGCAGGGGTCATTTACGGCCTCTTTTCAGTATCAGTTCACTCACCCACGGACGGCGTTTTTCCATCAGCAAAATCAAATTACGCATCCGCATCAGCGTCAGCGTGGCAGACAgg                           <  2:4472185/150‑1 (MQ=255)
                                                                                                                                cccATTGCGAAAAATTCATTCCAGGAAGCAAATGCAGGGGTCATTTACGGCCTCTTTTCAGTATCAGTTAACTCACCCACGGACGGCGTTTTTCCATCAGCAAAATCAAATTACGCATCCGCATCAGCGTCAGCGTGGCAGACAGGAGCa                       <  1:3130856/150‑1 (MQ=255)
                                                                                                                                 ccATTGCGACAAATTCATTCAAGGAAGCAAATGCCGGGGTCATGTACGGCCTCTTTTCAGTATCAGTGCACTCACCCACGGACGGCGGTTTTCCATCAGCACAATCATATTACGCATCCGCATCAGCGTCAGCGTGGCAGTCAGGCGCAg                      >  2:4925078/1‑150 (MQ=255)
                                                                                                                                 ccATTGCGAAAAATTCATTCCAGGAAGCAAATGCAGGGGTCATTTACGGCCTCTTTTCAGTATCAGTTCACTCACCCACGGACGGCGTTTTTCCATCAGCAAAATCAAATTACGCATCCGCATCAGCGTCAGCGTGGCAGACAGGAGCAg                      >  2:4252696/1‑150 (MQ=255)
                                                                                                                                 ccATTGCGAAAAATTCATTCAAGGAAGCAAATGCAGGGGTCATTTACGGCCTCTTTTCAGTATCAGTTCACTCACCCACGGACGGCGTTTTTCCATCAGCAAAATCAAATTACGCATCCGCATCAGCGTCAGCGTGGCAGACAGGAGCAg                      >  2:1943979/1‑150 (MQ=255)
                                                                                                                                  cATTGCGAAAAATTCATTCAAGGAAGCAAATGCAGGGGTCATTTACGGCCTCTTTTCAGTATCAGTTCACTCACCCACGGACGGCGTTTTTCCATCAGCAAAATCAAATTACGCATCCGCATCAGCGTCAGCGTGGCAGACAGGAGCAgg                     <  2:2878207/150‑1 (MQ=255)
                                                                                                                                   aTTGCGAAAAATTCATTCCAGGAAGCAAATGCAGGGGTCATTTACGGCCTCTTTTCAGTATCAGTTCACTCACCCACGGACGGCGTTTTTCCATCAGCAAAATCAAATTACGCATCCGCATCAGCGTCAGCGTGGCAGACAGGAGCAGGa                    <  2:570540/150‑1 (MQ=255)
                                                                                                                                   aTTGCGAAAAATTCATTCAAGGAAGCAAATGCAGGGGTCATTTACGGCCTCTTTTCAGTATCAGTTCACTCACCCACGGACGGCGTTTTTCCATCAGCAAAATCAAATTACGCATCCGCATCAGCGTCAGCGTGGCAGACAGGAGCAGGa                    >  2:1425569/1‑150 (MQ=255)
                                                                                                                                    ttGCGAAAAATTCATTCCAGGAAGCAAATGCAGGGGTCATTTACGGCCTCTTTTCAGTATCAGTTCACTCACCCACGGACGGCGTTTTTCCATCAGCAAAATCAAATTACGCATCCGCATCAGCGTCAGCGTGGCAGACAGGAGCAGGaa                   <  2:652452/150‑1 (MQ=255)
                                                                                                                                    ttGCGAAAAATTCATTCAAGGAAGCAAATGCAGGGGTCATTTACGGCCTCTTTTCAGTATCAGTTCACTCACCCACGGACGGCGTTTTTCCATCAGCAAAATCAAATTACGCATCCGCATCAGCGTCAGCGTGGCAGACAGGAGCAGGaa                   <  2:3010838/150‑1 (MQ=255)
                                                                                                                                        gAAAAATTCATTCCAGGAAGCAAATGCAGGGGTCATTTACGGCCTCTTTTCAGTATCAGTTCACTCACCCACGGACGGCGTTTTTCCATCAGCAAAATCAAATTACGCATCCGCATCAGCGTCAGCGTGGCAGACAGGAGCAGGAAGCCa               >  2:521649/1‑150 (MQ=255)
                                                                                                                                         aaaaaTTCATTCAAGGAAGCAAATGCAGGGGTCATTTACGGCCTCTTTTCAGTATCAGTTCACTCACCCACGGACGGCGTTTTTCCATCAGCAAAATCAAATTACGCATCCGCATCAGCGTCAGCGTGGCAGACAGGAGCAGGAAGCCaa              <  1:3326713/150‑1 (MQ=255)
                                                                                                                                          aaaaTTCATTCAAGGAAGCAAATGCAGGGGTCATTTACGGCCTCTTTTCAGTATCAGTTCACTCACCCACGGACGGCGTTTTTCCATCAGCAAAATCAAATTACGCATCCGCATCAGCGTCAGCGTGGCAGACAGGAGCAGGAAGCCaaa             <  1:5038261/150‑1 (MQ=255)
                                                                                                                                           aaaTTCATTCAAGGAAGCAAATGCAGGGGTCATTTACGGCCTCTTTTCAGTATCAGTTCACTCACCCACGGACGGCGTTTTTCCATCAGCAAAATCAAATTACGCATCCGCATCAGCGTCAGCGTGGCAGACAGGAGCAGGAAGCCaaaa            >  2:431760/1‑150 (MQ=255)
                                                                                                                                            aaTTCATTCCAGGAAGCAAATGCAGGGGTCATTTACGGCCTCTTTTCAGTATCAGTTCACTCACCCACGGACGGCGTTTTTCCATCAGCAAAATCAAATTACGCATCCGCATCAGCGTCAGCGTGGCAGACAGGAGCAGGAAGCCaaaaa           <  2:3810727/150‑1 (MQ=255)
                                                                                                                                             attcattcCAGGAAGCAAATGCAGGGGTCATTTACGGCCTCTTTTCAGTATCAGTTCACTCACCCACGGACGGCGTTTTTCCATCAGCAAAATCAAATTACGCATCCGCATCAGCGTCAGCGTGGCAGACAGGAGCAGGAAGCCAAAAAt          <  1:203163/150‑1 (MQ=255)
                                                                                                                                                 attcCAGGAAGCAAATGCAGGGGTCATTTACGGCCTCTTTTCAGTATCAGTTCACTCACCCACGGACGGCGTTTTTCCATCAGCAAAATCAAATTACGCATCCGCATCAGCGTCAGCGTGGCAGACAGGAGCAGGAAGCCAAAAATCGAc      >  2:675551/1‑150 (MQ=255)
                                                                                                                                                 attcCAGGAAGCAAATGCAGGGGTCATTTACGGCCTCTTTTCAGTATCAGTTCACTCACCCACGGACGGCGTTTTTCCATCAGCAAAATCAAATTACGCATCCGCATCAGCGTCAGCGTGGCAGACAGGAGCAGGAAGCCAAAAATCGAc      >  2:3253091/1‑150 (MQ=255)
                                                                                                                                                  ttcCAGGAAGCAAATGCAGGGGTCATTTACGGCCTCTTTTCAGTATCAGTTCACTCACCCACGGACGGCGTTTTTCCATCAGCAAAATCAAATTACGCATCCGCATCAGCGTCAGCGTGGCAGACAGGAGCAGGAAGCCAAAAATCGAcc     <  1:3192062/150‑1 (MQ=255)
                                                                                                                                                   tcCAGGAAGCAAATGCAGGGGTCATTTACGGCCTCTTTTCAGTATCAGTTCACTCACCCACGGACGGCGTTTTTCCATCAGCAAAATCAAATTACGCATCCGCATCAGCGTCAGCGTGGCAGACAGGAGCAGGAAGCCAAAAATCGACca    <  1:100488/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CCTCACGCGCCCGCTGTTGCGCCACGCCACGCAGAATTGCGCGATGTTGCATCACCGAGTGCACGACCAAAACCACCAGCGGAATAACGGTCATCACCACCGCCAGCCAGACAAAAAAGGCGTAACCGCCCATTGCGAAAAATTCATTCCAGGAAGCAAATGCAGGGGTCATTTACGGCCTCTTTTCAGTATCAGTTCACTCACCCACGGACGGCGTTTTTCCATCAGCAAAATCAAATTACGCATCCGCATCAGCGTCAGCGTGGCAGACAGGAGCAGGAAGCCAAAAATCGACCAGC  >  CP009273/2288894‑2289192

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: