Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 947,211 C→A 5.3% P600P (CCG→CCT pflB ← formate C‑acetyltransferase 1, anaerobic; pyruvate formate‑lyase 1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273947,2110CA 5.3% 1070.0 / 26.3 277P600P (CCG→CCTpflBformate C‑acetyltransferase 1, anaerobic; pyruvate formate‑lyase 1
Reads supporting (aligned to +/- strand):  ref base C (109/152);  new base A (5/11);  total (114/163)
Fisher's exact test for biased strand distribution p-value = 4.47e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 286 total aligned reads displayed.

TCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTCCGAACGGCGCGCCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGAACAGACTGAGTCGGGATAGCGTCACGGTAGGTGTGCAGTTTCTGAATTTTCTTCATGAAACGTTCTACCAGGTCAACAGCCAGGTCATCTACACGCGGATCATTGTTACCAAACTGCGGGTATTCGCCTTCGATTTCGAAGTCGATAGCCAGACCGTCTT  >  CP009273/947063‑947359
                                                                                                                                                    |                                                                                                                                                    
  agagagGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTCCGAACGGCGCGCCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGAACAGACTGAGTCgg                                                                                                                                                    <  1:1028495/149‑1 (MQ=255)
  agagagGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTCCGAACGGCGCGCCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGAACAGACTGAGTCgg                                                                                                                                                    >  1:4907436/1‑149 (MQ=255)
    agagGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTCCGAACGGCGCGCCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGAACAGACTGAGTAGGGa                                                                                                                                                  >  1:3898257/1‑149 (MQ=255)
    agagGCTACGGCACCTTTCTGGTCACGACCGTGCATCGGGATAGCACCCGGTCCGAACGGCGCGCCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACACCGTTAGAAGTGATGGTCTGAACAGACTGAGTCGGGa                                                                                                                                                  <  2:1562639/149‑1 (MQ=255)
     gagGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTCCGAACGGCGCGCCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGAACAGACTGAGTCGGGAt                                                                                                                                                 >  1:4410752/1‑149 (MQ=255)
     gagGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTCCGAACGGCGCGCCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGAACAGACTGAGTCGGGAt                                                                                                                                                 >  1:1593176/1‑149 (MQ=255)
     gagGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTCCGAACGGCGCGCCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGAACAGACTGAGTCGGGAt                                                                                                                                                 <  1:1767518/149‑1 (MQ=255)
           aCTGCACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTCCGAACGGCGCGCCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGAACAGACTGAGTCGGGATAGCGTc                                                                                                                                           <  1:2937022/149‑1 (MQ=255)
           aCTGCACCTGTCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTCCGACCGGCGCGCCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGAACAGACTGAGTCGGGATAGCGTc                                                                                                                                           <  2:1515749/149‑1 (MQ=255)
            cTGCACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTCCGAACGGCGCGCCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGAACAGACTGAGTCGGGATAGCGTCa                                                                                                                                          <  2:3808006/149‑1 (MQ=255)
             tGCACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTCCGAACGGCGCGCCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGAACAGACTGAGTCGGGATAGCGTCAc                                                                                                                                         <  1:447742/149‑1 (MQ=255)
              gCACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTCCGAACGGCGCGCCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGAACAGACTGAGTCGGGATAGCGTCACg                                                                                                                                        <  1:3578701/149‑1 (MQ=255)
                aCCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTCCGAACGGCGCGCCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGAACAGACTGAGTCGGGATAGCGTCACGGt                                                                                                                                      >  2:556388/1‑149 (MQ=255)
                  cTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTCCGAACGGCGCGCCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGAACAGACTGAGTCGGGATAGCGTCACGGTAg                                                                                                                                    <  2:479319/149‑1 (MQ=255)
                  cTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTCCGAACGGCGCGCCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGAACAGACTGAGTCGGGATAGCGTCACGGTAg                                                                                                                                    >  1:3115349/1‑149 (MQ=255)
                   tttCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTCCGAACGGCGCGCCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGAACAGACTGAGTCGGGATAGCGTCACGGTAgg                                                                                                                                   >  2:2965115/1‑149 (MQ=255)
                   tttCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTCCGAACGGCGCGCCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGAACAGACTGAGTCGGGATAGCGTCACGGTAgg                                                                                                                                   >  1:352569/1‑149 (MQ=255)
                    ttCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTCCGAACGGCGCGCCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGAACAGACTGAGTCGGGATAGCGTCACGGTAGgt                                                                                                                                  <  2:4899458/149‑1 (MQ=255)
                       tGGTCACGACCGTGCATCGGGTTAGCACCCGGTCCGAACGGCGCGCCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGAACAGACTGAGTCGGGATAGCGTCACGGTAGGtgtg                                                                                                                               <  2:3587292/149‑1 (MQ=255)
                        ggTCACGACCGTGCATCGGGTTAGCACCCGGTCCGAACGGCGCGCCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGAACAGACTGAGTCGGGATAGCGTCACGGTAGGTGTGc                                                                                                                              <  2:5129735/149‑1 (MQ=255)
                             cGACCGTGCATCGGGTTAGCACCCGGTCCGAACGGCGCGCCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGAACAGACTGAGTCGGGATAGCGTCACGGTAGGTGTGCAGttt                                                                                                                         <  2:4938753/149‑1 (MQ=255)
                              gACCGTGCATCGGGTTAGCACCCGGTCCGAACGGCGCGCCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGAACAGACTGAGTCGGGATAGCGTCACGGTAGGTGTGCAGTTTc                                                                                                                        >  2:1139670/1‑149 (MQ=255)
                              gACCGTGCATCGGGTTAGCACCCGGTCCGAACGGCGCGCCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGAACAGACTGAGTCGGGATAGCGTCACGGTAGGTGTGCAGTTTc                                                                                                                        <  1:742421/149‑1 (MQ=255)
                                  gTGCATCGGGTTAGCACCCGGTCCGAACGGCGCGCCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGAACAGACTGAGTCGGGATAGCGTCACGGTAGGTGTGCAGTTTCTGaa                                                                                                                    <  2:4310815/149‑1 (MQ=255)
                                  gTGCATCGGGTTAGCACCCGGTCCGAACGGCGCGCCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGAACAGACTGAGTCGGGATAGCGTCACGGTAGGTGTGCAGTTTCTGaa                                                                                                                    >  2:3942269/1‑149 (MQ=255)
                                  gTGCATCGGGTTAGCACCCGGTCCGAACGGCGCGCCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGAACAGACTGAGTCGGGATAGCGTCACGGTAGGTGTGCAGTTTCTGaa                                                                                                                    >  2:3264941/1‑149 (MQ=255)
                                  gTGCATCGGGTTAGCACCCGGTCCGAACGGCGCGCCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTGAGAAGTGATGGTCAGAACAGACTGAGTCGGGATAGCGTCACGGTAGGTGTGCAGTTTCTGaa                                                                                                                    >  2:2446258/1‑149 (MQ=255)
                                   tGCATCGGGTTAGCACCCGGTCCGAACGGCGCGCCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGAACAGACTGAGTCGGGATAGCGTCACGGTAGGTGTGCAGTTTCTGAAt                                                                                                                   <  2:1978378/149‑1 (MQ=255)
                                     cATCGGGTTAGCACCCGGTCCGAACGGCGCGCCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGAACAGACTGAGTCGGGATAGCGTCACGGTAGGTGTGCAGTTTCTGAAttt                                                                                                                 <  2:3289523/149‑1 (MQ=255)
                                     cATCGGGTTAGCACCCGGTCCGAACGGCGCGCCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGAACAGACTGAGTCGGGATAGCGTCACGGTAGGTGTGCAGTTTCTGAAttt                                                                                                                 <  2:30649/149‑1 (MQ=255)
                                     cATCGGGTTAGCACCCGGTCCGAACGGCGCGCCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGAACAGACTGAGTCGGGATAGCGTCACGGTAGGTGTGCAGTTTCTGAAttt                                                                                                                 >  1:1447231/1‑149 (MQ=255)
                                       tCGGGTTAGCACCCGGTCCGAACGGCGCGCCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGAACAGACTGAGTCGGGATAGCGTCACGGTAGGTGTGCAGTTTCTGAATTttc                                                                                                               <  2:338050/149‑1 (MQ=255)
                                        cGGGTTAGCACCCGGTCCGAACGGCGCGCCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGAACAGACTGAGTCGGGATAGCGTCACGGTAGGTGTGCAGTTTCTGAATTttct                                                                                                              <  1:1362066/149‑1 (MQ=255)
                                          ggTTAGCACCCGGTCCGAACGGCGCGCCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGAACAGACTGAGTCGGGATAGCGTCACGGTAGGTGTGCAGTTTCTGAATTttcttc                                                                                                            <  1:3939335/149‑1 (MQ=255)
                                            ttAGCACCCGGTCCGAACGGCGCGCCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGAACAGACTGAGTCGGGATAGCGTCACGGTAGGTGTGCAGTTTCTGAATTTTCTTCAt                                                                                                          <  2:2833164/149‑1 (MQ=255)
                                             tAGCACCCGGTCCGAACGGCGCGCCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGAACAGACTGAGTCGGGATAGCGTCACGGTAGGTGTGCAGTTTCTGAATTTTCTTCATg                                                                                                         <  1:3518792/149‑1 (MQ=255)
                                               gCACCCGGTCCGAACGGCGCGCCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGAACAGACTGAGTCGGGATAGCGTCACGGTAGGTGTGCAGTTTCTGAATTTTCTTCATGaa                                                                                                       >  1:3623383/1‑149 (MQ=255)
                                                   ccGGTCCGAACGGCGCGCCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGAACAGACTGAGTCGGGATAGCGTCACGGTAGGTGTGCAGTTTCTGAATTTTCTTCATGAAACGt                                                                                                   <  2:3615229/149‑1 (MQ=255)
                                                      gTCCGAACGGCGCGCCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGAACAGACTGAGTCGGGATAGCGTCACGGTAGGTGTGCAGTTTCTGAATTTTCTTCATGAAACGctct                                                                                                >  2:1230221/1‑149 (MQ=255)
                                                        ccGAACGGCGCGCCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGAACAGACTGAGTCGGGATAGCGTCACGGTAGGTGTGCAGTTTCTGAATTTTCTTCATGAAACGTTCTAc                                                                                              <  1:1008254/149‑1 (MQ=255)
                                                          gAACGGCGCGCCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGAACAGACTGAGTAGGGATAGCGTCACGGTAGGTGTGCAGTTTCTGAATTTTCTTCATGAAACGTTCTACCa                                                                                            >  1:4785759/1‑149 (MQ=255)
                                                           aaCGGCGCGCCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGAACAGACTGAGTCGGGATAGCGTCACGGTAGGTGTGCAGTTTCTGAATTTTCTTCATGAAACGTTCTACCAg                                                                                           >  2:1075835/1‑149 (MQ=255)
                                                             cGGCGCGCCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGAACAGACTGAGTCGGGATAGCGTCACGGTAGGTGTGCAGTTTCTGAATTTTCTTCATGAAACGTTCTACCAGGt                                                                                         <  1:1295251/149‑1 (MQ=255)
                                                                    ccAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGAACAGACTGAGTCGGGATAGCGTCACGGTAGGTGTGCAGTTTCTGAATTTTCTTCATGAAACGTTCTACCAGGTCAACAGc                                                                                  >  2:1263675/1‑149 (MQ=255)
                                                                       gCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGTCCGAACAGACTGAGTCGGGCTAGCGTCACGGTAGGTGTGCAGTGTCTGAATTTTCTTCATGAAACGTTCTACCAGGTCAACAGCCAg                                                                               >  1:4106260/1‑149 (MQ=255)
                                                                       gCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGAACAGACTGAGTCGGGATAGCGTCACGGTAGGTGTGCAGTTTCTGAATTTTCTTCATGAAACGTTCTACCAGGTCAACAGCCAg                                                                               <  2:1188814/149‑1 (MQ=255)
                                                                        cACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGAACAGACTGAGTCGGGATAGCGTCACGGTAGGTGTGCAGTTTCTGAATTTTCTTCATGAAACGTTCTACCAGGTCAACAGCCAgg                                                                              <  1:4286153/149‑1 (MQ=255)
                                                                              gacCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGAACAGACTGAGTCGGGATAGCGTCACGGTAGGTGTGCAGTTTCTGAATTTTCTTCATGAAACGTTCTACCAGGTCAACAGCCAGGTCATCt                                                                        <  2:1889288/149‑1 (MQ=255)
                                                                              gacCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGAACAGACTGAGTAGGGATAGCGTCACGGTAGGTGTGCAGTTTCTGAATTTTCTTCATGAAACGTTCTACCAGGTCAACAGCCAGGTCATCt                                                                        <  1:5248766/149‑1 (MQ=255)
                                                                                ccGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGAACAGACTGAGTAGGGATAGCGTCACGGTAGGTGTGCAGTTTCTGAATTTTCTTCATGAAACGTTCTACCAGGTCAACAGCCAGGTCATCTac                                                                      <  1:3227627/149‑1 (MQ=255)
                                                                                 cGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGAACAGACTGAGTCGGGATAGCGTCACGGTAGGTGTGCAGTTTCTGAATTTTCTTCATGAAACGTTCTACCAGGTCAACAGCCAGGTCATCTaca                                                                     <  2:3979002/149‑1 (MQ=255)
                                                                                   tCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGAACAGACTGAGTCGGGATAGCGTCACGGTAGGTGTGCAGTTTCTGAATTTTCTTCATGAAACGTTCTACCAGGTCAACAGCCAGGTCATCTACAcg                                                                   >  2:5005810/1‑149 (MQ=255)
                                                                                    cTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGAACAGACTGAGTCGGGATAGCGTCACGGTAGGTGTGCAGTTTCTGAATTTTCTTCATGAAACGTTCTACCAGGTCAACAGCCAGGTCATCTACAcgc                                                                  <  2:1674363/149‑1 (MQ=255)
                                                                                     tGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGAACAGACTGAGTCGGGATAGCGTCACGGTAGGTGTGCAGTTTCTGAATTTTCTTCATGAAACGTTCTACCAGGTCAACAGCCAGGTCATCTACAcgcg                                                                 <  1:3270148/149‑1 (MQ=255)
                                                                                        ggTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGAACAGACTGAGTCGGGATAGCGTCACGGTAGGTGTGCAGTTTCTGAATTTTCTTCATGAAACGTTCTACCAGGTCAACAGCCAGGTCATCTACACGCGGAt                                                              <  1:1779033/149‑1 (MQ=255)
                                                                                         gtgtTACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGAACAGACTGAGTCGGGATAGCGTCACGGTAGGTGTGCAGTTTCTGAATTTTCTTCATGAAACGTTCTACCAGGTCAACAGCCAGGTCATCTACACGCGGATc                                                             <  2:1358008/149‑1 (MQ=255)
                                                                                         gtgtTACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGAACAGACTGAGTCGGGATAGCGTCACGGTAGGTGTGCAGTTTCTGAATTTTCTTCATGAAACGTTCTACCAGGTCAACAGCCAGGTCATCTACACGCGGATc                                                             >  2:2094837/1‑149 (MQ=255)
                                                                                         gtgtTACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGAACAGACTGAGTAGGGATAGCGTCACGGTAGGTGTGCAGTTTCTGAATTTTCTTCATGAAACGTTCTACCAGGTCAACAGCCAGGTCATCTACACGCGGATc                                                             <  1:2361319/149‑1 (MQ=255)
                                                                                           gtTACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGAACAGACTGAGTCGGGATAGCGTCACGGTAGGTGTGCAGTTTCTGAATTTTCTTCATGAAACGTTCTACCAGGTCAACAGCCAGGTCATCTACACGCGGATCAt                                                           <  1:4989115/149‑1 (MQ=255)
                                                                                           gtTACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGAACAGACTGAGTCGGGATAGCGTCACGGTAGGTGTGCAGTTTCTGAATTTTCTTCATGAAACGTTCTACCAGGTCAACAGCCAGGTCATCTACACGCGGATCAt                                                           <  1:4146248/149‑1 (MQ=255)
                                                                                               cccGTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGAACAGACTGAGTAGGGATAGCGTCACGGTAGGTGTGCAGTTTCTGAATTTTCTTCATGAAACGTTCTACCAGGTCAACAGCCAGGTCATCTACACGCGGATCATTGtt                                                       >  2:4179969/1‑149 (MQ=255)
                                                                                                ccGTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGAACAGACTGAGTCGGGATAGCGTCACGGTAGGTGTGCAGTTTCTGAATTTTCTTCATGAAACGTTCTACCAGGTCAACAGCCAGGTCATCTACACGCGGATCATTGTTa                                                      >  2:5091515/1‑149 (MQ=255)
                                                                                                ccGTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGAACAGACTGAGTCGGGATAGCGTCACGGTAGGTGTGCAGTTTCTGAATTTTCTTCATGAAACGTTCTACCAGGTCAACAGCCAGGTCATCTACACGCGGATCATTGTTa                                                      >  1:2683389/1‑149 (MQ=255)
                                                                                                 cGTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGAACAGACTGAGTCGGGATAGCGTCACGGTAGGTGTGCAGTTTCTGAATTTTCTTCATGAAACGTTCTACCAGGTCAACAGCCAGGTCATCTACACGCGGATCATTGTTAc                                                     <  1:4051144/149‑1 (MQ=255)
                                                                                                  gTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGAACAGACTGAGTCGGGATAGCGTCACGGTAGGTGTGCAGTTTCTGAATTTTCTTCATGAAACGTTCTACCAGGTCAACAGCCAGGTCATCTACACGCGGATCATTGTTAcc                                                    >  2:1805715/1‑149 (MQ=255)
                                                                                                  gTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGAACAGACTGAGTCGGGATAGCGTCACGGTAGGTGTGCAGTTTCTGAATTTTCTTCATGAAACGTTCTACCAGGTCAACAGCCAGGTCATCTACACGCGGATCATTGTTAcc                                                    <  2:4259760/149‑1 (MQ=255)
                                                                                                      tCTTACCATACACACCGTTAGAAGTGATGGTCAGAACAGACTGAGTCGGGATAGCGTCACGGTAGGTGTGCAGTTTCTGAATGTTCTTCATGAAACGTTCTACCAGGTCAACAGCCAGGTCATCTACACGCGGATCATTGTTACCAAAc                                                <  2:4120847/149‑1 (MQ=255)
                                                                                                       cTTACCATACACAACGTTAGAAGTGATGGTCAGAACAGACTGAGTCGGGATAGCGTCACGGTAGGTGTGCAGTTTCTGAATTTTCTTCATGAAACGTTCTACCAGGTCAACAGCCAGGTCATCTACACGCGGATCATTGTTACCAAACt                                               >  1:2149720/1‑149 (MQ=255)
                                                                                                         tACCATACACAACGTTAGAAGTGATGGTCAGAACAGACTGAGTCGGGATAGCGTCACGGTAGGTGTGCAGTTTCTGAATTTTCTTCATGAAACGTTCTACCAGGTCAACAGCCAGGTCATCTACACGCGGATCATTGTTACCAAACTGc                                             >  2:270941/1‑149 (MQ=255)
                                                                                                             aTACACAACGTTAGAAGTGATGGTCAGAACAGACTGAGTCGGGATAGCGTCACGGTAGGTGTGCAGTTTCTGAATTTTCTTCATGAAACGTTCTACCAGGTCAACAGCCAGGTCATCTACACGCGGATCATTGTTACCAAACTGCGGGt                                         <  2:3061188/149‑1 (MQ=255)
                                                                                                             aTACACAACGTTAGAAGTGATGGTCAGAACAGACTGAGTCGGGATAGCGTCACGGTAGGTGTGCAGTTTCTGAATTTTCTTCATGAAACGTTCTACCAGGTCAACAGCCAGGTCATCTACACGCGGATCATTGTTACCAAACTGCGGGt                                         <  2:3228656/149‑1 (MQ=255)
                                                                                                              tACACAACGTTAGAAGTGATGGTCAGAACAGACTGAGTCGGGATAGCGTCACGGTAGGTGTGCAGTTTCTGAATTTTCTTCATGAAACGTTCTACCAGGTCAACAGCCAGGTCATCTACACGCGGATCATTGTTACCAAACTGCGGGTa                                        >  1:400952/1‑149 (MQ=255)
                                                                                                              tACACAACGTTAGAAGTGATGGTCAGAACAGACTGAGTCGGGATAGCGTCACGGTAGGTGTGCAGTTTCTGAATTTTCTTCATGAAACGTTCTACCAGGTCAACAGCCAGGTCATCTACACGCGGATCATTGTTACCAAACTGCGGGTa                                        >  2:2285784/1‑149 (MQ=255)
                                                                                                               acacAACGTTAGAAGTGATGGTCAGAACAGACTGAGTCGGGATAGCGTCACGGTAGGTGTGCAGTTTCTGAATTTTCTTCATGAAACGTTCTACCAGGTCAACAGCCAGGTCATCTACACGCGGATCATTGTTACCAAACTGCGGGTAt                                       >  2:4047077/1‑149 (MQ=255)
                                                                                                                 acaACGTTAGAAGTGATGGTCAGAACAGACTGAGTCGGGATAGCGTCACGGTAGGTGTGCAGTTTCTGAATTTTCTTCATGAAACGTTCTACCAGGTCAACAGCCAGGTCATCTACACGCGGATCATTGTTACCAAACTGCGGGTATTc                                     >  1:1721768/1‑149 (MQ=255)
                                                                                                                  caACGTTAGAAGTGATGGTCAGAACAGACTGAGTCGGGATAGCGTCACGGTAGGTGTGCAGTTTCTGAATTTTCTTCATGAAACGTTCTACCAGGTCAACAGCCAGGTCATCTACACGCGGATCATTGTTACCAAACTGCGGGTATTCg                                    >  2:2331842/1‑149 (MQ=255)
                                                                                                                    aCGTTAGAAGTGATGGTCAGAACAGACTGAGTCGGGATAGCGTCACGGTAGGTGTGCAGTTTCTGAATTTTCTTCATGAAACGTTCTACCAGGTCAACAGCCAGGTCATCTACACGCGGATCATTGTTACCAAACTGCGGGTATTCGcc                                  <  1:2965122/149‑1 (MQ=255)
                                                                                                                    aCGTTAGAAGTGATGGTCAGAACAGACTGAGTCGGGATAGCGTCACGGTAGGTGTGCAGTTTCTGAATTTTCTTCATGAAACGTTCTACCAGGTCAACAGCCAGGTCATCTACACGCGGATCATTGTTACCAAACTGCGGGTATTCGcc                                  <  1:1483358/149‑1 (MQ=255)
                                                                                                                     cGTTAGAAGTGATGGTCAGAACAGACTGAGTCGGGATAGCGTCACGGTAGGTGTGCAGTTTCTGAATTTTCTTCATGAAACGTTCTACCAGGTCAACAGCCAGGTCATCTACACGCGGATCATTGTTACCAAACTGCGGGTATTCGCCt                                 >  1:3840900/1‑149 (MQ=255)
                                                                                                                     cGTTAGAAGTGATGGTCAGAACAGACTGAGTCGGGATAGCGTCACGGTAGGTGTGCAGTTTCTGAATTTTCTTCATGAAACGTTCTACCAGGTCAACAGCCAGGTCATCTACACGCGGATCATTGTTACCAAACTGCGGGTATTCGCCt                                 <  1:2082377/149‑1 (MQ=255)
                                                                                                                      gTTAGAAGTGATGGTCAGAACAGACTGAGTCGGGATAGCGTCACGGTAGGTGTGCAGTTTCTGAATTTTCTTCATGAAACGTTCTACCAGGTCAACAGCCAGGTCATCTACACGCGGATCATTGTTACCAAACTGCGGGTATTCGCCtt                                <  1:3691714/149‑1 (MQ=255)
                                                                                                                       ttAGAAGTGATGGTCAGAACAGACTGAGTCGGGATAGCGTCACGGTAGGTGTGCAGTTTCTGAATTTTCTTCATGAAACGTTCTACCAGGTCAACAGCCAGGTCATCTACACGCGGATCATTGTTACCAAACTGCGGGTATTCGCCTTc                               >  1:1907850/1‑149 (MQ=255)
                                                                                                                          gAAGTGATGGTCAGAACAGACTGAGTCGGGATAGCGTCACGGTAGGTGTGCAGTTTCTGAATTTTCTTCATGAAACGTTCTACCAGGTCAACAGCCAGGTCATCTACACGCGGATCATTGTTACCAAACTGCGGGTATTCGCCTTCGAt                            <  2:5223212/149‑1 (MQ=255)
                                                                                                                                  ggTCAGAACAGACTGAGTCGGGATAGCGTCACGGTAGGTGTGCAGTTTCTGAATTTTCTTCATGAAACGTTCTACCAGGTCAACAGCCAGGTCATCTACACGCGGATCATTGTTACCAAACTGCGGGTATTCGCCTTCGATTTCGAAGt                    <  1:2486765/149‑1 (MQ=255)
                                                                                                                                      aGAACAGACTGAGTCGGGATAGCGTCACGGTAGGTGTGCAGTTTCTGAATTTTCTTCATGAAACGTTCTACCAGGTCAACAGCCAGGTCATCTACACGCGGATCATTGTTACCAAACTGCGGGTATTCGCCTTCGATTTCGAAGTCGAt                <  2:139305/149‑1 (MQ=255)
                                                                                                                                       gAACAGACTGAGTCGGGATAGCGTCACGGTAGGTGTGCAGTTTCTGAATTTTCTTCATGAAACGTTCTACCAGGTCAACAGCCAGGTCATCTACACGCGGATCATTGTTACCAAACTGCGGGTATTCGCCTTCGATTTCGAAGTCGATa               >  2:3482936/1‑149 (MQ=255)
                                                                                                                                       gAACAGACTGAGTCGGGATAGCGTCACGGTAGGTGTGCAGTTTCTGAATTTTCTTCATGAAACGTTCTACCAGGTCAACAGCCAGGTCATCTACACGCGGATCATTGTTACCAAACTGCGGGTATTCGCCTTCGATTTCGAAGTCGATa               >  2:255042/1‑149 (MQ=255)
                                                                                                                                        aaCAGACTGAGTAGGGATAGCGTCACGGTAGGTGTGCAGTTTCTGAATTTTCTTCATGAAACGTTCTACCAGGTCAACAGCCAGGTCATCTACACGCGGATCATTGTTACCAAACTGCGGGTATTCGCCTTCGATTTCGAAGTCGATAg              <  2:1982606/149‑1 (MQ=255)
                                                                                                                                         aCAGACTGAGTCGGGATAGCGTCACGGTAGGTGTGCAGTTTCTGAATTTTCTTCATGAAACGTTCTACCAGGTCAACAGCCAGGTCATCTACACGCGGATCATTGTTACCAAACTGCGGGTATTCGCCTTCGATTTCGAAGTCGATAGc             >  2:5172199/1‑149 (MQ=255)
                                                                                                                                         aCAGACTGAGTCGGGATAGCGTCACGGTAGGTGTGCAGTTTCTGAATTTTCTTCATGAAACGTTCTACCAGGTCAACAGCCAGGTCATCTACACGCGGATCATTGTTACCAAACTGCGGGTATTCGCCTTCGATTTCGAAGTCGATAGc             >  2:2993295/1‑149 (MQ=255)
                                                                                                                                          cAGACTGAGTCGGGATAGCGTCACGGTAGGTGTGCAGTTTCTGAATTTTCTTCATGAAACGTTCTACCAGGTCAACAGCCAGGTCATCTACACGCGGATCATTGTTACCAAACTGCGGGTATTCGCCTTCGATTTCGAAGTCGATAGcc            <  2:2189704/149‑1 (MQ=255)
                                                                                                                                           aGACTGAGTCGGGATAGCGTCACGGTAGGTGTGCAGTTTCTGAATTTTCTTCATGAAACGTTCTACCAGGTCAACAGCCAGGTCATCTACACGCGGATCATTGTTACCAAACTGCGGGTATTCGCCTTCGATTTCGAAGTCGATAGCCa           >  2:4530961/1‑149 (MQ=255)
                                                                                                                                             aCTGAGTCGGGATAGCGTCACGGTAGGTGTGCAGTTTCTGAATTTTCTTCATGAAACGTTCTACCAGGTCAACAGCCAGGTCATCTACACGCGGATCATTGTTACCAAACTGCGGGTATTCGCCTTCGATTTCGAAGTCGATAGCCAGa         <  1:1958294/149‑1 (MQ=255)
                                                                                                                                               tGAGTCGGGATAGCGTCACGGTAGGTGTGCAGTTTCTGAATTTTCTTCATGAAACGTTCTACCAGGTCAACAGCCAGGTCATCTACACGCGGATCATTGTTACCAAACTGCGGGTATTCGCCTTCGATTTCGAAGTCGATAGCCAGAcc       <  1:2880144/149‑1 (MQ=255)
                                                                                                                                               tGAGTCGGGATAGCGTCACGGTAGGTGTGCAGTTTCTGAATTTTCTTCATGAAACGTTCTACCAGGTCAACAGCCAGGTCATCTACACGCGGATCATTGTTACCAAACTGCGGGTATTCGCCTTCGATTTCGAAGTCGATAGCCAGAcc       >  2:2941895/1‑149 (MQ=255)
                                                                                                                                               tGAGTCGGGATAGCGTCACGGTAGGTGTGCAGTTTCTGAATTTTCTTCATGAAACGTTCTACCAGGTCAACAGCCAGGTCATCTACACGCGGATCATTGTTACCAAACTGCGGGTATTCGCCTTCGATTTCGAAGTCGATAGCCAGAcc       >  1:132718/1‑149 (MQ=255)
                                                                                                                                               tGAGTAGGGATAGCGTCACGGTAGGTGTGCAGTTTCTGAATTTTCTTCATGAAACGTTCTACCAGGTCAACAGCCAGGTCATCTACACGCGGATCATTGTTACCAAACTGCGGGTATTCGCCTTCGATTTCGAAGTCGATAGCCAGAcc       <  1:2218378/149‑1 (MQ=255)
                                                                                                                                                gAGTCGGGATAGCGTCACGGTAGGTGTGCAGTTTCTGAATTTTCTTCATGAAACGTTCTACCAGGTCAACAGCCAGGTCATCTACACGCGGATCATTGTTACCAAACTGCGGGTATTCGCCTTCGATTTCGAAGTCGATAGCCAGACCg      <  2:3442248/149‑1 (MQ=255)
                                                                                                                                                    cGGGATAGCGTCACGGTAGGTGTGCAGTTTCTGAATTTTCTTCATGAAACGTTCTACCGGGTCAACAGCCAGGTCATCTACACGCGGATCATTGTTACCAAACTGCGGGTATTCGCCTTCGATTTCGAAGTCGATAGCCAGACCGTCtt  <  1:1527155/149‑1 (MQ=255)
                                                                                                                                                    cGGGATAGCGTCACGGTAGGTGTGCAGTTTCTGAATTTTCTTCATGAAACGTTCTACCAGGTCAACAGCCAGGTCATCTACACGCGGATCATTGTTACCAAACTGCGGGTATTCGCCTTCGATTTCGAAGTCGATAGCCAGACCGTCtt  <  1:1128970/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
TCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTCCGAACGGCGCGCCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGAACAGACTGAGTCGGGATAGCGTCACGGTAGGTGTGCAGTTTCTGAATTTTCTTCATGAAACGTTCTACCAGGTCAACAGCCAGGTCATCTACACGCGGATCATTGTTACCAAACTGCGGGTATTCGCCTTCGATTTCGAAGTCGATAGCCAGACCGTCTT  >  CP009273/947063‑947359

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: