Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,050,591 C→A 13.8% A393A (GCG→GCT)  serA ← D‑3‑phosphoglycerate dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,050,5910CA13.8% 1199.6 / 127.6 369A393A (GCG→GCT) serAD‑3‑phosphoglycerate dehydrogenase
Reads supporting (aligned to +/- strand):  ref base C (168/150);  new base A (22/29);  total (190/179)
Fisher's exact test for biased strand distribution p-value = 2.28e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.82e-01

Only 100 of 378 total aligned reads displayed.

CGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACGTCTTCGTCGGCTTCAATATCAATAACCACATAACCCATCTGGGCGGAAGTTTGCAGATATTGCGCGGCGATGTTGACGCCCTGCTCGGCGAAGATTTTGTTCAGCGCAGTTAGCACGCCCGGACGGTTTTCGT  >  CP009273/3050442‑3050739
                                                                                                                                                     |                                                                                                                                                    
 gCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAgcg                                                                                                                                                     >  1:5559442/1‑150 (MQ=255)
   aaaGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCt                                                                                                                                                   >  2:317716/1‑150 (MQ=255)
   aaaGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCt                                                                                                                                                   <  1:2868066/150‑1 (MQ=255)
   aaaGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCt                                                                                                                                                   >  2:4819592/1‑150 (MQ=255)
    aaGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCtt                                                                                                                                                  <  2:3445816/150‑1 (MQ=255)
     aGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCttt                                                                                                                                                 <  2:2900190/150‑1 (MQ=255)
      gTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCtttt                                                                                                                                                >  2:3567732/1‑150 (MQ=255)
       tAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCttttt                                                                                                                                               <  2:954630/150‑1 (MQ=255)
            tGCGCAACGCGGCACCGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGACGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCa                                                                                                                                          <  2:4876535/150‑1 (MQ=255)
            tGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCa                                                                                                                                          <  2:2388337/150‑1 (MQ=255)
              cgcACCGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAAc                                                                                                                                        <  1:2153792/150‑1 (MQ=255)
              cgcAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAAc                                                                                                                                        <  2:5068268/150‑1 (MQ=255)
               gcAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACg                                                                                                                                       >  1:3563647/1‑150 (MQ=255)
                cAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACGt                                                                                                                                      >  2:1960535/1‑150 (MQ=255)
                cAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACGt                                                                                                                                      >  2:2799244/1‑150 (MQ=255)
                cAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGAGCTTTTTCGGCAACGt                                                                                                                                      <  2:17423/150‑1 (MQ=255)
                    gcgGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACGTCTtc                                                                                                                                  <  1:3753276/150‑1 (MQ=255)
                     cgGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACGTCTtcg                                                                                                                                 >  1:4969428/1‑150 (MQ=255)
                        cAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACGTCTtcgtcg                                                                                                                              <  1:2731997/150‑1 (MQ=255)
                         aaCGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCACCGTCTTCGTCgg                                                                                                                             >  2:1989715/1‑150 (MQ=255)
                           cGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACGTCTTCGTCGGCt                                                                                                                           <  2:3816071/150‑1 (MQ=255)
                             gtgtGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACGTCTTCGTCGGCTtc                                                                                                                         <  2:4775993/150‑1 (MQ=255)
                             gtgtGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACGTCTTCGTCGGCTtc                                                                                                                         <  2:1361224/150‑1 (MQ=255)
                                 ggAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACGTCTTCGTCGGCTtcaata                                                                                                                     >  1:808522/1‑150 (MQ=255)
                                 ggAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACGTCTTCGTCGGCTtcaata                                                                                                                     >  2:220899/1‑150 (MQ=255)
                                  gagaAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACGTCTTCGTCGGCTtcaatat                                                                                                                    >  1:2560842/1‑150 (MQ=255)
                                   agaAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACGTCTTCGTCGGCTtcaatatc                                                                                                                   >  2:1772342/1‑150 (MQ=255)
                                    gaAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACGTCTTCGTCGGCTtcaatatca                                                                                                                  <  2:4438707/150‑1 (MQ=255)
                                     aaGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGAGCTTTTTCGGCAACGTCTTCGTCGGCTtcaatatcaa                                                                                                                 >  1:2658542/1‑150 (MQ=255)
                                      aGGGATAAAAAAACGGGCAAGTCATTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGTCAACGTCTTCGTCGGCTtcaatatcaaa                                                                                                                >  2:3385237/1‑149 (MQ=255)
                                         gATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACGTCTTCGTCGGCTTCAATATCAATAAc                                                                                                             >  1:3981791/1‑150 (MQ=255)
                                          aTAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACGTCTTCGTCGGCTTCAATATCAATAAcc                                                                                                            <  1:309221/150‑1 (MQ=255)
                                          aTAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACGTCTTCGTCGGCTTCAATATCAATAAcc                                                                                                            >  2:783353/1‑150 (MQ=255)
                                            aaaaaaaCGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGAGCTTTTTCGGCAACGTCTTCGTCGGCTTCAATATCAATAACcac                                                                                                          <  2:5359611/150‑1 (MQ=255)
                                                aaaCGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACGTCTTCGTCGGCTTCAATATCAATAACCACATaa                                                                                                      >  1:1797336/1‑150 (MQ=255)
                                                 aaCGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACGTCTTCGTCGGCTTCAATATCAATAACCACATAAc                                                                                                     <  1:2018826/150‑1 (MQ=255)
                                                   cGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACGTCTTCGTCGGCTTCAATATCAATAACCACATAAccc                                                                                                   <  1:3611870/150‑1 (MQ=255)
                                                    ggtcAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACGTCTTCGTCGGCTTCAATATCAATAACCACATAACCCa                                                                                                  <  2:4370775/147‑1 (MQ=255)
                                                     ggCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACGTCTTCGTCGGCTTCAATATCAATAACCACATAACCCAt                                                                                                 >  1:5152041/1‑150 (MQ=255)
                                                         aGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACGTCTTCGTCGGCTTCAATATCAATAACCACATAACCCATCTgg                                                                                             >  2:2357907/1‑150 (MQ=255)
                                                             aGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGAGCTTTTTCGGCAACGTCTTCGTCGGCTTCAATATCAATAACCACATAACCCATCTGGGCgg                                                                                         >  2:4516057/1‑150 (MQ=255)
                                                              gTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACGTCTTCGTCGGCTTCAATATCAATAACCACATAACCCATCTGGGCGGa                                                                                        <  1:1722356/150‑1 (MQ=255)
                                                                     gCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACGTCTTCGTCGGCTTCAATATCAATAACCACATAACCCATCTGGGCGGAAGTTTGc                                                                                 <  1:2331054/150‑1 (MQ=255)
                                                                       ccGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACGTCTTCGTCGGCTTCAATATCAATAACCACATAACCCATCTGGGCGGAAGTTTGCAg                                                                               <  2:1116938/150‑1 (MQ=255)
                                                                        cGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACGTCTTCGTCGGCTTCAATATCAATAACCACATAACCCATCTGGGCGGAAGTTTGCAGa                                                                              <  1:2093874/150‑1 (MQ=255)
                                                                        cGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACGTCTTCGTCGGCTTCAATATCAATAACCACATAACCCATCTGGGCGGAAGTTTGCAGa                                                                              >  1:5376477/1‑150 (MQ=255)
                                                                            gATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGAGCTTTTTCGGCAACGTCTTCGTCGGCTTCAATATCAATAACCACATAACCCATCTGGGCGGAAGTTTGCAGATAtt                                                                          >  2:2859479/1‑150 (MQ=255)
                                                                             aTTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACGTCTTCGTCGGCTTCAATATCAATAACCACATAACCCATCTGGGCGGAAGTTTGCAGATATTg                                                                         >  1:1081604/1‑150 (MQ=255)
                                                                              ttttCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACGTCTTCGTCGGCTTCAATATCAATAACCACATAACCCATCTGGGCGGAAGTTTGCAGATATTgc                                                                        <  1:1575457/150‑1 (MQ=255)
                                                                                  cAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACGTCTTCGTCGGCTTCAATATCAATAACCACATAACCCATCTGGGCGGAAGTTTGCAGATATTGCgcgg                                                                    <  2:824658/150‑1 (MQ=255)
                                                                                   agagAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACGTCTTCGTCGGCTTCAATATCAATAACCACATAACCCATCTGGGCGGAAGTTTGCAGATATTGCgcggc                                                                   >  2:1010363/1‑150 (MQ=255)
                                                                                   agagAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACGTCTTCGTCGGCTTCAATATCAATAACCACATAACCCATCTGGGCGGAAGTTTGCAGATATTGCgcggc                                                                   >  2:1678572/1‑150 (MQ=255)
                                                                                   agagAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACGTCTTCGTCGGCTTCAATATCAATAACCACATAACCCATCTGGGCGGAAGTTTGCAGATATTGCgcggc                                                                   >  1:4889722/1‑150 (MQ=255)
                                                                                      gaAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGAGCTTTTTCGGCAACGTCTTCGTCGGCTTCAATATCAATAACCACATAACCCATCTGGGCGGAAGTTTGCAGATATTGCGCGGCGAt                                                                <  1:1920155/150‑1 (MQ=255)
                                                                                       aaGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACGTCTTCGTCGGCTTCAATATCAATAACCACATAACCCATCTGGGCGGAAGTTTGCAGATATTGCGCGGCGATg                                                               <  1:5397550/150‑1 (MQ=255)
                                                                                         ggggAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACGTCTTCGTCGGCTTCAATATCAATAACCACATAACCCATCTGGGCGGAAGTTTGCAGATATTGCGCGGCGATGtt                                                             <  1:2214101/150‑1 (MQ=255)
                                                                                         ggggAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACGTCTTCGTCGGCTTCAATATCAATAACCACATAACCCATCTGGGCGGAAGTTTGCAGATATTGCGCGGCGATGtt                                                             <  1:44320/150‑1 (MQ=255)
                                                                                           ggAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCAGTGCCTGCAGCGCTTTTTCGGCAACGTCTTCTTCGGCTTCAATATCAATAACCACATAACCCATCTGGGCGGAAGTTTGCAGATATTGCGCGGCGATGTTGa                                                           <  2:5447359/150‑1 (MQ=255)
                                                                                           ggAATTAGAACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACGTCTTCGTCGGCTTCAATATCAATAACCACATAACCCATCTGGGCGGAAGTTTGCAGATATTGCGCGGCGATGTTGa                                                           <  1:5291301/150‑1 (MQ=255)
                                                                                            gAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACGTCTTCGTCGGCTTCAATATCAATAACCACATAACCCATCTGGGCGGAAGTTTGCAGATATTGCGCGGCGATGTTGAc                                                          <  1:3805943/150‑1 (MQ=255)
                                                                                            gAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGAGCTTTTTCGGCAACGTCTTCGTCGGCTTCAATATCAATAACCACATAACCCATCTGGGCGGAAGTTTGCAGATATTGCGCGGCGATGTTGAc                                                          >  1:4804202/1‑150 (MQ=255)
                                                                                             aaTTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACGTCTTCGTCGGCTTCAATATCAATAACCACATAACCCATCTGGGCGGAAGTTTGCAGATATTGCGCGGCGATGTTGACg                                                         <  1:5586355/150‑1 (MQ=255)
                                                                                             aaTTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACGTCTTCGTCGGCTTCAATATCAATAACCACATAACCCATCTGGGCGGAAGTTTGCAGATATTGCGCGGCGATGTTGACg                                                         >  1:1053771/1‑150 (MQ=255)
                                                                                              aTTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCACCGTCTTCGTCGGCTTCCATATCAATAACCACATAACCCATCTGGGCGGAAGTTTGCAGATATTGCGCGGCGATGTTGACGc                                                        >  2:3883273/1‑150 (MQ=255)
                                                                                                 aGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACGTCTTCGTCGGCTTCAATATCAATAACCACATAACCCATCTGGGCGGAAGTTTGCAGATATTGCGCGGCGATGTTGACGCCCt                                                     >  1:4831843/1‑150 (MQ=255)
                                                                                                 aGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACGTCTTCGTCGGCTTCAATATCAATAACCACATAACCCATCTGGGCGGAAGTTTGCAGATATTGCGCGGCGATGTTGACGCCCt                                                     <  1:4063892/150‑1 (MQ=255)
                                                                                                  gTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACGTCTTCGTCGGCTTCAATATCAATAACCACATAACCCATCTGGGCGGAAGTTTGCAGATATTGCGCGGCGATGTCGACGCCCTg                                                    >  1:2430644/1‑150 (MQ=255)
                                                                                                     cagcagACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACGTCTTCGTCGGCTTCAATATCAATAACCACATACCCCATCTGGGCGGAAGTTTGCAGATATTGCGCGGCGATGTTGACGCCCTGCTc                                                 <  2:662159/150‑1 (MQ=255)
                                                                                                     cagcagACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACGTCTTCGTCGGCTTCAATATCAATAACCACATAACCCATCTGGGCGGAAGTTTGCAGATATTGCGCGGCGATGTTGACGCCCTGCTc                                                 <  2:2065201/150‑1 (MQ=255)
                                                                                                     cagcagACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACGTCTTCGTCGGCTTCAATATCAATAACCACATAACCCATCTGGGCGGAAGTTTGCAGATATTGCGCGGCGATGTTGACGCCCTGCTc                                                 <  1:1459221/150‑1 (MQ=255)
                                                                                                      agcagACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACGTCTTCGTCGGCTTCAATATCAATAACCACATAACCCATCTGGGCGGAAGTTTGCAGATATTGCGCGGCGATGTTGACGCCCTGCTCg                                                <  1:5098143/150‑1 (MQ=255)
                                                                                                              ggCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACGTCTTCGTCGGCTTCAATATCAATAACCACATAACCCATCTGGGCGGAAGTTTGCAGATATTGCGCGGCGATGTTGACGCCCTGCTCGGCGAAGAt                                        >  1:3239418/1‑150 (MQ=255)
                                                                                                              ggCGCGAATGGTACCCGGAATAGCTTTCAGTGCCTGCAGCGCTTTTTCGGCAACGTCTTCGTCGGCTTCAATATCAATAACCACATAACCCATCTGGGCGGAAGTTTGCAGATATTGCGCGGCGATGTTGACGCCCTGCTCGGCGAAGAt                                        <  1:3032305/150‑1 (MQ=255)
                                                                                                               gcgcGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGAGCTTTTTCGGCAACGTCTTCGTCGGCTTCAATATCAATAACCACATAACCCATCTGGGCGGAAGTTTGCAGATATTGCGCGGCGATGTTGACGCCCTGCTCGGCGAAGAtt                                       <  2:5602393/150‑1 (MQ=255)
                                                                                                                cgcgAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACGTCTTCGTCGGCTTCAATATCAATAACCACATAACCCATCTGGGCGGAAGTTTGCAGATATTGCGCGGCGATGTTGACGCCCTGCTCGGCGAAGAttt                                      <  1:3797764/150‑1 (MQ=255)
                                                                                                                  cgAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACGTCTTCGTCGGCTTCAATATCAATAACCACATAACCCATCTGGGCGGAAGTTTGCAGATATTGCGCGGCGATGTTGACGCCCTGCTCGGCGAAGATTTTg                                    <  2:2688114/150‑1 (MQ=255)
                                                                                                                   gAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACGTCTTCGTCGGCTTCAATATCAATAACCACATAACCCATCTGGGCGGAAGTTTGCAGATATTGCGCGGCGATGTTGACGCCCTGCTCGGCGAAGATTTTGt                                   <  2:5405406/150‑1 (MQ=255)
                                                                                                                    aaTGGTACCCGGAATAGCTTTCATTGCCTGCAGAGCTTTTTCGGCAACGTCTTCGTCGGCTTCAATATCCATAACCACATAACCCATCTGGGCGGAAGTTTGCAGATATTGCGCGGCGATGTTGACGCCCTGCTCGGCGAAGATTTTGtt                                  >  2:3045634/1‑150 (MQ=255)
                                                                                                                     aTGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACGTCTTCGTCGGCTTCAATATCAATAACCACATAACCCATCTGGGCGGAAGTTTGCAGATATTGCGCGGCGATGTTGACGCCCTGCTCGGCGAAGATTTTGTTc                                 <  1:4550314/150‑1 (MQ=255)
                                                                                                                           cccGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACGTCTTCGTCGGCTTCAATATCAATAACCACATAACCCATCTGGGCGGAAGTTTGCAGATATTGCGCGGCGATGTTGACGCCCTGCTCGGCGAAGATTTTGTTCAGCGCa                           <  1:5123788/150‑1 (MQ=255)
                                                                                                                           cccGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACGTCTTCGTCGGCTTCAATATCAATAACCACATAACCCATCTGGGCGGAAGTTTGCAGATATTGCGCGGCGATGTTGACGCCCTGCTCGGCGAAGATTTTGTTCAGCGCa                           <  1:667984/150‑1 (MQ=255)
                                                                                                                           cccGGAATAGCTTTCATTGCCTGCAGAGCTTTTTCGGCAACGTCTTCGTCGGCTTCAATATCAATAACCACATAACCCATCTGGGCGGAAGTTTGCAGATATTGCGCGGCGATGTTGACGCCCTGCTCGGCGAAGATTTTGTTCAGCGCa                           <  2:2531818/150‑1 (MQ=255)
                                                                                                                             cGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACGTCTTCGTCGGCTTCAATATCAATAACCACATAACCCATCTGGGCGGAAGTTTGCAGATATTGCGCGGCGATGTTGACGCCCTGCTCGGCGAAGATTTTGTTCAGCGCAGt                         >  2:5672744/1‑150 (MQ=255)
                                                                                                                              ggAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACGTCTTCGTCGGCTTCAATATCAATAACCACATAACCCATCTGGGCGGAAGTTTGCAGATATTGCGCGGCGATGTTGACGCCCTGCTCGGCGAAGATTTTGTTCAGCGCAGtt                        <  1:4326209/150‑1 (MQ=255)
                                                                                                                               gAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACGTCTTCGTCGGCTTCAATATCAATAACCACATAACCCATCTGGGCGGAAGTTTGCAGATATTGCGCGGCGATGTTGACGCCCTGCTCGGCGAAGATTTTGTTCAGCGCAGTTa                       <  2:2653728/150‑1 (MQ=255)
                                                                                                                                  tAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACGTCTTCGTCGGCTTCAATATCAATAACCACATAACCCATCTGGGCGGAAGTTTGCAGATATTGCGCGGCGATGTTGACGCCCTGCTCGGCGAAGATTTTGTTCAGCGCAGTTAGCa                    >  2:407567/1‑150 (MQ=255)
                                                                                                                                  tAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACGTCTTCGTCGGCTTCAATATCAATAACCACATAACCCATCTGGGCGGAAGTTTGCAGATATTGCGCGGCGATGTTGACGCCCTGCTCGGCGAAGATTTTGTTCAGCGCAGTTAGCa                    >  1:1863320/1‑150 (MQ=255)
                                                                                                                                   aGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACGTCTTCGTCGGCTTCAATATCAATAACCACATAACCCATCTGGGCGGAAGTTTGCAGATATTGCGCGGCGATGTTGACGCCCTGCTCGGCGAAGATTTTGTTCAGCGCAGTTAGCAc                   >  2:3948882/1‑150 (MQ=255)
                                                                                                                                   aGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACGTCTTCGTCGGCTTCAATATCAATAACCACATAACCCATCTGGGCGGAAGTTTGCAGATATTGCGCGGCGATGTTGACGCCCTGCTCGGCGAAGATTTTGTTCAGCGCAGTTAGCAc                   >  2:418555/1‑150 (MQ=255)
                                                                                                                                    gCTTTCATTGCCTGCAGCGCTTTTTCGGCAACGTCTTCGTCGGCTTCAATATCAATAACCACATAACCCATCTGGGCGGAAGTTTGCAGATATTGCGCGGCGATGTTGACGCCCTGCTCGGCGAAGATTTTGTTCAGCGCAGTTAGCACg                  >  2:1255741/1‑150 (MQ=255)
                                                                                                                                          aTTGCCTGCAGCGCTTTTTCGGCAACGTCTTCGTCGGCTTCAATATCAATAACCACATAACCCATCTGGGCGGAAGTTTGCAGATATTGCGCGGCGATGTTGACGCCCTGCTCGGCGAAGATTTTGTTCAGCGCAGTTAGCACGCCCGGa            <  1:1940542/150‑1 (MQ=255)
                                                                                                                                            tGCCTGCAGAGCTTTTTCGGCAACGTCTTCGTCGGCTTCAATATCAATAACCACATAACCCATCTGGGCGGAAGTTTGCAGATATTGCGCGGCGATGTTGACGCCCTGCTCGGCGAAGATTTTGTTCAGCGCAGTTAGCACGCCAGGACg          <  1:2859479/150‑1 (MQ=255)
                                                                                                                                             gCCTGCAGCGCTTTTTCGGCAACGTCTTCGTCGGCTTCAATATCAATAACCACATAACCCATCTGGGCGGAAGTTTGCAGATATTGCGCGGCGATGTTGACGCCCTGCTCGGCGAAGATTTTGTTCAGCGCAGTTAGCACGCCCGGACgg         >  2:1847864/1‑150 (MQ=255)
                                                                                                                                             gCCTGCAGCGC‑TTTTCGGCAACGTCTTCGTCGGCTTCAATATCAAAAACAACATAACCCATCTGGGCGGAAGTTTGCAGATATTGCGCGGCGATGTTGACGCCCTGCTCGGCGAAGATTTTGTTCAGCGCAGTAATCACGCCAGGACgga        >  2:748629/1‑149 (MQ=255)
                                                                                                                                               cTGCAGCGCTTTTTCGGCAACGTCTTCGTCGGCTTCAATATCAATAACCACATAACCCATCTGGGCGGAAGTTTGCAGATATTGCGCGGCGATGTTGACGCCCTGCTCGGCGAAGATTTTGTTCAGCGCAGTTAGCACGCCCGGACGGtt       >  1:3393440/1‑150 (MQ=255)
                                                                                                                                                 gCAGCGCTTTTTCGGCAACGTCTTCGTCGGCTTCAATATCAATAACCACATAACCCATCTGGGCGGAAGTTTGCAGATATTGCGCGGCGATGTTGACGCCCTGCTCGGCGAAGATTTTGTTCAGCGCAGTTAGCACGCCCGGACGGtttt     <  1:1214281/150‑1 (MQ=255)
                                                                                                                                                  cAGCGCTTTTTCGGCAACGTCTTCGTCGGCTTCAATATCAATAACCACATAACCCATCTGGGCGGAAGTTTGCAGATATTGCGCGGCGATGTTGACGCCCTGCTCGGCGAAGATTTTGTTCAGCGCAGTTAGCACGCCCGGACGGTTTTc    <  2:3283145/150‑1 (MQ=255)
                                                                                                                                                   aGCGCTTTTTCGGCAACGTCTTCGTCGGCTTCAATATCAATAACCACATAACCCATCTGGGCGGAAGTTTGCAGATATTGCGCGGCGATGTTGACGCCCTGCTCGGCGAAGATTTTGTTCAGCGCAGTTAGCACGCCCGGACGGTTTTCg   >  2:4981680/1‑150 (MQ=255)
                                                                                                                                                   aGCGCTTTTTCGGCAACGTCTTCGTCGGCTTCAATATCAATAACCACATAACCCAACTGGGCGGAAGTTTGCAGATATTGCGCGGCGATGTTGACTCCCTGCTCGGCGAAGATTTTGTTCAGCGCAGATAGCACGCCCGGACGGTTTTCg   >  2:1821807/1‑150 (MQ=255)
                                                                                                                                                    gcgcTTTTTCGGCAACGTCTTCGTCGGCTTCAATATCAATAACCACATAACCCATCTGGGCGGAAGTTTGCAGATATTGCGCGGCGATGTTGACGCCCTGCTCGGCGAAGATTTTGTTCAGCGCAGTTAGCACGCCCGGACGGTTTTCGt  <  1:1300003/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
CGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACGTCTTCGTCGGCTTCAATATCAATAACCACATAACCCATCTGGGCGGAAGTTTGCAGATATTGCGCGGCGATGTTGACGCCCTGCTCGGCGAAGATTTTGTTCAGCGCAGTTAGCACGCCCGGACGGTTTTCGT  >  CP009273/3050442‑3050739

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: ‑