Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,050,591 C→A 9.6% 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,5910CA 9.6% 866.8 / 51.6 249A393A (GCG→GCTserAD‑3‑phosphoglycerate dehydrogenase
Reads supporting (aligned to +/- strand):  ref base C (104/121);  new base A (11/13);  total (115/134)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.79e-01

Only 100 of 252 total aligned reads displayed.

GCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACGTCTTCGTCGGCTTCAATATCAATAACCACATAACCCATCTGGGCGGAAGTTTGCAGATATTGCGCGGCGATGTTGACGCCCTGCTCGGCGAAGATTTTGTTCAGCGCAGTTAGCACGCCCGGACGGTTTTCGT  >  CP009273/3050443‑3050739
                                                                                                                                                    |                                                                                                                                                    
gCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAgcg                                                                                                                                                     >  2:405077/1‑150 (MQ=255)
  aaaGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCt                                                                                                                                                   >  1:2506537/1‑150 (MQ=255)
  aaaGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCt                                                                                                                                                   >  1:5203411/1‑150 (MQ=255)
   aaGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCtt                                                                                                                                                  <  1:1381315/150‑1 (MQ=255)
   aaGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCtt                                                                                                                                                  >  2:3952420/1‑150 (MQ=255)
    aGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCttt                                                                                                                                                 >  1:4866410/1‑150 (MQ=255)
      tAGAATGCGCAACGCGGCAACGGTGTGGGGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGACTGGTACCCGGACTAGCTTTCATTGCCTGCAGCGCttttt                                                                                                                                               >  1:4996939/1‑150 (MQ=255)
      tAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCttttt                                                                                                                                               <  2:3621118/150‑1 (MQ=255)
      tAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGAGCttttt                                                                                                                                               <  2:1520923/150‑1 (MQ=255)
       aGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTc                                                                                                                                              >  2:539627/1‑150 (MQ=255)
        gAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCg                                                                                                                                             <  1:3931806/150‑1 (MQ=255)
             cgcAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGAGCTTTTTCGGCAAc                                                                                                                                        <  1:3328713/150‑1 (MQ=255)
              gcAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACg                                                                                                                                       <  2:261123/150‑1 (MQ=255)
                 aCGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACGTCt                                                                                                                                    <  1:1918173/150‑1 (MQ=255)
                  cgcgGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACGTCtt                                                                                                                                   <  1:2391643/150‑1 (MQ=255)
                   gcgGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACGTCTtc                                                                                                                                  <  1:787872/150‑1 (MQ=255)
                   gcgGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACGTCTtc                                                                                                                                  <  1:4720063/150‑1 (MQ=255)
                     ggCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACGTCTtcgt                                                                                                                                <  1:1213364/150‑1 (MQ=255)
                      gCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGACTGGTACCCGGAATAGCTTTCATTGCCTGCAGAGCTTTTTCGGCAACGTCTtcgtc                                                                                                                               >  2:3383718/1‑150 (MQ=255)
                       cAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACGTCTtcgtcg                                                                                                                              >  1:1832241/1‑150 (MQ=255)
                       cAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACGTCTtcgtcg                                                                                                                              <  2:1066072/150‑1 (MQ=255)
                       cAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGAGCTTTTTCGGCAACGTCTtcgtcg                                                                                                                              >  2:4239478/1‑150 (MQ=255)
                          tggTGTGGAGCATGGATTACTAATCGGGCAAGACAGTGTCCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTTCAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACGTCTTCGTCGGCt                                                                                                                           <  2:2836402/149‑1 (MQ=255)
                          cGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACGTCTTCGTCGGCt                                                                                                                           <  1:1082158/150‑1 (MQ=255)
                            gtgtGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACGTCTTCGTCGGCTtc                                                                                                                         >  1:924491/1‑150 (MQ=255)
                             tgtgGAGAAGGGATAAAAAAACGGGCAAGTCATTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACGTCTTCGTCGGCTtca                                                                                                                        <  2:2012049/150‑1 (MQ=255)
                             tgtgGAGAAGGGATAAAAAAACGGGCAAGTCAGGGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACGTCTTCGTCGGCTtca                                                                                                                        <  2:1715874/150‑1 (MQ=255)
                                     aGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACGTCTTCGTCGGCTtcaatatcaat                                                                                                                <  1:633612/150‑1 (MQ=255)
                                      gggATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACGTCTTCGTCGGCTtcaatatcaata                                                                                                               >  2:377853/1‑150 (MQ=255)
                                      gggATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACGTCTTCGTCGGCTtcaatatcaata                                                                                                               <  2:4738149/150‑1 (MQ=255)
                                      gggATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACGTCTTCGTCGGCTtcaatatcaata                                                                                                               <  1:4341642/150‑1 (MQ=255)
                                      gggATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACGTCTTCGTCGGCTtcaatatcaata                                                                                                               <  2:4738266/150‑1 (MQ=255)
                                        gATAAAAAAACGGGCAAGTCATTGACCTGCACGTTTATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGAGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACGTCTTCGTCGGCTTCAATATCAATAAc                                                                                                             <  1:2790895/150‑1 (MQ=255)
                                         aTAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACGTCTTCGTCGGCTTCAATATCAATAAcc                                                                                                            <  1:3056545/150‑1 (MQ=255)
                                           aaaaaaaCGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACGTCTTCGTCGGCTTCAATATCAATAACcac                                                                                                          <  2:165450/150‑1 (MQ=255)
                                            aaaaaaCGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGAGCTTTTTCGGCAACGTCTTCGTCGGCTTCAATATCAATAACcaca                                                                                                         >  2:92496/1‑150 (MQ=255)
                                             aaaaaCGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACGTCTTCGTCGGCTTCACTATCAATAACcacaa                                                                                                        >  2:2752243/1‑149 (MQ=255)
                                             aaaaaCGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACGTCTTCGTCGGCTTCAATATCAATAACCACAt                                                                                                        >  2:3286300/1‑150 (MQ=255)
                                               aaaCGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACGTCTTCGTCGGCTTCAATATCAATAACCACATaa                                                                                                      <  1:5082265/150‑1 (MQ=255)
                                                aaCGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACGTCTTCGTCGGCTTCAATATCAATAACCACATAAc                                                                                                     <  2:227678/150‑1 (MQ=255)
                                                aaCGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACGTCTTCGTCGGCTTCAATATCAATAACCACATAAc                                                                                                     <  2:655850/150‑1 (MQ=255)
                                                   gggCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTGTCATTGCCTGCAGCGCTTTTTCGGCACCGTCTTCGTCGGCTTCAATCTCAATAACCACATAACCCa                                                                                                  >  2:2659486/1‑150 (MQ=255)
                                                     gCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGGAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACGTCTTCGTCGGCTTCAATATCAATAACCACATAACCCATc                                                                                                <  2:3065472/150‑1 (MQ=255)
                                                        aGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACGTCTTCGTCGGCTTCAATATCAATAACCACATAACCCATCTgg                                                                                             <  1:3602833/150‑1 (MQ=255)
                                                        aGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACGTCTTCGTCGGCTTCAATATCAATAACCACATAACCCATCTgg                                                                                             >  1:546002/1‑150 (MQ=255)
                                                        aGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCAAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACGTCTTCGTCGGCTTCAATATCAATAACCACATAACCCATCTgg                                                                                             <  1:3431961/150‑1 (MQ=255)
                                                           cAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCATCGTCTTCGTCGGCTTCAATATCAATAACCACATAACCCATCTGGGCg                                                                                          <  2:1576312/150‑1 (MQ=255)
                                                           cAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACGTCTTCGTCGGCTTCAATATCAATAACCACATAACCCATCTGGGCg                                                                                          >  2:2090445/1‑150 (MQ=255)
                                                            aGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACGTCTTCGTCGGCTTCAATATCAATAACCACATAACCCATCTGGGCgg                                                                                         <  2:5185961/150‑1 (MQ=255)
                                                             gTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACGTCTTCGTCGGCTTCAATATCAATAACCACATAACCCATCTGGGCGGa                                                                                        >  2:2564293/1‑150 (MQ=255)
                                                              tGACCTGCCCGTTGATTTTCAGTGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCTCTTTTTCGGCATCGTATTAGTCGGCTTCAATATCAATAACCACATAACCCATCTGGGCGGaa                                                                                       <  2:3243751/150‑1 (MQ=255)
                                                               gACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACGTCTTCGTCGGCTTCAATATCAATAACCACATAACCCATCTGGGCGGAAg                                                                                      >  1:3091966/1‑150 (MQ=255)
                                                               gACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACGTCTTCGTCGGCTTCAATATCAATAACCACATAACCCATCTGGGCGGAAg                                                                                      >  1:3188574/1‑150 (MQ=255)
                                                                  cTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACGTCTTCGTCGGCTTCAATATCAATAACCACATAACCCATCTGGGCGGAAGttt                                                                                   >  1:4041054/1‑150 (MQ=255)
                                                                   tGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACGTCTTCGTCGGCTTCAATATCAATAACCACATAACCCATCTGGGCGGAAGTTTg                                                                                  <  1:683293/150‑1 (MQ=255)
                                                                       cGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACGTCTTCGTCGGCTTCAATATCAATAACCACATAACCCATCTGGGCGGAAGTTTGCAGa                                                                              >  1:174037/1‑150 (MQ=255)
                                                                        gTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACGTCTTCGTCGGCTTCAATATCAATAACCACATAACCCATCTGGGCGGAAGTTTGCAGat                                                                             >  1:4196032/1‑150 (MQ=255)
                                                                        gTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACGTCTTCGTCGGCTTCAATATCAATAACCACATAACCCATCTGGGCGGAAGTTTGCAGat                                                                             >  2:1349690/1‑150 (MQ=255)
                                                                             ttttCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACGTCTTCGTCGGCTTCAATATCAATAACCACATAACCCATCTGGGCGGAAGTTTGCAGATATTgc                                                                        >  1:3451786/1‑150 (MQ=255)
                                                                              tttCAGAGAAGGGGAATTAGTACCTCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACGTCTTCGTCGGCTTCAATATCAATAACCACATAACCCATCTGGGCGGAAGTTTGCAGATATTgcg                                                                       <  1:1660370/150‑1 (MQ=255)
                                                                                   gagaAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGACTGCAGCGATTTTTCTGCAACGTCTTCGTCGGCTTCAATATCAATAACCACATAACCCATCTGGGCGGAAGTTTGCAGATATTGCgctgcg                                                                  >  1:1482425/1‑150 (MQ=255)
                                                                                    agaAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACGTCTTCGTCGGCTTCAATATCAATAACCACATAACCCATCTGGGCGGAAGTTTGCAGATATTGCGCGGCGa                                                                 >  1:997882/1‑150 (MQ=255)
                                                                                     gaAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACGTCTTCGTCGGCTTCAATATCAATAACCACATAACCCATCTGGGCGGAAGTTTGCAGATATTGCGCGGCGAt                                                                <  2:11519/150‑1 (MQ=255)
                                                                                       aGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACGTCTTCGTCGGCTTCAATATCAATAACCACATAACCCATCTGGGCGGAAGTTTGCAGATATTGCGCGGCGATGt                                                              >  1:4404832/1‑150 (MQ=255)
                                                                                        ggggAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGAGCTTTTTCGGCAACGTCTTCGTCGGCTTCAATATCAATAACCACATAACCCATCTGGGCGGAAGTTTGCAGATATTGCGCGGCGATGtt                                                             <  2:4050365/150‑1 (MQ=255)
                                                                                          ggAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACGTCTTCGTCGGCTTCAATATCAATAACCACATAACCCATCTGGGCGGAAGTTTGCAGATATTGCGCGGCGATGTTGa                                                           >  1:4692568/1‑150 (MQ=255)
                                                                                            aaTTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACGTCTTCGTCGGCTTCAATATCAATAACCACATAACCCATCTGGGCGGAAGTTTGCAGATATTGCGCGGCGATGTTGACg                                                         <  1:159481/150‑1 (MQ=255)
                                                                                              ttAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACGTCTTCGTCGGCTTCAATATCAATAACCACATAACCCATCTGGGCGGAAGTTTGCAGATATTGCGCGGCGATGTTGACGcc                                                       <  1:2080691/150‑1 (MQ=255)
                                                                                                aGTACAGCAGACGGGCTCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACGTCTTCGTCGGCTTCAATATCAATAACCACATAACCCATCTGGGCGGAAGTTTGCAGATATTGCGCGGCGATGTTGACGCCCt                                                     >  1:2298864/1‑150 (MQ=255)
                                                                                                  tACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACGTCTTCGTCGGCTTCAATATCAATAACCACATAACCCATCTGGGCGGAAGTTTGCAGATATTGCGCGGCGATGTTGACGCCCTGc                                                   >  2:2942254/1‑150 (MQ=255)
                                                                                                   aCAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACGTCTTCGTCGGCTTCAATATCAATAACCACATAACCCATCTGGGCGGAAGTTTGCAGATATTGCGCGGCGATGTTGACGCCCTGCt                                                  >  2:1814516/1‑150 (MQ=255)
                                                                                                   aCAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACGTCTTCGTCGGCTTCAATATCAATAACCACATAACCCATCTGGGCGGAAGTTTGCAGATATTGCGCGGCGATGTTGACGCCCTGCt                                                  <  2:3756278/150‑1 (MQ=255)
                                                                                                     agcagACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACGTCTTCGTCGGCTTCAATATCAATAACCACATAACCCATCTGGGCGGAAGTTTGCAGATATTGCGCGGCGATGTTGACGCCCTGCTCg                                                <  2:2170851/150‑1 (MQ=255)
                                                                                                            gggCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACGTCTTCGTCGGCTTCAATATCAATAACCACATAACCCATCTGGGCGGAAGTTTGCAGATATTGCGCGGCGATGTTGACGCCCTGCTCGGCGAAGa                                         <  2:409/150‑1 (MQ=255)
                                                                                                              gcgcGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACGTCTTCGTCGGCTTCAATATCAATAACCACATAACCCATCTGGGCGGAAGTTTGCAGATATTGCGCGGCGATGTTGACGCCCTGCTCGGCGAAGAtt                                       <  2:4917983/150‑1 (MQ=255)
                                                                                                              gcgcGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGAGCTTTTTCGGCAACGTCTTCGTCGGCTTCAATATCAATAACCACATAACCCATCTGGGCGGAAGTTTGCAGATATTGCGCGGCGATGTTGACGCCCTGCTCGGCGAAGAtt                                       <  1:668977/150‑1 (MQ=255)
                                                                                                                gcgAATGGTACCCGGAATAGCTGTCATTGCCTGCAGAGCTTTTTCGGCAACGTCTTCGTCGGCTTCAATATCAATAACCACATAACCCATCTGGGCGGAAGTTTGCAGATATTGCGCGGCGATGTTGACGCCCTGCTCGGCGAAGAtttt                                     <  2:2811059/150‑1 (MQ=255)
                                                                                                                 cgAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACGTCTTCGTCGGCTTCAATATCAATAACCACATAACCCATCTGGGCGGAAGTTTGCAGATATTGCGCGGCGATGTTGACGCCCTGCTCGGCGAAGATTTTg                                    <  2:4404832/150‑1 (MQ=255)
                                                                                                                 cgAATGGTACCCGGAATAGCTTTCATTGCCTGCAGAGCTTTTTCGGCAACGTCTTCGTCGGCTTCAATATCAATAACCACATAACCCATCTGGGCGGAAGTTTGCAGATATTGCGCGGCGATGTTGACGCCCTGCTCGGCGAAGATTTTg                                    <  2:756520/150‑1 (MQ=255)
                                                                                                                   aaTGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACGTCTTCGTCGGCTTCAATATCAATAACCACATAACCCATCTGGGCGGAAGTTTGCAGATATTGCGCGGCGATGTTGACGCCCTGCTCGGCGAAGATTTTGtt                                  <  2:846840/150‑1 (MQ=255)
                                                                                                                      ggTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACGTCTTCGTCGGCTTCAATATCAATAACCACATAACCCATCTGGGCGGAAGTTTGCAGATATTGCGCGGCGATGTTGACGCCCTGCTCGGCGAAGATTTTGTTCAg                               <  2:2015055/150‑1 (MQ=255)
                                                                                                                      ggTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACGTCTTCGTCGGCTTCAATATCAATAACCACATAACCCATCTGGGCGGAAGTTTGCAGATATTGCGCGGCGATGTTGACGCCCTGCTCGGCGAAGATTTTGTTCAg                               >  1:1943570/1‑150 (MQ=255)
                                                                                                                      ggTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACGTCTTCGTCGGCTTCAATATCAATAACCACATAACCCATCTGGGCGGAAGTTTGCAGATATTGCGCGGCGATGTTGACGCCCTGCTCGGCGAAGATTTTGTTCAg                               <  1:1545362/150‑1 (MQ=255)
                                                                                                                          cccGGAATAGCTTTCATTGCCTGCAGAGCTTTTTCGGCAACGTCTTCGTCGGCTTCAATATCAATAACCACATAACCCATCTGGGCGGAAGTTTGCAGATATTGCGCGGCGATGTTGACGCCCTGCTCGGCGAAGATTTTGTTCAGCGCa                           >  1:2588763/1‑150 (MQ=255)
                                                                                                                           ccGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACGTCTTCGTCGGCTTCAATATCAATAACCACATAACCCATCTGGGCGGAAGTTTGCAGATATTGCGCGGCGATGTTGACGCCCTGCTCGGCGAAGATTTTGTTCAGCGCAg                          <  2:1300448/150‑1 (MQ=255)
                                                                                                                            cGGAATAGCTTTCAGTGCCTGCAGCGCTTTTTCGGCAACGTCTTCGTCGGCTTCAATATCCAGAACCACATAACCCATCTGGGCGGAAGTTTGCAGATATTGCGCGGCGATGTTGACGCCCTGCTCGGCGAAGATTTTGTTCAGCGCAGt                         <  2:2303066/150‑1 (MQ=255)
                                                                                                                              gAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACGTCTTCGTCGGCTTCAATATCAATAACCACATAACCCATCTGGGCGGAAGTTTGCAGATATTGCGCGGCGATGTTGACGCCCTGCTCGGCGAAGATTTTGTTCAGCGCAGTTa                       <  1:5132842/150‑1 (MQ=255)
                                                                                                                               aaTAGCTTTCATTGCCTGCATCGCTTTTTCTGCAACGTCTTCGTCGTCTTCAATATCAATAACCACATAACCCATCTGGGCGGAAGTTTGCAGATATTGCGCGGCGATGTTGACGCCCTGCTCGGCGAAGATTTTGTTCAGCGCAGTTAg                      <  2:664676/150‑1 (MQ=255)
                                                                                                                               aaTAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACGTCTTCGTCGGCTTCAATATCAATAACCACATAACCCATCTGGGCGGAAGTTTGCAGATATTGCGCGGCGATGTTGACGCCCTGCTCGGCGAAGATTTTGTTCAGCGCAGTTAg                      >  1:4622977/1‑150 (MQ=255)
                                                                                                                               aaTAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACGTCTTCGTCGGCTTCAATATCAATAACCACATAACCCATCTGGGCGGAAGTTTGCAGATATTGCGCGGCGATGTTGACGCCCTGCTCGGCGAAGATTTTGTTCAGCGCAGTTAg                      <  1:4218706/150‑1 (MQ=255)
                                                                                                                                  aGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACGTCTTCGTCGGCTTCAATATCAATAACCACATAACCCATCTGGGCGGAAGTTTGCAGATATTGCGCGGCGATGTTGACGCCCTGCTCGGCGAAGATTTTGTTCAGCGCAGTTAGCAc                   <  1:3357144/150‑1 (MQ=255)
                                                                                                                                     tttCATTGCCTGCAGCGCTTTTTCGGCAACGTCTTCGTCGGCTTCAATATCAATAACCACATAACCCATCTGGGCGGAAGTTTGCAGATATTGCGCGGCGATGTTGACGCCCTGCTCGGCGAAGATTTTGTTCAGCGCAGTTAGCACGcc                <  2:4513975/150‑1 (MQ=255)
                                                                                                                                      ttCATTGCCTGCAGCGCTTTTTCGGCAACGTCTTCGTCGGCTTCAATATCAATAACCACATAACCCATCTGGGCGGACGTTTGCAGATATTGCGCGGCGATGTTGACGCCCTGCTCGGCGAAGATTTTGTTCAGCGCAGTTAGCACGccc               >  2:3206296/1‑150 (MQ=255)
                                                                                                                                      ttCATTGCCTGCAGCGCTTTTTCGGCAACGTCTTCGTCGGCTTCAATATCAATAACCACATAACCCATCTGGGCGGAAGTTTGCAGATATTGCGCGGCGATGTTGACGCCCTGCTCGGCGAAGATTTTGTTCAGCGCAGTTAGCACGccc               >  1:3260243/1‑150 (MQ=255)
                                                                                                                                       tCATTGCCTGCAGCGCTTTTTCGGCAACGTCTTCGTCGGCTTCAATATCAATAACCACATAACCCATCTGGGCGGAAGTTTGCAGATATTGCGCGGCGATGTTGACGCCCTGCTCGGCGAAGATTTTGTTCAGCGCAGTTAGCACGCCCg              <  1:1116074/150‑1 (MQ=255)
                                                                                                                                        cATTGCCTGCAGCGCTTTTTCGGCAACGTCTTCGTCGGCTTCAATATCAATAACCACATAACCCATCTGGGCGGAAGTTTGCAGATATTGCGCGGCGATGTTGACGCCCTGCTCGGCGAAGATTTTGTTCAGCGCAGTTAGCACGCCCgg             >  1:2110703/1‑150 (MQ=255)
                                                                                                                                         aTTGCCTGCAGAGCTTTTTCGGCAACGTCTTCGTCGGCTTCAATATCAATAACCACATAACCCATCTGGGCGGAAGTTTGCAGATATTGCGCGGCGATGTTGACGCCCTGCTCGGCGAAGATTTTGTTCAGCGCAGTTAGCACGCCCGGa            <  2:491106/150‑1 (MQ=255)
                                                                                                                                           tGCCTGCAGCGCTTTTTCGGCAACGTCTTCGTCGGCTTCAATATCAATAACCACATAACCCATCTGGGCGGAAGTTTGCAGATATTGCGCGGCGATGTTGACGCCCTGCTCGGCGAAGATTTTGTTCAGCGCAGTTAGCACGCCCGGACg          >  2:1183566/1‑150 (MQ=255)
                                                                                                                                           tGCCTGCAGCGCTTTTTCGGCAACGTCTTCGTCGGCTTCAATATCAATAACCACATAACCCATCTGGGCGGAAGTTTGCAGATATTGCGCGGCGATGTTGACGCCCTGCTCGGCGAAGATTTTGTTCAGCGCAGTTAGCACGCCCGGACg          <  1:377853/150‑1 (MQ=255)
                                                                                                                                                  aGCGCTTTTTCGGCAACGTCTTCGTCGGCTTCAATATCAATAACCACATAACCCATCTGGGCGGAAGTTTGCAGATATTGCGCGGCGATGTTGACGCCCTGCTCGGCGAAGATTTTGTTCAGCGCAGTTAGCACGCCCGGACGGTTTTCg   >  2:4210522/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
GCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACGTCTTCGTCGGCTTCAATATCAATAACCACATAACCCATCTGGGCGGAAGTTTGCAGATATTGCGCGGCGATGTTGACGCCCTGCTCGGCGAAGATTTTGTTCAGCGCAGTTAGCACGCCCGGACGGTTTTCGT  >  CP009273/3050443‑3050739

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