Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,245,537 G→T 100% N286K (AAC→AAA plsB ← glycerol‑3‑phosphate O‑acyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,245,5370GT98.7% 761.1 / 0.3 239N286K (AAC→AAAplsBglycerol‑3‑phosphate O‑acyltransferase
Reads supporting (aligned to +/- strand):  ref base G (2/1);  new base T (100/136);  total (102/137)
Fisher's exact test for biased strand distribution p-value = 5.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.20e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

Only 100 of 242 total aligned reads displayed.

AGGTTGATCCCGGCGGCGATATGCGGCGGCACCAGCCCCTGGTGATACAGCACGTAAGAAAGCAGCAGGTAGTCCATGTGACTGCGGTGGCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCGTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGAAATTCGCCGCAATCTCTTCCATCAGTGCAATCGCGTTCTGCTGCGCTTTTTCATGGGAGATT  >  CP009273/4245391‑4245681
                                                                                                                                                  |                                                                                                                                                
aGGTTGATCCCGGCGGCGATATGCGGCGGCACCAGCCCCTGGTGATACAGCACGTAAGAAAGCAGCAGGTAGTCCATGTGACTGCGGTGGCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCttt                                                                                                                                                >  2:426585/1‑149 (MQ=255)
aGGTTGATCCCGGCGGCGATATGCGGCGGCACCAGCCCCTGGTGATACAGCACGTAAGAAAGCAGCAGGTAGTCCATGTGACTGCGGTGGCAAGGCACATATACCAGCTCATGACAGTCGTGGGCCAGCTGGCGAACGCGCTCAGCttt                                                                                                                                                >  1:3044555/1‑149 (MQ=255)
 ggTTGATCCCGGCGGCGATATGCGGCGGCACCAGCCCCTGGTGATACAGCACGTAAGAAAGCAGCAGGTAGTCCATGTGACTGCGGTGGCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTa                                                                                                                                               >  2:3009586/1‑149 (MQ=255)
     gATCCCGGCGGCGATATGCGGCGGCACCAGCCCCTGGTGATACAGCACGTAAGAAAGCAGCAGGTAGTCCATGTGACTGCGGTGGCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGa                                                                                                                                           <  2:1742003/149‑1 (MQ=255)
        cccGGCGGCGATATGCGGCGGCACCAGCCCCTGGTGATACAGCACGTAAGAAAGCAGCAGGTAGTCCATGTGACTGCGGTGGCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGt                                                                                                                                        <  2:149646/149‑1 (MQ=255)
        cccGGCGGCGATATGCGGCGGCACCAGCCCCTGGTGATACAGCACGTAAGAAAGCAGCAGGTAGTCCATGTGACTGCGGTGGCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGt                                                                                                                                        <  2:151684/149‑1 (MQ=255)
        cccGGCGGCGATATGCGGCGGCACCAGCCCCTGGTGATACAGCACGTAAGAAAGCAGCAGGTAGTCCATGTGACTGCGGTGGCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGt                                                                                                                                        >  2:3044590/1‑149 (MQ=255)
            gcggcgATATGCGGCGGCACCAGCCCCTGGTGATACAGCACGTAAGAAAGCAGCAGGTAGTCCATGTGACTGCGGTGGCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGAt                                                                                                                                    <  1:1601213/149‑1 (MQ=255)
            gcggcgATATGCGGCGGCACCAGCCCCTGGTGATACAGCACGTAAGAAAGCAGCAGGTAGTCCATGTGACTGCGGTGGCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGAt                                                                                                                                    <  1:3640872/149‑1 (MQ=255)
                cgATATGCGGCGGCACCAGCCCCTGGTGATACAGCACGTAAGAAAGCAGCAGGTAGTCCATGTGACTGCGGTGGCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGccc                                                                                                                                <  1:38789/149‑1 (MQ=255)
                 gATATGCGGCGGCACCAGACCCTGGTGATACAGCACGTAAGAAAGCAGCAGGTAGTCCATGTGACTGCGGTGGCAAGGCACATATACCAGCTCATGAACGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCt                                                                                                                               <  1:3472176/149‑1 (MQ=255)
                  atatGCGGCGGCACCAGCCCCTGGTGATACAGCACGTAAGAAAGCAGCAGGTAGTCCATGTGACTGCGGTGGCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTg                                                                                                                              >  1:2279879/1‑149 (MQ=255)
                       cggcggCACCAGCCCCTGGTGATACAGCACGTAAGAAAGCAGCAGGTAGTCCATGTGACTGCGGTGGCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTaaa                                                                                                                         >  2:3319010/1‑149 (MQ=255)
                       cggcggCACCAGCCCCTGGTGATACAGCACGTAAGAAAGCAGCAGGTAGTCCATGTGACTGCGGTGGCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTaaa                                                                                                                         <  2:1007766/149‑1 (MQ=255)
                          cggcACCAGCCCCTGGTGATACAGCACGTAAGAAAGCAGCAGGTAGTCCATGTGACTGCGGTGGCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTc                                                                                                                      >  2:501172/1‑149 (MQ=255)
                             cACCAGCCCCTGGTGATACAGCACGTAAGAAAGCAGCAGGTAGTCCATGTGACTGCGGTGGCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGt                                                                                                                   <  2:852211/149‑1 (MQ=255)
                             cACCAGCCCCTGGTGATACAGCACGTAAGAAAGCAGCAGGTAGTCCATGTGACTGCGGTGGCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGt                                                                                                                   <  1:856034/149‑1 (MQ=255)
                              aCCAGCCCCTGGTGATACAGCACGTAAGAAAGCAGCAGGTAGTCCATGTGACTGCGGTGGCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGtt                                                                                                                  <  1:1069225/149‑1 (MQ=255)
                               ccAGCCCCTGGTGATACAGCACGTAAGAAAGCAGCAGGTAGTCCATGTGACTGCGGTGGCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTc                                                                                                                 >  2:3615602/1‑149 (MQ=255)
                               ccAGCCCCTGGTGATACAGCACGTAAGAAAGCAGCAGGTAGTCCATGTGACTGCGGTGGCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTc                                                                                                                 <  2:3970737/149‑1 (MQ=255)
                                 aGCCCCTGGTGATACAGCACGTAAGAAAGCAGCAGGTAGTCCATGTGACTGCGGTGGCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCGTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCa                                                                                                               >  1:3233546/1‑149 (MQ=255)
                                  gCCCCTGGTGATACAGCACGTAAGAAAGCAGCAGGTAGTCCATGTGACTGCGGTGGCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAg                                                                                                              >  1:3384620/1‑149 (MQ=255)
                                  gCCCCTGGTGATACAGCACGTAAGAAAGCAGCAGGTAGTCCATGTGACTGCGGTGGCAAGGCACATATACCAGATCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAg                                                                                                              <  1:1020043/149‑1 (MQ=255)
                                   ccccTGGTGATACAGCACGTAAGAAAGCAGCAGGTAGTCCATGTGACTGCGGTGGCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCGGGTAAAGTCGGTTCCAgg                                                                                                             >  2:3893623/1‑149 (MQ=255)
                                     ccTGGTGATACAGCACGTAAGAAAGCAGCAGGTAGTCCATGTGACTGCGGTGGCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTg                                                                                                           <  2:158267/149‑1 (MQ=255)
                                       tGGTGATACAGCACGTAAGAAAGCAGCAGGTAGTCCATGTGACTGCGGTGGCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGaa                                                                                                         <  2:2755535/149‑1 (MQ=255)
                                        ggTGATACAGCACGTAAGAAAGCAGCAGGTAGTCCATGTGACTGCGGTGGCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAg                                                                                                        >  2:1298421/1‑149 (MQ=255)
                                        ggTGATACAGCACGTAAGAAAGCAGCAGGTAGTCCATGTGACTGCGGTGGCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAg                                                                                                        >  1:3467454/1‑149 (MQ=255)
                                          tGATACAGCACGTAAGAAAGCAGCAGGTAGTCCATGTGACTGCGGTGGCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGcc                                                                                                      <  2:242577/149‑1 (MQ=255)
                                          tGATACAGCACGTAAGAAAGCAGCAGGTAGTCCATGTGACTGCGGTGGCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGcc                                                                                                      >  1:3830586/1‑149 (MQ=255)
                                           gATACAGCACGTAAGAAAGCAGCAGGTAGTCCATGTGACTGCGGTGGCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGccc                                                                                                     >  1:3518864/1‑149 (MQ=255)
                                            aTACAGCACGTAAGAAAGCAGCAGGTAGTCCATGTGACTGCGGTGGCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCa                                                                                                    <  2:3853206/149‑1 (MQ=255)
                                             tACAGCACGTAAGAAAGCAGCAGGTAGTCCATGTGACTGCGGTGGCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAg                                                                                                   >  1:10368/1‑149 (MQ=255)
                                             tACAGCACGTAAGAAAGCAGCAGGTAGTCCATGTGACTGCGGTGGCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAg                                                                                                   >  1:2484638/1‑149 (MQ=255)
                                              aCAGCACGTAAGAAAGCAGCAGGTAGTCCATGTGACTGCGGTGGCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGa                                                                                                  >  2:2589464/1‑149 (MQ=255)
                                               cAGCACGTAAGAAAGCAGCAGGTAGTCCATGTGACTGCGGTGGCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGaa                                                                                                 >  2:1251719/1‑149 (MQ=255)
                                                aGCACGTAAGAAAGCAGCAGGTAGTCCATGTGACTGCGGTGGCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAAt                                                                                                <  1:3488500/149‑1 (MQ=255)
                                                 gCACGTAAGAAAGCAGCAGGTAGTCCATGTGACTGCGGTGGCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATa                                                                                               <  1:182058/149‑1 (MQ=255)
                                                 gCACGTAAGAAAGCAGCAGGTAGTCCATGTGACTGCGGTGGCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATa                                                                                               <  2:2215799/149‑1 (MQ=255)
                                                   aCGTAAGAAAGCAGCAGGTAGTCCATGTGACTGCGGTGGCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACg                                                                                             >  1:2074250/1‑149 (MQ=255)
                                                     gTAAGAAAGCAGCAGGTAGTCCATGTGACTGCGGTGGCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGgt                                                                                           <  1:3494559/149‑1 (MQ=255)
                                                     gTAAGAAAGCAGCAGGTAGTCCATGTGACTGCGGTGGCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGgt                                                                                           <  1:2550801/149‑1 (MQ=255)
                                                      tAAGAAAGCAGCAGGTAGTCCATGTGACTGCGGTGGCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGgtc                                                                                          <  1:3074682/149‑1 (MQ=255)
                                                         gAAAGCAGCAGGTAGTCCATGTGACTGCGGTGGCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGgtcagt                                                                                       >  2:2287916/1‑149 (MQ=255)
                                                          aaaGCAGCAGGTAGTCCATGTGACTGCGGTGGCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGgtcagtc                                                                                      <  1:1631104/149‑1 (MQ=255)
                                                          aaaGCAGCAGGTAGTCCATGTGACTGCGGTGGCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGgtcagtc                                                                                      >  2:3186347/1‑149 (MQ=255)
                                                           aaGCAGCAGGTAGTCCATGTGACTGCGGTGGCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGgtcagtca                                                                                     <  1:438623/149‑1 (MQ=255)
                                                            agcagcAGGTAGTCCATGTGACTGCGGTGGCAAGGCACATATACGAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGtcagtcag                                                                                    <  1:3439755/149‑1 (MQ=255)
                                                               agcagGTAGTCCATGTGACTGCGGTGGCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCg                                                                                 >  1:1687452/1‑149 (MQ=255)
                                                                gcagGTAGTCCATGTGACTGCGGTGGCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGa                                                                                <  1:2437816/149‑1 (MQ=255)
                                                                gcagGTAGTCCATGTGACTGCGGTGGCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGa                                                                                >  1:4026346/1‑149 (MQ=255)
                                                                  agGTAGTCCATGTGACTGCGGTGGCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAat                                                                              <  1:3637167/149‑1 (MQ=255)
                                                                        tCCATGTGACTGCGGTGGCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCAtctc                                                                        <  1:1692432/149‑1 (MQ=255)
                                                                        tCCATGTGACTGCGGTGGCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCAtctc                                                                        <  2:3425147/149‑1 (MQ=255)
                                                                                  tGCGGTGGCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGaaa                                                              >  2:739216/1‑149 (MQ=255)
                                                                                   gCGGTGGCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGAAAt                                                             <  1:1965672/149‑1 (MQ=255)
                                                                                   gCGGTGGCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGAAAt                                                             <  1:3976070/149‑1 (MQ=255)
                                                                                    cGGTGGCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGAAAtt                                                            <  1:2661658/149‑1 (MQ=255)
                                                                                    cGGTGGCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGAAAtt                                                            <  1:3108746/149‑1 (MQ=255)
                                                                                       tGGCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGAAATTcgc                                                         <  1:2010853/149‑1 (MQ=255)
                                                                                       tGGCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGAAATTcgc                                                         <  1:2901177/149‑1 (MQ=255)
                                                                                        ggCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGAAATTcgcc                                                        >  1:90709/1‑149 (MQ=255)
                                                                                        ggCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGAAATTcgcc                                                        <  2:1381302/149‑1 (MQ=255)
                                                                                        ggCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGAAATTcgcc                                                        <  2:3477105/149‑1 (MQ=255)
                                                                                        ggCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGAAATTcgcc                                                        >  1:397230/1‑149 (MQ=255)
                                                                                             ggCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGAAATTCGCCGCAAt                                                   >  2:3106427/1‑149 (MQ=255)
                                                                                                   tataCCAGATCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATTCGGTGAGTCAGGCGATTCATCTCGTAAGAGACATTCGCCGCAATGTCTTc                                             >  2:3240357/1‑149 (MQ=255)
                                                                                                       ccAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGAAATTCGCCGCAATCTCTTCCATc                                         <  2:1579700/149‑1 (MQ=255)
                                                                                                       ccAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGAAATTCGCCGCAATCTCTTCCATc                                         <  1:1784653/149‑1 (MQ=255)
                                                                                                       ccAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGAAATTCGCCGCAATCTCTTCCATc                                         <  1:3319010/149‑1 (MQ=255)
                                                                                                          gCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGAAATTCGCCGCAATCTCTTCCATCAGt                                      >  2:403816/1‑149 (MQ=255)
                                                                                                          gCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGAAATTCGCCGCAATCTCTTCCATCAGt                                      >  1:271341/1‑149 (MQ=255)
                                                                                                          gCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGAAATTCGCCGCAATCTCTTCCATCAGt                                      <  1:3916855/149‑1 (MQ=255)
                                                                                                           cTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGAAATTCGCCGCAATCTCTTCCATCAGTg                                     >  2:451092/1‑149 (MQ=255)
                                                                                                            tCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGAAATTCGCCGCAATCTCTTCCATCAGTGc                                    >  1:519653/1‑149 (MQ=255)
                                                                                                                gACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGACATTCGCCGCAATCTCTTCCATCAGTGCAATc                                >  2:380271/1‑149 (MQ=255)
                                                                                                                    gtcgtGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGAAATTCGCCGCAATCTCTTCCATCAGTGCAATCGCGt                            >  1:1114896/1‑149 (MQ=255)
                                                                                                                     tcgtGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGAAATTCGCCGCAATCTCTTCCATCAGTGCAATCGCGtt                           <  1:735025/149‑1 (MQ=255)
                                                                                                                      cgtGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGAAATTCGCCGCAATCTCTTCCATCAGTGCAATCGCGTTc                          <  1:1565190/149‑1 (MQ=255)
                                                                                                                        tGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGAAATTCGCCGCAATCTCTTCCATCAGTGCAATCGCGTTctg                        >  1:3399070/1‑149 (MQ=255)
                                                                                                                        tGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGAAATTCGCCGCAATATCTTCCATCAGTGCAATCGCGTTctg                        <  2:394309/149‑1 (MQ=255)
                                                                                                                         gggCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGAAATTCGCCGCAATCTCTTCCATCAGTGCAATCGCGTTctgc                       <  1:193491/149‑1 (MQ=255)
                                                                                                                         gggCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGAAATTCGCCGCAATCTCTTCCATCAGTGCAATCGCGTTctgc                       >  2:818921/1‑149 (MQ=255)
                                                                                                                            ccAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGAAATTCGCCGCAATCTCTTCCATCAGTGCAATCGCGTTCtgctgc                    <  2:1787895/149‑1 (MQ=255)
                                                                                                                             cAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGAAATTCGCCGCAATCTCTTCCATCAGTGCAATCGCGTTCTGCTgcg                   >  1:2612360/1‑149 (MQ=255)
                                                                                                                                cTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGAAATTCGCCGCAATCTCTTCCATCAGTGCAATCGCGTTCTGCTGCGCtt                <  1:1873613/149‑1 (MQ=255)
                                                                                                                                cTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGAAATTCGCCGCAATCTCTTCCATCAGTGCAATCGCGTTCTGCTGCGCtt                <  1:1298421/149‑1 (MQ=255)
                                                                                                                                 tGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGAAATTCGCCGCAATCTCTTCCATCAGTGCAATCGCGTTCTGCTGCGCttt               <  2:1938208/149‑1 (MQ=255)
                                                                                                                                   gCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGAAATTCGCCGCAATCTCTTCCATCAGTGCAATCGCGTTCTGCTGCGCttttt             <  2:2687440/149‑1 (MQ=255)
                                                                                                                                    cGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGAAATTCGCCGCAATCTCTTCCATCAGTGCAATCGCGTTCTGCTGCGCTTTTTc            <  2:3916379/149‑1 (MQ=255)
                                                                                                                                    cGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGAAATTCGCCGCAATCTCTTCCATCAGTGCAATCGCGTTCTGCTGCGCTTTTTc            <  1:3146863/149‑1 (MQ=255)
                                                                                                                                     gAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGAAATTCGCCGCAATCTCTTCCATCAGTGCAATCGCGTTCTGCTGCGCTTTTTCa           <  2:2150228/149‑1 (MQ=255)
                                                                                                                                      aaCGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGAAATTCGCCGCAATCTCTTCCATCAGTGCAATCGCGTTCTGCTGCGCTTTTTCAt          >  2:1967625/1‑149 (MQ=255)
                                                                                                                                       aCGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGAAATTCGCCGCAATCTCTTCCATCAGTGCAATCGCGTTCTGCTGCGCTTTTTCATg         <  2:3006792/149‑1 (MQ=255)
                                                                                                                                          cgcTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGAAATTCGCCGCAATCTCTTCCATCAGTGCAATCGCGTTCTGCTGCGCTTTTTCATGGga      <  2:517932/149‑1 (MQ=255)
                                                                                                                                          cgcTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGAAATTCGCCGCAATCTCTTCCATCAGTGCAATCGCGTTCTGCTGCGCTTTTTCATGGga      <  2:682487/149‑1 (MQ=255)
                                                                                                                                           gcTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGAAATTCGCCGCAATCTCTTCCATCAGTGCAATCGCGTTCTGCTGCGCTTTTTCATGGgag     >  2:1029950/1‑149 (MQ=255)
                                                                                                                                            cTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGAAATTCGCCGCAATCTCTTCCATCAGTGCAATCGCGTTCTGCTGCGCTTTTTCATGGgaga    <  1:4042692/149‑1 (MQ=255)
                                                                                                                                             tCAGCTTTATGGACGTTGATGCCCTGGAAAAGTCGGTTCCAGGTGAAGCCCAGAAAACGGTCAGTCAGGCGAATCATCTCGTAAGAGAAATTCGCCGCAATCTCTTCCATCAGTGCAATCGCGTTCTGCTGCGCTTTTTCATGGGAGAt   >  2:2823809/1‑149 (MQ=255)
                                                                                                                                              cAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGAAATTCGCCGCAATCTCTTCCATCAGTGCAATCGCGTTCTGCTGCGCTTTTTCATGGGAGAtt  <  2:4110380/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                
AGGTTGATCCCGGCGGCGATATGCGGCGGCACCAGCCCCTGGTGATACAGCACGTAAGAAAGCAGCAGGTAGTCCATGTGACTGCGGTGGCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCGTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGAAATTCGCCGCAATCTCTTCCATCAGTGCAATCGCGTTCTGCTGCGCTTTTTCATGGGAGATT  >  CP009273/4245391‑4245681

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: