Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,822,975 C→A 47.4% A339A (GCG→GCT astD ← succinylglutamic semialdehyde dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,822,9750CA47.4% 162.0 / inf 244A339A (GCG→GCTastDsuccinylglutamic semialdehyde dehydrogenase
Reads supporting (aligned to +/- strand):  ref base C (55/73);  new base A (72/44);  total (127/117)
Fisher's exact test for biased strand distribution p-value = 3.23e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 247 total aligned reads displayed.

AAAGTATCATAACGCCAGACGCGCAATAACGGTCCGAACACCTCTTCATCTGGTACGCCAGCAACGCCTGTCATTTCAATGATCCCCGGCGTCAGCAACGATGTCCCTGCTTGTAATAAGCGCGGCGCAAGCAGGGGTCGTCCGCCCATCGCTTCCAGTTGCTGCCATGCAGTAACCACCTGCTGTGCGGCCTGTTCAGAAATCAGCCCGCCAATAAACGGCTGCGGTTCGTCATCCCAGTTGCCCGGCGTTAATCGCTGGCTGACGGCAACCAGACGAGCAAGAAACGCATCGCCCTG  >  CP009273/1822826‑1823124
                                                                                                                                                     |                                                                                                                                                     
aaaGTATCATAACGCCAGACGCGCAATAACGGTCCGAACACCTCTTCATCTGGTACGCCAGCAACGCCTGTCATTTCAATGATCCCCGGCGTCAGCAACGATGTCCCTGCTTGTAATAAGCGCGGCGCAAGCAGGGGTCGTCCGCCCATc                                                                                                                                                       <  2:2615207/150‑1 (MQ=255)
   gTATCATAACGCCAGACGCGCAATAACGGTCCGAACACCTCTTCATCTGGTACGCCAGCAACGCCTGTCATTTCAATGATCCCCGGCGTCAGCAACGATGTCCCTGCTTGTAATAAGCGCGGCGCAAGCAGGGGTCGTCCGCCCATCGCt                                                                                                                                                    <  2:1385980/150‑1 (MQ=255)
     aTCATAACGCCAGACGCGCAATAACGGTCCGAACACCTCTTCATCTGGTACGCCAGCAACGCCTGTCATTTCAATGATCCCCGGCGTCAGCAACGATGTCCCTGCTTGTAATAAGCGCGGCGCAAGCAGGGGTCGTCCGCCCATAGCTTc                                                                                                                                                  <  1:2300059/150‑1 (MQ=255)
      tCATAACGCCAGACGCGCAATAACGGTCCGAACACCTCTTCATCTGGTACGCCAGCAACGCCTGTCATTTCAATGATCCCCGGCGTCAGCAACGATGTCCCTGCTTGTAATAAGCGCGGCGCAAGCAGGGGTCGTCCGCCCATAGCTTcc                                                                                                                                                 <  1:1243439/150‑1 (MQ=255)
        aTAACGCCAGACGCGCAATAACGGTCCGAACACCTCTTCATCTGGTACGCCAGCAACGCCGGTCATTTCAATGATCCCCGGCGTCAGCAACGATGTCCCTGCTTGTAATAAGCGCGGCGCAAGCAGGGGTCGTCCGCCCATAGCTTCCAg                                                                                                                                               >  2:3529782/1‑150 (MQ=255)
         tAACGCCAGACGCGCAATAACGGTCCGAACACCTCTTCATCTGGTACGCCAGCAACGCCTGTCATTTCAATGATCCCCGGCGTCAGCAACGATGTCCCTGCTTGTAATAAGCGCGGCGCAAGCAGGGGTCGTCCGCCCATCGCTTCCAGt                                                                                                                                              >  2:124462/1‑150 (MQ=255)
         tAACGCCAGACGCGCAATAACGGTCCGAACACCTCTTCATCTGGTACGCCAGCAACGCCTGTCATTTCAATGATCCCCGGCGTCAGCAACGATGTCCCTGCTTGTAATAAGCGCGGCGCAAGCAGGGGTCGTCCGCCCATAGCTTCCAGt                                                                                                                                              >  2:388536/1‑150 (MQ=255)
             gCCAGACGCGCAATAACGGTCCGAACACCTCTTCATCTGGTACGCCAGCAACGCCTGTCATTTCAATGATCCCCGGCGTCAGCAACGATGTCCCTGCTTGTAATAAGCGCGGCGCAAGCAGGGGTCGTCCGCCCATAGCTTCCAGTtgct                                                                                                                                          >  1:4112660/1‑150 (MQ=255)
              ccAGACGCGCAATAACGGTCCGAACACCTCTTCATCTGGTACGCCAGCAACGCCTGTCATTTCAATGATCCCCGGCGTCAGCAACGATGTCCCTGCTTGTAATAAGCGCGGCGCAAGCAGGGGTCGTCCGCCCATCGCTTCCAGTtgctg                                                                                                                                         <  2:4080277/150‑1 (MQ=255)
                  aCGCGCAATAACGGTCCGAACACCTCTTCATCTGGTACGCCAGCAACGCCTGTCATTTCAATGATCCCCGGCGTCAGCAACGATGTCCCTGCTTGTAATAAGCGCGGCGCAAGCAGGGGTCGTCCGCCCATAGCTTCCAGTTGCTGCCAt                                                                                                                                     <  2:1054094/150‑1 (MQ=255)
                       cAATAACGGTCCGAACACCTCTTCATCTGGTACGCCAGCAACGCCTGTCATTTCAATGATCCCCGGCGTCAGCAACGATGTCCCTGCTTGTAATAAGCGCGGCGCAAGCAGGGGTCGTCCGCCCATCGCTTCCAGTTGCTGCCATGCAGt                                                                                                                                >  2:3431120/1‑150 (MQ=255)
                        aaTAACGGTCCGAACACCTCTTCATCTGGTACGCCAGCAACGCCTGTCATTTCAATGATCCCCGGCGTCAGCAACGATGTCCCTGCTTGTAATAAGCGCGGCGCAAGCAGGGGTCGTCCGCCCATCGCTTCCAGTTGCTGCCATGCAGTa                                                                                                                               >  2:3785173/1‑150 (MQ=255)
                        aaTAACGGTCCGAACACCTCTTCATCTGGTACGCCAGCAACGCCTGTCATTTCAATGATCCCCGGCGTCAGCAACGATGTCCCTGCTTGTAATAAGCGCGGCGCAAGCAGGGGTCGTCCGCCCATAGCTTCCAGTTGCTGCCATGCAGTa                                                                                                                               >  1:3377040/1‑150 (MQ=255)
                        aaTAACGGTCCGAACACCTCTTCATCTGGTACGCCAGCAACGCCTGTCATTTCAATGATCCCCGGCGTCAGCAACGATGTCCCTGCTTGTAATAAGCGCGGCGCAAGCAGGGGTCGTCCGCCCATAGCTTCCAGTTGCTGCCATGCAGTa                                                                                                                               <  1:3754402/150‑1 (MQ=255)
                         aTAACGGTCCGAACACCTCTTCATCTGGTACGCCAGCAACGCCTGTCATTTCAATGATCCCCGGCGTCAGCAACGATGTCCCTGCTTGTAATAAGCGCGGCGCAAGCAGGGGTCGTCCGCCCATCGCTTCCAGTTGCTGCCATGCAGTaa                                                                                                                              <  1:4055282/150‑1 (MQ=255)
                         aTAACGGTCCGAACACCTCTTCATCTGGTACGCCAGCAACGCCTGTCATTTCAATGATCCCCGGCGTCAGCAACGATGTCCCTGCTTGTAATAAGCGCGGCGCAAGCAGGGGTCGTCCGCCCATAGCTTCCAGTTGCTGCCATGCAGTaa                                                                                                                              >  1:330920/1‑150 (MQ=255)
                          tAACGGTCCGAACACCTCTTCATCTGGTACGCCAGCAACGCCTGTCATTTCAATGATCCCCGGCGTCAGCAACGATGTCCCTGCTTGTAATAAGCGCGGCGCAAGCAGGGGTCGTCCGCCCATCGCTTCCAGTTGCTGCCATGCAGTAac                                                                                                                             >  1:3575373/1‑150 (MQ=255)
                          tAACGGTCCGAACACCTCTTCATCTGGTACGCCAGCAACGCCTGTCATTTCAATGATCCCCGGCGTCAGCAACGATGTCCCTGCTTGTAATAAGCGCGGCACAAGCAGGGGTCGTCCGCCCATCGCTTCCAGTTGCTGCCATGCAGTAac                                                                                                                             <  2:2343442/150‑1 (MQ=255)
                           aaCGGTCCGAACACCTCTTCATCTGGTACGCCAGCAACGCCTGTCATTTCAATGATCCCCGGCGTCAGCAACGATGTCCCTGCTTGTAATAAGCGCGGCGCAAGCAGGGGTCGTCCGCCCATAGCTTCCAGTTGCTGCCATGCAGTAacc                                                                                                                            <  2:3756302/150‑1 (MQ=255)
                              ggTCCGAACACCTCTTCATCTGGTACGCCAGCAACGCCTGTCATTTCAATGATCCCCGGCGTCAGCAACGATGTCCCTGCTTGTAATAAGCGCGGCGCAAGCAGGGGTCGTCCGCCCATCGCTTCCAGTTGCTGCCATGCAGTAaccacc                                                                                                                         <  1:1550169/150‑1 (MQ=255)
                              ggTCCGAACACCTCTTCATCTGGTACGCCAGCAACGCCTGTCATTTCAATGATCCCCGGCGTCAGCAACGATGTCCCTGCTTGTAATAAGCGCGGCGCAAGCAGGGGTCGTCCGCCCATAGCTTCCAGTTGCTGCCATGCAGTAaccacc                                                                                                                         >  1:2981261/1‑150 (MQ=255)
                                 ccGAACACCTCTTCATCTGGTACGCCAGCAACGCCTGTCATTTCAATGATCCCCGGCGTCAGCAACGATGTCCCTGCTTGTAATAAGCGCGGCGCAAGCAGGGGTCGTCCGCCCATAGCTTCCAGTTGCTGCCATGCAGTAACCACctgc                                                                                                                      >  1:3268540/1‑150 (MQ=255)
                                    aaCACCTCTTCATCTGGTACGCCAGCAACGCCTGTCATTTCAATGATCCCCGGCGTCAGCAACGATGTCCCTGCTTGTAATAAGCGCGGCGCAAGCAGGGGTCGTCCGCCCATCGCTTCCAGTTGCTGCCATGCAGTAACCACCTGCtgt                                                                                                                   <  2:969156/150‑1 (MQ=255)
                                    aaCACCTCTTCATCTGGTACGCCAGCAACGCCTGTCATTTCAATGATCCCCGGCGTCAGCAACGATGTCCCTGCTTGTAATAAGCGCGGCGCAAGCAGGGGTCGTCCGCCCATAGCTTCCAGTTGCTGCCATGCAGTAACCACCTGCtgt                                                                                                                   >  2:3372441/1‑150 (MQ=255)
                                       acCTCTTCATCTGGTACGCCAGCAACGCCTGTCATTTCAATGATCCCCGGCGTCAGCAACGATGTCCCTGCTTGTAATAAGCGCGGCGCAAGCAGGGGTCGTCCGCCCATCGCTTCCAGTTGCTGCCATGCAGTAACCACCTGCTGTGCg                                                                                                                <  2:1852253/150‑1 (MQ=255)
                                       acCTCTTCATCTGGTACGCCAGCAACGCCTGTCATTTCAATGATCCCCGGCGTCAGCAACGATGTCCCTGCTTGTAATAAGCGCGGCGCAAGCAGGGGTCGTCCGCCCATAGCTTCCAGTTGCTGCCATGCAGTAACCACCTGCTGTGCg                                                                                                                >  2:1672437/1‑150 (MQ=255)
                                            ttCATCTGGTACGCCAGCAACGCCTGTCATTTCAATGATCCCCGGCGTCAGCAACGATGTCCCTGCTTGTAATAAGCGCGGCGCAAGCAGGGGTCGTCCGCCCATCGCTTCCAGTTGCTGCCATGCAGTAACCACCTGCTGTGCGGCCTg                                                                                                           <  1:4088852/150‑1 (MQ=255)
                                            ttCATCTGGTACGCCAGCAACGCCTGTCATTTCAATGATCCCCGGCGTCAGCAACGATGTCCCTGCTTGTAATAAGCGCGGCGCAAGCAGGGGTCGTCCGCCCATAGCTTCCAGTTGCTGCCATGCAGTAACCACCTGCTGTGCGGCCTg                                                                                                           >  1:2165146/1‑150 (MQ=255)
                                              cATCTGGTATGCCAGCAACGCCTGTCATTTCAATGATCCCCGGCGTCAGCAACGATGTCCCTGCTTGTAATAAGCGCGGCGCAAGCAGGGGTCGTCCGCCCATAGCTTCCAGTTGCTGCCATGCAGTAACCACCTGCTGTGCGGCCTGtt                                                                                                         >  2:2576444/1‑150 (MQ=255)
                                                 cTGGTACGCCAGCAACGCCTGTCATTTCAATGATCCCCGGCGTCAGCAACGATGTCCCTGCTTGTAATAAGCGCGGCGCAAGCAGGGGTCGTCCGCCCATAGCTTCCAGTTGCTGCCATGCAGTAACCACCTGCTGTGCGGCCTGTTCAg                                                                                                      >  1:8205/1‑150 (MQ=255)
                                                 cTGGTACGCCAGCAACGCCTGTCATTTCAATGATCCCCGGCGTCAGCAACGATGTCCCTGCTTGTAATAAGCGCGGCGCAAGCAGGGGTCGTCCGCCCATAGCTTCCAGTTGCTGCCATGCAGTAACCACCTGCTGTGCGGCCTGTTCAg                                                                                                      >  1:1575375/1‑150 (MQ=255)
                                                  tGGTACGCCAGCAACGCCTGTCATTTCAATGATCCCCGGCGTCAGCAACGATGTCCCTGCTTGTAATAAGCGCGGCGCAAGCAGGGGTCGTCCGCCCATAGCTTCCAGTTGCTGCCATGCAGTAACCACCTGCTGTGCGGCCTGTTCAGa                                                                                                     <  1:3488829/150‑1 (MQ=255)
                                                     tACGCCAGCAACGCCTGTCATTTCAATGATCCCCGGCGTCAGCAACGATGTCCCTGCTTGTAATAAGCGCGGCGCAAGCAGGGGTCGTCCGCCCATAGCTTCCAGTTGCTGCCATGCAGTAACCACCTGCTGTGCGGCCTGTTCAGAAAt                                                                                                  >  1:393698/1‑150 (MQ=255)
                                                     tACGCCAGCAACGCCTGTCATTTCAATGATCCCCGGCGTCAGCAACGATGTCCCTGCTTGTAATAAGCGCGGCGCAAGCAGGGGTCGTCCGCCCATAGCTTCCAGTTGCTGCCATGCAGTAACCACCTGCTGTGCGGCCTGTTCAGAAAt                                                                                                  >  1:393603/1‑150 (MQ=255)
                                                          cAGCAACGCCTGTCATTTCAATGATCCCCGGCGTCAGCAACGATGTCCCTGCTTGTAATAAGCGCGGCGCAAGCAGGGGTCGACCGCCCATAGCTTCCAGTTGCTGCCATGCAGTAACCACCAGCTGTGCGGCCTGTTCAGAAATCAGcc                                                                                             >  1:923095/1‑150 (MQ=255)
                                                           aGCAACGCCTGTCATTTCAATGATCCCCGGCGTCAGCAACGATGTCCCTGCTTGTAATAAGCGCGGCGCAAGCAGGGGTCGTCCGCCCATAGCTTCCAGTTGCTGCCATGCAGTAACCACCTGCTGTGCGGCCTGTTCAGAAATCAGccc                                                                                            <  2:3924541/150‑1 (MQ=255)
                                                             cAACGCCTGTCATTTCAATGATCCCCGGCGTCAGCAACGATGTCCCTGCTTGTAATAAGCGCGGCGCAAGCAGGGGTCGTCCGCCCATCGCTTCCAGTTGCTGCCATGCAGTAACCACCTGCTGTGCGGCCTGTTCAGAAATCAGCCCGc                                                                                          >  2:3227161/1‑150 (MQ=255)
                                                               aCGCCTGTCATTTCAATGATCCCCGGCGTCAGCAACGATGTCCCTGCTTGTAATAAGCGCGGCGCAAGCAGGGGTCGTCCGCCCATCGCTTCCAGTTGCTGCCATGCAGTAACCACCTGCTGTGCGGCCTGTTCAGAAATCAGCCCGCCa                                                                                        >  1:804690/1‑150 (MQ=255)
                                                                 gCCTGTCATTTCAATGATCCCCGGCGTCAGCAACGATGTCCCTGCTTGTAATAAGCGCGGCGCAAGCAGGGGTCGTCCGCCCATCGCTTCCAGTTGCTGCCATGCAGTAACCACCTGCTGTGCGGCCTGTTCAGAAATCAGCCCGCCAAt                                                                                      >  1:1301515/1‑150 (MQ=255)
                                                                 gCCTGTCATTTCAATGATCCCCGGCGTCAGCAACGATGTCCCTGCTTGTAATAAGCGCGGCGCAAGCAGGGGTCGTCCGCCCATCGCTTCCAGTTGCTGCCATGCAGTAACCACCTGCTGTGCGGCCTGTTCAGAAATCAGCCCGCCAAt                                                                                      >  1:1985428/1‑150 (MQ=255)
                                                                 gCCTGTCATTTCAATGATCCCCGGCGTCAGCAACGATGTCCCTGCTTGTAATAAGCGCGGCGCAAGCAGGGGTCGTCCGCCCATCGCTTCCAGTTGCTGCCATGCAGTAACCACCTGCTGTGCGGCCTGTTCAGAAATCAGCCCGCCAAt                                                                                      <  1:3278883/150‑1 (MQ=255)
                                                                  ccTGTCATTTCAATGATCCCCGGCGTCAGCAACGATGTCCCTGCTTGTAATAAGCGCGGCGCAAGCAGGGGTCGTCCGCCCATAGCTTCCAGTTGCTGCCATGCAGTAACCACCTGCTGTGCGGCCTGTTCAGAAATCAGCCCGCCAATa                                                                                     >  2:932962/1‑150 (MQ=255)
                                                                   cTGTCATTTCAATGATCCCCGGCGTCAGCAACGATGTCCCTGCTTGTAATAAGCGCGGCGCAAGCAGGGGTCGTCCGCCCATCGCTTCCAGTTGCTGCCATGCAGTAACCACCTGCTGTGCGGCCTGTTCAGAAATCAGCCCGCCAATaa                                                                                    <  1:3046844/150‑1 (MQ=255)
                                                                     gTCATTTCAATGATCCCCGGCGTCAGCAACGATGTCCCTGCTTGTAATAAGCGCGGCGCAAGCAGGGGTCGTCCGCCCATCGCTTCCAGTTGCTGCCATGCAGTAACCACCTGCTGTGCGGCCTGTTCAGAAATCAGCCCGCCAATAAAc                                                                                  >  1:2664165/1‑150 (MQ=255)
                                                                           tCAATGATCCCCGGCGTCAGCAACGATGTCCCTGCTTGTAATAAGCGCGGCGCAAGCAGGGGTCGTCCGCCCATCGCTTCCAGTTGCTGCCATGCAGTAACCACCTGCTGTGCGGCCTGTTCAGAAATCAGCCCGCCAATAAACGGCTGc                                                                            >  1:3693268/1‑150 (MQ=255)
                                                                                gATCCCCGGCGTCAGCAACGATGTCCCTGCTTGTAATAAGCGCGGCGCAAGCAGGGGTCGTCCGCCCATAGCTTCCAGTTGCTGCCATGCAGTAACCACCTGCTGTGCGGCCTGTTCAGAAATCAGCCCGCCAATAAACGGCTGCGGTTc                                                                       >  1:2813951/1‑150 (MQ=255)
                                                                                gATCCCCGGCGTCAGCAACGATGTCCCTGCTTGTAATAAGCGCGGCGCAAGCAGGGGTCGTCCGCCCATAGCTTCCAGTTGCTGCCATGCAGTAACCACCTGCTGTGCGGCCTGTTCAGAAATCAGCCCGCCAATAAACGGCTGCGGTTc                                                                       >  1:2813931/1‑150 (MQ=255)
                                                                                 aTCCCCGGTGTCAGCAACGATGTCCCTGCTTGTAATAAGCGCGGCGCAAGCAGGGGTCGTCCGCCCATAGCTTCCAGTTGCTGCCATGCAGTAACCACCTGCTGTGCGGCCTGTTCAGAAATCAGCCCGCCAATAAACGGCTGCGGTTCg                                                                      >  2:623255/1‑150 (MQ=255)
                                                                                      cGGCGTCAGCAACGATGTCCCTGCTTGTAATAAGCGCGGCGCAAGCAGGGGTCGTCCGCCCATCGCTTCCAGTTGCTGCCATGCAGTAACCACCTGCTGTGCGGCCTGTTCAGAAATCAGCCCGCCAATAAACGGCTGCGGTTCGTCATc                                                                 >  2:2212598/1‑150 (MQ=255)
                                                                                            cAGCAACGATGTCCCTGCTTGTAATAAGCGCGGCGCAAGCAGGGGTCGTCCGCCCATCGCTTCCAGTTGCTGCCATGCAGTAACCACCTGCTGTGCGGCCTGTTCAGAAATCAGCCCGCCAATAAACGGCTGCGGTTCGTCATCCCAGtt                                                           <  1:3542382/150‑1 (MQ=255)
                                                                                             aGCAACGATGTCCCTGCTTGTAATAAGCGCGGCGCAAGCAGGGGTCGTCCGCCCATAGCTTCCAGTTGCTGCCATGCAGTAACCACCTGCTGTGCGGCCTGTTCAGAAATCAGCCCGCCAATAAACGGCTGCGGTTCGTCATCCCAGTTg                                                          >  2:1792001/1‑150 (MQ=255)
                                                                                             aGCAACGATGTCCCTGCTTGTAATAAGCGCGGCGCAAGCAGGGGTCGTCCGCCCATAGCTTCCAGTTGCTGCCATGCAGTAACCACCTGCTGTGCGGCCTGTTCAGAAATCAGCCCGCCAATAAACGGCTGCGGTTCGTCATCCCAGTTg                                                          >  1:556885/1‑150 (MQ=255)
                                                                                               cAACGATGTCCCTGCTTGTAATAAGCGCGGCGCAAGCAGGGGTCGTCCGCCCATCGCTTCCAGTTGCTGCCATGCAGTAACCACCTGCTGTGCGGCCTGTTCAGAAATCAGCCCGCCAATAAACGGCTGCGGTTCGTCATCCCAGTTGcc                                                        >  1:705302/1‑150 (MQ=255)
                                                                                                 aCGATGTCCCTGCTTGTAATAAGCGCGGCGCAAGCAGGGGTCGTCCGCCCATAGCTTCCAGTTGCTGCCATGCAGTAACCACCTGCTGTGCGGCCTGTTCAGAAATCAGCCCGCCAATAAACGGCTGCGGTTCGTCATCCCAGTTGCCCg                                                      <  1:2627680/150‑1 (MQ=255)
                                                                                                   gATGTCCCTGCTTGTAATAAGCGCGGCGCAAGCAGGGGTCGTCCGCCCATCGCTTCCAGTTGCTGCCATGCAGTAACCACCTGCTGTGCGGCCTGTTCAGAAATCAGCCCGCCAATAAACGGCTGCGGTTCGTCATCCCAGTTGCCCGGc                                                    >  2:2847133/1‑150 (MQ=255)
                                                                                                   gATGTCCCTGCTTGTAATAAGCGCGGCGCAAGCAGGGGTCGTCCGCCCATAGCTTCCAGTTGCTGCCATGCAGTAACCACCTGCTGTGCGGCCTGTTCAGAAATCAGCCCGCCAATAAACGGCTGCGGTTCGTCATCCCAGTTGCCCGGc                                                    >  2:3629757/1‑150 (MQ=255)
                                                                                                    aTGTCCCTGCTTGTAATAAGCGCGGCGCAAGCAGGGGTCGTCCGCCCATAGCTTCCAGTTGCTGCCATGCAGTAACCACCTGCTGTGCGGCCTGTTCAGAAATCAGCCCGCCAATAAACGGCTGCGGTTCGTCATCCCAGTTGCCCGGCg                                                   >  2:840378/1‑150 (MQ=255)
                                                                                                     tGTCCCTGCTTGTAATAAGCGCGGCGCAAGCAGGGGTCGTCCGCCCATCGCTTCCAGTTGCTGCCATGCAGTAACCACCTGCTGTGCGGCCTGTTCAGAAATCAGCCCGCCAATAAACGGCTGCGGTTCGTCATCCCAGTTGCCCGGCGt                                                  >  2:3317422/1‑150 (MQ=255)
                                                                                                       tCCCTGCTTGTAATAAGCGCGGCGCAAGCAGGGGTCGTCCGCCCATCGCTTCCAGTTGCTGCCATGCAGTAACCACCTGCTGTGCGGCCTGTTCAGAAATCAGCCCGCCAATAAACGGCTGCGGTTCGTCATCCCAGTTGCCCGGCGTTa                                                <  1:831746/150‑1 (MQ=255)
                                                                                                          cTGCTTGTAATAAGCGCGGCGCAAGCAGGGGTCGTCCGCCCATCGCTTCCAGTTGCTGCCATGCAGTAACCACCTGCTGTGCGGCCTGTTCAGAAATCAGCCCGCCAATAAACGGCTGCGGTTCGTCATCCCAGTTGCCCGGCGTTAATc                                             <  1:1453477/150‑1 (MQ=255)
                                                                                                          cTGCTTGTAATAAGCGCGGCGCAAGCAGGGGTCGTCCGCCCATAGCTTCCAGTTGCTGCCATGCAGTAACCACCTGCTGTGCGGCCTGTTCAGAAATCAGCCCGCCAATAAACGGCTGCGGTTCGTCATCCCAGTTGCCCGGCGTTAATc                                             >  1:1400837/1‑150 (MQ=255)
                                                                                                             cTTGTAATAAGCGCGGCGCAAGCAGGGGTCGTCCGCCCATCGCTTCCAGTTGCTGCCATGCAGTAACCACCTGCTGTGCGGCCTGTTCAGAAATCAGCCCGCCAATAAACGGCTGCGGTTCGTCATCCCAGTTGCCCGGCGTTAATCgct                                          <  1:1018912/150‑1 (MQ=255)
                                                                                                             cTTGTAATAAGCGCGGCGCAAGCAGGGGTCGTCCGCCCATCGCTTCCAGTTGCTGCCATGCAGTAACCACCTGCTGTGCGGCCTGTTCAGAAATCAGCCCGCCAATAAACGGCTGCGGTTCGTCATCCCAGTTGCCCGGCGTTAATCgct                                          >  2:792757/1‑150 (MQ=255)
                                                                                                             cTTGTAATAAGCGCGGCGCAAGCAGGGGTCGTCCGCCCATCGCTTCCAGTTGCTGCCATGCAGTAACCACCTGCTGTGCGGCCTGTTCAGAAATCAGCCCGCCAATAAACGGCTGCGGTTCGTCATCCCAGTTGCCCGGCGTTAATCgct                                          <  1:1737000/150‑1 (MQ=255)
                                                                                                             cTTGTAATAAGCGCGGCGCAAGCAGGGGTCGTCCGCCCATAGCTTCCAGTTGCTGCCATGCAGTAACCACCTGCTGTGCGGCCTGTTCAGAAATCAGCCCGCCAATAAACGGCTGCGGTTCGTCATCCCAGTTGCCCGGCGTTAATCgct                                          <  2:3633343/150‑1 (MQ=255)
                                                                                                             cTTGTAATAAGCGCGGCGCAAGCAGGGGTCGTCCGCCCATAGCTTCCAGTTGCTGCCATGCAGTAACCACCTGCTGTGCGGCCTGTTCAGAAATCAGCCCGCCAATAAACGGCTGCGGTTCGTCATCCCAGTTGCCCGGCGTTAATCgct                                          <  1:1659186/150‑1 (MQ=255)
                                                                                                             cTTGTAATAAGCGCGGCGCAAGCAGGGGTCGTCCGCCCATAGCTTCCAGTTGCTGCCATGCAGTAACCACCTGCTGTGCGGCCTGTTCAGAAATCAGCCCGCCAATAAACGGCTGCGGTTCGTCATCCCAGTTGCCCGGCGTTAATCgct                                          >  1:2636823/1‑150 (MQ=255)
                                                                                                               tGTAATAAGCGCGGCGCAAGCAGGGGTCGTCCGCCCATCGCTTCCAGTTGCTGCCATGCAGTAACCACCTGCTGTGCGGCCTGTTCAGAAATCAGCCCGCCAATAAACGGCTGCGGTTCGTCATCCCAGTTGCCCGGCGTTAATCgctgg                                        >  1:3998212/1‑150 (MQ=255)
                                                                                                                 taataaGCGCGGCGCAAGCAGGGGTCGTCCGCCCATCGCTTCCAGTTGCTGCCATGCAGTAACCACCTGCTGTGCGGCCTGTTCAGAAATCAGCCCGCCAATAAACGGCTGCGGTTCGTCATCCCAGTTGCCCGGCGTTAATCgctggct                                      <  2:2682941/150‑1 (MQ=255)
                                                                                                                 taataaGCGCGGCGCAAGCAGGGGTCGTCCGCCCATAGCTTCCAGTTGCTGCCATGCAGTAACCACCTGCTGTGCGGCCTGTTCAGAAATCAGCCCGCCAATAAACGGCTGCGGTTCGTCATCCCAGTTGCCCGGCGTTAATCgctggct                                      <  1:2249784/150‑1 (MQ=255)
                                                                                                                  aataaGCGCGGCGCAAGCAGGGGTCGTCCGCCCATCGCTTCCAGTTGCTGCCATGCAGTAACCACCTGCTGTGCGGCCTGTTCAGAAATCAGCCCGCCAATAACCGGCTGCGGTTCGTCATCCCAGTTGCCCGGCGTTAATCgctggctg                                     <  2:1784025/150‑1 (MQ=255)
                                                                                                                  aataaGCGCGGCGCAAGCAGGGGTCGTCCGCCCATAGCTTCCAGTTGCTGCCATGCAGTAACCACCTGCTGTGCGGCCTGTTCAGAAATCAGCCCGCCAATAAACGGCTGCGGTTCGTCATCCCAGTTGCCCGGCGTTAATCgctggctg                                     <  1:3545608/150‑1 (MQ=255)
                                                                                                                   ataaGCGCGGCGCAAGCAGGGGTCGTCCGCCCATCGCTTCCAGTTGCTGCCATGCAGTAACCACCTGCTGTGCGGCCTGTTCAGAAATCAGCCCGCCAATAAACGGCTGCGGTTCGTCATCCCAGTTGCCCGGCGTTAATCGCTGGCTGa                                    <  2:2484136/150‑1 (MQ=255)
                                                                                                                   ataaGCGCGGCGCAAGCAGGGGTCGTCCGCCCATAGCTTCCAGTTGCTGCCATGCAGTAACCACCTGCTGTGCGGCCTGTTCAGAAATCAGCCCGCCAATAAACGGCTGCGGTTCGTCATCCCAGTTGCCCGGCGTTAATCGCTGGCTGa                                    <  2:1522965/150‑1 (MQ=255)
                                                                                                                    taaGCGCGGCGCAAGCAGGGGTCGTCCGCCCATCGCTTCCAGTTGCTGCCATGCAGTAACCACCTGCTGTGCGGCCTGTTCAGAAATCAGCCCGCCAATAAACGGCTGCGGTTCGTCATCCCAGTTGCCCGGCGTTAATCGCTGGCTGAc                                   >  2:3996808/1‑150 (MQ=255)
                                                                                                                    taaGCGCGGCGCAAGCAGGGGTCGTCCGCCCATAGCTTCCAGTTGCTGCCATGCAGTAACCACCTGCTGTGCGGCCTGTTCAGAAATCAGCCCGCCAATAAACGGCTGCGGTTCGTCATCCCAGTTGCCCGGCGTTAATCGCTGGCTGAc                                   <  2:1018216/150‑1 (MQ=255)
                                                                                                                     aaGCGCGGCGCAAGCAGGGGTCGTCCGCCCATCGCTTCCAGTTGCTGCCATGCAGTAACCACCTGCTGTGCGGCCTGTTCAGAAATCAGCCCGCCAATAAACGGCTGCGGTTCGTCATCCCAGTTGCCCGGCGTTAATCGCTGGCTGACg                                  >  1:2451744/1‑150 (MQ=255)
                                                                                                                     aaGCGCGGCGCAAGCAGGGGTCGTCCGCCCATAGCTTCCAGTTGCTGCCATGCAGTAACCACCTGCTGTGCGGCCTGTTCAGAAATCAGCCCGCCAATAAACGGCTGCGGTTCGTCATCCCAGTTGCCCGGCGTTAATCGCTGGCTGACg                                  <  2:2665770/150‑1 (MQ=255)
                                                                                                                        cgcgGCGCAAGCAGGGGTCGTCCGCCCATCGCTTCCAGTTGCTGCCATGCAGTAACCACCTGCTGTGCGGCCTGTTCAGAAATCAGCCCGCCAATAAACGGCTGCGGTTCGTCATCCCAGTTGCCCGGCGTTAATCGCTGGCTGACGGCa                               <  1:1340664/150‑1 (MQ=255)
                                                                                                                           ggcgCAAGCAGGGGTCGTCCGCCCATCGCTTCCAGTTGCTGCCATGCAGTAACCACCTGCTGTGCGGCCTGTTCAGAAATCAGCCCGCCAATAAACGGCTGCGGTTCGTCATCCCAGTTGCCCGGCGTTAATCGCTGGCTGACGGCAAcc                            <  1:2254093/150‑1 (MQ=255)
                                                                                                                             cgcAAGCAGGGGTCGTCCGCCCATCGCTTCCAGTTGCTGCCATGCAGTAACCACCTGCTGTGCGGCCTGTTCAGAAATCAGCCCGCCAATAAACGGCTGCGGTTCGTCATCCCAGTTGCCCGGCGTTAATCGCTGGCTGACGGCAACCAg                          <  1:472834/150‑1 (MQ=255)
                                                                                                                                aaGCAGGGGTCGTCCGCCCATCGCTTCCAGTTGCTGCCATGCAGTAACCACCTGCTGTGCGGCCTGTTCAGAAATCAGCCCGCCAATAAACGGCTGCGGTTCGTCATCCCAGTTGCCCGGCGTTAATCGCTGGCTGACGGCAACCAGACg                       >  1:342307/1‑150 (MQ=255)
                                                                                                                                aaGCAGGGGTCGTCCGCCCATCGCTTCCAGTTGCTGCCATGCAGTAACCACCTGCTGTGCGGCCTGTTCAGAAATCAGCCCGCCAATAAACGGCTGCGGTTCGTCATCCCAGTTGCCCGGCGTTAATCGCTGGCTGACGGCAACCAGACg                       >  2:211147/1‑150 (MQ=255)
                                                                                                                                aaGCAGGGGTCGTCCGCCCATCGCTTCCAGTTGCTGCCATGCAGTAACCACCTGCTGTGCGGCCTGTTCAGAAATCAGCCCGCCAATAAACGGCTGCGGTTCGTCATCCCAGTTGCCCGGCGTTAATCGCTGGCTGACGGCAACCAGACg                       <  2:109364/150‑1 (MQ=255)
                                                                                                                                aaGCAGGGGTCGTCCGCCCATAGCTTCCAGTTGCTGCCATGCAGTAACCACCTGCTGTGCGGCCTGTTCAGAAATCAGCCCGCCAATAAACGGCTGCGGTTCGTCATCCCAGTTGCCCGGCGTTAATCGCTGGCTGACGGCAACCAGACg                       <  2:850503/150‑1 (MQ=255)
                                                                                                                                 cgCAGGGGTCGTCCGCCCATAGCTTCCAGTTGCTGCCATGCTGTAACCACCTGCTGTGCGGCCTGTTCAGAAATCAGCCCGCCAATAAACGGCTGCGGTTCGTCATCCCAGTTGCCCGGCGTTAATCGCTGGCTGACGGCAACCTGACGa                      <  1:1766611/149‑1 (MQ=255)
                                                                                                                                  gCAGGGGTCGTCCGCCCATAGCTTCCAGTTGCTGCCATGCAGTAACCACCTGCTGTGCGGCCTGTTCAGAAATCAGCCCGCCAATAAACGGCTGCGGTTCGTCATCCCAGTTGCCCGGCGTTAATCGCTGGCTGACGGCAACCAGACGAg                     <  1:2572415/150‑1 (MQ=255)
                                                                                                                                    aGGGGTCGTCCGCCCATAGCTTCCAGTTGCTGCCATGCAGTAACCACCTGCTGTGCGGCCTGTTCAGAAATCAGCCCGCCAATAAACGGCTGCGGTTCGTCATCCCAGTTGCCCGGCGTTAATCGCTGGCTGACGGCAACCAGACGAGCa                   >  1:605695/1‑150 (MQ=255)
                                                                                                                                     ggggTCGTCCGCCCATAGCTTCCAGTTGCTGCCATGCAGTAACCACCTGCTGTGCGGCCTGTTCAGAAATCAGCCCGCCAATAAACGGCTGCGGTTCGTCATCCCAGTTGCCCGGCGTTAATCGCTGGCTGACGGCAACCAGACGAGCaa                  >  2:2749418/1‑150 (MQ=255)
                                                                                                                                       ggTCGTCCGCCCATCGCTTCCAGTTGCTGCCATGCAGTAACCACCTGCTGTGCGGCCTGTTCAGAAATCAGCCCGCCAATAAACGGCTGCGGTTCGTCATCCCAGTTGCCCGGCGTTAATCGCTGGCTGACGGCAACCAGACGAGCAAGa                <  1:1220430/150‑1 (MQ=255)
                                                                                                                                         tcgtcCGCCCATCGCTTCCAGTTGCTGCCATGCAGTAACCACCTGCTGTGCGGCCTGTTCAGAAATCAGCCCGCCAATAAACGGCTGCGGTTCGTCATCCCAGTTGCCCGGCGTTAATCGCTGGCTGACGGCAACCAGACGAGCAAGaaa              <  2:889379/150‑1 (MQ=255)
                                                                                                                                          cgtcCGCCCATAGCTTCCAGTTGCTGCCATGCAGTAACCACCTGCTGTGCGGCCTGTTCAGAAATCAGCCCGCCAATAAACGGCTGCGGTTCGTCATCCCAGTTGCCCGGCGTTAATCGCTGGCTGACGGCAACCAGACGAGCAAGAAAc             <  2:2927495/150‑1 (MQ=255)
                                                                                                                                          cgtcCGCCCATAGCTTCCAGTTGCTGCCATGCAGTAACCACCTGCTGTGCGGCCTGTTCAGAAATCAGCCCGCCAATAAACGGCTGCGGTTCGTCATCCCAGTTGCCCGGCGTTAATCGCTGGCTGACGGCAACCAGACGAGCAAGAAAc             <  2:2151348/150‑1 (MQ=255)
                                                                                                                                            tcCGCCCATCGCTTCCAGTTGCTGCCATGCAGTAACCACCTGCTGTGCGGCCTGTTCAGAAATCAGCCCGCCAATAAACGGCTGCGGTTCGTCATCCCAGTTGCCCGGCGTTAATCGCTGGCTGACGGCAACCAGACGAGCAAGAAACGc           >  1:2729007/1‑150 (MQ=255)
                                                                                                                                               gCCCATCGCTTCCAGTTGCTGCCATGCAGTAACCACCTGCTGTGCGGCCTGTTCAGAAATCAGCCCGCCAATAAACGGCTGCGGTTCGTCATCCCAGTTGCCCGGCGTTAATCGCTGGCTGACGGCAACCAGACGAGCAAGAAACGCATc        <  1:933596/150‑1 (MQ=255)
                                                                                                                                               gCCCATAGCTTCCAGTTGCTGCCATGCAGTAACCACCTGCTGTGCGGCCTGTTCAGAAATCAGCCCGCCAATAAACGGCTGCGGTTCGTCATCCCAGTTGCCCGGCGTTAATCGCTGGCTGACGGCAACCAGACGAGCAAGAAACGCATc        <  1:1145218/150‑1 (MQ=255)
                                                                                                                                                cccATAGCTTCCAGTTGCTGCCATGCAGTAACCACCTGCTGTGCGGCCTGTTCAGAAATCAGCCCGCCAATAAACGGCTGCGGTTCGTCATCCCAGTTGCCCGGCGTTAATCGCTGGCTGACGGCAACCAGACGAGCAAGAAACGCATCg       >  2:3617466/1‑150 (MQ=255)
                                                                                                                                                   aTCGCTTCCAGTTGCTGCCATGCAGTAACCACCTGCTGTGCGGCCTGTTCAGAAATCAGCCCGCCAATAAACGGCTGCGGTTCGTCATCCCAGTTGCCCGGCGTTAATCGCTGGCTGACGGCAACCAGACGAGCAAGAAACGCATCGccc    <  2:553949/150‑1 (MQ=255)
                                                                                                                                                    tCGCTTCCAGTTGCTGCCATGCAGTAACCACCTGCTGTGCGGCCTGTTCAGAAATCAGCCCGCCAATAAACGGCTGCGGTTCGTCATCCCAGTTGCCCGGCGTTAATCGCTGGCTGACGGCAACCAGACGAGCAAGAAACGCATCGCCCt   <  2:1871066/150‑1 (MQ=255)
                                                                                                                                                     cGCTTCCAGTTGCTGCCATGCAGTAACCACCTGCTGTGCGGCCTGTTCAGAAATCAGCCCGCCAATAAACGGCTGCGGTTCGTCATCCCAGTTGCCCGGCGTTAATCGCTGGCTGACGGCAACCAGACGAGCAAGAAACGCATCGCCCTg  >  2:1668247/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
AAAGTATCATAACGCCAGACGCGCAATAACGGTCCGAACACCTCTTCATCTGGTACGCCAGCAACGCCTGTCATTTCAATGATCCCCGGCGTCAGCAACGATGTCCCTGCTTGTAATAAGCGCGGCGCAAGCAGGGGTCGTCCGCCCATCGCTTCCAGTTGCTGCCATGCAGTAACCACCTGCTGTGCGGCCTGTTCAGAAATCAGCCCGCCAATAAACGGCTGCGGTTCGTCATCCCAGTTGCCCGGCGTTAATCGCTGGCTGACGGCAACCAGACGAGCAAGAAACGCATCGCCCTG  >  CP009273/1822826‑1823124

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: