Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 631,020 G→T 30.2% Q336K (CAG→AAG)  ybdN ← PAPS reductase‑like domain protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273631,0200GT30.2% 606.0 / 231.1 317Q336K (CAG→AAG) ybdNPAPS reductase‑like domain protein
Reads supporting (aligned to +/- strand):  ref base G (100/121);  new base T (35/61);  total (135/182)
Fisher's exact test for biased strand distribution p-value = 1.74e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

Only 100 of 319 total aligned reads displayed.

ATAGTTCTTCGATTTTGTGGGGCTAAATGATAATGCCCGACACCAGTAATCGTTATTGAGTAAAACTTTGCATATCCGCCGCCAGGAGGGGATATCTTTCGCGCCAATGTCCCCTTGCTGGGTTTGTGGAACCTCGATGCCTTTTTTCTGATACCAGTGCAAATAAATGGCGATTTTATTGCGGTAATGCTCAGCTGTTTTTTCCGGCATGCTATTGAGCAGCAATAATGCATATTCTTGCCAGTCTAAATGTTCTGGCTTGAGGATTTTCCGATGCCCATAGAAATGATTGTCATGT  >  CP009273/630871‑631168
                                                                                                                                                     |                                                                                                                                                    
  aGTTCTTCGATTTTGTGGGGCTAAATGATAATGCCCGACACCAGTAATCGTTATTGAGTAAAACTTTGCATATCCGCCGCCAGGAGGGGATATCTTTCGCGCCAATGTCCCCTTGCTGGGTTTGTGGAACCTCGATGCCTTTTTTCTGAt                                                                                                                                                    <  1:120680/150‑1 (MQ=255)
    ttcttcGATTTTGTGGGGCTAAATGATAATGCCCGACACCAGTAATCGTTATTGAGTAAAACTTTGCATATCCGCCGCCAGGAGGGGATATCTTTCGCGCCAATGTCCCCTTGCTGGGTTTGTGGAACCTCGATGCCTTTTTTCTGATAc                                                                                                                                                  <  2:4988579/150‑1 (MQ=255)
      cttcGATTTTGTGGGGCTAAATGATAATGCCCGACACCAGTAATCGTTATTGAGTAAAACTTTGCATATCCGCCGCCAGGAGGGGATATCTTTCGCGCCAATGTCCCCTTGCTGGGTTTGTGGAACCTCGATGCCTTTTTTCTTATACCa                                                                                                                                                >  1:5089691/1‑150 (MQ=255)
      cttcGATTTTGTGGGGCTAAATGATAATGCCCGACACCAGTAATCGTTATTGAGTAAAACTTTGCATATCCGCCGCCAGGAGGGGATATCTTTCGCGCCAATGTCCCCTTGCTGGGTTTGTGGAACCTCGATGCCTTTTTTCTGATACCa                                                                                                                                                >  1:3372933/1‑150 (MQ=255)
         cGATTTTGTGGGGCTAAATGATAATGCCCGACACCAGTAATCGTTATTGAGTAAAACTTTGCATATCCGCCGCCAGGAGGGGATATCTTTCGCGCCAATGTCCCCTTGCTGGGTTTGTGGAACCTCGATGCCTTTTTTCTGATACCAGTg                                                                                                                                             <  2:1753216/150‑1 (MQ=255)
         cGATTTTGTGGGGCTAAATGATAATGCCCGACACCAGTAATCGTTATTGAGTAAAACTTTGCATATCCGCCGCCAGGAGGGGATATCTTTCGCGCCAATGTCCCCTTGCTGGGTTTGTGGAACCTCGATGCCTTTTTTCTGATACCAGTg                                                                                                                                             <  1:3913206/150‑1 (MQ=255)
          gATTTTGTGGGGCTAAATGATAATGCCCGACACCAGTAATCGTTATTGAGTAAAACTTTGCATATCCGCCGCCAGGAGGGGATATCTTTCGCGCCAATGTCCCCTTGCTGGGTTTGTGGAACCTAGATGCCTTTTTTCTTATACCAGTGc                                                                                                                                            >  1:2427549/1‑150 (MQ=255)
              ttGTGGGGCTAAATGATAATGCCCGACACCAGTAATCGTTATTGAGTAAAACTTTGCATATCCGCCGCCAGGAGGGGATATCTTTCGCGCCAATGTCCCCTTGCTGGGTTTGTGGAACCTCGATGCCTTTTTTCTTATACCAGTGCaaat                                                                                                                                        <  2:1111882/150‑1 (MQ=255)
              gtgtgGGGCTAAATGATAATGCCCGACACCAGTAATCGTTATTGAGTAAAACTTTGCATATCCGCCGCCAGGAGGGGATATCTTTCGCGCCAATGTCCCCTTGCTGGGTTTGTGGAACCTCGATGCCTTTTTTCTTATACCAGTGCaaat                                                                                                                                        <  2:4940330/149‑1 (MQ=255)
                  ggggCTAAATGATAATGCCCGACACCAGTAATCGTTATTGAGTAAAACTTTGCATATCCGCCGCCAGGAGGGGATATCTTTCGCGCCAATGTCCCCTTGCTGGGTTTGTGGAACCTCGATGCCTTTTTTCTGATACCAGTGCaaataaat                                                                                                                                    <  1:1636125/150‑1 (MQ=255)
                   gggCTAAATGATAATGCCCGACACCAGTAATCGTTATTGAGTAAAACTTTGCATATCCGCCGCCAGGAGGGGATATCTTTCGCGCCAATGTCCCCTTGCTGGGTTTGTGGAACCTCGATGCCTTTTTTCTGATACCAGTGCAAATAAATg                                                                                                                                   >  2:3190651/1‑150 (MQ=255)
                   gggCTAAATGATAATGCCCGACACCAGTAATCGTTATTGAGTAAAACTTTGCATATCCGCCGCCAGGAGGGGATATCTTTCGCGCCAATGTCCCCTTGCTGGGTTTGTGGAACCTCGATGCCTTTTTTCTGATACCAGTGCAAATAAATg                                                                                                                                   <  1:975124/150‑1 (MQ=255)
                       tAAATGATAATGCCCGACACCAGTAATCGTTATTGAGTAAAACTTTGCATATCCGCCGCCAGGAGGGGATATCTTTCGCGCCAATGTCCCCTTGCTGGGTTTGTGGAACCTCGATGCCTTTTTTCTTATACCAGTGCAAATAAATGGCGa                                                                                                                               >  2:3500359/1‑150 (MQ=255)
                       tAAATGATAATGCCCGACACCAGTAATCGTTATTGAGTAAAACTTTGCATATCCGCCGCCAGGAGGGGATATCTTTCGCGCCAATGTCCCCTTGCTGGGTTTGTGGAACCTCGATGCCTTTTTTCTTATACCAGTGCAAATAAATGGCGa                                                                                                                               >  2:826413/1‑150 (MQ=255)
                        aaaTGATAATGCCCGACACCAGTAATCGTTATTGAGTAAAACTTTGCATATCCGCCGCCAGGAGGGGATATCTTTCGCGCCAATGTCCCCTTGCTGGGTTTGTGGAACCTCGATGCCTTTTTTCTGATACCAGTGCAAATAAATGGCGAt                                                                                                                              <  2:2135451/150‑1 (MQ=255)
                           tGATAATGCCCGACACCAGTAATCGTTATTGAGTAAAACTTTGCATATCCGCCGCCAGGAGGGGATATCTTTCGCGCCAATGTCCCCTTGCTGGGTTTGTGGAACCTCGATGCCTTTTTTCTGATACCAGTGCAAATAAATGGCGAtttt                                                                                                                           <  2:4084158/150‑1 (MQ=255)
                           tGATAATGCCCGACACCAGTAATCGTTATTGAGTAAAACTTTGCATATCCGCCGCCAGGAGGGGATATCTTTCGCGCCAATGTCCCCTTGCTGGGTTTGTGGAACCTCGATGCCTTTTTTCTGATACCAGTGCAAATAAATGGCGAtttt                                                                                                                           <  2:1421894/150‑1 (MQ=255)
                            gATAATGCCCGACACCAGTAATCGTTATTGAGTAAAACTTTGCATATCCGCCGCCAGGAGGGGATATCTTTCGCGCCAATGTCCCCTTGCTGGGTTTGTGGAACCTCGATGCCTTTTTTCTGATACCAGTGCAAATAAATGGCGATTTTa                                                                                                                          <  2:4932425/150‑1 (MQ=255)
                             aTAATGCCCGACACCAGTAATCGTTATTGAGTAAAACTTTGCATATCCGCCGCCAGGAGGGGATATCTTTCGCGCCAATGTCCCCTTGCTGGGTTTGTGGAACCTCGATGCCTTTTTTCTGATACCAGTGCAAATAAATGGCGATTTTAt                                                                                                                         <  1:1285939/150‑1 (MQ=255)
                              tAATGCCCGACACCAGTAATCGTTATTGAGTAAAACTTTGCATATCCGCCGCCAGGAGGGGATATCTTTCGCGCCAATGTCCCCTTGCTGGGTTTGTGGAACCTCGATGCCTTTTTTCTGATACCAGTGCAAATAAATGGCGATTTTAtt                                                                                                                        >  2:692088/1‑150 (MQ=255)
                               aaTGCCCGACACCAGTAATCGTTATTGAGTAAAACTTTGCATATCCGCCGCCAGGAGGGGATATCTTTCGCGCCAATGTCCCCTTGCTGGGTTTGTGGAACCTCGATGCCTTTTTTCTGATACCAGTGCAAATAAATGGCGATTTTATTg                                                                                                                       >  2:3784617/1‑150 (MQ=255)
                                aTGCCCGACACNAGTAATCGTTATTGAGTAAAACTTTGCATATCCGCCGCCAGGAGGGGATATCTTTCGCGCCAATGTCCCCTTGCTGGGTTTGTGGAACCTCGATGCCTTTTTTCTGATACCAGTGCAAATAAATGGCGATTTTATTGc                                                                                                                      >  1:2761037/1‑150 (MQ=255)
                                 tGCCCGACACCAGTAATCGTTATTGAGTAAAACTTTGCATATCCGCCGCCAGGAGGGGATATCTTTCGCGCCAATGTCCCCTTGCTGGGTTTGTGGAACCTCGATGCCTTTTTTCTGATACCAGTGCAAATAAATGGCGATTTTATTGCg                                                                                                                     >  2:1575026/1‑150 (MQ=255)
                                      gACACCAGTAATCGTTATTGAGTAAAACTTTGCATATCCGCCGCCAGGAGGGGATATCTTTCGCGCCAATGTCCCCTTGCTGGGTTTGTGGAACCTCGATGCCTTTTTTCTGATACCAGTGCAAATAAATGGCGATTTTATTGCGGTAAt                                                                                                                >  2:4591284/1‑150 (MQ=255)
                                      gACACCAGTAATCGTTATTGAGTAAAACTTTGCATATCCGCCGCCAGGAGGGGATATCTTTCGCGCCAATGTCCCCTTGCTGGGTTTGTGGAACCTCGATGCCTTTTTTCTGATACCAGTGCAAATAAATGGCGATTTTATTGCGGTAAt                                                                                                                <  2:1200639/150‑1 (MQ=255)
                                       acacCAGTAATCGTTATTGAGTAAAACTTTGCATATCCGCCGCCAGGAGGGGATATCTTTCGCGCCAATGTCCCCTTGCTGGGTTTGTGGAACCTCGATGCCTTTTTTCTGATACCAGTGCAAATAAATGGCGATTTTATTGCGGTAATg                                                                                                               <  2:4480275/150‑1 (MQ=255)
                                            aGTAATCGTTATTGAGTAAAACTTGGCATATCCGCCGCCAGGAGGGGATATCTTTCGCGCCAATGTCCCCTTGCTGGGTTTGTGGAACCTCGAGGCCTTTTTTCTTATACCAGTGCAAATAAATGGCGATTTTATTGCGGTAATGCTCAg                                                                                                          <  2:54991/150‑1 (MQ=255)
                                             gTAATCGTTATTGAGTAAAACTTTGCATATCCGCCGCCAGGAGGGGATATCTTTCGCGCCAATGTCCCCTTGCTGGGTTTGTGGAACCTCGATGCCTTTTTTCTTATACCAGTGCAAATAAATGGCGATTTTATTGCGGTAATGCTCAGc                                                                                                         <  1:624536/150‑1 (MQ=255)
                                               aaTCGTTATTGAGTAAAACTTTGCATATCCGCCGCCAGGAGGGGATATCTTTCGCGCCAATGTCCCCTTGCTGGGTTTGTGGAACCTCGATGCCTTTTTTCTTATACCAGTGCAAATAAATGGCGATTTTATTGCGGTAATGCTCAGCTg                                                                                                       >  2:3569862/1‑150 (MQ=255)
                                                   gTTATTGAGTAAAACTTTGCATATCCGCCGCCAGGAGGGGATATCTTTCGCGCCAATGTCCCCTTGCTGGGTTTGTGGAACCTCGATGCCTTTTTTCTGATACCAGTGCAAATAAATGGCGATTTTATTGCGGTAATGCTCAGCTGtttt                                                                                                   >  1:4676630/1‑150 (MQ=255)
                                                    ttATTGAGTAAAACTTTGCATATCCGCCGCCAGGAGGGGATATCTTTCGCGCCAATGTCCCCTTGCTGGGTTTGTGGAACCTCGATGCCTTTTTTCTGATACCAGTGCAAATAAATGGCGATTTTATTGCGGTAATGCTCAGCTGttttt                                                                                                  <  2:2812948/150‑1 (MQ=255)
                                                    ttATTGAGTAAAACTTTGCATATCCGCCGCCAGGAGGGGATATCTTTCGCGCCAATGTCCCCTTGCTGGGTTTGTGGAACCTCGATGCCTTTTTTCTGATACCAGTGCAAATAAATGGCGATTTTATTGCGGTAATGCTCAGCTGttttt                                                                                                  <  1:2066003/150‑1 (MQ=255)
                                                     tATTGAGTAAAACTTTGCATATCCGCCGCCAGGAGGGGATATCTTTCGCGCCAATGTCCCCTTGCTGGGTTTGTGGAACCTCGATGCCTTTTTTCTTATACCAGTGCAAATAAATGGCGATTTTATTGCGGTAATGCTCAGCTGtttttt                                                                                                 >  2:4534736/1‑150 (MQ=255)
                                                      aTTGAGTAAAACTTTGCATATCCGCCGCCAGGAGGGGATATCTTTCGCGCCAATGTCCCCTTGCTGGGTTTGTGGAACCTCGATGCCTTTTTTCTGATACCAGTGCAAATAAATGGCGATTTTATTGCGGTAATGCTCAGCTGTTTTTTc                                                                                                >  1:3318800/1‑150 (MQ=255)
                                                       ttGAGTAAAACTTTGCATATCCGCCGCCAGGAGGGGATATCTTTCGCGCCAATGTCCCCTTGCTGGGTTTGTGGAACCTCGATGCCTTTTTTCTGATACCAGTGCAAATAAATGGCGATTTTATTGCGGTAATGCTCAGCTGTTTTTTcc                                                                                               <  2:3052653/150‑1 (MQ=255)
                                                        tGAGTAAAACTTTGCATATCCGCCGCCAGGAGGGGATATCTTTCGCGCCAATGTCCCCTTGCTGGGTTTGTGGAACCTCGATGCCTTTTTTCTGATACCAGTGCAAATAAATGGCGATTTTATTGCGGTAATGCTCAGCTGTTTTTTCCg                                                                                              >  2:4273740/1‑150 (MQ=255)
                                                          aGTAAAACTTTGCATATCCGCCGCCAGGAGGGGATATCTTTCGCGCCAATGTCCCCTTGCTGGGTTTGTGGAACCTCGATGCCTTTTTTCTTATACCAGTGCAAATAAATGGCGATTTTATTGCGGTAATGCTCAGCTGTTTTTTCCGGc                                                                                            >  1:4653060/1‑150 (MQ=255)
                                                           gTAAAACTTTGCATATCCGCCGCCAGGAGGGGATATCTTTCGCGCCAATGTCCCCTTGCTGGGTTTGTGGAACCTCGATGCCTTTTTTCTGATACCAGTGCAAATAAATGGCGATTTTATTGCGGTAATGCTCAGCTGTTTTTTCCGGCa                                                                                           <  1:3457655/150‑1 (MQ=255)
                                                             aaaaCTTTGCATATCCGCCGCCAGGAGGGGATATCTTTCGCGCCAATGTCCCCTTGCTGGGTTTGTGGAACCTCGATGCCTTTTTTCTGATACCAGTGCAAATAAATGGCGATTTTATTGCGGTAATGCTCAGCTGTTTTTTCCGGCATg                                                                                         >  2:3974504/1‑150 (MQ=255)
                                                             aaaaCTTTGCATATCCGCCGCCAGGAGGGGATATCTTTCGCGCCAATGTCCCCTTGCTGGGTTTGTGGAACCTCGATGCCTTTTTTCTGATACCAGTGCAAATAAATGGCGATTTTATTGCGGTAATGCTCAGCTGTTTTTTCCGGCATg                                                                                         >  1:4792932/1‑150 (MQ=255)
                                                                    tGCATATCCGCCGCCAGGAGGGGATATCTTTCGCGCCAATGTCCCCTTGCTGGGTTTGTGGAACCTCGATGCCTTTTTTCTTATACCAGTGCAAATAAATGGCGATTTTATTGCGGTAATGCTCAGCTGTTTTTTCCGGCATGCTATTGa                                                                                  >  1:3016595/1‑150 (MQ=255)
                                                                    tGCATATCCGCCGCCAGGAGGGGATATCTTTCGCGCCAATGTCCCCTTGCTGGGTTTGTGGAACCTCGATGCCTTTTTTCTGATACCAGTGCAAATAAATGGCGATTTTATTGCGGTAATGCTCAGCTGTTTTTTCCGGCATGCTATTGa                                                                                  >  2:566315/1‑150 (MQ=255)
                                                                     gCATATCCGCCGCCAGGAGGGGATATCTTTCGCGCCAATGTCCCCTTGCTGGGTTTGTGGAACCTCGATGCCTTTTTTCTTATACCAGTGCAAATAAATGGCGATTTTATTGCGGTAATGCTCAGCTGTTTTTTCCGGCATGCTATTGag                                                                                 <  1:1888726/150‑1 (MQ=255)
                                                                       atatCCGCCGCCAGGAGGGGATATCTTTCGCGCCAATTTCCCCTTGCTGGGTTTGTGGAACCTCGATGCCTTTTTTCTTATACCAGTGCAAATAAATGGCGATTTTATTGCGGTAATGCTCAGCTGTTTTTTCCGGCATGCTATTGagca                                                                               <  1:5196686/150‑1 (MQ=255)
                                                                        tatCCGCCGCCAGGAGGGGATATCTTTCGCGCCAATGTCCCCTTGCTGGGTTTGTGGAACCTCGATGCCTTTTTTCTTATACCAGTGCAAATAAATGGCGATTTTATTGCGGTAATGCTCAGCTGTTTTTTCCGGCATGCTATTGagcag                                                                              >  2:2179728/1‑150 (MQ=255)
                                                                        tatCCGCCGCCAGGAGGGGATATCTTTCGCGCCAATGTCCCCTTGCTGGGTTTGTGGAACCTCGATGCCTTTTTTCTGATACCAGTGCAAATAAATGGCGATTTTATTGCGGTAATGCTCAGCTGTTTTTTCCGGCATGCTATTGagcag                                                                              >  1:4943702/1‑150 (MQ=255)
                                                                        tatCCGCCGCCAGGAGGGGATATCTTTCGCGCCAATGTCCCCTTGCTGGGTTTGTGGAACCTCGATGCCTTTTTTCTGATACCAGTGCAAATAAATGGCGATTTTATTGCGGTAATGCTCAGCTGTTTTTTCCGGCATGCTATTGagcag                                                                              >  2:4405342/1‑150 (MQ=255)
                                                                         atCCGCCGCCAGGAGGGGATATCTTTCGCGCCAATGTCCCCTTGCTGGGTTTGTGGAACCTCGATGCCTTTTTTCTGATACCAGTGCAAATAAATGGCGATTTTATTGCGGTAATGCTCAGCTGTTTTTTCCGGCATGCTATTGagcagc                                                                             <  2:238103/150‑1 (MQ=255)
                                                                          tCCGCCGCCAGGAGGGGATATCTTTCGCGCCAATGTCCCCTTGCTGGGTTTGTGGAACCTCGATGCCTTTTTTCTTATACCAGTGCAAATAAATGGCGATTTTATTGCGGTAATGCTCAGCTGTTTTTTCCGGCATGCTATTGAgcagca                                                                            >  1:1744831/1‑150 (MQ=255)
                                                                            cgccgcCAGGAGGGGATATCTTTCGCGCCAATGTCCCCTTGCTGGGTTTGTGGAACCTCGATGCCTTTTTTCTTATACCAGTGCAAATAAATGGCGATTTTATTGCGGTAATGCTCAGCTGTTTTTTCCGGCATGCTATTGAGCAGCaat                                                                          >  2:2709868/1‑150 (MQ=255)
                                                                            cgccgcCAGGAGGGGATATCTTTCGCGCCAATGTCCCCTTGCTGGGTTTGTGGAACCTCGATGCCTTTTTTCTGATACCAGTGCAAATAAATGGCGATTTTATTGCGGTAATGCTCAGCTGTTTTTTCCGGCATGCTATTGAGCAGCaat                                                                          <  1:134156/150‑1 (MQ=255)
                                                                               cgccAGGAGGGGATATCTTTCGCGCCAATGTCCCCTTGCTGGGTTTGTGGAACCTCGATGCCTTTTTTCTTATACCAGTGCAAATAAATGGCGATTTTATTGCGGTAATGCTCAGCTGTTTTTTCCGGCATGCTATTGAGCAGCaataat                                                                       <  2:4293858/150‑1 (MQ=255)
                                                                                gccAGGAGGGGATATCTTTCGCGCCAATGTCCCCTTGCTGGGTTTGTGGAACCTCGATGCCTTTTTTCTGATACCAGTGCAAATAAATGGCGATTTTATTGCGGTAATGCTCAGCTGTTTTTTCCGGCATGCTATTGAGCAGCAATAATg                                                                      <  2:2517423/150‑1 (MQ=255)
                                                                                 ccAGGAGGGGATATCTTTCGCGCCAATGTCCCCTTGCTGGGTTTGTGGAACCTCGATGCCTTTTTTCTGATACCAGTGCAAATAAATGGCGATTTTATTGCGGTAATGCTCAGCTGTTTTTTCCGGCATGCTATTGAGCAGCAATAATGc                                                                     >  1:2830793/1‑150 (MQ=255)
                                                                                  cAGGAGGGGATATCTTTCGCGCCAATGTCCCCTTGCTGGGTTTGTGGAACCTCGATGCCTTTTTTCTGATACCAGTGCAAATAAATGGCGATTTTATTGCGGTAATGCTCAGCTGTTTTTTCCGGCATGCTATTGAGCAGCAATAATGCa                                                                    >  1:4537498/1‑150 (MQ=255)
                                                                                  cAGGAGGGGATATCTTTCGCGCCAATGTCCCCTTGCTGGGTTTGTGGAACCTCGATGCCTTGTTTCTGATACCAGTGCAAATAAATGGCGATTTTATTGCGGAAATGCTCAGCTGTTTTTTCCGGCATGCTATTGAGCAGCAATAATGCa                                                                    >  1:3504267/1‑150 (MQ=255)
                                                                                      aggGGATATCTTTCGCGCCAATGTCCCCTTGCTGGGTTTGTGGAACCTCGATGCCTTTTTTCTGATACCAGTGCAAATAAATGGCGATTTTATTGCGGTAATGCTCAGCTGTTTTTTCCGGCATGCTATTGAGCAGCAATAATGCATAtt                                                                <  2:779589/150‑1 (MQ=255)
                                                                                      aggGGATATCTTTCGCGCCAATGTCCCCTTGCTGGGTTTGTGGAACCTCGATGCCTTTTTTCTGATACCAGTGCAAATAAATGGCGATTTTATTGCGGTAATGCTCAGCTGTTTTTTCCGGCATGCTATTGAGCAGCAATAATGCATAtt                                                                <  2:2072524/150‑1 (MQ=255)
                                                                                        gggATATCTTTCGCGCCAATGGCCCCTTGCTGGGTTTGTGGAACCTCGATGCCTTTTTTCTGATACCAGTGCAAATAAATGGCGATTTTATTGCGGTAATGCTCAGCTGTTTTTTCCGGCATGCTATTGAGCAGCAATAATGCATATTCt                                                              <  1:1081973/150‑1 (MQ=255)
                                                                                         ggATATCTTTCGCGCCAATGTCCCCGTGCTGGGTTTGTGGAACCTCGATGCCTTTTTTCTGATACCAGTGCAAATCAATGGCGATTTTATTGCGGTAATGCTCAGATTTTTTTTCCGGCATGCTATTGAGCTGCAATAATGCATATTCtt                                                             <  1:2247558/150‑1 (MQ=255)
                                                                                          gATATCTTTCGCGCCAATGTCCCCTTGCTGGGTTTGTGGAACCTCGATGCCTTTTTTCTTATACCAGTGCAAATAAATGGCGATTTTATTGCGGTAATGCTCAGCTGTTTTTTCCGGCATGCTATTGAGCAGCAATAATGCATATTCTTg                                                            >  2:2220291/1‑150 (MQ=255)
                                                                                          gATATCTTTCGCGCCAATGTCCCCTTGCTGGGTTTGTGGAACCTCGATGCCTTTTTTCTGATACCAGTGCAAATAAATGGCGATTTTATTGCGGTAATGCTCAGCTGTTTTTTCCGGCATGCTATTGAGCAGCAATAATGCATATTCTTg                                                            >  1:670697/1‑150 (MQ=255)
                                                                                          gATATCTTTCGCGCCAATGTCCCCTTGCTGGGTTTGTGGAACCTCGATGCCTTTTTTCTGATACCAGTGCAAATAAATGGCGATTTTATTGCGGTAATGCTCAGCTGTTTTTTCCGGCATGCTATTGAGCAGCAATAATGCATATTCTTg                                                            <  2:2461369/150‑1 (MQ=255)
                                                                                            tatCTTTCGCGCCAATGTCCCCTTGCTGGGTTTGTGGAACCTCGATGCCTTTTTTCTGATACCAGTTCAAATAAATGGCGATTTTATTGCGGTAATGCTCAGCTGTTTTTTCCGGCATGCTATTGAGCAGCAATAATGCATATTCTTGcc                                                          <  1:749488/150‑1 (MQ=255)
                                                                                            tatCTTTCGCGCCAATGTCCCCTTGCTGGGTTTGTGGAACCTCGATGCCTTTTTTCTGATACCAGTGCAAATAAATGGCGATTTTATTGCGGTAATGCTCAGCTGTTTTTTCCGGCATGCTATTGAGCAGCAATAATGCATATTCTTGcc                                                          <  2:1655145/150‑1 (MQ=255)
                                                                                             atCTTTCGCGCCAATGTCCCCTTGCTGGGTTTGTGGAACCTCGATGCCTTTTTTCTTATACCAGTGCAAATAAATGGCGATTTTATTGCGGTAATGCTCAGCTGTTTTTTCCGGCATGCTATTGAGCAGCAATAATGCATATTCTTGCCa                                                         >  1:4408500/1‑150 (MQ=255)
                                                                                             atCTTTCGCGCCAATGTCCCCTTGCTGGGTTTGTGGAACCTCGATGCCTTTTTTCTGATACCAGTGCAAATAAATGGCGATTTTATTGCGGTAATGCTCAGCTGTTTTTTCCGGCATGCTATTGAGCAGCAATAATGCATATTCTTGCCa                                                         >  1:906543/1‑150 (MQ=255)
                                                                                              tCTTTCGCGCCAATGTCCCCTTGCTGGGTTTGTGGAACCTCGATGCCTTTTTTCTGATACCAGTGCAAATAAATGGCGATTTTATTGCGGTAATGCTCAGCTGTTTTTTCCGGCATGCTATTGAGCAGCAATAATGCATATTCTTGCCAg                                                        <  2:1508482/150‑1 (MQ=255)
                                                                                                     cgcCAATGTCCCCTTGCTGGGTTTGTGGAACCTCGATGCCTTTTTTCTTATACCAGTGCAAATAAATGGCGATTTTATTGCGGTAATGCTCAGCTGTTTTTTCCGGCATGCTATTGAGCAGCAATAATGCATATTCTTGCCAGTCTAAAt                                                 <  2:1876309/150‑1 (MQ=255)
                                                                                                         aaTGTCCCCTTGCTGGGTTTGTGGAACCTCGATGCCTTTTTTCTTATACCAGTGCAAATAAATGGCGATTTTATTGCGGTAATGCTCAGCTGTTTTTTCCGGCATGCTATTGAGCAGCAATAATGCATATTCTTGCCAGTCTAAATGTTc                                             >  1:4314239/1‑150 (MQ=255)
                                                                                                          aTGTCCCCTTTCTTGGTTTGTTGAAACTCGATGCATTTTTTATTATACCAGTGCAAATAAATGGCGATTTTATTGCGGTAATGCTAATCTGTTTTTTCCGGCTTGCTATTGAGCAGCAATAATGCATATTCTTGCCAGTCTAAATGTTCt                                            <  1:2703011/150‑1 (MQ=255)
                                                                                                                ccTTGCTGGGTTTGTGGAACCTCGATGCCTTTTTTCTGATACCAGTGCAAATAAATGGCGATTTTATTGCGGTAATGCTCAGCTGTTTTTTCCGGCATGCTATTGAGCAGCAATAATGCATATTCTTGCCAGTCTAAATGTTCTGGCTTg                                      <  2:4391054/150‑1 (MQ=255)
                                                                                                                ccTTGCTGGGTTTGTGGAACCTCGATGCCTTTTTTCTGATACCAGTGCAAATAAATGGCGATTTTATTGCGGTAATGCTCAGCTGTTTTTTCCGGCATGCTATTGAGCAGCAATAATGCATATTCTTGCCAGTCTAAATGTTCTGGCTTg                                      <  1:1570057/150‑1 (MQ=255)
                                                                                                                  ttGCTGGGTTTGTGGAACCTCGATGCCTTTTTTCTTATACCAGTGCAAATAAATGGCGATTTTATTGCGGTAATGCTCAGCTGTTTTTTCCGGCATGCTATTGAGCAGCAATAATGCATATTCTTGCCAGTCTAAATGTTCTGGCTTGAg                                    <  2:3537594/150‑1 (MQ=255)
                                                                                                                   tGCTGGGTTTGTGGAACCTCGATGCCTTTTTTCTGATACCAGTGCAAATAAATGGCGATTTTATTGCGGTAATGCTCAGCTGTTTTTTCCGGCATGCTATTGAGCAGCAATAATGCATATTCTTGCCAGTCTAAATGTTCTGGCTTGAgg                                   <  2:1165881/150‑1 (MQ=255)
                                                                                                                       gggTTTGTGGAACCTCGATGCCTTTTTTCTTATACCAGTGCAAATAAATGGCGATTTTATTGCGGTAATGCTCAGCTGTTTTTTCCGGCATGCTATTGAGCAGCAATAATGCATATTCTTGCCAGTCTAAATGTTCTGGCTTGAGGAttt                               >  2:5205393/1‑150 (MQ=255)
                                                                                                                          tttGTGGAACCTCGATGCCTTTTTTCTTATACCAGTGCAAATAAATGGCGATTTTATTGCGGTAATGCTCAGCTGTTTTTTCCGGCATGCTATTGAGCAGCAATAATGCATATTCTTGCCAGTCTAAATGTTCTGGCTTGAGGATTTTcc                            <  1:4112683/150‑1 (MQ=255)
                                                                                                                            tgtgGAACCTCGATGCCTTTTTTCTGATACCAGTGCAAATAAATGGCGATTTTATTGCGGTAATGCTCAGCTGTTTTTTCCGGCATGCTATTGAGCAGCAATAATGCATATTCTTGCCAGTCTAAATGTTCTGGCTTGAGGATTTTCCGa                          <  1:1114324/150‑1 (MQ=255)
                                                                                                                              tgGAACCTCGATTCCTTTTTTCTGATACCAGTGCAAATAAATGGCGATTTTATTTCGGTAATGCTCAGATGTTTTTTCCGGCATGCTATTGAGCAGCAATAATGCATATTCTTGCCAGTCTAAATTTTCTGTCTTTAGGATTTTACGATg                        >  1:2491305/1‑150 (MQ=255)
                                                                                                                              tgGAACCTCGATGCCTTTTTTCTTATACCAGTGCAAATAAATGGCGATTTTATTGCTGTAATGCTCAGCTGTTTTTTCCGGCATGCTATTGAGCAGCAATAATGCATATTCTTGCCAGTCTAAATGTTCTGGCTTGAGGATTTTCCGATg                        <  1:2183825/150‑1 (MQ=255)
                                                                                                                              tgGAACCTCGATGCCTTTTTTCTGATACCAGTGCAAATAAATGGCGATTTTATTGCGGTAATGCTCAGCTGTTTTTTCCGGCATGCTATTGAGCAGCAATAATGCATATTCTTGCCAGTCTAAATGTTCTGGCTTGAGGATTTTCCGATg                        >  1:3986759/1‑150 (MQ=255)
                                                                                                                              tgGAACCTAGATGCCTTTTTTCTTATACCAGTGCAAATAAATGGCGATTTTATTGCTGTAATGCTCAGCTGTTTTTTCCGGCATGCTATTGAGCATCAATAATGCATATTCTTGCCAGTATAAATGTTCTGGCTTGAGGATTTTCCGATg                        <  1:4038055/150‑1 (MQ=255)
                                                                                                                                 aaCCTCGATGCCTTTTTTCTTATACCAGTGCAAATAAATGGCGATTTTATTGCGGTAATGCTCAGCTGTTTTTTCCGGCATGCTATTGAGCAGCAATAATGCATATTCTTGCCAGTCTAAATGTTCTGGCTTGAGGATTTTCCGATGccc                     >  1:3744575/1‑150 (MQ=255)
                                                                                                                                    cTCGATGCCTTTTTTCTGATACCAGTGCAAATAAATGGCGATTTTATTGCGGTAATGCTCAGCTGTTTTTTCCGGCATGCTATTGAGCAGCAATAATGCATATTCTTGCCAGTCTAAATGTTCTGGCTTGAGGATTTTCCGATGCCCATa                  >  1:1062822/1‑150 (MQ=255)
                                                                                                                                     tCGATGCCTTTTTTCTTATACCAGTGCAAATAAATGGCGATTTTATTGCGGTAATGCTCAGCTGTTTTTTCCGGCATGCTATTGAGCAGCAATAATGCATATTCTTGCCAGTCTAAATGTTCTGGCTTGAGGATTTTCCGATGCCCATAg                 <  2:4667890/150‑1 (MQ=255)
                                                                                                                                     tCGATGCCTTTTTTCTTATACCAGTGCAAATAAATGGCGATTTTATTGCGGTAATGCTCAGCTGTTTTTTCCGGCATGCTATTGAGCAGCAATAATGCATATTCTTGCCAGTCTAAATGTTCTGGCTTGAGGATTTTCCGATGCCCATAg                 <  1:5003872/150‑1 (MQ=255)
                                                                                                                                     tCGATGCCTTTTTTCTGATACCAGTGCAAATAAATGGCGATTTTATTGCGGTAATGCTCAGCTGTTTTTTCCGGCATGCTATTGAGCAGCAATAATGCATATTCTTGCCAGTCTAAATGTTCTGGCTTGAGGATTTTCCGATGCCCATAg                 <  1:3190651/150‑1 (MQ=255)
                                                                                                                                     tCGATGCCTTTTTTCTGATACCAGTGCAAATAAATGGCGATTTTATTGCGGTAATGCTCAGCTGTTTTTTCCGGCATGCTATTGAGCAGCAATAATGCATATTCTTGCCAGTCTAAATGTTCTGGCTTGAGGATTTTCCGATGCCCATAg                 <  1:3289496/150‑1 (MQ=255)
                                                                                                                                      cGATGCCTTTTTTCTGATACCAGTGCAAATAAATGTCGATTTTATTGCGGTAATGCGCAGCTGTTTTTTCCGGCATGCTATTGAGCAGCAATAATGCATATTCTTGCCAGTCTAAATGTTCTGGCTTGAGGATTTTCCGATGCCAATAGa                >  1:2591926/1‑150 (MQ=255)
                                                                                                                                       gATGCCTTTTTTCTGATACCAGTGCAAATAAATGGCGATTTTATTGCGGTAATGCTCAGCTGTTTTTTCCGGCATGCTATTGAGCAGCAATAATGCATATTCTTGCCAGTCTAAATGTTCTGGCTTGAGGATTTTCCGATGCCCATAGaa               >  1:3139409/1‑150 (MQ=255)
                                                                                                                                       gATGCCTTTTTTCTGATACCAGTGCAAATAAATGGCGATTTTATTGCGGTAATGCTCAGCTGTTTTTTCCGGCATGCTATTGAGCAGCAATAATGCATATTCTTGCCAGTCTAAATGTTCTGGCTTGAGGATTTTCCGATGCCCATAGaa               >  1:843463/1‑150 (MQ=255)
                                                                                                                                          gCCTTTTTTCTGATACCAGTGCAAATAAATGGCGATTTTATTGCGGTAATGCTCAGCTGTTTTTTCCGGCATGCTATTGAGCAGCAATAATGCATATTCTTGCCAGTCTAAATGTTCTGGCTTGAGGATTTTCCGATGCCCATAGAAATg            <  1:463017/150‑1 (MQ=255)
                                                                                                                                            cTTTTTTCTGATACCAGTGCAAATAAATGGCGATTTTATTGCGGTAATGCTCAGCTGTTTTTTCCGGCATGCTATTGAGCAGCAATAATGCATATTCTTGCCAGTCTAAATGTTCTGGCTTGAGGATTTTCCGATGCCCATAGAAATGAt          >  1:4159181/1‑150 (MQ=255)
                                                                                                                                            cTTTTTTCTGATACCAGTGCAAATAAATGGCGATTTTATTGCGGTAATGCTCAGCTGTTTTTTCCGGCATGCTATTGAGCAGCAATAATGCATATTCTTGCCAGTCTAAATGTTCTGGCTTGAGGATTTTCCGATGCCCATAGAAATGAt          >  1:1496383/1‑150 (MQ=255)
                                                                                                                                             ttttttCTGATACCAGTGCAAATAAATGGCGATTTTATTGCGGTAATGCTCAGCTGTTTTTTCCGGCATGCTATTGAGCAGCAATAATGCATATTCTTGCCAGTCTAAATGTTCTGGCTTGAGGATTTTCCGATGCCCATAGAAATGAtt         <  1:2877305/150‑1 (MQ=255)
                                                                                                                                             ttttttCTGATACCAGTGCAAATAAATGGCGATTTTATTGCGGTAATGCTCAGCTGTTTTTTCCGGCATGCTATTGAGCAGCAATAATGCATATTCTTGCCAGTCTAAATGTTCTGGCTTGAGGATTTTCCGATGCCCATAGAAATGAtt         >  2:3695449/1‑150 (MQ=255)
                                                                                                                                                tttCTTATACCAGTGCAAATAAATGGCGATTTTATTGCGGTAATGCTCAGCTGTTTTTTCCGGCATGCTATTGAGCAGCAATAATGCATATTCTTGCCAGTCTAAATGTTCTGGCTTGAGGATTTTCCGATGCCCATAGAAATGATTGTc      <  1:3826886/150‑1 (MQ=255)
                                                                                                                                                tttCTGATACCAGTGCAAATAAATGGCGATTTTATTGCGGTAATGCTCAGCTGTTTTTTCCGGCATGCTATTGAGCAGCAATAATGCATATTCTTGCCAGTCTAAATGTTCTGGCTTGAGGATTTTCCGATGCCCATAGAAATGATTGTc      >  2:291155/1‑150 (MQ=255)
                                                                                                                                                tttCTGATACCAGTGCAAATAAATGGCGATTTTATTGCGGTAATGCTCAGCTGTTTTTTCCGGCATGCTATTGAGCAGCAATAATGCATATTCTTGCCAGTCTAAATGTTCTGGCTTGAGGATTTTCCGATGCCCATAGAAATGATTGTc      >  1:1515636/1‑150 (MQ=255)
                                                                                                                                                    tGATACCAGTGCAAATAAATGGCGATTTTATTGCGGTAATGCTCAGCTGTTTTTTCCGGCATGCTATTGAGCAGCAATAATGCATATTCTTGCCAGTCTAAATGTTCTGGCTTGAGGATTTTCCGATGCCCATAGAAATGATTGTCATGt  >  1:866253/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
ATAGTTCTTCGATTTTGTGGGGCTAAATGATAATGCCCGACACCAGTAATCGTTATTGAGTAAAACTTTGCATATCCGCCGCCAGGAGGGGATATCTTTCGCGCCAATGTCCCCTTGCTGGGTTTGTGGAACCTCGATGCCTTTTTTCTGATACCAGTGCAAATAAATGGCGATTTTATTGCGGTAATGCTCAGCTGTTTTTTCCGGCATGCTATTGAGCAGCAATAATGCATATTCTTGCCAGTCTAAATGTTCTGGCTTGAGGATTTTCCGATGCCCATAGAAATGATTGTCATGT  >  CP009273/630871‑631168

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: