Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,058,769 C→T 88.1% I203I (ATC→ATT scpC → propionyl‑CoA:succinate CoA transferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,058,7690CT88.1% 510.3 / 54.0 177I203I (ATC→ATTscpCpropionyl‑CoA:succinate CoA transferase
Reads supporting (aligned to +/- strand):  ref base C (7/14);  new base T (66/90);  total (73/104)
Fisher's exact test for biased strand distribution p-value = 4.87e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

Only 100 of 177 total aligned reads displayed.

AGAAAGTGATCATTGAACTCAATCACTATCACGATCCGCGCGTTGCAGAACTGGCGGATATTGTGATTCCTGGCGCGCCACCGCGGCGCAATAGCGTGTCGATCTTCCATGCAATGGATCGCGTCGGTACCCGCTATGTGCAAATCGATCCGAAAAAGATTGTCGCCGTCGTGGAAACCAACTTGCCCGACGCCGGTAATATGCTGGATAAGCAAAATCCCATGTGCCAGCAGATTGCCGATAACGTGGTCACGTTCTTATTGCAGGAAATGGCGCATGGGCGTATTCCG  >  CP009273/3058624‑3058913
                                                                                                                                                 |                                                                                                                                                
agaaAGTGATCATTGAACTCAATCACTATCACGATCCGCGCGTTGCAGAACTGGCGGATATTGTGATTCCTGGCGCGCCACCGCGGCGCAATAGCGTGTCGATCTTCCATGCAATGGATCGCGTCGGTACCCGCTATGTGCAAATTGa                                                                                                                                                <  1:1352611/148‑1 (MQ=255)
 gaaAGTGATCATTGAACTCAATCACTATCACGATCCGCGCGTTGCAGAACTGGCGGATATTGTGATTCCTGGCGCGCCACCGCGGCGCAATAGCGTGTCGATCTTCCATGCAATGGATCGCGTCGGTACCCGCTATGTGCAAATTGAt                                                                                                                                               <  1:2783529/148‑1 (MQ=255)
 gaaAGTGATCATTGAACTCAATCACTATCACGATCCGCGCGTTGCAGAACTGGCGGATATTGTGATTCCTGGCGCGCCACCGCGGCGCAATAGCGTGTCGATCTTCCATGCAATGGATCGCGTCGGTACCCGCTATGTGCAAATTGAt                                                                                                                                               <  1:2860351/148‑1 (MQ=255)
 gaaAGTGATCATTGAACTCAATCACTATCACGATCCGCGCGTTGCAGAACTGGCGGATATTGTGATTCCTGGCGCGCCACCGCGGCGCAATAGCGTGTCGATCTTCCATGCAATGGATCGCGTCGGTACCCGCTATGTGCAAATTGAt                                                                                                                                               <  2:3186057/148‑1 (MQ=255)
  aaaGTGATCATTGAACTCAATCACTATCACGATCCGCGCGTTGCAGAACTGGCGGATATTGTGATTCCTGGCGCGCCACCGCGGCGCAATAGCGTGTCGATCTTCCATGCAATGGATCGCGTCGGTACCCGCTATGTGCAAATTGATc                                                                                                                                              >  1:1605372/1‑148 (MQ=255)
         tCATTGAACTCAATCACTATCACGATCCGCGCGTTGCAGAACTGGCGGATATTGTGATTTCTGGCGCGCCACCGCGGCGCAATAGCGTGTCGATCTTCCATGCAATGGATCGCGTCGGTACCCGCTATGTGCAAATCGATCCGaaaaa                                                                                                                                       <  1:1379916/148‑1 (MQ=255)
         tCATTGAACTCAATCACTATCACGATCCGCGCGTTGCAGAACTGGCGGATATTGTGATTCCTGGCGCGCCACCGCGGCGCAATAGCGTGTCGATCTTCCATGCAATGGATCGCGTCGGTACCCGCTATGTGCAAATTGATCCGaaaaa                                                                                                                                       >  1:2933472/1‑148 (MQ=255)
          cATTGAACTCAATCACTATCACGATCCGCGCGTTGCAGAACTGGCGGATATTGTGATTTCTGGCGCGCCACCGCGGCGCAATAGCGTGTCGATCTTCCATGCAATGGATCGCGTCGGTACCCGCTATGTGCAAATCGATCCGAAAAAg                                                                                                                                      <  2:2958154/148‑1 (MQ=255)
          cATTGAACTCAATCACTATCACGATCCGCGCGTTGCAGAACTGGCGGATATTGTGATTCCTGGCGCGCCACCGCGGCGCAATAGCGTGTCGATCTTCCATGCAATGGATCGCGTCGGTACCCGCTATGTGCAAATTGATCCGAAAAAg                                                                                                                                      <  2:1484064/148‑1 (MQ=255)
           aTTGAACTCAATCACTATCACGATCCGCGCGTTGCAGAACTGGCGGATATTGTGATTTCTGGCGCGCCACCGCGGCGCAATAGCGTGTCGATCTTCCATGCAATGGATCGCGTCGGTACCCGCTATGTGCAAATCGATCCGAAAAAGa                                                                                                                                     <  1:1340734/148‑1 (MQ=255)
           aTTGAACTCAATCACTATCACGATCCGCGCGTTGCAGAACTGGCGGATATTGTGATTCCTGGCGCGCCACCGCGGCGCAATAGCGTGTCGATCTTCCATGCAATGGATCGCGTCGGTACCCGCTATGTGCAAATTGATCCGAAAAAGa                                                                                                                                     <  1:2694412/148‑1 (MQ=255)
           aTTGAACTCAATCACTATCACGATCCGCGCGTTGCAGAACTGGCGGATATTGTGATTCCTGGCGCGCCACCGCGGCGCAATAGCGTGTCGATCTTCCATGCAATGGATCGCGTCGGTACCCGCTATGTGCAAATTGATCCGAAAAAGa                                                                                                                                     <  2:1778541/148‑1 (MQ=255)
            ttGAACTCAATCACTATCACGATCCGCGCGTTGCAGAACTGGCGGATATTGTGATTCCTGGCGCGCCACCGCGGCGCAATAGCGTGTCGATCTTCCATGCAATGGATCGCGTCGGTACCCGCTATGTGCAAATTGATCCGAAAAAGAt                                                                                                                                    >  2:184112/1‑148 (MQ=255)
             tGAACTCAATCACTATCACGATCCGCGCGTTGCAGAACTGGCGGATATTGTGATTTCTGGCGCGCCACCGCGGCGCAATAGCGTGTCGATCTTCCATGCAATGGATCGCGTCGGTACCCGCTATGTGCAAATCGATCCGAAAAAGAtt                                                                                                                                   <  1:1282595/148‑1 (MQ=255)
             tGAACTCAATCACTATCACGATCCGCGCGTTGCAGAACTGGCGGATATTGTGATTCCTGGCGCGCCACCGCGGCGCAATAGCGTGTCGATCTTCCATGCAATGGATCGCGTCGGTACCCGCTATGTGCAAATTGATCCGAAAAAGAtt                                                                                                                                   <  2:3294893/148‑1 (MQ=255)
             tGAACTCAATCACTATCACGATCCGCGCGTTGCAGAACTGGCGGATATTGTGATTCCTGGCGCGCCACCGCGGCGCAATAGCGTGTCGATCTTCCATGCAATGGATCGCGTCGGTACCCGCTATGTGCAAATTGATCCGAAAAAGAtt                                                                                                                                   <  2:2013907/148‑1 (MQ=255)
             tGAACTCAATCACTATCACGATCCGCGCGTTGCAGAACTGGCGGATATTGTGATTCCTGGCGCGCCACCGCGGCGCAATAGCGTGTCGATCTTCCATGCAATGGATCGCGTCGGTACCCGCTATGTGCAAATTGATCCGAAAAAGAtt                                                                                                                                   <  2:2109620/148‑1 (MQ=255)
             tGAACTCAATCACTATCACGATCCGCGCGTTGCAGAACTGGCGGATATTGTGATTCCTGGCGCGCCACCGCGGCGCAATAGCGTGTCGATCTTCCATGCAATGGATCGCGTCGGTAACCGCTATGTGCAAATTGATCCGAAAAAGAtt                                                                                                                                   <  1:119531/148‑1 (MQ=255)
               aaCTCAATCACTATCACGATCCGCGCGTTGCAGAACTGGCGGATATTGTGATTCCTGGCGCGCCACCGCGGCGCAATAGCGTGTCGATCTTCCATGCAATGGATCGCGTCGGTACCCGCTATGTGCAAATTGATCCGAAAAAGATTGt                                                                                                                                 >  1:173433/1‑148 (MQ=255)
                  tCAATCACTATCACGATCCGCGCGTTGCAGAACTGGCGGATATTGTGATTCCTGGCGCGCCACCGCGGCGCAATAGCGTGTCGATCTTCCATGCAATGGATCGCGTCGGTACCCGCTATGTGCAAATTGATCCGAAAAAGATTGTCGc                                                                                                                              >  1:1193546/1‑148 (MQ=255)
                  tCAATCACTATCACGATCCGCGCGTTGCAGAACTGGCGGATATTGTGATTCCTGGCGCGCCACCGCGGCGCAATAGCGTGTCGATCTTACATGCAATGGATCGCGTCGGTACCCGCTATGTGCAAATTGATCCGAAAAAGATTGTCGc                                                                                                                              <  1:1099285/148‑1 (MQ=255)
                    aaTCACTATCACGATCCGCGCGTTGCAGAACTGGCGGATATTGTGATTTCTGGCGCGCCACCGCGGCGCAATAGCGTGTCGATCTTCCATGCAATGGATCGCGTCGGTACCCGCTATGTGCAAATCGATCCGAAAAAGATTGTCGCcg                                                                                                                            >  1:949211/1‑148 (MQ=255)
                       cACTATCACGATCCGCGCGTTGCAGAACTGGCGGATATTGTGATTTCTGGCGCGCCACCGCGGCGCAATAGCGTGTCGATCTTCCATGCAATGGATCGCGTCGGTACCCGCTATGTGCAAATCGATCCGAAAAAGATTGTCGCcgtcg                                                                                                                         >  1:3232596/1‑148 (MQ=255)
                       cACTATCACGATCCGCGCGTTGCAGAACTGGCGGATATTGTGATTCCTGGCGCGCCACCGCGGCGCAATAGCGTGTCGATCTTCCATGCAATGGATCGCGTCGGTACCCGCTATGTGCAAATTGATCCGAAAAAGATTGTCGCcgtcg                                                                                                                         <  1:1659349/148‑1 (MQ=255)
                              aCGATCCGCGCGTTGCAGAACTGGCGGATATTGTGATTTCTGGCGCGCCACCGCGGCGCAATAGCGTGTCGATCTTCCATGCAATGGATCGCGTCGGTACCCGCTATGTGCAAATCGATCCGAAAAAGATTGTCGCCGTCGTGGAAAc                                                                                                                  <  2:2272147/148‑1 (MQ=255)
                              aCGATCCGCGCGTTGCAGAACTGGCGGATATTGTGATTCCTGGCGCGCCACCGCGGCGCAATAGCGTGTCGATCTTCCATGCAATGGATCGCGTCGGTACCCGCTATGTGCAAATTGATCCGAAAAAGATTGTCGCCGTCGTGGAAAc                                                                                                                  <  2:270863/148‑1 (MQ=255)
                               cGATCCGCGCGTTGCAGAACTGGCGGATATTGTGATTCCTGGCGCGCCACCGCGGCGCAATAGCGTGTCGATCTTCCATGCAATGGATCGCGTCGGTACCCGCTATGTGCAAATTGATCCGAAAAAGATTGTCGCCGTCGTGGAAAcc                                                                                                                 >  2:1617899/1‑148 (MQ=255)
                                gATCCGCGCGTTGCAGAACTGGCGGATATTGTGATTCCTGGCGCGCCACCGCGGCGCAATAGCGTGTCGATCTTCCATGCAATGGATCGCGTCGGTACCCGCGATGTGCAAATTGATCCGAAAAAGATTGTCGCCGTCGTGGAAACCa                                                                                                                <  2:2357455/148‑1 (MQ=255)
                                 aTCCGCGCGTTGCAGAACTGGCGGATATTGTGATTCCTGGCGCGCCACCGCGGCGCAATAGCGTGTCGATCTTCCATGCAATGGATCGCGTCGGTACCCGCTATGTGCAAATTGATCCGAAAAAGATTGTCGCCGTCGTGGAAACCaa                                                                                                               <  2:2933472/148‑1 (MQ=255)
                                 aTCCGCGCGTTGCAGAACTGGCGGATATTGTGATTCCTGGCGCGCCACCGCGGCGCAATAGCGTGTCGATCTTCCATGCAATGGATCGCGTCGGTACCCGCTATGTGCAAATTGATCCGAAAAAGATTGTCGCCGTCGTGGAAACCaa                                                                                                               <  1:571574/148‑1 (MQ=255)
                                 aTCCGCGCGTTGCAGAACTGGCGGATATTGTGATTCCTGGCGCGCCACCGCGGCGCAATAGCGTGTCGATCTTCCATGCAATGGATCGCGTCGGTACCCGCTATGTGCAAATTGATCCGAAAAAGATTGTCGCCGTCGTGGAAACCaa                                                                                                               <  1:184112/148‑1 (MQ=255)
                                  tCCGCGCGTTGCAGAACTGGCGGATATTGTGATTCCTGGCGCGCCACCGCGGCGCAATAGCGTGTCGATCTTCCATGCAATGGATCGCGTCGGTACCCGCTATGTGCAAATTGATCCGAAAAAGATTGTCGCCGTCGTGGAAACCAAc                                                                                                              >  2:682254/1‑148 (MQ=255)
                                       gcgTTGCAGAACTGGCGGATATTGTGATTCCTGGCGCGCCACCGCGGCGCAATAGCGTGTCGATCTTCCATGCAATGGATCGCGTCGGTACCCGCTATGTGCAAATTGATCCGAAAAAGATTGTCGCCGTCGTGGAAACCAACTTGcc                                                                                                         <  1:406780/148‑1 (MQ=255)
                                           tGCAGAACTGGCGGATATTGTGATTCCTGGCGCGCCACCGCGGCGCAATAGCGTGTCGATCTTCCATGCAATGGATCGCGTCGGTACCCGCTATGTGCAAATTGATCCGAAAAAGATTGTCGCCGTCGTGGAAACCAACTTGCCCGAc                                                                                                     >  1:2552647/1‑148 (MQ=255)
                                             cAGAACTGGCGGATATTGTGATTCCTGGCGCGCCACCGCGGCGCAATAGCGTGTCGATCTTCCATGCAATGGATCGCGTCGGTACCCGCTATGTGCAAATTGATCCGAAAAAGATTGTCGCCGTCGTGGAAACCAACTTGCCCGACGc                                                                                                   <  1:1817183/148‑1 (MQ=255)
                                                  cTGGCGGATATTGTGATTCCTGGCGCGCCACCGCGGCGCAATAGCGTGTCGATCTTCCATGCAATGGATCGCGTCGGTACCCGCTATGTGCAAATTGATCCGAAAAAGATTGTCGCCGTCGTGGAAACCAACTTGCCCGACGCCGGTa                                                                                              <  1:2199339/148‑1 (MQ=255)
                                                  cTGGCGGATATTGTGATTCCTGGCGCGCCACCGCGGCGCAATAGCGTGTCGATCTTCCATGCAATGGATCGCGTCGGTACCCGCTATGTGCAAATTGATCCGAAAAAGATTGTCGCCGTCGTGGAAACCAACTTGCCCGACGCCGGTa                                                                                              <  1:2949365/148‑1 (MQ=255)
                                                  cTGGCGGATATTGTGAGTCCTGGCGCGCCACCGCGGCGCAATAGCGTGTCGATCTTCCATGCAATGGATCGCGTCGGTACCCGCTATGTGCAAATTGATCCGAAAAAGATTGTCGCCGTCGTGGAAACCAACTTGCCCGACGCCGGTa                                                                                              <  1:1345549/148‑1 (MQ=255)
                                                   tGGCGGATATTGTGATTCCTGGCGCGCCACCGCGGCGCAATAGCGTGTCGATCTTCCATGCAATGGATCGCGTCGGTACCCGCTATGTGCAAATTGATCCGAAAAAGATTGTCGCCGTCGTGGAAACCAACTTGCCCGACGCCGGTaa                                                                                             <  2:1513635/148‑1 (MQ=255)
                                                   tGGCGGATATTGTGATTCCTGGCGCGCCACCGCGGCGCAATAGCGTGTCGATCTTCCATGCAATGGATCGCGTCGGTACCCGCTATGTGCAAATTGATCCGAAAAAGATTGTCGCCGTCGTGGAAACCAACTTGCCCGACGCCGGTaa                                                                                             <  2:769830/148‑1 (MQ=255)
                                                    ggCGGATATTGTGATTCCTGGCGCGCCACCGCGGCGCAATAGCGTGTCGATCTTCCATGCAATGGATCGCGTCGGTACCCGCTATGTGCAAATTGATCCGAAAAAGATTGTCGCCGTCGTGGAAACCAACTTGCCCGACGCCGGTAat                                                                                            <  2:744779/148‑1 (MQ=255)
                                                    ggCGGATATTGTGATTCCTGGCGCGCCACCGCGGCGCAATAGCGTGTCGATCTTCCATGCAATGGATCGCGTCGGTACCCGCTATGTGCAAATTGATCCGAAAAAGATTGTCGCCGTCGTGGAAACCAACTTGCCCGACGCCGGTAat                                                                                            <  2:782293/148‑1 (MQ=255)
                                                         atatTGTGATTCCTGGCGCGCCACCGCGGCGCAATAGCGTGTCGATCTTCCATGCAATGGATCGCGTCGGTACCCGCTATGTGCAAATTGATCCGAAAAAGATTGTCGCCGTCGTGGAAACCAACTTGCCCGACGCCGGTAATATGCt                                                                                       >  1:853907/1‑148 (MQ=255)
                                                         atatTGTGATTCCTGGCGCGCCACCGCGGCGCAATAGCGTGTCGATCTTCCATGCAATGGATCGCGTCGGTACCCGCTATGTGCAAATTGATCCGAAAAAGATTGTCGCCGTCGTGGAAACCAACTTGCCCGACGCCGGTAATATGCt                                                                                       >  1:1383562/1‑148 (MQ=255)
                                                          tatTGTGATTCCTGGCGCGCCACCGCGGCGCAATAGCGTGTCGATCTTCCATGCAATGGATCGCGTCGGTACCCGCTATGTGCAAATTGATCCGAAAAAGATTGTCGCCGTCGTGGAAACCAACTTGCCCGACGCCGGTAATATGCTg                                                                                      <  1:2887550/148‑1 (MQ=255)
                                                           atTGTGATTCCTGGCGCGCCACCGCGGCGCAATAGCGTGTCGATCTTCCATGCAATGGATCGCGTCGGTACCCGCTATGTGCAAATTGATCCGAAAAAGATTGTCGCCGTCGTGGAAACCAACTTGCCCGACGCCGGTAATATGCTgg                                                                                     >  2:3112712/1‑148 (MQ=255)
                                                            ttGTGATTCCTGGCGCGCCACCGCGGCGCAATAGCGTGTCGATCTTCCATGCAATGGATCGCGTCGGTACCCGCTATGTGCAAATTGATCCGAAAAAGATTGTCGCCGTCGTGGAAACCAACTTGCCCGACGCCGGTAATATGCTGGa                                                                                    <  1:1497575/148‑1 (MQ=255)
                                                            ttGTGAGTCCTGGCGCGCCACCGCGGCGCAATAGCGTGTCGATCTTCCATGCAATGGATCGCGTCGGTACCCGCTATGTGCAAATTGATCCGAAAAAGATTGTCGCCGTCGTGGAAACCAACTTGCCCGACGCCGGTAATATGCTGGa                                                                                    <  1:1497544/148‑1 (MQ=255)
                                                              gtgATTCCTGGCGCGCCACCGCGGCGCAATAGCGTGTCGATCTTCCATGCAATGGATCGCGTCGGTACCCGCTATGTGCAAATTGATCCGAAAAAGATTGTCGCCGTCGTGGAAACCAACTTGCCCGACGCCGGTAATATGCTGGATa                                                                                  <  2:1136763/148‑1 (MQ=255)
                                                                    ccTGGCGCGCCACCGCGGCGCAATAGCGTGTCGATCTTCCATGCAATGGATCGCGTCGGTACCCGCTATGTGCAAATTGATCCGAAAAAGATTGTCGCCGTCGTGGAAACCAACTTGCCCGACGCCGGTAATATGCTGGATAAGCaaa                                                                            <  1:1885304/148‑1 (MQ=255)
                                                                     cTGGCGCGCCACCGCGGCGCAATAGCGTGTCGATCTTCCATGCAATGGATCGCGTCGGTACCCGCTATGTGCAAATTGATCCGAAAAAGATTGTCGCCGTCGTGGAAACCAACTTGCCCGACGCCGGTAATATGCTGGATAAGCaaaa                                                                           <  1:2146220/148‑1 (MQ=255)
                                                                       ggCGCGCCACCGCGGCGCAATAGCGTGTCGATCTTCCATGCAATGGATCGCGTCGGTACCCGCTATGTGCAAATTGATCCGAAAAAGATTGTCGCCGTCGTGGAAACCAACTTGCCCGACGCCGGTAATATGCTGGATAAGCAAAATc                                                                         >  2:2564437/1‑148 (MQ=255)
                                                                        gcgcgcCACCGCGGCGCAATAGCGTGTCGATCTTCCATGCAATGGATCGCGTCGGTACCCGCTATGTGCAAATTGATCCGAAAAAGATTGTCGCCGTCGTGGAAACCAACTTGCCCGACGCCGGTAATATGCTGGATAAGCAAAATcc                                                                        <  2:2361732/148‑1 (MQ=255)
                                                                        gcgcgcCACCGCGGCGCAATAGCGTGTCGATCTTCCATGCAATGGATCGCGTCGGTACCCGCTATGTGCAAATTGATCCGAAAAAGATTGTCGCCGTCGTGGAAACCAACTTGCCCGACGCCGGTAATATGCTGGATAAGCAAAATcc                                                                        >  1:2841917/1‑148 (MQ=255)
                                                                        gcgcgcCACCGCGGCGCAATAGCGTGTCGATCTTCCATGCAATGGATCGCGTCGGTACCCGCTATGTGCAAATTGATCCGAAAAAGATTGTCGCCGTCGTGGAAACCAACTTGCCCGACGCCGGTAATATGCTGGATAAGCAAAATcc                                                                        >  1:1002156/1‑148 (MQ=255)
                                                                         cgcgcCACCGCGGCGCAATAGCGTGTCGATCTTCCATGCAATGGATCGCGTCGGTACCCGCTATGTGCAAATTGATCCGAAAAAGATTGTCGCCGTCGTGGAAACCAACTTGCCCGACGCCGGTAATATGCTGGATAAGCAAAATccc                                                                       <  2:534953/148‑1 (MQ=255)
                                                                         cgcgcCACCGCGGCGCAATAGCGTGTCGATCTTCCATGCAATGGATCGCGTCGGTACCCGCTATGTGCAAATTGATCCGAAAAAGATTGTCGCCGTCGTGGAAACCAACTTGCCCGACGCCGGTAATATGCTGGATAAGCAAAATccc                                                                       <  2:535692/148‑1 (MQ=255)
                                                                            gcCACCGCGGCGCAATAGCGTGTCGATCTTCCATGCAATGGATCGCGTCGGTACCCGCTATGTGCAAATTGATCCGAAAAAGATTGTCGCCGTCGTGGAAACCAACTTGCCCGACGCCGGTAATATGCTGGATAAGCAAAATCCCAtg                                                                    <  1:3113568/148‑1 (MQ=255)
                                                                               aCCGCGGCGCAATAGCGTGTCGATCTTCCATGCAATGGATCGCGTCGGTACCCGCTATGTGCAAATTGATCCGAAAAAGATTGTCGCCGTCGTGGAAACCAACTTGCCCGACGCCGGTAATATGCTGGATAAGCAAAATCCCATGTGc                                                                 <  2:2293595/148‑1 (MQ=255)
                                                                                 cgcgGCGCAATAGCGTGTCGATCTTCCATGCAATGGATCGCGTCGGTACCCGCTATGTGCAAATTGATCCGAAAAAGATTGTCGCCGTCGTGGAAACCAACTTGCCCGACGCCGGTAATATGCTGGATAAGCAAAATCCCATGTGCca                                                               >  1:1866888/1‑148 (MQ=255)
                                                                                  gcggcgCAATAGCGTGTCGATCTTCCATGCAATGGATCGCGTCGGTACCCGCTATGTGCAAATTGATCCGAAAAAGATTGTCGCCGTCGTGGAAACCAACTTGCCCGACGCCGGTAATATGCTGGATAAGCAAAATCCCATGTGCcag                                                              <  2:673302/148‑1 (MQ=255)
                                                                                      cgcAATAGCGTGTCGATCTTCCATGCAATGGATCGCGTCGGTACCCGCTATGTGCAAATTGATCCGAAAAAGATTAGCGCCGTCGTGGAAACCAACTTGCCCGACGCCGGTAATATGCTGGATAAGCAAAATCCCATGTGCCAGCAGa                                                          <  2:861307/148‑1 (MQ=255)
                                                                                       gcAATAGCGTGTCGATCTTCCATGCAATGGATCGCGTCGGTACCCGCTATGTGCAAATCGATCCGAAAAAGATTGTCGCCGTCGTGGAAACCAACTTGCCCGACGCCGGTAATATGCTGGATAAGCAAAATCCCATGTGCCAGCAGAt                                                         <  2:1856171/148‑1 (MQ=255)
                                                                                             gCGTGTCGATCTTCCATGCAATGGATCGCGTCGGTACCCGCTATGTGCAAATTGATCCGAAAAAGATTGTCGCCGTCGTGGAAACCAACTTGCCCGACGCCGGTAATATGCTGGATAAGCAAAATCCCATGTGCCAGCAGATTGCCGa                                                   >  2:3150397/1‑148 (MQ=255)
                                                                                             gCGTGTCGATCTTCCATGCAATGGATCGCGTCGGTACCCGCTATGTGCAAATTGATCCGAAAAAGATTGTCGCCGTCGTGGAAACCAACTTGCCCGACGCCGGTAATATGCTGGATAAGCAAAATCCCATGTGCCAGCAGATTGCCGa                                                   <  1:2049449/148‑1 (MQ=255)
                                                                                               gtgtCGATCTTCCATGCAATGGATCGCGTCGGTACCCGCTATGTGCAAATCGATCCGAAAAAGATTGTCGCCGTCGTGGAAACCAACTTGCCCGACGCCGGTAATATGCTGGATAAGCAAAATCCCATGTGCCAGCAGATTGCCGATa                                                 <  2:949211/148‑1 (MQ=255)
                                                                                                tgtCGATCTTCCATGCAATGGATCGCGTCGGTACCCGCTATGTGCAAATTGATCCGAAAAAGATTGTCGCCGTCGTGGAAACCAACTTGCCCGACGCCGGTAATATGCTGGATAAGCAAAATCCCATGTGCCAGCAGATTGCCGATaa                                                >  1:1030461/1‑148 (MQ=255)
                                                                                                tgtCGATCTTCCATGCAATGGATCGCGTCGGTACCCGCTATGTGCAAATTGATCCGAAAAAGATTGTCGCCGTCGTGGAAACCAACTTGCCCGACGCCGGTAATATGCTGGATAAGCAAAATCCCATGTGCCAGCAGATTGCCGATaa                                                <  2:1118579/148‑1 (MQ=255)
                                                                                                 gtCGATCTTCCATGCAATGGATCGCGTCGGTACCCGCTATGTGCAAATTGATCCGAAAAAGATTGTCGCCGTCGTGGAAACCAACTTGCCCGACGCCGGTAATATGCTGGATAAGCAAAATCCCATCTGCCAGCAGATTGCCGATAAc                                               >  2:2109052/1‑148 (MQ=255)
                                                                                                  tCGATCTTCCATGCAATGGATCGCGTCGGTACCCGCTATGTGCAAATTGATCCGAAAAAGATTGTCGCCGTCGTGGAAACCAACTTGCCCGACGCCGGTAATATGCTGGATAAGCAAAATCCCATGTGCCAGCAGATTGCCGATAACg                                              <  2:83843/148‑1 (MQ=255)
                                                                                                  tCGATCTTCCATGCAATGGATCGCGTCGGTACCCGCTATGTGCAAATTGATCCGAAAAAGATTGTCGCCGTCGTGGAAACCAACTTGCCCGACGCCGGTAATATGCTGGATAAGCAAAATCCCATGTGCCAGCAGATTGCCGATAACg                                              <  2:277709/148‑1 (MQ=255)
                                                                                                    gATCTTCCATGCAATGGATCGCGTCGGTACCCGCTATGTGCAAATCGATCCGAAAAAGATTGTCGCCGTCGTGGAAACCAACTTGCCCGACGCCGGTAATATGCTGGATAAGCAAAATCCCATGTGCCAGCAGATTGCCGATAACGTg                                            <  1:490458/148‑1 (MQ=255)
                                                                                                    gATCTTCCATGCAATGGATCGCGTCGGTACCCGCTATGTGCAAATCGATCCGAAAAAGATTGTCGCCGTCGTGGAAACCAACTTGCCCGACGCCGGTAATATGCTGGATAAGCAAAATCCCATGTGCCAGCAGATTGCCGATAACGTg                                            <  1:490859/148‑1 (MQ=255)
                                                                                                     aTCTTCCATGCAATGGATCGCGTCGGTACCCGCTATGTGCAAATTGATCCGAAAAAGATTGTCGCCGTCGTGGAAACCAACTTGCCCGACGCCGGTAATATGCTGGATAAGCAAAATCCCATGTGCCAGCAGATTGCCGATAACGTgg                                           >  2:388370/1‑148 (MQ=255)
                                                                                                      tCTTCCATGCAATGGATCGCGTCGGTACCCGCTATGTGCAAATTGATCCGAAAAAGATTGTCGCCGTCGTGGAAACCAACTTGCCCGACGCCGGTAATATGCTGGATAAGCAAAATCCCATGTGCCAGCAGATTGCCGATAACGTGGt                                          >  2:542556/1‑148 (MQ=255)
                                                                                                           cATGCAATGGATCGCGTCGGTACCCGCTATGTGCAAATTGATCCGAAAAAGATTGTCGCCGTCGTGGAAACCAACTTGCCCGACGCCGGTAATATGCTGGATAAGCAAAATCCCATGTGCCAGCAGATTGCCGATAACGTGGTCACGt                                     >  1:2765006/1‑148 (MQ=255)
                                                                                                                aaTGGATCGCGTCGGTACCCGCTATGTGCAAATTGATCCGAAAAAGATTGTCGCCGTCGTGGAAACCAACTTGCCCGACGCCGGTAATATGCTGGATAAGCAAAATCCCATGTGCCAGCAGATTGCCGATAACGTGGTCACGTTCTTa                                <  2:1092592/148‑1 (MQ=255)
                                                                                                                aaTGGATCGCGTCGGTACCCGCTATGTGCAAATTGATCCGAAAAAGATTGTCGCCGTCGTGGAAACCAACTTGCCCGACGCCGGTAATATGCTGGATAAGCAAAATCCCATGTGCCAGCAGATTGCCGATAACGTGGTCACGTTCTTa                                <  2:2507194/148‑1 (MQ=255)
                                                                                                                aaTGGATCGCGTCGGTACCCGCTATGTGCAAATTGATCCGAAAAAGATTGTCGCCGTCGTGGAAACCAACTTGCCCGACGCCGGTAATATGCTGGATAAGCAAAATCCCATGTGCCAGCAGATTGCCGATAACGTGGTCACGTTCTTa                                <  2:1749895/148‑1 (MQ=255)
                                                                                                                aaTGGATCGCGTCGGTACCCGCTATGTGCAAATTGATCCGAAAAAGATTGTCGCCGTCGTGGAAACCAACTTGCCCGACGCCGGTAATATGCTGGATAAGCAAAATCCCATGTGCCAGCAGATTGCCGATAACGTGGTCACGTTCTTa                                <  2:1956778/148‑1 (MQ=255)
                                                                                                                aaTGGATCGCGTCGGTACCCGCTATGTGCAAATTGATCCGAAAAAGATTGTCGCCGTCGTGGAAACCAACTTGCCCGACGCCGGTAATATGCTGGATAAGCAAAATCCCATGTGCCAGCAGATTGCCGATAACGTGGTCACGTTCTTa                                <  1:304661/148‑1 (MQ=255)
                                                                                                                 aTGGATCGCGTCGGTACCCGCTATGTGCAAATTGATCCGAAAAAGATTGTCGCCGTCGTGGAAACCAACTTGCCCGACGCCGGTAATATGCTGGATAAGCAAAATCCCATGTGCCAGCAGATTGCCGATAACGTGGTCACGTTCTTAt                               >  1:1811023/1‑148 (MQ=255)
                                                                                                                  tGGATCGCGTCGGTACCCGCTATGTGCAAATTGATCCGAAAAAGATTGTCGCCGTCGTGGAAACCAACTTGCCCGACGCCGGTAATATGCTGGATAAGCAAAATCCCATGTGCCAGCAGATTGCCGATAACGTGGTCACGTTCTTAtt                              >  2:301730/1‑148 (MQ=255)
                                                                                                                     aTCGCGTCGGTACCCGCTATGTGCAAATTGATCCGAAAAAGATTGTCGCCGTCGTGGAAACCAACTTGCCCGACGCCGGTAATATGCTGGATAAGCAAAATCCCATGTGCCAGCAGATTGCCGATAACGTGGTCACGTTCTTATTGCa                           <  2:2894641/148‑1 (MQ=255)
                                                                                                                      tCGCGTCGGTACCCGCTATGTGCAAATTGATCCGAAAAAGATTGTCGCCGTCGTGGAAACCAACTTGCCCGACGCCGGTAATATGCTGGATAAGCAAAATCCCATGTGCCAGCAGATTGCCGATAACGTGGTCACGTTCTTATTGCAg                          >  2:715466/1‑148 (MQ=255)
                                                                                                                           tCGGTACCCGCTATGTGCAAATTGATCCGAAAAAGATTGTCGCCGTCGTGGAAACCAACTTGCCCGACGCCGGTAATATGCTGGATAAGCAAAATCCCATGTGCCAGCAGATTGCCGATAACGTGGTCACGTTCTTATTGCAGGAAAt                     <  2:552794/148‑1 (MQ=255)
                                                                                                                             ggTACCCGCTATGTGCAAATTGATCCGAAAAAGATTGTCGCCGTCGTGGAAACCAACTTGCCCGACGCCGGTAATATGCTGGATAAGCAAAATCCCATGTGCCAGCAGATTGCCGATAACGTGGTCACGTTCTTATTGCAGGAAATgg                   <  1:1568363/148‑1 (MQ=255)
                                                                                                                               tACCCGCTATGTGCAAATTGATCCGAAAAAGATTGTCGCCGTCGTGGAAACCAACTTGCCCGACGCCGGTAATATGCTGGATAAGCAAAATCCCATGTGCCAGCAGATTGCCGATAACGTGGTCACGTTCTTATTGCAGGAAATGgcg                 <  2:3022783/148‑1 (MQ=255)
                                                                                                                                 cccGCTATGTGCAAATTGATCCGAAAAAGATTGTCGCCGTCGTGGAAACCAACTTGCCCGACGCCGGTAATATGCTGGATAAGCAAAATCCCATGTGCCAGCAGATTGCCGATAACGTGGTCACGTTCTTATTGCAGGAAATGGCGCa               >  2:134683/1‑148 (MQ=255)
                                                                                                                                  ccGCTATGTGCAAATTGATCCGAAAAAGATTGTCGCCGTCGTGGAAACCAACTTGCCCGACGCCGGTAATATGCTGGATAAGCAAAATCCCATGTGCCAGCAGATTGCCGATAACGTGGTCACGTTCTTATTGCAGGAAATGGCGCAt              >  2:1957709/1‑148 (MQ=255)
                                                                                                                                  ccGCTATGTGCAAATTGATCCGAAAAAGATTGTCGCCGTCGTGGAAACCAACTTGCCCGACGCCGGTAATATGCTGGATAAGCAAAATCCCATGTGCCAGCAGATTGCCGATAACGTGGTCACGTTCTTATTGCAGGAAATGGCGCAt              <  1:842823/148‑1 (MQ=255)
                                                                                                                                     cTATGTGCAAATTGATCCGAAAAAGATTGTCGCCGTCGTGGAAACCAACTTGCCCGACGCCGGTAATATGCTGGATAAGCAAAATCCCATGTGCCAGCAGATTGCCGATAACGTGGTCACGTTCTTATTGCAGGAAATGGCGCATggg           >  2:3246826/1‑148 (MQ=255)
                                                                                                                                     cTATGTGCAAATTGATCCGAAAAAGATTGTCGCCGTCGTGGAAACCAACTTGCCCGACGCCGGTAATATGCTGGATAAGCAAAATCCCATGTGCCAGCAGATTGCCGATAACGTGGTCACGTTCTTATTGCAGGAAATGGCGCATggg           >  1:1272600/1‑148 (MQ=255)
                                                                                                                                     cTATGTGCAAATCGATCCGAAAAAGATTGTCGCCGTCGTGGAAACCAACTTGCCCGACGCCGGTAATATGCTGGATAAGCAAAATCCCATGTGCCAGCAGATTGCCGATAACGTGGTCACGTTCTTATTGCAGGAAATGGCGCATggg           >  1:1168299/1‑148 (MQ=255)
                                                                                                                                       aTGTGCAAATTGATCCGAAAAAGATTGTCGCCGTCGTGGAAACCAACTTGCCCGACGCCGGTAATATGCTGGATAAGCAAAATCCCATGTGCCAGCAGATTGCCGATAACGTGGTCACGTTCTTATTGCAGGAAATGGCGCATGGGCg         >  2:3012137/1‑148 (MQ=255)
                                                                                                                                        tgtgCAAATTGATCCGAAAAAGATTGTCGCCGTCGTGGAAACCAACTTGCCCGACGCCGGTAATATGCTGGATAAGCAAAATCCCATGTGCCAGCAGATTGCCGATAACGTGGTCACGTTCTTATTGCAGGAAATGGCGCATGGGCGt        <  1:2371675/148‑1 (MQ=255)
                                                                                                                                           gCAAATTGATCCGAAAAAGATTGTCGCCGTCGTGGAAACCAACTTGCCCGACGCCGGTAATATGCTGGATAAGCAAAATCCCATGTGCCAGCAGATTGCCGATAACGTGGTCACGTTCTTATTGCAGGAAATGGCGCATGGGCGTAtt     <  1:3307350/148‑1 (MQ=255)
                                                                                                                                            cAAATTGATCCGAAAAAGATTGTCGCCGTCGTGGAAACCAACTTGCCCGACGCCGGTAATATGCTGGATAAGCAAAATCCCATGTGCCAGCAGATTGCCGATAACGTGGTCACGTTCTTATTGCAGGAAATGGCGCATGGGCGTATTc    >  1:2271018/1‑148 (MQ=255)
                                                                                                                                            cAAATTGATCCGAAAAAGATTGTCGCCGTCGTGGAAACCAACTTGCCCGACGCCGGTAATATGCTGGATAAGCAAAATCCCATGTGCCAGCAGATTGCCGATAACGTGGTCACGTTCTTATTGCAGGAAATGGCGCATGGGCGTATTc    <  2:1193546/148‑1 (MQ=255)
                                                                                                                                              aaTTGATCCGAAAAAGATTGTCGCCGTCGTGGAAACCAACTTGCCCGACGCCGGTAATATGCTGGATAAGCAAAATCCCATGTGCCAGCAGATTGCCGATAACGTGGTCACGTTCTTATTGCAGGAAATGGCGCATGGGCGTATTccg  >  1:674052/1‑148 (MQ=255)
                                                                                                                                                 |                                                                                                                                                
AGAAAGTGATCATTGAACTCAATCACTATCACGATCCGCGCGTTGCAGAACTGGCGGATATTGTGATTCCTGGCGCGCCACCGCGGCGCAATAGCGTGTCGATCTTCCATGCAATGGATCGCGTCGGTACCCGCTATGTGCAAATCGATCCGAAAAAGATTGTCGCCGTCGTGGAAACCAACTTGCCCGACGCCGGTAATATGCTGGATAAGCAAAATCCCATGTGCCAGCAGATTGCCGATAACGTGGTCACGTTCTTATTGCAGGAAATGGCGCATGGGCGTATTCCG  >  CP009273/3058624‑3058913

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: