Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,245,568 T→C 6.0% N276S (AAC→AGC)  plsB ← glycerol‑3‑phosphate O‑acyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,245,5680TC 6.0% 230.5 / 3.2 83N276S (AAC→AGC) plsBglycerol‑3‑phosphate O‑acyltransferase
Reads supporting (aligned to +/- strand):  ref base T (49/29);  new base C (2/3);  total (51/32)
Fisher's exact test for biased strand distribution p-value = 3.69e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.01e-01

TGGTGATACAGCACGTAAGAAAGCAGCAGGTAGTCCATGTGACTGCGGTGGCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCGTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGAAATTCGCCGCAATCTCTTCCATCAGTGCAATCGCGTTCTGCTGCGCTTTTTCATGGGAGATTTTTTTGCTGCGCGCTTCATCTTCTACCGCTTTGG  >  CP009273/4245430‑4245715
                                                                                                                                          |                                                                                                                                                   
tGGTGATACCGCACGTAAGAAAGCAGCAGGTAGTCCATGTGACTGCGGTGGCAAGGCACATATACCAGCCCATGACCGTCTTGGGCCAGCTGGCGAACGCGCTCATCGTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCa                                                                                                                                                 <  2:565390/143‑1 (MQ=255)
     aTACAGCACGTAAGAACGCAGAAGGTAGCCCATGTGACTGCGGTGGCAAGGCACCTATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCGTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCa                                                                                                                                      <  1:450820/149‑1 (MQ=255)
     aTACAGCACGTAAGAAAGCAGCAGGTAGTCCATGTGACTGCGGTGGCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTAGCGAACGCGCTCAGCGTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCCGGTGAAGcccc                                                                                                                                      >  1:878912/1‑148 (MQ=255)
       aCAGCACGTAAGAAAGCAGCAGGTAGTCCATGTGACTGCGGTGGCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCGTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAg                                                                                                                                          >  2:236546/1‑143 (MQ=255)
       aCAGCACGTAAGAAAGCAGCAGGTAGTCCATGTGACTGCGGTGGCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCGTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGa                                                                                                                                    <  1:421053/149‑1 (MQ=255)
        cAGCACGTAAGAAAGCAGCAGGTAGTCCATGTGACTGCGGTGGCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCGGTATGGACGTTGATGCCCTGGGAAAGTCGGTTCCAGGTGGAGc                                                                                                                                         >  2:263985/1‑143 (MQ=255)
          gCACGTAAGAAAGCAGCAGGTAGTCCATGTGACTGCGGTGGCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGTACGTGCTCAACGTTATGGACGTTGATGCCCTGGTAAAGTCGGCTCCCGGTGAAGACCAGAATa                                                                                                                                 >  1:626371/1‑149 (MQ=255)
          gCACGTAAGAAAGCAGCAGGTAGTCCATGTGACTGCGGTGGCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCGTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGccc                                                                                                                                       >  2:455936/1‑143 (MQ=255)
          gCACGTAAGAAAGCAGCAGGTAGTCCATGTGACTGCGGTGGCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCGTTATGGACGTTGATGCCCTGGTAAAGTCGGCTCCAGGTGAAGCCCAGAATa                                                                                                                                 >  1:317826/1‑149 (MQ=255)
                aaGAAAGCAGCAGGTAGTCCATGTGACTGCGGTGGCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCGTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATa                                                                                                                                 <  2:216034/143‑1 (MQ=255)
                  gAAAGCAGCAGGTAGTCCATGTGACTGCGGTGGCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCGTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCCGGTGGAGCCCGGGATACGgtcagt                                                                                                                         >  1:17653/1‑149 (MQ=255)
                  gAAAGCAGCAGGTAGTCCATGTGACTGCGGTGGCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCGTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACg                                                                                                                               <  2:437457/143‑1 (MQ=255)
                  gAAAGCAGCAGGTAGTCCATGTGACTGCAGTGGCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCGTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGgtcagt                                                                                                                         >  1:693524/1‑149 (MQ=255)
                   aaaGCAGCAGGTAGTCCATGTGACTGCGGTGGCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCGTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGGGAAGCCCAGAATAAgg                                                                                                                              >  2:282832/1‑143 (MQ=255)
                   aaaGCAGCAGGTAGTCCATGTGACTGCGGTGGCAAGGCACATATACCAGATCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCGTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGGGGAGCCCAGAATACGtacagtc                                                                                                                        >  1:397303/1‑149 (MQ=255)
                          cagGTAGTCCATGTGACTGCGGTGGCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCGTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGgtcagtca                                                                                                                       >  2:23225/1‑143 (MQ=255)
                           agGTAGTCCATGTGACTGCGGTGGCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCGTTATGGACGTTGATGCCCTGGTAAAGTCGGCTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAat                                                                                                                <  1:434490/149‑1 (MQ=255)
                            ggTAGTCCATGTGACTGCGGTGGCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCGTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAgg                                                                                                                     <  2:993400/143‑1 (MQ=255)
                                  ccATGTGACTGCGGTGGCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCGTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCg                                                                                                         <  1:197560/149‑1 (MQ=255)
                                           tGCGGTGGCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCGTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGaaa                                                                                                >  1:577792/1‑149 (MQ=255)
                                           tGCGGGGGCAAGGCACATATACCAGCCCCTGACCGACGTGCGCCAGCTGGCGAACGCGCACAGCGTTATGGACGTTGATGCCCCGGTAAAGTCGGTTCCAGGTGAAACCCAGAATACGGTCAGTCAGGCGAATCATCTCGTaa                                                                                                      >  2:325494/1‑143 (MQ=255)
                                              ggTGGCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCGTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCACGTAAGAGAAAATc                                                                                             >  1:1118096/1‑149 (MQ=255)
                                                tGGCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCGTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGAAATTcgc                                                                                           >  1:275740/1‑149 (MQ=255)
                                                  gCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCGTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGAAAt                                                                                               >  2:316143/1‑143 (MQ=255)
                                                  gCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCGTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGAAATTcgccg                                                                                         >  1:639465/1‑149 (MQ=255)
                                                   cAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCGTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTACGAGAAAtt                                                                                              >  2:131310/1‑143 (MQ=255)
                                                   cAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCGTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGGGAAGCCCAGAATCAGGTCAGTCAGGCGAATCATCTCGTAAGAGAAAtt                                                                                              >  2:954057/1‑143 (MQ=255)
                                                    aacgCACATATACCAGCTCCTGCCCGTCCTGGGCCAGCTGGCGAACGCGCTCAGCGTTATGGACGTTGATGCCCTGTGAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGCCAGTCAGGCGAATCATCTCGTAAGAGAAATTCGCCGCa                                                                                       <  1:347234/146‑1 (MQ=255)
                                                    aaGGCACATATACCAGCTCATGACCGTCGTGGGCCGGCTGGCGAACGCGCTCAGCGTTATGGACGTTGATGCCCTGGTAAAGTTGGTTCCAGGTGAAGCCCAGAATACGGTCATTCAGGCGAATCATCTCGTAAGAGAAATTCGCCGCa                                                                                       >  1:222155/1‑149 (MQ=255)
                                                     aGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCGTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGGGAAATTcg                                                                                            >  2:132179/1‑143 (MQ=255)
                                                            tataCCAGCTCATGGCCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCGTTAGGGACGGTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGAAATTCGCCGCAATCTCTTc                                                                               <  1:1144896/149‑1 (MQ=255)
                                                                ccAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCGTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGAAATTCGCCGCAATCTCTTCCATc                                                                           >  1:672064/1‑149 (MQ=255)
                                                                 cAGCTCCTGCCCGGCGTGGGCCAGCTGGCGAACGCGCTCAGCGTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGAAATTCGCCGCAATCTCTTCCATCa                                                                          <  1:758804/149‑1 (MQ=255)
                                                                 cAGCTCATGACCGTCGTGGGCCAGCCGGCGAACGCGCTCAGCGTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGAAATTCGCCGCAATCTCTTCCATCa                                                                          <  1:668806/149‑1 (MQ=255)
                                                                 cAGCTCATGACCGGCGTGGGCCAGCTGGCGAACGCGCTCAGCGTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGAAATTCGCCGCAATCTCtt                                                                                <  2:625575/143‑1 (MQ=255)
                                                                 cAGCCCATCACCTCCGTGGGCCAGCTGGCGAACGCGCTCAGCGTTATGGACGTTGATGCCCTGGTAAAGGCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGAAATTCGCCGCAATCTCtt                                                                                <  2:142581/143‑1 (MQ=255)
                                                                    cTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCGTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGAAATTCGCCGCAATCTCTTCCATCAGTg                                                                       >  1:1130410/1‑149 (MQ=255)
                                                                    cTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCGTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGAAATTCGCCGCAATCTCTTCCATCAGTg                                                                       <  1:125266/149‑1 (MQ=255)
                                                                     tCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCGTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGAAATTCGCCGCAATCTCTTCCAt                                                                            >  2:886318/1‑143 (MQ=255)
                                                                       aTGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCGTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGAAATTCGCCGCAATCTCTTCCATCAGTGCaa                                                                    >  1:899972/1‑149 (MQ=255)
                                                                                   ggCCAGCTGGCGAACGCGCTCAGCGTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGAAATTCGCCGCAATCTCTTCCATCAGTGCAATCGCGTTctgct                                                        >  1:374044/1‑149 (MQ=255)
                                                                                    gCCAGCTGGCGAACGCGCTCAGCGTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCCGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGCAATTCGCCGCAATCTCTTCCAACAGTGCAATCGCGtt                                                             >  2:75733/1‑143 (MQ=255)
                                                                                    gCCAGCTGGCGAACGCGCTCAGCGTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGAAATTCGCCGCAATCTCTTCCATCAGTGCAATCGCGtt                                                             >  2:75731/1‑143 (MQ=255)
                                                                                         cTGGCGAACGCGCTCAGCGTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCCGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGAAATTCGCCGCAATCTCTTCCATCAGTGCAATCACGGTCTGCTGCGCtt                                                  >  1:264663/1‑149 (MQ=255)
                                                                                            gCGAACGCGCTCAGCGTTATGGCCGGTGATGCCCTGGTAAAGGCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGAAATTCGCCGCAATCTCTTCCATCAGTGCAATCGCGTTCTGCTGCGCttttt                                               <  1:48026/149‑1 (MQ=255)
                                                                                             cGAACGCGCTCAGCGTTATGGACGTTGCTGCACTTGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGAAATTCGCCGCAATCTCTTCCATCAGTGCAATCGCGGTCTGCTGCGCTTTTTc                                              >  1:636790/1‑149 (MQ=255)
                                                                                             cGAACGCGCTCAGCGTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGAGATTCGCCGCAATCTCTTCCATCAGTGCAATCGCGTTCTCCTGCGCTTTTTc                                              >  1:379800/1‑149 (MQ=255)
                                                                                              gAACGCGCTCAGCGTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGAAATTCGCCGCAATCTCTTCCCTCAGTGCAATCGCGTTCTGCTGCGCTTTGTCa                                             >  1:1106270/1‑149 (MQ=255)
                                                                                               aaCGCGCTCAGCGTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGAAATTCGCCGCAATCTCTTCCATCAGTGCAATCGCGTTCTGCTGCGCttttaaat                                            >  1:292712/1‑145 (MQ=255)
                                                                                                aCGCGCTCAGCGTTATGGACGTCGATGACCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGAAATTCGCCGCAATCTCTTCCATCAGTGCAATCGCGTTCTGCTGCGCTTTTTCATg                                           <  1:1294987/149‑1 (MQ=255)
                                                                                                 cgcgcgGAGCGTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGGGACTCATCTCGTAAGCGAAATTCGCCGCAATCTCTTCCATCAGTGCAATCGCGTTCTGCTGCGCTTTTTCATgg                                          <  1:505407/149‑1 (MQ=255)
                                                                                                 cgcgCTCAGCGTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGAAATTCGCCGCAATCTCTTCCATCAGTGCAATCGCGTTCTGCTGCGCTTTTTCATgg                                          <  1:621236/149‑1 (MQ=255)
                                                                                                  gcgcgcAGCGTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCAGCTCGTAAGAGAAATTCGCCGCAATCTCTTCCATCAGTGCAATCGCGTTCTGCTGCGCttttt                                               <  2:768953/143‑1 (MQ=255)
                                                                                                      tCAGCGTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGAAATTCGCCGCAATCTCTTCCATCAGTGCAATCGCGTTCTGCTGCGCTTTTTCATg                                           >  2:989278/1‑143 (MQ=255)
                                                                                                       cAGCGTTATGGACGTTGATGCCATGGTAAAGTCGGTTCCAGGTGAAGCCCAGCTTTCCTCCAGTTGGGCGAATCATCTCGTAAGAGAAATTCGCCGCAATCTCTTCCATCAGTGCAATCGCGTTCTGCTGCGCTTTTTCATGGGAGAtt                                    <  1:624729/149‑1 (MQ=255)
                                                                                                          cGTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGAAATTCGCCGCAATCTCTTCCATCAGTGCAATCGCGTTCTGCTGCGCTTTTTCATGGGAGAttttt                                 >  1:800455/1‑149 (MQ=255)
                                                                                                          cGTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGAAATTCGCCGCAATCTCTTCCATCAGTGCAATCGCGTTCTGCTGCGCTTTTTCATGGGAGAttttt                                 <  1:530623/149‑1 (MQ=255)
                                                                                                          cGTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGAAATTCGCCGCAATCTCTTCCATCAGTGCAATCGCGTTCTGCTGCGCCTTTTCATGGgag                                       >  2:784299/1‑143 (MQ=255)
                                                                                                            ttATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGAAATTCGCCGCAATCTCTTCCATCAGTGCAATCGCGTTCTGCTGCGCTTTTTCATGGGAGAttttttt                               <  1:52504/149‑1 (MQ=255)
                                                                                                            tgagggACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGAAATTCGCCGCAATCTCTTCCATCAGTGCAATCGCGTTCTGCTGCGCTTTTTCATGGGAGAttttttt                               <  1:52843/145‑1 (MQ=255)
                                                                                                              aTGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGAAATTCGCCGCAATCTCTTCCATCAGTGCAATCGCGTTCTGCTGCGCTTTTTCATGGGAGATTTTTTtgc                             >  1:283705/1‑149 (MQ=255)
                                                                                                                ggACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGAAATTCGCCGCAATCTCTTCCATCAGTGCAATCGCGTTCTGCTGCGCTTTTTCATGGGAGATTTTTTtgctg                           >  1:1032205/1‑149 (MQ=255)
                                                                                                                   cGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGAAATTCGCCGCAATCTCTTCCATCAGTGCAATCGCGTTCTGCTGCGCTTTTTCATGGGAGATTTTTTTGCTgcgc                        <  1:1047103/149‑1 (MQ=255)
                                                                                                                       gATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGAAATTCGCCGCAATCTCTTCCATCAGTGCAATCGCGTTCTGCTGCGCTTTTTCATGGGAGATTTTTTTGCTGCGCGCtt                    <  1:1298292/149‑1 (MQ=255)
                                                                                                                       gATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGAAATTCGCCGCAATCTCTTCCATCAGTGCAATCGCGTTCTGCTGCGCTTTTTCATGGGAGATTTTTTTGCTGCGCGCtt                    >  1:366855/1‑149 (MQ=255)
                                                                                                                        aTGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGAAATTCGCCGCAATCTCTTCCATCAGTGCAATCGCGTTCTGCTGCGCTTTTTCATGGGAGATTTTTTTGCTGCGCGCTTc                   <  1:1103474/149‑1 (MQ=255)
                                                                                                                        aTGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGAAATTCGCCGCAATCTCTTCCATCAGTGCAATCGCGTTCTGCTGCGCTTTTTCATGGGAGAATTTTTTGCTgcg                         >  2:447197/1‑143 (MQ=255)
                                                                                                                        aTGCCCTGGTAAAGTCGGCTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGAAATTCGCCGCAATCTCTTCCATCAGTGCAATCGCGTTCTGCTGCGCTTTTTCATGGGAGATTTTTTTGCTgcg                         <  2:535839/143‑1 (MQ=255)
                                                                                                                        aTGCCCTGGTAAAGTCGGCTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGAAATTCGCCGCAATCTCTTCCATCAGTGCAATCGCGTTCTGCTGCGCTTTTTCATGGGAGATTTTTTTGCTgcg                         <  2:277410/143‑1 (MQ=255)
                                                                                                                         tGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGAAATTCGCCGCAATCTCTTCCATCAGTGCAATCGCGTTCTGCTGCGCTTTTTCATGGGAGATTTTTTTGCTgcgc                        >  2:121745/1‑143 (MQ=255)
                                                                                                                         tGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGAAATTCGCCGCAATCTCTTCCATCAGTGCAATCGCGTTCTGCTGCGCTTTTTCATGGGAGATTTTTTTGCTgcgc                        <  2:988981/143‑1 (MQ=255)
                                                                                                                         tGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGAAATTCGCCGCAATCTCTTCCATCAGTGCAATCGCGTTCTGCTGCGCTTTTTCATGGGAGATTTTTTTGCTGCGCGCTTCa                  >  1:1123602/1‑149 (MQ=255)
                                                                                                                          gCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGAAATTCGCCGCAATCCCTTCCATCAGTGCAATCGCGTTCTGCTGCGCCTTTTCATGGGAGATTTTTTTGCTgcgcg                       >  2:207010/1‑143 (MQ=255)
                                                                                                                               ggTAAAGTCGGTTCCAGTTGAAGCCCAGAATACGGTCAGTCAGGCGAATCAGCTCGCAAGAGAAATTCGCCGCAATCTCTTCCATCAGTGCAATCGCGTTCTGCTGCGCTTTTTCATGGGAGATTTTTTTGCTGCGCGCTTCa                  <  2:279626/143‑1 (MQ=255)
                                                                                                                                gTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGAAATTCGCCGCAATCTCTTCCATCAGTGCAATCGCGTTCTGCTGCGCTTTTTCATGGGAGATTTTTTTGCTGCGCGCTTCAt                 >  2:620935/1‑143 (MQ=255)
                                                                                                                                gTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGAAATTCGCCGCAATCTCTTCCATCAGTGCAATCGCGTTCTGCTGCGCTTTTTCATGGGAGATTTTTTTGCTGCGCGCTTCATCTTCTa           >  1:355980/1‑149 (MQ=255)
                                                                                                                                gTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGAAATTCGCCGCAAGCTCTTCCATCAGTGCAATCGCGTTCTGCTGCGCTTTTTCATGGGAGATTTTTTTGCTGCGCGCTTCATCTTCTa           >  1:335821/1‑149 (MQ=255)
                                                                                                                                 tAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTTAGTCAGGCGAATCATCTCGTAAGAGAAATTCGCCGCAATCTCTTCCATCAGTGCAATCGCGTTCTGCTGCGCTTTTTCATGGGAGATTTTTTTGCTGCGCGCTTCATCTTCTAc          <  1:271092/149‑1 (MQ=255)
                                                                                                                                    aGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATGTCGTGAGAGAAATTCGCCGCAATCTCTTCCATCAGTGCAATCGCGTTCTGCTGCGCTTTTTCATGGGAGATTTTTTTGCTGCGCGCTTCCtcttc             >  2:561441/1‑143 (MQ=255)
                                                                                                                                     gTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGAAATTCGCCGCAATCTCTTCCATCAGTGCAATCGCGTTCTGCTGCGCTTTTTCATGGGAGATTTTTTTGCTGCGCGCTTCATCTTCTACCGCt      >  1:371827/1‑149 (MQ=255)
                                                                                                                                         gTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGAAATTCGCCGCAATCTCTTCCATCAGTGCAATCGCTTTCTGCTGCGCTTTTTCATGGGAGATTTTTTTGCTGCGCGCTTCATCCTCTACCGCTTTgg  >  1:216097/1‑149 (MQ=255)
                                                                                                                                         gTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGAAATTCGCCGCAATCTCTTCCATCAGTGCAATCGCGTTCTGCTGCGCTTTTTCATGGGAGATTTTTTTGCTGCGCGCTTCATCTTCTACCg        >  2:409820/1‑143 (MQ=255)
                                                                                                                                         gTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGAAATTCGCCGCAATCTCTTCCATCAGTGCAATCGCGTTCTGCGGCGCTTTATCATGGGAGATTTTTTTGCTGCGCGCTTCCTCTTCTACCGCTTTgg  >  1:466095/1‑149 (MQ=255)
                                                                                                                                          |                                                                                                                                                   
TGGTGATACAGCACGTAAGAAAGCAGCAGGTAGTCCATGTGACTGCGGTGGCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCGTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGAAATTCGCCGCAATCTCTTCCATCAGTGCAATCGCGTTCTGCTGCGCTTTTTCATGGGAGATTTTTTTGCTGCGCGCTTCATCTTCTACCGCTTTGG  >  CP009273/4245430‑4245715

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: