Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 37,692 G→T 8.7% Q45K (CAG→AAG)  caiC ← putative crotonobetaine/carnitine‑CoA ligase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP00927337,6920GT 8.7% 951.0 / 23.4 307Q45K (CAG→AAG) caiCputative crotonobetaine/carnitine‑CoA ligase
Reads supporting (aligned to +/- strand):  ref base G (164/102);  new base T (2/38);  total (166/141)
Fisher's exact test for biased strand distribution p-value = 1.65e-12
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 311 total aligned reads displayed.

GGCGTTAATCGGCACCATAATCGCGCCAATTTTTGCCAGCCCGAACCAGCAAAAGATAAATTCCGGGCAGTTGTCGAGATGTAGTGCAACCTTGTCGCCTTTGCGAATCCCCAGCGTATAAAACAGGTTTGCCGTGCGGTTAATCTCCTGATTTAACTCAAGATAACTATACCGGTTAACGACTCCGCCGCTGGATTCACAAATCAGCGCCGTTTTATGACCGTAAACGTCCGCAAGATCGTCCCACATTTGACGTAGATGTTGTCCGCCAATGATATCCATTGCACCTCTATCCATTT  >  CP009273/37543‑37841
                                                                                                                                                       |                                                                                                                                                     
ggcgTTAATCGGCACCATAATCGCGCCAATTTTTGCCAGCCCGAACCAGCAAAAGATAAATTCCGGGCAGTTGTCGAGATGTAGTGCAACCTTGTCGCCTTTGCGAATCCCCAGCGTATAAAACAGGTTTGCCGTGCGGTTAATCTCCTg                                                                                                                                                       <  2:2869482/150‑1 (MQ=255)
   gTTAATCGGCACCATAATCGCGCCAATTTTTGCCAGCCCGAACCAGCAAAAGATAAATTCCGGGCAGTTGTCGAGATGTAGTGCAACCTTGTCGCCTTTGCGAATCCCCAGCGTATAAAACAGGTTTGCCGTGCGGTTAATCTCCTGAtt                                                                                                                                                    >  1:3732835/1‑150 (MQ=255)
     tAATCGGCACCATAATCGCGCCAATTTTTGCCAGCCCGAACCAGCAAAAGATAAATTCCGGGCAGTTGTCGAGATGTAGTGCAACCTTGTCGCCTTTGCGAATCCCCAGCGTATAAAACAGGTTTGCCGTGCGGTTAATCTCCTGATTTa                                                                                                                                                  >  1:2745211/1‑150 (MQ=255)
         cGGCACCATAATCGCGCCAATTTTTGCCAGCCCGAACCAGCAAAAGATAAATTCCGGGCAGTTGTCGAGATGTAGTGCAACCTTGTCGCCTTTGCGAATCCCCAGCGTATAAAACAGGTTTGCCGTGCGGTTAATCTCCTGATTTAACTc                                                                                                                                              >  2:3389400/1‑150 (MQ=255)
         cGGCACCATAATCGCGCCAATTTTTGCCAGCCCGAACCAGCAAAAGATAAATTCCGGGCAGTTGTCGAGATGTAGTGCAACCTTGTCGCCTTTGCGAATCCCCAGCGTATAAAACAGGTTTGCCGTGCGGTTAATCTCCTGATTGAACTc                                                                                                                                              >  2:3091109/1‑150 (MQ=255)
          ggCACCATAATCGCGCCAATTTTTGCCAGCCCGAACCAGCAAAAGATAAATTCCGGGCAGTTGTCGAGATGTAGTGCAACCTTGTCGCCTTTGCGAATCCCCAGCGTATAAAACAGGTTTGCCGTGCGGTTAATCTCCTGATTTAACTCa                                                                                                                                             <  2:4362466/150‑1 (MQ=255)
            cACCATAATCGCGCCAATTTTTGCCAGCCCGAACCAGCAAAAGATAAATTCCGGGCAGTTGTCGAGATGTAGTGCAACCTTGTCGCCTTTGCGAATCCCCAGCGTATAAAACAGGTTTGCCGTGCGGTTAATCTCCTGATTTAACTCAAg                                                                                                                                           <  1:2928124/150‑1 (MQ=255)
            cACCATAATCGCGCCAATTTTTGCCAGCCCGAACCAGCAAAAGATAAATTCCGGGCAGTTGTCGAGATGTAGTGCAACCTTGTCGCCTTTGCGAATCCCCAGCGTATAAAACAGGTTTGCCGTGCGGTTAATCTCCTGATTTAACTCAAg                                                                                                                                           >  1:1814088/1‑150 (MQ=255)
             aCCATAATCGCGCCAATTTTTGCCAGCCCGAACCAGCAAAAGATAAATTCCGGGCAGTTGTCGAGATGTAGTGCAACCTTGTCGCCTTTGCGAATCCCCAGCGTATAAAACAGGTTTGCCGTGCGGTTAATCTCCTGATTTAACTCAAGa                                                                                                                                          >  2:3003148/1‑150 (MQ=255)
              ccATAATCGCGCCAATTTTTGCCAGCCCGAACCAGCAAAAGATAAATTCCGGGCAGTTGTCGAGATGTAGTGCAACCTTGTCGCCTTTGCGAATCCCCAGCGTATAAAACAGGTTTGCCGTGCGGTTAATCTCCTGATTTAACTCAAGAt                                                                                                                                         >  2:2601700/1‑150 (MQ=255)
               cATAATCGCGCCAATTTTTGCCAGCCCGAACCAGCAAAAGATAAATTCCGGGCAGTTGTCGAGATGTAGTGCAACCTTGTCGCCTTTGCGAATCCCCAGCGTATAAAACAGGTTTGCCGTGCGGTTAATCTCCTGATTTAACTCAAGATa                                                                                                                                        <  1:3124192/150‑1 (MQ=255)
                 tAATCGCGCCAATTTTTGCCAGCCCGAACCAGCAAAAGATAAATTCCGGGCAGTTGTCGAGATGTAGTGCAACCTTGTCGCCTTTGCGAATCCCCAGCGTATAAAACAGGTTTGCCGTGCGGTTAATCTCCTGATTTAACTCAAGATAAc                                                                                                                                      <  1:2144961/150‑1 (MQ=255)
                  aaTCGCGCCAATTTTTGCCAGCCCGAACCAGCAAAAGATAAATTCCGGGCAGTTGTCGAGATGTAGTGCAACCTTGTCGCCTTTGCGAATCCCCAGCGTATAAAACAGGTTTGCCGTGCGGTTAATCTCCTGATTTAACTCAAGATAACt                                                                                                                                     <  1:703016/150‑1 (MQ=255)
                    tCGCGCCAATTTTTGCCAGCCCGAACCAGCAAAAGATAAATTCCGGGCAGTTGTCGAGATGTAGTGCAACCTTGTCGCCTTTGCGAATCCCCAGCGTATAAAACAGGTTTGCCGTGCGGTTAATCTCCTGATTTAACTCAAGATAACtat                                                                                                                                   <  2:639249/150‑1 (MQ=255)
                     cgcgCCAATTTTTGCCAGCCCGAACCAGCAAAAGATAAATTCCGGGCAGTTGTCGAGATGTAGTGCAACCTTGTCGCCTTTGCGAATCCCCAGCGTATAAAACAGGTTTGCCGTGCGGTTAATCTCCTGATTTAACTCAAGATAACtata                                                                                                                                  <  1:355914/150‑1 (MQ=255)
                     cgcgCCAATTTTTGCCAGCCCGAACCAGCAAAAGATAAATTCCGGGCAGTTGTCGAGATGTAGTGCAACCTTGTCGCCTTTGCGAATCCCCAGCGTATAAAACAGGTTTGCCGTGCGGTTAATCTCCTGATTTAACTCAAGATAACtata                                                                                                                                  <  1:3676695/150‑1 (MQ=255)
                        gcCAATTTTTGCCAGCCCGAACCAGCAAAAGATAAATTCCGGGCAGTTGTCGAGATGTAGTGCAACCTTGTCGCCTTTGCGAATCCCCAGCGTATAAAACAGGTTTGCCGTGCGGTTAATCTCCTGATTTAACTCAAGATAACTATACCg                                                                                                                               >  1:2417233/1‑150 (MQ=255)
                        gcCAATTTTTGCCAGCCCGAACCAGCAAAAGATAAATTCCGGGCAGTTGTCGAGATGTAGTGCAACCTTGTCGCCTTTGCGAATCCCCAGCGTATAAAACAGGTTTGCCGTGCGGTTAATCTCCTGATTTAACTCAAGATAACTATACCg                                                                                                                               >  2:866492/1‑150 (MQ=255)
                        gcCAATTTTTGCCAGCCCGAACCAGCAAAAGATAAAGTCCGGGCAGTTGTCGAGATGTAGTGCAACCTTGTCGCCTTTGCGAATCCCCAGCGTATAAAACAGGTTTGCCGTGCGGTTAATCTCCTGATTTAACTCAAGATAACTATACCg                                                                                                                               <  2:505958/150‑1 (MQ=255)
                          cAATTTTTGCCAGCCCGAACCAGCAAAAGATAAATTCCGGGCAGTTGTCGAGATGTAGTGCAACCTTGTCGCCTTTGCGAATCCCCAGCGTATAAAACAGGTTTGCCGTGCGGTTAATCTCCTGATTTAACTCAAGATAACTATACCGGt                                                                                                                             >  2:2421779/1‑150 (MQ=255)
                           aaTTTTTGCCAGCCCGAACCAGCAAAAGATAAATTCCGGGCAGTTGTCGAGATGTAGTGCAACCTTGTCGCCTTTGCGAATCCCCAGCGTATAAAACAGGTTTGCCGTGCGGTTAATCTCCTGATTTAACTCAAGATCACTCTACCGGtt                                                                                                                            >  2:3203076/1‑150 (MQ=255)
                            aTTTTTGCCAGCCCGAACCAGCAAAAGATAAATTCCGGGCAGTTGTCGAGATGTAGTGCAACCTTGTCGCCTTTGCGAATCCCCAGCGTATAAAACAGGTTTGCCGTGCGGTTAATCTCCTGATTTAACTCAAGATAACTATACCGGTTa                                                                                                                           >  2:1761406/1‑150 (MQ=255)
                            aTTTTTGCCAGCCCGAACCAGCAAAAGATAAATTCCGGGCAGTTGTCGAGATGTAGTGCAACCTTGTCGCCTTTGCGAATCCCCAGCGTATAAAACAGGTTTGCCGTGCGGTTAATCTCCTGATTTAACTCAAGATAACTATACCGGTTa                                                                                                                           >  2:2899079/1‑150 (MQ=255)
                                ttGCCAGCCCGAACCAGCAAAAGATAAATTCCGGGCAGTTGTCGAGATGTAGTGCAACCTTGTCGCCTTTGCGAATCCCCAGCGTATAAAACAGGTTTGCCGTGCGGTTAATCTCCTGATTTAACTCAAGATAACTATACCGGTTAACGa                                                                                                                       >  1:1012999/1‑150 (MQ=255)
                                  gCCAGCCCGAACCAGCAAAAGATAAATTCCGGGCAGTTGTCGAGATGTAGTGCAACCTTGTCGCCTTTGCGAATCCCCAGCGTATAAAACAGGTTTGCCGTGCGGTTAATCTCCTGATTTAACTCAAGATAACTATACCGGTTAACGACt                                                                                                                     >  1:1306583/1‑150 (MQ=255)
                                    cAGCCCGAACCAGCAAAAGATAAATTCCGGGCAGTTGTCGAGATGTAGTGCAACCTTGTCGCCTTTGCGAATCCCCAGCGTATAAAACAGGTTTGCCGTGCGGTTAATCTCCTGATTTAACTCAAGATAACTATACCGGTTAACGACTcc                                                                                                                   >  2:4258162/1‑150 (MQ=255)
                                    cAGCCCGAACCAGCAAAAGATAAATTCCGGGCAGTTGTCGAGATGTAGTGCAACCTTGTCGCCTTTGCGAATCCCCAGCGTATAAAACAGGTTTGCCGTGCGGTTAATCTCCTGATTTAACTCAAGATAACTATACCGGTTAACGACTcc                                                                                                                   >  2:2589190/1‑150 (MQ=255)
                                    cAGCCCGAACCAGCAAAAGATAAATTCCGGGCAGTTGTCGAGATGTAGTGCAACCTTGTCGCCTTTGCGAATCCCCAGCGTATAAAACAGGTTTGCCGTGCGGTTAATCTCCTGATTTAACTCAAGATAACTATACCGGTTAACGACTcc                                                                                                                   >  2:1500258/1‑150 (MQ=255)
                                          gAACCAGCAAAAGATAAATTCCGGGCAGTTGTCGAGATGTAGTGCAACCTTGTCGCCTTTGCGAATCCCCAGCGTATAAAACAGGTTTGCCGTGCGGTTAATCTCCTGATTTAACTCAAGATAACTATACCGGTTAACGACTCCGCCGCt                                                                                                             >  2:2122702/1‑150 (MQ=255)
                                             ccAGCAAAAGATAAATTCCGGGCAGTTGTCGAGATGTAGTGCAACCTTGTCGCCTTTGCGAATCCCCAGCGTATAAAACAGGTTTGCCGTGCGGTTAATCTCCTGATTTAACTCAAGATAACTATACCGGTTAACGACTCCGCCGCTGGa                                                                                                          >  1:1343685/1‑150 (MQ=255)
                                             ccAGAAAAAGATAAATTCCGGGCTGTTGTCGAGATGTAGTGCAACCTTGTCGCCTTTGCGAATCCCCAGCTTATAAAACTGGTTTGCCGTGCGGTTAATCTCCTGATTTAACTCAAGATAACTATACCGGTTAACGACTCCGCCGCTGGa                                                                                                          <  2:1098203/150‑1 (MQ=255)
                                              cAGCAAATGATAAATGCCGGGCAGTTGTCGAGATGTAGTGAAACCTTGTCGCCTTTGCGAATCCCCAGCGTATAAAACTGGTTTGCCGTGCGGTTAATCTCCTGATTTAACTCAATATAACTATACCGGTTAACGACTCCGCCGCTGGAt                                                                                                         <  2:4306907/150‑1 (MQ=255)
                                               agaaaaaGATAAATTCCGGGCAGTTGGCGAGATGTAGTGCAACCTTGTCGCCTTTGCGAAACCCCAGCGTATAAAATAGGTTTGCCGTGCGGTTAATCTCCTGATTTAACTCAAGATAACTATACCGGTTAACGACTCCGCCGCTGGAtt                                                                                                        <  2:3270715/147‑1 (MQ=255)
                                               aGCAAAAGATAAATTCCGGGCAGTTGTCGAGATGTAGTGCAACCTTGTCGCCTTTGCGAATCCCCAGCGTATAAAACAGGTTTGCCGTGCGGTTAATCTCCTGATTTAACTCAAGATAACTATACCGGTTAACGACTCCGCCGCTGGAtt                                                                                                        >  2:3411322/1‑150 (MQ=255)
                                                gaaaaaTATAAATTCCGGGCAGTGGTCGATTTGTTGTGAAACCTTGTCGCCTTTGCTAATCCCCCGCTTATAAAACTGGTTTGCCGTGGGTTTAATCTCATTATTTAACTCAATATAAATATACCGGTAAACGACTCCGCCGCTGGATTc                                                                                                       <  1:4004165/148‑1 (MQ=255)
                                                  aaaaGATAAATTCCGGGCAGTTGTCGAGATGTAGTGCAACCTTGTCGCCTTTGCGAATCCCCAGCGTATAAAACAGGTTTGCCGTGCGGTTAATCTCCTGATTTAACTCAAGATAACTATACCGGTTAACGACTCCGCCGCTGGATTcac                                                                                                     >  1:84055/1‑150 (MQ=255)
                                                  aaaaGATAAATTCCGGGCAGTTGTCGAGATGTAGTGCAACCTTGTCGCCTTTGCGAATCCCCAGCGTATAAAACAGGTTTGCCGTGCGGTTAATCTCCTGATTTAACTCAAGATAACTACACCGGTTAACGACTCCGCCGCTGGATTcac                                                                                                     >  1:1720663/1‑150 (MQ=255)
                                                       ttAAATTCCGGGCAGTTGTCGAGATGTAGTGCAACCTTGTCGCCTTTGCGAATCCCCAGCGTATAAAACAGGTTTGCCGTGCGGTTAATCTCCTGATTTAACTCAAGATAACTATACCGGTTAACGACTCCGCCGCTGGATTCACAAATc                                                                                                <  2:2413676/149‑1 (MQ=255)
                                                       aTAAATTCCGGGCTGTTGTCGAGATGTAGTGCAACCTTGTCGCCTTTGCGAAGCCCACGCGTATAAAACAGGTTTGCCGTGCGTTTAATCTCCTTATTTAACTCAAGATAAATATACCGGTTAACGACTCCGCCGCTGGATTCACAAATc                                                                                                <  2:3599581/150‑1 (MQ=255)
                                                       aTAAATTCCGGGCAGTTGTCGAGATGTAGTGCAACCTTGTCGCCTTTGCGAATCCCCAGCGTATAAAACAGGTTTGCCGTGCGGTTAATCTCCTGATTTAACTCAAGATAACTATACCGGTTAACGACTCCGCCGCTGGATTCACAAATc                                                                                                >  1:1988766/1‑150 (MQ=255)
                                                           agtCCGGGCAGTTGTCGAGATGTAGTGCAACCTTGTCGCCTTTGCGAATCCCCAGCGTATAAAACAGGTTTGCCGTGCGGTTAATCTCCTGATTTAACTCAAGATAACTATACCGGTTAACGACTCCGCCGCTGGATTCACAAATCAgcg                                                                                            <  2:2215309/148‑1 (MQ=255)
                                                           aTTCCGGGCAGTTGTCGAGATGTAGTGCAACCTTGTCGCCTTTGCGAATCCCCCGCGTATAAAACAGGTTTGCCGTGCGGTTAATCTCCTGATTTAACTCAAGATAACTATACCGGTTAACGACTCCGCCGCTGGATTCACAAATCAgcg                                                                                            <  2:1363121/150‑1 (MQ=255)
                                                           aTTCCGGGCAGTTGTCGAGATGTAGTGCAACCTTGTCGCCTTTGCGAATCCCCAGCGTATAAAACAGGTTTGCCGTGCGGTTAATCTCCTGATTTAACTCAAGATAACTATACCGGTTAACGACTCCGCCGCTGGATTCACAAATCAgcg                                                                                            >  1:4370631/1‑150 (MQ=255)
                                                             tCCGGGCAGTTGTCGTGATGTAGTGCAACCTTGTCGCCTTTGCGAATCCCCAGCGTATAAAACAGGTTTGCCGTGCGGTTAATCTCCTGATTTAACTCAAGATAACTATACCGGTTAACGACTCCGCCGCTGGATTCACAAATCAGCGcc                                                                                          <  2:368998/150‑1 (MQ=255)
                                                                       tGTCGAGATGTAGTGCAACCTTGTCGCCTTTGCGAATCCCCAGCGTATAAAACAGGTTTGCCGTGCGGTTAATCTCCTGATTTAACTCAAGATAACTATACCGGTTAACGACTCCGCCGCTGGATTCACAAATCAGCGCCGTTTTATGAc                                                                                >  1:3821410/1‑150 (MQ=255)
                                                                       tGTCGAGATGTAGTGCAACCTTGTCGCCTTTGCGAATCCCCAGCGTATAAAACAGGTTTGCCGTGCGGTTAATCTCCTGATTTAACTCAAGATAACTATACCGGTTAACGACTCCGCCGCTGGATTCACAAATCAGCGCCGTTTTATGAc                                                                                >  1:4528769/1‑150 (MQ=255)
                                                                       tGTCGAGATGTAGTGCAACCTTGTCGCCTTTGCGAATCCCCAGCGTATAAAACAGGTTTGCCGTGCGGTTAATCTCCTGATTTAACTCAAGATAACTATACCGGTTAACGACTCCGCCGCTGGATTCACAAATCAGCGCCGTTTTATGAc                                                                                <  2:1646280/150‑1 (MQ=255)
                                                                        gTCGAGATGTAGTGCAACCTTGTCGCCTTTGCGAATCCCCAGCGTATAAAACAGGTTTGCCGTGCGGTTAATCTCCTGATTTAACTCAAGATAACTATCCCGGTTAACGACTCCGCCGCTGGATTCACAAATCAGCGCCGTTTTATGAcc                                                                               <  2:1037146/150‑1 (MQ=255)
                                                                           gtgaTGTATTGAAACCTTTTCGATTTTGCGAATCCCCAGCTTATAAAACAGTTTTGCCGTGCGGTTAATATCGTTTTTTATCTATATATTACTATACCGGTTAACGACTCCGCCGCTGGATTCACAAATCAGCGCCGTTTTATGACCGTa                                                                            <  2:1794018/148‑1 (MQ=255)
                                                                           gagaTGTAGTGCAACCTTGTCGCCTTTGCGAATCCCCAGCTTATAAAACTGGTTTGCCGTGCGGTTAATCTCCTTATTTAACTCAATATAACTATACCGGTTAACGACTCCGCCGCTGGATTCACAAATCAGCGCCGTTTTATGACCGTa                                                                            <  2:2286328/150‑1 (MQ=255)
                                                                              aTGTAGTGCAACCTTGTCGCCTTTGCGAATCCCCAGCGTATAAAACAGGTTTGCCGTGCGGTTAATCTCCTGATTTAACTCAAGATAACTATACCGGTTAACGACTCCGCCGCTGGATTCACAAATCAGCGCCGTTTTATGACCGTAAAc                                                                         >  1:2099276/1‑150 (MQ=255)
                                                                                      cAACCTTGTCGCCTTTGCGAATCCCCAGCGTATAAAACAGGTTTGCCGTGCGGTTAATCTCCTGATTTAACTCAAGATAACTATACCGGTTAACGACTCCGCCGCTGGATTCACAAATCAGCGCCGTTTTATGACCGTAAACGTCCGCaa                                                                 >  2:589009/1‑150 (MQ=255)
                                                                                       aaCCTTGTCGCCTTTGCGAATCCCCAGCGTATAAAACAGGTTTGCCGTGCGGTTAATCTCCTGATTGAACTCACGATAACTGTCCCGGTTAACGACTCCGCCGCTGGCTTCACACATCAGCGCCGTTTTCTGTCCGTACACGTCCGCatg                                                                >  2:1916045/1‑148 (MQ=255)
                                                                                        agcTTGTGGCTTTTGCTAACCCCCTGCGTATAAAACTGGTTTGCGTTGGGGTTAATCTCCTTTTTTAACTAATTATAATTATCCGGTTTAACGACTCCGCCGCTGGATTCACAAATCAGCGCCGTTTTATGACCGTAAACGTCCGCAAGa                                                               <  1:2549320/148‑1 (MQ=255)
                                                                                            tGTCGCCTTTGCGAATCCCCAGCGTATAAAACAGGTTTGCCGTGCGGTTAATCTCCTGATTTAACTCAAGATAACTATACCGGTTAACGACTCCGCCGCTGGATTCACAAATCAGCGCCGTTTTATGACCGTAAACGTCCGCAAGATCGt                                                           >  2:4108670/1‑150 (MQ=255)
                                                                                             gTCGCTTTTGCCAACCCCCTGCTAAAAAACCTGTTTGGCCGTGGGTTAATTCACCTTATTTAAATAAATATAACTATACGGGTTAACGACTCCGCCGCTGGATTCACAAATCAGCGCCTTTTTATGACCGTAAACGTCCGCAAGATCGTc                                                          <  1:645872/150‑1 (MQ=255)
                                                                                               cGCTTTTGCGAACCCCAAGCTTAAAAAACAGTTTTGCGGTGCGGTTAATCTCTTTTTTTAACTCAATATACATATCCCGTTTCACGACTCCGCCGCTGGATTCACAAATCAGCGCCGTTTTATGACCGTAAACGTCCGCAAGATCGTccc                                                        <  2:1179826/150‑1 (MQ=255)
                                                                                                gCCTTTGCGAATCCCCAGCGTATAAAACAGGTTTGCCGTGCGGTTAATCTCCTGATTTAACTCAAGATAACTATACCGGTTAACGACTCCGCCGCTGGATTCACAAATCAGCGCCGTTTTATGACCGTAAACGTCCGCAAGATCGTCCca                                                       >  2:184158/1‑150 (MQ=255)
                                                                                                 ccTTTGCGAATCCCCAGCGTATAAAACAGGTTTGCCGTGCGGTTAATCTCCTGATTTAACTCAAGATAACTATACCGGTTAACGACTCCGCCGCTGGATTCACAAATCAGCGCCGTTTTATGACCGTCAACGTCCGCAAGATAGTCCcac                                                      >  1:3273593/1‑150 (MQ=255)
                                                                                                  cTTTGCGAATCCCCAGCGTATAAAACCGGTTTGCCGTGCGGTTAATCTCCTGATTTAACTCAAGATAAATATCCCGGTTAACGACTCCGCCGCTGGATTCACAAATCAGCGCCGTTTTATGACCGTAAACGTCCGCAAGATCGTCCcaca                                                     <  2:4372237/150‑1 (MQ=255)
                                                                                                  cTTTGCGAATCCCCAGCGTATAAAACAGGTTTGCCGTGCGGTTAATCTCCTGATTTAACTCAAGATAACTATACCGGTTAACGACTCCGCCGCTGGATTCACAAATCAGCGCCGTTTTATGACCGTAAACGTCCGCAAGATCGTCCcaca                                                     >  2:1272870/1‑150 (MQ=255)
                                                                                                   tttGCGAATCCCCAGCGTATAAAACAGGTTTGCCGTGCGGTTAATCTCCTGATTTAACTCAAGATAACTATACCGGTTAACGACTCCGCCGCTGGATTCACAAATCGGCGCCGTTTTATGACCGTAAACGTCCGCCAGATCGTCCCACAt                                                    >  1:991577/1‑150 (MQ=255)
                                                                                                     tGCGAATCCCCAGCGTATAAAACAGGTTTGCCGTGCGGTTAATCTCCTGATTTAACTCAAGATAACTATACCGGTTAACGACTCCGCCGCTGGATTCACAAATCAGCGCCGTTTTATGACCGTAAACGTCCGCAAGATCGTCCCACAttt                                                  >  2:3711094/1‑150 (MQ=255)
                                                                                                         aaTCCCCAGCGTATAAAACAGGTTTGCCGTGCGGTTAATCTCCTGATTTAACTCAAGATAACTATACCGGTTAACGACTCCGCCGCTGGATTCACAAATCAGCGCCGTGTTATGACCGTAAACGTCCGCAAGATCGTCCCACAGTTGACg                                              >  1:188050/1‑150 (MQ=255)
                                                                                                            ccccGGGGTATAAACCTGTTTGGCGGGGGGTTTAATGTCGTTTTTTAACTCATTATAACTATCCGGGTTAACGACTCCGCCGCTGGATTCACAAATCAGCGCCGTTTTATGACCGTAAACGTCCGCAAGATCGTCCCACATTTGACGTAg                                           <  1:4296283/150‑1 (MQ=255)
                                                                                                               cAGCGTATAAAACAGGTTTGCCGTGCGGTTAATCTCCTGATTTAACTCAAGATAACTATACCGGTTAACGACTCCGCCGCTGGATTCACAAATCAGCGCCGTTTTATGACCGTAAACGTCCGCAAGATCGTCCCACATTTGACGTAGAtg                                        >  1:3512168/1‑150 (MQ=255)
                                                                                                                 gcttaaaaaaCTGTTTTGCGGTGGGTTTAAGCACCTGATTTAAATAAATATAAATATACCGTTAAACGACTCCGCGGCTGGATTCACAAATCAGCGCCGTTTTATGACCGTAAACGTCCGCAAGATCGTCCCACATTTGACGTAGAtgtt                                      <  1:1626068/144‑1 (MQ=255)
                                                                                                                  cGTATAAAACAGGTTTGCCGTGCGGTTAATCTCCTGATTTAACTCAAGATAACTATACCGGTTAACGACTCCGCCGCTGGATTCACAAATCAGCGCCGTTTTATGACCGTAAACGTCCGCAAGATCGTCCCACATTTGACGTCGAtgttg                                     >  2:1817526/1‑150 (MQ=255)
                                                                                                                   gTATAAAACAGGTTTGCCGTGCGGTTAATCTCCTGATTTAACTCAAGATAACTATACCGGTTAACGACTCCGCCGCTGGATTCACAAATCAGCGCCGTTTTATGACCGTAAACGTCCGCAAGATCGTCCCACATTTGACGTAGAtgttgt                                    >  2:2831877/1‑150 (MQ=255)
                                                                                                                    tataAAACAGGTTTGCCGTGCGGTTAATCTCCTGATTTAACTCAAGATAACTATACCGGTTAACGACTCCGCCGCTGGATTCACAAATCAGCGCCGTTTTATGACCGTAAACGTCCGCAAGATCGTCCCACATTTGACGTAGATGTTGTc                                   >  1:916804/1‑150 (MQ=255)
                                                                                                                      taAAACTGGTTTGCGGTGCGGTTAATCTCATGATTTAACTCAAGATAACTATACGGTTTAACGACTCCGCCGCTGGATTCACAAATCAGCGCCGTTTTATGACCGTAAACGTCCGCAAGATCGTCCCACATTTGACGTAGATGTTGTCCg                                 <  1:3003148/150‑1 (MQ=255)
                                                                                                                      taAAACAGGTTTGCCGTGCGGTTAATCTCCTGATTTAACTCAAGATAACTATACCGGTTAACGACTCCGCCGCTGGATTCACAAATCAGCGCCGTTTTATGACCGTAAACGTCCGCAAGATCGTCCCACATTTGACGTAGATGTTGTCCg                                 >  2:4141950/1‑150 (MQ=255)
                                                                                                                       aaaaCAGGTTTGCCGTGCGGTTAATCTCCTGATTTAACTCAAGATAACTATACCGGTTAACGACTCCGCCGCTGGATTCACAAATCAGCGCCGTTTTATGACCGTAAACGTCCGCAAGATCGTCCCACATTTGACGTAGATGTTGTCCGc                                >  1:2809805/1‑150 (MQ=255)
                                                                                                                           ctggTTTGCCGTGCGGTTAATCTCCTGATTTAACTCAAGATAACTATACCGGTTAACGACTCCGCCGCTGGATTCACAAATCAGCGCCGTTTTATGACCGTAAACGTCCGCAAGATCGTCCCACATTTGACGTAGATGTTGTCCGCCAAt                            <  2:3685037/148‑1 (MQ=255)
                                                                                                                            tgttttGCCGTGGGTTAAATCTCTTTATTTAACTCATTATAAATATCCCGGTAAACGACTCCGCCGCTGGATTCACAAATCTGCGCCGTTTTATGACCGTAAACGTCCGCAAGATCGTCCCACATTTGACGTAGATGTTGTCCGCCAATg                           <  2:6749/147‑1 (MQ=255)
                                                                                                                              gTTTGCCGTGCGGTTAATCTCCTGATTTAACTCAAGATAACTATACCGGTTAACGACTCCGCCGCTGGATTCACAAATCAGCGCCGTTTTATGACCGTAAACGTCCGCAAGATCGTCCCACATTTGACGTAGATGTTGTCCGCCAATGat                         <  2:2907638/150‑1 (MQ=255)
                                                                                                                                ttGCCGTGCGGTTAATCTCCTGATTTAACTCAAGATAACTATACCGGTTAACGACTCCGCCGCTGGATTCACAAATCAGCGCCGTTTTATGACCGTAAACGTCCGCAAGATCGTCCCACATTTGACGTAGATGTTGTCCGCCAATGatat                       >  2:1220046/1‑150 (MQ=255)
                                                                                                                                   gCCGTGCGTTTAATCTCCTGTTTTAACTCAAGATAACTCTACCGGTTAACGACTCCGCCGCTGGATTCACAAATCAGCGCCGTTTTATGACCGTAAACGTCCGCAAGATCGTCCCACATTTGACGTAGATGTTGTCCGCCAATGATATcc                     <  2:2276341/150‑1 (MQ=255)
                                                                                                                                     cGTGCGGTTAATCTCCTGATTTAACTCAAGATAACTATACCGGTTAACGACTCCGCCGCTGGATTCACAAATCAGCGCCGTTTTATGACCGTAAACGTCCGCAAGATCGTCCCACATTTGACGTAGATGTTGTCCGCCAATGATATCCag                   >  2:3430052/1‑149 (MQ=255)
                                                                                                                                      ttGCGTTTAATCTCCTTATTTAACTAATGATAAATATACGGGTTAACGACTCCGCCGCTGGATTCACAAATCAGCGCCGTTTTATGACCGTAAACGTCCGCAAGATCGTCCCACATTTGACGTAGATGTTGTCCGCCAATGATATCCAtt                  <  1:1761406/149‑1 (MQ=255)
                                                                                                                                      gTGCGGTTATTCTCTTTTTTTAACTAATGATAACTATACGGGTTCACGACTCCGCCGGTGGATTCACAAATCAGCGCCGTTTTATGACCGTAAACGTCCGCAAGATCGTCCCACATTTGACGTAGATGTTGTCCGCCAATGATATCCAtt                  <  1:2339498/150‑1 (MQ=255)
                                                                                                                                      gTGCGGTTAATCTCCTTATTTAACTCAAGATAACTATACCGGTTAACGACTCCGCCGCTGGATTCACAAATCAGCGCCGTTTTATGACCGTAAACGTCCGCAAGATCGTCCCACATTTGACGTAGATGTTGTCCGCCAATGATATCCAtt                  <  2:2363131/150‑1 (MQ=255)
                                                                                                                                      gTGCGGTTAATCTCCTGATTTAACTCAAGATAACTATACCGGTTAACGACTCCGCCGCTGGATTCACAAATCAGCGCCGTTTTATGACCGTAAACGTCCGCAAGATCGTCCCACATTTGACGTAGATGTTGTCCGCCAATGATATCCAtt                  >  1:1674502/1‑150 (MQ=255)
                                                                                                                                       tGCGGTTAATCTCCTGATTTAACTAAATATAACTATACCGGTAAACGACTCCGCCGCTGGATTCACAAATCAGCGCCGTTTTATGACCGTAAACGTCCGCAAGATCGTCCCACATTTGACGTAGATGTTGTCCGCCAATGATATCCATTg                 <  2:1988094/150‑1 (MQ=255)
                                                                                                                                       ggCGTTTAATCTCCTTATTTAACTCAAGATAAATATCCCGGTTAACGACTCCGCCGCTGGATTCACAAATCAGCGCCGTTTTATGACCGTAAACGTCCGCAAGATCGTCCCACATTTGACGTAGATGTTGTCCGCCAATGATATCCATTg                 <  2:103526/149‑1 (MQ=255)
                                                                                                                                         cGGTTAATCTCCTGATTTAACTCAAGATAACTATACCGGTTAACGACTCCGCCGCTGGATTCACAAATCAGCGCCGTTTTATGACCGTAAACGTCCGCAAGATCGTCCCACATTTGACGTAGATGTTGTCCGCCAATGATATCCATTGCa               >  2:1807018/1‑150 (MQ=255)
                                                                                                                                           gTTAATCTCCTGATTTAACTCAAGATAACTATACCGGTTAACGACTCCGCCGCTGGATTCACAAATCAGCGCCGTTTTATGACCGTAAACGTCCGCAAGATCGTCCCACATTTGACGTAGATGTTGTCCGCCAATGATATCCATTGCAcc             >  2:3970742/1‑150 (MQ=255)
                                                                                                                                           gTTAATCTCCTGATTTAACTCAAGATAACTATACCGGTTAACGACTCCGCCGCTGGATTCACAAATCAGCGCCGTTTTATGACCGTAAACGTCCGCAAGATCGTCCCACATTTGACGTAGATGTTGTCCGCCAATGATATCCATTGCAcc             >  2:4216056/1‑150 (MQ=255)
                                                                                                                                           gTTAATCTCCTGATTTAACTCAAGATAACTATACCGGTTAACGACTCCGCCGCTGGATTCACAAATCAGCGCCGTTTTATGACCGTAAACGTCCGCAAGATCGTCCCACATTTGACGTAGATGTTGTCCGCCAATGATATCCATTGCAcc             <  2:1590306/150‑1 (MQ=255)
                                                                                                                                           gTTAATCTCCTGATTTAACTCAAGATAACTATACCGGTTAACGACTCCGCCGCTGGATTCACAAATCAGCGCCGTTTTATGACCGTAAACGTCCGCAAGATCGTCCCACATTTGACGTAGATGTTGTCCGCCAATGATATCCATTGCAcc             >  2:2644370/1‑150 (MQ=255)
                                                                                                                                           gTTAATCTCCTGATTTAACTCAAGATAACTATACCGGTTAACGACTCCGCCGCTGGATTCACAAATCAGCGCCGTTTTATGACCGTAAACGTCCGCAAGATCGTCCCACATTTGACGTAGATGTTGTCCGCCAATGATATCCATTGCAcc             >  2:3477000/1‑150 (MQ=255)
                                                                                                                                                 tctcCTGATTTAACTCAAGATAACTATACCGGTTAACGACTCCGCCGCTGGATTCACAAATCAGCGCCGTTTTATGACCGTAAACGTCCGCAAGATCGTCCCACATTTGACGTAGATGTTTTCCGCCAATTATATCCAATGCACCTCTAt        >  1:1584097/1‑150 (MQ=255)
                                                                                                                                                 tctcCTGATTTAACTCAAGATAACTATACCGGTTAACGACTCCGCCGCTGGATTCACAAATCAGCGCCGTTTTATGACCGTAAACGTCCGCAAGATCGTCCCACATTTGACGTAGATGTTGTCCGCCAATGATATCCATTGCACCTCTAt        <  2:4477875/150‑1 (MQ=255)
                                                                                                                                                 tctcCTGATTTAACTCAAGATAACTATACCGGTTAACGACTCCGCCGCTGGATTCACAAATCAGCGCCGTTTTATGACCGTAAACGTCCGCAAGATCGTCCCACATTTGACGTAGATGTTGTCCGCCAATGATATCCATTGCACCTCTAt        >  2:2365784/1‑150 (MQ=255)
                                                                                                                                                  ctcCTGATTTAACTCAAGATAACTATACCGGTTAACGACTCCGCCGCTGGATTCACAAATCAGCGCCGTTTTATTACCGTAAACGTCCGCAAAATCGTCCCACATTTGACGTAGATGTTTTACGCCAATGATATCCATTGCACCACTAta       >  1:1887910/1‑149 (MQ=255)
                                                                                                                                                   tcCTGATTTAACTCAATATAACTATACCGTTTAACGACTCCGCCGCTGGATTCACAAATCAGCGCCGTTTTATGACCGTAAACGTCCGCAAGATCGTCCCACATTTGACGTAGATGTTGTCCGCCAATGATATCCATTGCACCTCTATcc      <  2:3325118/150‑1 (MQ=255)
                                                                                                                                                    ccTGATTTAACTCAAGATAACTATACCGGTTAACGACTCCGCCGCTGGATTCACAAATCAGCGCCGTTTTATGACCGTAAACGTCCGCAAGATCGTCCCACATTTGACGTAGATGTTGTCCGCCAATGATATCCATTGCACCTCTATCCa     >  2:2341555/1‑150 (MQ=255)
                                                                                                                                                     cTGATTTAACTCAAGATAACTATACCGGTTAACGACTCCGCCGCTGGATTCACAAATCAGCGCCGTTTTATGACCGTAAACGTCCGCAAGATCGTCCCACATTTGACGTAGATGTTGTCCGCCAATGATATCCATTGCACCTCTATCCAt    >  2:2944564/1‑150 (MQ=255)
                                                                                                                                                     cTGATTTAACTCAAGATAACTATACCGGTTAACGACTCCGCCGCTGGATTCACAAATCAGCGCCGTTTTATGACCGTAAACGTCCGCAAGATCGTCCCACATTTGACGTAGATGTTGTCCGCCAATGATATCCATTGCACCTCTATCCAt    >  2:760654/1‑150 (MQ=255)
                                                                                                                                                       gATTTAACTCAAGATAACTATACCGGTTAACGACTCCGCCGCTGGATTCACAAATCAGCGCCGTTTTATGACCGTAAACGTCCGCAAGATCGTCCCACATTTGACGTAGATGTTGTCCGCCAATGATATCCATTTCACCTCTATCCAttt  >  1:222959/1‑150 (MQ=255)
                                                                                                                                                       |                                                                                                                                                     
GGCGTTAATCGGCACCATAATCGCGCCAATTTTTGCCAGCCCGAACCAGCAAAAGATAAATTCCGGGCAGTTGTCGAGATGTAGTGCAACCTTGTCGCCTTTGCGAATCCCCAGCGTATAAAACAGGTTTGCCGTGCGGTTAATCTCCTGATTTAACTCAAGATAACTATACCGGTTAACGACTCCGCCGCTGGATTCACAAATCAGCGCCGTTTTATGACCGTAAACGTCCGCAAGATCGTCCCACATTTGACGTAGATGTTGTCCGCCAATGATATCCATTGCACCTCTATCCATTT  >  CP009273/37543‑37841

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: