Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,086,004 C→A 86.7% A42A (GCC→GCA hisC → histidinol‑phosphate aminotransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,086,0040CA86.7% 644.8 / 97.0 267A42A (GCC→GCAhisChistidinol‑phosphate aminotransferase
Reads supporting (aligned to +/- strand):  ref base C (19/17);  new base A (97/134);  total (116/151)
Fisher's exact test for biased strand distribution p-value = 2.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 267 total aligned reads displayed.

AACGCCCTTAAGGAGCAAGCATGAGCACCGTGACTATTACCGATTTAGCGCGTGAAAACGTCCGCAACCTGACGCCGTATCAGTCGGCGCGTCGTCTGGGCGGTAACGGCGATGTCTGGCTGAACGCCAACGAATACCCCACTGCCGTGGAGTTTCAGCTTACTCAGCAAACGCTCAACCGCTACCCGGAATGCCAGCCGAAAGCGGTGATTGAAAATTACGCGCAATATGCAGGCGTAAAACCGGAGCAGGTGCTGGTCAGCCGTGGCGCGGACGAAGGTATTGAACTGC  >  CP009273/2085859‑2086149
                                                                                                                                                 |                                                                                                                                                 
 aCGCCCTTAAGGAGCAAGCATGAGCACCGTGACTATTACCGATTTAGCGCGTGAAAACGTCCGCAACCTGACGCCGTATCAGTCGGCGCGTCGTCTGGGCGGTAACGGCGATGTCTGGCTGAACGCCAACGAATACCCCACTGCAGTgg                                                                                                                                               <  2:3843152/149‑1 (MQ=255)
  cGCCCTTAAGGAGCAAGCATGAGCACCGTGACTATTACCGATTTAGCGCGTGAAAACGTCCGCAACCTGACGCCGTATCAGTCGGCGCGTCGTCTGGGCGGTAACGGCGATGTCTGGCTGAACGCCAACGAATACCCCACTGCAGTGGa                                                                                                                                              <  1:4644121/149‑1 (MQ=255)
   gCCCTTAAGGAGCAAGCATGAGCACCGTGACTATTACCGATTTAGCGCGTGAAAACGTCCGCAACCTGACGCCGTATCAGTCGGCGCGTCGTCTGGGCGGTAACGGCGATGTCTGGCTGAACGCCAACGAATACCCCACTGCAGTGGAg                                                                                                                                             >  2:83174/1‑149 (MQ=255)
   gCCCTTAAGGAGCAAGCATGAGCACCGTGACTATTACCGATTTAGCGCGTGAAAACGTCCGCAACCTGACGCCGTATCAGTCGGCGCGTCGTCTGGGCGGTAACGGCGATGTCTGGCTGAACGCCAACGAATACCCCACTGCAGTGGAg                                                                                                                                             >  2:4664937/1‑149 (MQ=255)
   gCCCTTAAGGAGCAAGCATGAGCACCGTGACTATTACCGATTTAGCGCGTGAAAACGTCCGCAACCTGACGCCGTATCAGTCGGCGCGTCGTCTGGGCGGTAACGGCGATGTCTGGCTGAACGCCAACGAATACCCCACTGCAGTGGAg                                                                                                                                             >  2:470592/1‑149 (MQ=255)
    cccTTAAGGAGCAAGCATGAGCACCGTGACTATTACCGATTTAGCGCGTGAAAACGTCCGCAACCTGACGCCGTATCAGTCGGCGCGTCGTCTGGGCGGTAACGGCGATGTCTGGCTGAACGCCAACGAATACCCCACTGCCGTGGAGt                                                                                                                                            >  2:3615225/1‑149 (MQ=255)
    cccTTAAGGAGCAAGCATGAGCACCGTGACTATTACCGATTTAGCGCGTGAAAACGTCCGCAACCTGACGCCGTATCAGTCGGCGCGTCGTCTGGGCGGTAACGGCGATGTCTGGCTGAACGCCAACGAATACCCCACTGCCGTGGAGt                                                                                                                                            >  2:739519/1‑149 (MQ=255)
     ccTTAAGGAGCAAGCATGAGCACCGTGACTATTACCGATTTAGCGCGTGAAAACGTCCGCAACCTGACGCCGTATCAGTCGGCGCGTCGTCTGGGCGGTAACGGCGATGTCTGGCTGAACGCCAACGAATACCCCACTGCCGTGGAGtt                                                                                                                                           <  2:4468274/149‑1 (MQ=255)
     ccTTAAGGAGCAAGCATGAGCACCGTGACTATTACCGATTTAGCGCGTGAAAACGTCCGCAACCTGACGCCGTATCAGTCGGCGCGTCGTCTGGGCGGTAACGGCGATGTCTGGCTGAACGCCAACGAATACCCCACTGCAGTGGAGtt                                                                                                                                           >  2:4532755/1‑149 (MQ=255)
       ttAAGGAGCAAGCATGAGCACCGTGACTATTACCGATTTAGCGCGTGAAAACGTCCGCAACCTGACGCCGTATCAGTCGGCGCGTCGTCTGGGCGGTAACGGCGATGTCTGGCTGAACGCCAACGAATACCCCACTGTCGTGGAGTTTc                                                                                                                                         <  1:1006269/149‑1 (MQ=255)
       ttAAGGAGCAAGCATGAGCACCGTGACTATTACCGATTTAGCGCGTGAAAACGTCCGCAACCTGACGCCGTATCAGTCGGCGCGTCGTCTGGGCGGTAACGGCGATGTCTGGCTGAACGCCAACGAATACCCCACTGCAGTGGAGTTTc                                                                                                                                         <  1:2921459/149‑1 (MQ=255)
        tAAGGAGCAAGCATGAGCACCGTGACTATTACCGATTTAGCGCGTGAAAACGTCCGCAACCTGACGCCGTATCAGTCGGCGCGTCGTCTGGGCGGTAACGGCGATGTCTGGCTGAACGCCAACGAATACCCCACTGCAGTGGAGTTTCa                                                                                                                                        <  2:1214785/149‑1 (MQ=255)
         aaGGAGCAAGCATGAGCACCGTGACTATTACCGATTTAGCGCGTGAAAACGTCCGCAACCTGACGCCGTATCAGTCGGCGCGTCGTCTGGGCGGTAACGGCGATGTCTGGCTGAACGCCAACGAATACCCCACTGCCGTGGAGTTTCAg                                                                                                                                       >  2:1475773/1‑149 (MQ=255)
         aaGGAGCAAGCATGAGCACCGTGACTATTACCGATTTAGCGCGTGAAAACGTCCGCAACCTGACGCCGTATCAGTCGGCGCGTCGTCTGGGCGGTAACGGCGATGTCTGGCTGAACGCCAACGAATACCCCACTGCAGTGGAGTTTCAg                                                                                                                                       <  2:4420091/149‑1 (MQ=255)
          aGGAGCAAGCATGAGCACCGTGACTATTACCGATTTAGCGCGTGAAAACGTCCGCAACCTGACGCCGTATCAGTCGGCGCGTCGTCTGGGCGGTAACGGCGATGTCTGGCTGAACGCCAACGAATACCCCACTGCAGTGGAGTTTCAGc                                                                                                                                      >  2:941759/1‑149 (MQ=255)
           ggAGCAAGCATGAGCACCGTGACTATTACCGATTTAGCGCGTGAAAACGTCCGCAACCTGACGCCGTATCAGTCGGCGCGTCGTCTGGGCGGTAACGGCGATGTCTGGCTGAACGCCAACGAATACCCCACTGCCGTGGAGTTTCAGCt                                                                                                                                     >  2:1236972/1‑149 (MQ=255)
              gcaagcaTGAGCACCGTGACTATTACCGATTTAGCGCGTGAAAACGTCCGCAACCTGACGCCGTATCAGTCGGCGCGTCGTCTGGGCGGTAACGGCGATGTCTGGCTGAACGCCAACGAATACCCCACTGCAGTGGAGTTTCAGCTTAc                                                                                                                                  >  2:3378879/1‑149 (MQ=255)
                aagcaTGAGCACCGTGACTATTACCGATTTAGCGCGTGAAAACGTCCGCAACCTGACGCCGTATCAGTCGGCGCGTCGTCTGGGCGGTAACGGCGATGTCTGGCTGAACGCCAACGAATACCCCACTGCAGTGGAGTTTCAGCTTACTc                                                                                                                                >  1:451786/1‑149 (MQ=255)
                 agcgTGAGCACCGTGACTATTACCGATTTAGCGCGTGAAAACGTCCGCAACCTGACGCCGTATCAGTCGGCGCGTCGTCTGGGCGGTAACGGCGATGTCTGGCTGAACGCCAACGAATACCCCACTGCAGTGGAGTTTCAGCTTACTCa                                                                                                                               >  1:1931129/1‑149 (MQ=255)
                 agcaTGAGCACCGTGACTATTACCGCTTTAGCGCGTGAAAACGTCCGCAACCTGACGCCGTATCAGTCGGCGCGTCGTCTGGGCGGTAACGGCGATGTCTGGCTGAACGCCAACGAATACCCCACTGCAGTGGAGTTTCAGCTTACTCa                                                                                                                               <  2:785519/149‑1 (MQ=255)
                 agcaTGAGCACCGTGACTATTACCGATTTAGCGCGTGAAAACGTCCGCAACCTGACGCCGTATCAGTCGGCGCGTCGTCTGGGCGGTAACGGCGATGTCTGGCTGAACGCCAACGAATACCCCACTGCCGTGGAGTTTCAGCTTACTCa                                                                                                                               <  2:2106994/149‑1 (MQ=255)
                   caTGAGCACGGTGACTATTACCGATTTAGCGCGTGAAAACGTCCGCAACCTGACGCCGTATCAGTCGGCGCGTCGTCTGGGCGGTAACGGCGATGTCTGGCTGAACGCCAACGAATACCCCACTGCAGTGGAGTTTCAGCTTACTCAGc                                                                                                                             <  1:2925687/149‑1 (MQ=255)
                    aTGAGCACCGTGACTATTACCGATTTAGCGCGTGAAAACGTCCGCAACCTGACGCCGTATCAGTCGGCGCGTCGTCTGGGCGGTAACGGCGATGTCTGGCTGAACGCCAACGAATACCCCACTGCCGTGGAGTTTCAGCTTACTCAGCa                                                                                                                            >  1:38963/1‑149 (MQ=255)
                    aTGAGCACCGTGACTATTACCGATTTAGCGCGTGAAAACGTCCGCAACCTGACGCCGTATCAGTCGGCGCGTCGTCTGGGCGGTAACGGCGATGTCTGGCTGAACGCCAACGAATACCCCACTGCAGTGGAGTTTCAGCTTACTCAGCa                                                                                                                            >  1:2420530/1‑149 (MQ=255)
                      gAGCACCGTGACTATTACCGATTTAGCGCGTGAAAACGTCCGCAACCTGACGCCGTATCAGTCGGCGCGTCGTCTGGGCGGTAACGGCGATGTCTGGCTGAACGCCAACGAATACCCCACTGCCGTGGAGTTTCAGCTTACTCAGCaaa                                                                                                                          >  2:1769405/1‑149 (MQ=255)
                       aGCACCGTGACTATTACCGATTTAGCGCGTGAAAACGTCCGCAACCTGACGCCGTATCAGTCGGCGCGTCGTCTGGGCGGTAACGGCGATGTCTGGCTGAACGCCAACGAATACCCCACTGCAGTGGAGTTTCAGCTTACTCAGCAAAc                                                                                                                         <  2:4297153/149‑1 (MQ=255)
                        gCACCGTGACTATTACCGATTTAGCGCGTGAAAACGTCCGCAACCTGACGCCGTATCAGTCGGCGCGTCGTCTGGGCGGTAACGGCGATGTCTGGCTGAACGCCAACGAATACCCCACTGCAGTGGAGTTTCAGCTTACTCAGCAAACg                                                                                                                        >  1:1538603/1‑149 (MQ=255)
                          aCCGTGACTATTACCGATTTAGCGCGTGAAAACGTCCGCAACCTGACGCCGTATCAGTCGGCGCGTCGTCTGGGCGGTAACGGCGATGTCTGGCTGAACGCCAACGAATACCCCACTGCCGTGGAGTTTCAGCTTACTCAGCAAACGCt                                                                                                                      >  1:1484591/1‑149 (MQ=255)
                          aCCGTGACTATTACCGATTTAGCGCGTGAAAACGTCCGCAACCTGACGCCGTATCAGTCGGCGCGTCGTCTGGGCGGTAACGGCGATGTCTGGCTGAACGCCAACGAATACCCCACTGCAGTGGAGTTTCAGCTTACTCAGCAAACGCt                                                                                                                      <  2:620261/149‑1 (MQ=255)
                          aCCGTGACTATTACCGATTTAGCGCGTGAAAACGTCCGCAACCTGACGCCGTATCAGTCGGCGCGTCGTCTGGGCGGTAACGGCGATGTCTGGCTGAACGCCAACGAATACCCCACTGCAGTGGAGTTTCAGCTTACTCAGCAAACGCt                                                                                                                      >  2:4752554/1‑149 (MQ=255)
                           ccGTGACTATTACCGATTTAGCGCGTGAAAACGTCCGCAACCTGACGCCGTATCAGTCGGCGCGTCGTCTGGGCGGTAACGGCGATGTCTGGCTGAACGCCAACGAATACCCCACTGCAGTGGAGTTTCAGCTTACTCAGCAAACGCTc                                                                                                                     <  1:244678/149‑1 (MQ=255)
                               gACTATTACCGATTTAGCGCGTGAAAACGTCCGCAACCTGACGCCGTATCAGTCGGCGCGTCGTCTGGGCGGTAACGGCGATGTCTGGCTGAACGCCAACGAATACCCCACTGCAGTGGAGTTTCAGCTTACTCAGCAATCGCTCAAcc                                                                                                                 <  2:2354819/149‑1 (MQ=255)
                               gACTATTACCGATTTAGCGCGTGAAAACGTCCGCAACCTGACGCCGTATCAGTCGGCGCGTCGTCTGGGCGGTAACGGCGATGTCTGGCTGAACGCCAACGAATACCCCACTGCAGTGGAGTTTCAGCTTACTCAGCAAACGCTCAAcc                                                                                                                 >  2:541691/1‑149 (MQ=255)
                                 cTATTACCGATTTAGCGCGTGAAAACGTCCGCAACCTGACGCCGTATCAGTCGGCGCGTCGTCTGGGCGGTAACGGCGATGTCTGGCTGAACGCCAACGAATACCCCACTGCAGTGGAGTTTCAGCTTACTCAGCAAACGCTCAACCGc                                                                                                               <  2:4240175/149‑1 (MQ=255)
                                      aCCGATTTAGCGCGTGAAAACGTCCGCAACCTGACGCCGTATCAGTCGGCGCGTCGTCTGGGCGGTAACGGCGATGTCTGGCTGAACGCCAACGAATACCCCACTGCAGTGGAGTTTCAGCTTACTCAGCAAACGCTCAACCGCTAccc                                                                                                          >  1:3334506/1‑149 (MQ=255)
                                          aTTTAGCGCGTGAAAACGTCCGCAACCTGACGCCGTATCAGTCGGCGCGTCGTCTGGGCGGTAACGGCGATGTCTGGCTGAACGCCAACGAATACCCCACTGCAGTGGAGTTTCAGCTTACTCAGCAAACGCTCAACCGCTACCCGGaa                                                                                                      <  2:4075420/149‑1 (MQ=255)
                                           tttAGCGCGTGAAAACGTCCGCAACCTGACGCCGTATCAGTCGGCGCGTCGTCTGGGCGGTAACGGCGATGTCTGGCTGAACGCCAACGAATACCCCACTGCAGTGGAGTTTCAGCTTACTCAGCAAACGCTCAACCGCTACCCGGAAt                                                                                                     >  2:532693/1‑149 (MQ=255)
                                               gcgcGTGAAAACGTCCGCAACCTGACGCCGTATCAGTCGGCGCGTCGTCTGGGCGGTAACGGCGATGTCTGGCTGAACGCCAACGAATACCCCACTGCAGTGGAGTTTCAGCTTACTCAGCAAACGCTCAACCGCTACCCGGAATGCCa                                                                                                 <  2:1558970/149‑1 (MQ=255)
                                               gcgcGTGAAAACGTCCGCAACCTGACGCCGTATCAGTCGGCGCGTCGTCTGGGCGGTAACGGCGATGTCTGGCTGAACGCCAACGAATACCCCACTGCAGTGGAGTTTCAGCTTACTCAGCAAACGCTCAACCGCTACCCGGAATGCCa                                                                                                 <  2:2815772/149‑1 (MQ=255)
                                                  cgTGAAAACGTCCGCAACCTGACGCCGTATCAGTCGGCGCGTCGTCTGGGCGGTAACGGCGATGTCTGGCTGAACGCCAACGAATACCCCACTGCAGTGGAGTTTCAGCTTACTCAGCAAACGCTCAACCGCTACCCGGAATGCCAGcc                                                                                              <  1:2876455/149‑1 (MQ=255)
                                                     gAAAACGTCCGCAACCTGACGCCGTATCAGTCGGCGCGTCGTCTGGGCGGTAACGGCGATGTCTGGCTGAACGCCAACGAATACCCCACTGCAGTGGAGTTTCAGCTTACTCAGCAAACGCTCAACCGCTACCCGGAATGCCAGCCGaa                                                                                           >  2:3050573/1‑149 (MQ=255)
                                                       aaaCGTCCGCAACCTGACGCCGTATCAGTCGGCGCGTCGTCTGGGCGGTAACGGCGATGTCTGGCTGAACGCCAACGAATACCCCACTGCAGTGGAGTTTCAGCTTACTCAGCAAACGCTCAACCGCTACCCGGAATGCCAGCCGAAAg                                                                                         >  1:2095674/1‑149 (MQ=255)
                                                       aaaCGTCCGCAACCTGACGCCGTATCAGTCGGCGCGTCGTCTGGGCGGTAACGGCGATGTCTGGCTGAACGCCAACGAATACCCCACTGCAGTGGAGTTTCAGCTTACTCAGCAAACGCTCAACCGCTACCCGGAATGCCAGCCGAAAg                                                                                         >  2:4204722/1‑149 (MQ=255)
                                                        aaCGTCCGCAACCTGACGCCGTATCAGTCGGCGCGTCGTCTGGGCGGTAACGGCGATGTCTGGCTGAACGCCAACGAATACCCCACTGCAGTGGAGTTTCAGCTTACTCAGCAAACGCTCAACCGCTACCCGGAATGCCAGCCGAAAGc                                                                                        <  1:4270312/149‑1 (MQ=255)
                                                        aaCGTCCGCAACCTGACGCCGTATCAGTCGGCGCGTCGTCTGGGCGGTAACGGCGATGTCTGGCTGAACGCCAACGAATACCCCACTGCAGTGGAGTTTCAGCTTACTCAGCAAACGCTCAACCGCTACCCGGAATGCCAGCCGAAAGc                                                                                        <  1:3585351/149‑1 (MQ=255)
                                                           gTCCGCAACCTGACGCCGTATCAGTCGGCGCGTCGTCTGGGCGGTAACGGCGATGTCTGGCTGAACGCCAACGAATACCCCACTGCCGTGGAGTTTCAGCTTACTCAGCAAACGCTCAACCGCTACCCGGAATGCCAGCCGAAAGCGGt                                                                                     <  2:2870542/149‑1 (MQ=255)
                                                           gTCCGCAACCTGACGCCGTATCAGTCGGCGCGTCGTCTGGGCGGTAACGGCGATGTCTGGCTGAACGCCAACGAATACCCCACTGCAGTGGAGTTTCAGCTTACTCAGCAAACGCTCAACCGCTACCCGGAATGCCAGCCGAAAGCGGt                                                                                     >  1:1139312/1‑149 (MQ=255)
                                                           gTCCGCAACCTGACGCCGTATCAGTCGGCGCGTCGTCTGGGCGGTAACGGCGATGTCTGGCTGAACGCCAACGAATACCCCACTGCAGTGGAGTTTCAGCTTACTCAGCAAACGCTCAACCGCTACCCGGAATGCCAGCCGAAAGCGGt                                                                                     >  1:1877064/1‑149 (MQ=255)
                                                             ccGCAACCTGACGCCGTATCAGTCGGCGCGTCGTCTGGGCGGTAACGGCGATGTCTGGCTGAACGCCAACGAATACCCCACTGCAGTGGAGTTTCAGCTTACTCAGCAAACGCTCAACCGCTACCCGGAATGCCAGCCGAAAGCGGTGa                                                                                   <  1:162300/149‑1 (MQ=255)
                                                              cGCAACCTGACGCCGTATCAGTCGGCGCGTCGTCTGGGCGGTAACGGCGATGTCTGGCTGAACGCCAACGAATACCCCACTGCAGTGGAGTTTCAGCTTACTCAGCAAACGCTCAACCGCTACCCGGAATGCCAGCCGAAAGCGGTGAt                                                                                  <  1:3394746/149‑1 (MQ=255)
                                                                cAACCTGACGCCGTATCAGTCGGCGCGTCGTCTGGGCGGTAACGGCGATGTCTGGCTGAACGCCAACGAATACCCCACTGCAGTGGAGTTTCAGCTTACTCAGCAAACGCTCAACCGCTACCCGGAATGCCAGCCGAAAGCGGTGATTg                                                                                >  2:1488833/1‑149 (MQ=255)
                                                                cAACCTGACGCCGTATCAGTCGGCGCGTCGTCTGGGCGGTAACGGCGATGTCTGGCTGAACGCCAACGAATACCCCACTGCAGTGGAGTTTCAGCTTACTCAGCAAACGCTCAACCGCTACCCGGAATGCCAGCCGAAAGCGGTGATTg                                                                                <  1:1708445/149‑1 (MQ=255)
                                                                cAACCTGACGCCGTATCAGTCGGCGCGTCGTCTGGGCGGTAACGGCGATGTCTGGCTGAACGCCAACGAATACCCCACTGCAGTGGAGTTTCAGCTTACTCAGCAAACGCTCAACCGCTACCCGGAATGCCAGCCGAAAGCGGTGATTg                                                                                <  2:2578389/149‑1 (MQ=255)
                                                                  aCCTGACGCCGTATCAGTCGGCGCGTCGTCTGGGCGGTAACGGCGATGTCTGGCTGAACGCCAACGAATACCCCACTGCAGTGGAGTTTCAGCTTACTCAGCAAACGCTCAACCGCTACCCGGAATGCCAGCCGAAAGCGGTGATTGaa                                                                              >  2:295911/1‑149 (MQ=255)
                                                                   ccTGACGCCGTATCAGTCGGCGCGTCGTCTGGGCGGTAACGGCGATGTCTGGCTGAACGCCAACGAATACCCCACTGCAGTGGAGTTTCAGCTTACTCAGCAAACGCTCAACCGCTACCCGGAATGCCAGCCGAAAGCGGTGATTGaaa                                                                             <  1:3630055/149‑1 (MQ=255)
                                                                      gACGCCGTATCAGTCGGCGCGTCGTCTGGGCGGTAACGGCGATGTCTGGCTGAACGCCAACGAATACCCCACTGCAGTGGAGTTTCAGCTTACTCAGCAAACGCTCAACCGCTACCCGGAATGCCAGCCGAAAGCGGTGATTGAAAAtt                                                                          <  2:2137394/149‑1 (MQ=255)
                                                                        cGCCGTATCAGTCGGCGCGTCGTCTGGGCGGTAACGGCGATGTCTGGCTGAACGCCAACGAATACCCCACTGCAGTGGAGTTTCAGCTTACTCAGCAAACGCTCAACCGCTACCCGGAATGCCAGCCGAAAGCGGTGATTGAAAATTAc                                                                        <  1:3778664/149‑1 (MQ=255)
                                                                           cGTATCAGTCGGCGCGTCGTCTGGGCGGTAACGGCGATGTCTGGCTGAACGCCAACGAATACCCCACTGCAGTGGAGTTTCAGCTTACTCAGCAAACGCTCAACCGCTACCCGGAATGCCAGCCGAAAGCGGTGATTGAAAATTAcgcg                                                                     <  1:4206041/149‑1 (MQ=255)
                                                                            gTATCAGTCGGCGCGTCGTCTGGGCGGTAACGGCGATGTCTGGCTGAACGCCAACGAATACCCCACTGCCGTGGAGTTTCAGCTTACTCAGCAAACGCTCAACCGCTACCCGGAATGCCAGCCGAAAGCGGTGATTGAAAATTACgcgc                                                                    <  2:2497742/149‑1 (MQ=255)
                                                                                cAGTCGGCGCGTCGTCTGGGCGGTAACGGCGATGTCTGGCTGAACGCCAACGAATACCCCACTGCAGTGGAGTTTCAGCTTACTCAGCAAACGCTCAACCGCTACCCGGAATGCCAGCCGAAAGCGGTGATTGAAAATTACGCGCAata                                                                <  2:2671351/149‑1 (MQ=255)
                                                                                    cGGCGCGTCGTCTGGGCGGTAACGGCGATGTCTGTCTGAACGCCAACGAATACCCCACTGCAGTGGAGTTTCAGCTTACTCAGCAAACGCTCAACCGCTACCCGGAATGCCAGCCGAAAGCGGTGATTGAAAATTACGCGCAATATGCa                                                            <  2:3514031/149‑1 (MQ=255)
                                                                                     ggCGCGTCGTCTGGGCGGTAACGGCGATGTCTGGCTGAACGCCAACGAATACCCCACTGCAGTGGAGTTTCAGCTTACTCAGCAAACGCTCAACCGCTACCCGGAATGCCAGCCGAAAGCGGTGATTGAAAATTACGCGCAATATGCAg                                                           <  1:1784972/149‑1 (MQ=255)
                                                                                     ggCGCGTCGTCTGGGCGGTAACGGCGATGTCTGGCTGAACGCCAACGAATACCCCACTGCAGTGGAGTTTCAGCTTACTCAGCAAACGCTCAACCGCTACCCGGAATGCCAGCCGAAAGCGGTGATTGAAAATTACGCGCAATATGCAg                                                           <  2:3236978/149‑1 (MQ=255)
                                                                                       cgcgTCGTCTGGGCGGTAACGGCGATGTCTGGCTGAACGCCAACGAATACCCCACTGCAGTGGAGTTTCAGCTTACTCAGCAAACGCTCAACCGCTACCCGGAATGCCAGCCGAAAGCGGTGATTGAAAATTACGCGCAATATGCAGGc                                                         <  1:318338/149‑1 (MQ=255)
                                                                                        gcgTCGTCTGGGCGGTAACGGCGATGTCTGGCTGAACGCCAACGAATACCCCACTGCAGTGGAGTTTCAGCTTACTCAGCAAACGCTCAACCGCTACCCGGAATGCCAGCCGAAAGCGGTGATTGAAAATTACGCGCAATATGCAGGCg                                                        >  2:3943838/1‑149 (MQ=255)
                                                                                           tcgtcTGGGCGGTAACGGCGATGTCTGGCTGAACGCCAACGAATACCCCACTGCAGTGGAGTTTCAGCTTACTCAGCAAACGCTCAACCGCTACCCGGAATGCCAGCCGAAAGCGGTGATTGAAAATTACGCGCAATATGCAGGCGTaa                                                     >  1:4479549/1‑149 (MQ=255)
                                                                                              tcTGGGCGGTAACGGCGATGTCTGGCTGAACGCCAACGAATACCCCACTGCAGTGGAGTTTCAGCTTACTCAGCAAACGCTCAACCGCTACCCGGAATGCCAGCCGAAAGCGGTGATTGAAAATTACGCGCAATATGCAGGCGTAAAAc                                                  >  2:1418204/1‑149 (MQ=255)
                                                                                                    cGGTAACGGCGATGTCTGGCTGAACGCCAACGAATACCCCACTGCAGTGGAGTTTCAGCTTACTCAGCAAACGCTCAACCGCTACCCGGAATGCCAGCCGAAAGCGGGGATTGAAAATTACGCGCAATATGCAGGCGTAAAACCGGAGc                                            >  1:781911/1‑149 (MQ=255)
                                                                                                         aCGGCGATGTCTGGCTGAACGCCAACGAATACCCCACTGCCGTGGAGTTTCAGCTTACTCAGCAAACGCTCAACCGCTACCCGGAATGCCAGCCGAAAGCGGTGATTGAAAATTAAGCGCAATATGCAGGCGTAAAACCGGAGCAGGTg                                       >  1:593283/1‑149 (MQ=255)
                                                                                                         aCGGCGATGTCTGGCTGAACGCCAACGAATACCCCACTGCAGTGGAGTTTCAGCTTACTCAGCAAACGCTCAACCGCTACCCGGAATGCCAGCCGAACGCGGTGATTGAAAATTACGCGCAATATGCAGGCGTAAAACCGGAGCAGGTg                                       >  2:795273/1‑149 (MQ=255)
                                                                                                         aCGGCGATGTCTGGCTGAACGCCAACGAATACCCCACTGCAGTGGAGTTTCAGCTTACTCAGCAAACGCTCAACCGCTACCCGGAATGCCAGCCGAAAGCGGTGATTGAAAATTACGCGCAATATGCAGGCGTAAAACCGGAGCAGGTg                                       >  1:2279039/1‑149 (MQ=255)
                                                                                                          cGGCGATGTCTGGCTGAACGCCAACGAATACCCCACTGCCGTGGAGTTTCAGCTTACTCAGCAAACGCTCAACCGCTACCCGGAATGCCAGCCGAAAGCGGTGATTGAAAATTACGCGCAATATGCAGGCGTAAAACCGGAGCAGGTGc                                      >  1:623916/1‑149 (MQ=255)
                                                                                                            gCGATGTCTGGCTGAACGCCAACGAATACCCCACTGCCGTGGAGTTTCAGCTTACTCAGCAAACGCTCAACCGCTACCCGGAATGCCAGCCGAAAGCGGTGATTGAAAATTACGCGCAATATGCAGGCGTAAAACCGGAGCAGGTGCTg                                    >  2:4682117/1‑149 (MQ=255)
                                                                                                            gCGATGTCTGGCTGAACGCCAACGAATACCCCACTGCCGTGGAGTTTCAGCTTACTCAGCAAACGCTCAACCGCTACCCGGAATGCCAGCCGAAAGCGGTGATTGAAAATTACGCGCAATATGCAGGCGTAAAACCGGAGCAGGTGCTg                                    >  2:2298534/1‑149 (MQ=255)
                                                                                                            gCGATGTCTGGCTGAACGCCAACGAATACCCCACTGCAGTGGAGTTTCAGCTTACTCAGCAAACGCTCAACCGCTACCCGGAATGCCAGCCGAAAGCGGTGATTGAAAATTACGCGCACTATGCAGGCGTAAAACCGGAGCAGGTGCTg                                    >  1:1850065/1‑149 (MQ=255)
                                                                                                               aTGTCTGGCTGAACGCCAACGAATACCCCACTGCAGTGGAGTTTCAGCTTACTCAGCAAACGCTCAACCGCTACCCGGAATGCCAGCCGAAAGCGGTGATTGAAAATTACGCGCAATATGCAGGCGTAAAACCGGAGCAGGTGCTGGTc                                 >  1:2222063/1‑149 (MQ=255)
                                                                                                                cgTCTGGCTGAACGCCAACGAATACCCCACTGCAGTGGAGTTTCAGCTTACTCAGCAAACGCTCAACCGCTACCCGGAATGCCAGCCGAAAGCGGTGATTGAAAATTACGCGCAATATGCAGGCGTAAAACCGGAGCAGGTGCTGGTCa                                <  2:4776717/148‑1 (MQ=255)
                                                                                                                 gTCTGGCTGAACGCCAACGAATACCCCACTGCAGTGGAGTTTCAGCTTACTCAGCAAACGCTCAACCGCTACCCGGAATGCCAGCCGAAAGCGGTGATTGAAAATTACGCGCAATATGCAGGCGTAAAACCGGAGCAGGTGCTGGTCAg                               >  2:4128802/1‑149 (MQ=255)
                                                                                                                       cTGAACGCCAACGAATACCCCACTGCAGTGGAGTTTCAGCTTACTCAGCAAACGCTCAACCGCTACCCGGAATGCCAGCCGAAAGCGGTGATTGAAAATTACGCGCAATATGCAGGCGTAAAACCGGAGCAGGTGCTGGTCAGCCGTgg                         <  2:4364077/149‑1 (MQ=255)
                                                                                                                        tGAACGCCAACGAATACCCCACTGCAGTGGAGTTTCAGCTTACTCAGCAAACGCTCAACCGCTACCCGGAATGCCAGCCGAAAGCGGTGATTGAAAATTACGCGCAATATGCAGGCGTAAAACCGGAGCAGGTGCTGGTCAGCCGTGgc                        >  2:3449990/1‑149 (MQ=255)
                                                                                                                         gAACGCCAACGAATACCCCACTGCAGTGGAGTTTCAGCTTACTCAGCAAACGCTCAACCGCTACCCGGAATGCCAGCCGAAAGCGGTGATTGAAAATTACGCGCAATATGCAGGCGTAAAACCGGAGCAGGTGCTGGTCAGCCGTGgcg                       >  2:3876620/1‑149 (MQ=255)
                                                                                                                         gAACGCCAACGAATACCCCACTGCAGTGGAGTTTCAGCTTACTCAGCAAACGCTCAACCGCTACCCGGAATGCCAGCCGAAAGCGGTGATTGAAAATTACGCGCAATATGCAGGCGTAAAACCGGAGCAGGTGCTGGTCAGCCGTGgcg                       >  1:2205617/1‑149 (MQ=255)
                                                                                                                           aCGCCAACGAATACCCCACTGCAGTGGAGTTTCAGCTTACTCAGCAAACGCTCAACCGCTACCCGGAATGCCAGCCGAAAGCGGTGATTGAAAATTACGCGCAATATGCAGGCGTAAAACCGGAGCAGGTGCTGGTCAGCCGTGGcgcg                     >  1:435295/1‑149 (MQ=255)
                                                                                                                           aCGCCAACGAATACCCCACTGCAGTGGAGTTTCAGCTTACTCAGCAAACGCTCAACCGCTACCCGGAATGCCAGCCGAAAGCGGTGATTGAAAATTACGCGCAATATGCAGGCGTAAAACCGGAGCAGGTGCTGGTCAGCCGTGGcgcg                     >  2:186744/1‑149 (MQ=255)
                                                                                                                            cGCCAACGAATACCCCACTGCAGTGGAGTTTCAGCTTACTCAGCAAACGCTCAACCGCTACCCGGAATGCCAGCCGAAAGCGGTGATTGAAAATTACGCGCAATATGCAGGCGTAAAACCGGAGCAGGTGCTGGTCAGCCGTGGCGCgg                    <  2:3566567/149‑1 (MQ=255)
                                                                                                                             gCCAACGAATACCCCACTGCAGTGGAGTTTCAGCTTACTCAGCAAACGCTCAACCGCTACCCGGAATGCCAGCCGAAAGCGGTGATTGAAAATTACGCGCAATATGCAGGCGTAAAACCGGAGCAGGTGCTGGTCAGCCGTGGCGCGGa                   <  1:749726/149‑1 (MQ=255)
                                                                                                                              ccAACGAATACCCCACTGCAGTGGAGTTTCAGCTTACTCAGCAAACGCTCAACCGCTACCCGGAATGCCAGCCGAAAGCGGTGATTGAAAATTACGCGCAATATGCAGGCGTAAAACCGGAGCAGGTGCTGGTCAGCCGTGGCGCGGAc                  >  1:2559504/1‑149 (MQ=255)
                                                                                                                                 aCGAATACCCCACTGCAGTGGAGTTTCAGCTTACTCAGCAAACGCTCAACCGCTACCCGGAATGCCAGCCGAAAGCGGTGATTGAAAATTACGCGCAATATGCAGGCGTAAAACCGGAGCAGGTGCTGGTCAGCCGTGGCGCGGACGaa               <  2:2021736/149‑1 (MQ=255)
                                                                                                                                   gAATACCCCACTGCAGTGGAGTTTCAGCTTACTCAGCAAACGCTCAACCGCTACCCGGAATGCCAGCCGAAAGCGGTGATTGAAAATTACGCGCAATATGCAGGCGTAAAACCGGAGCAGGTGCTGGTCAGCCGTGGCGCGGACGAAgg             >  2:3070010/1‑149 (MQ=255)
                                                                                                                                     aTACCCCACTGCAGTGGAGTTTCAGCTTACTCAGCAAACGCTCAACCGCTACCCGGAATGCCAGCCGAAAGCGGTGATTGAAAATTACGCGCAATATGCAGGCGTAAAACCGGAGCAGGTGCTGGTCAGCCGTGGCGCGGACGAAGGtc           >  2:4287853/1‑148 (MQ=255)
                                                                                                                                     aTACCCCACTGCAGTGGAGTTTCAGCTTACTCAGCAAACGCTCAACCGCTACCCGGAATGCCAGCCGAAAGCGGTGATTGAAAATTACGCGCAATATGCAGGCGTAAAACCGGAGCAGGTGCTGGTCAGCCGTGGCGCGGACGAAGGTa           >  2:2289831/1‑149 (MQ=255)
                                                                                                                                     aTACCCCACTGCAGTGGAGTTTCAGCTTACTCAGCAAACGCTCAACCGCTACCCGGAATGCCAGCCGAAAGCGGTGATTGAAAATTACGCGCAATATGCAGGCGTAAAACCGGAGCAGGTGCTGGTCAGCCGTGGCGCGGACGAAGGTa           <  1:541412/149‑1 (MQ=255)
                                                                                                                                      tACCCCACTGCAGTGGAGTTTCAGCTTACTCAGCAAACGCTCAACCGCTACCCGGAATGCCAGCCGAAAGCGGTGATTGAAAATTACGCGCAATATGCAGGCGTAAAACCGGAGCAGGTGCTGGTCAGCCGTGGCGCGGACGAAGGTAt          >  1:894520/1‑149 (MQ=255)
                                                                                                                                       aCCCCACTGCAGTGGAGTTTCAGCTTACTCAGCAAACGCTCAACCGCTACCCGGAATGCCAGCCGAAAGCGGTGATTGAAAATTACGCGCAATATGCAGGCGTAAAACCGGAGCAGGTGCTGGTCAGCCGTGGCGCGGACGAAGGTAtt         <  1:661987/149‑1 (MQ=255)
                                                                                                                                       aCCCCACTGCAGTGGAGTTTCAGCTTACTCAGCAAACGCTCAACCGCTACCCGGAATGCCAGCCGAAAGCGGTGATTGAAAATTACGCGCAATATGCAGGCGTAAAACCGGAGCAGGTGCTGGTCAGCCGTGGCGCGGACGAAGGTAtt         <  2:1141794/149‑1 (MQ=255)
                                                                                                                                         cccACTGCAGTGGAGTTTCAGCTTACTCAGCAAACGCTCAACCGCTACCCGGAATGCCAGCCGAAAGCGGTGATTGAAAATTACGCGCAATATGCAGGCGTAAAACCGGAGCAGGTGCTGGTCAGCCGTGGCGCGGACGAAGGTATTGa       >  1:118371/1‑149 (MQ=255)
                                                                                                                                          ccACTGCAGTGGAGTTTCAGCTTACTCAGCAAACGCTCAACCGCTACCCGGAATGCCAGCCGAAAGCGGTGATTGAAAATTACGCGCAATATGCAGGCGTAAAACCGGAGCAGGTGCTGGTCAGCCGTGGCGCGGACGAAGGTATTGaa      <  1:4073792/149‑1 (MQ=255)
                                                                                                                                          ccACTGCAGTGGAGTTTCAGCTTACTCAGCAAACGCTCAACCGCTACCCGGAATGCCAGCCGAAAGCGGTGATTGAAAATTACGCGCAATATGCAGGCGTAAAACCGGAGCAGGTGCTGGTCAGCCGTGGCGCGGACGAAGGTATTGaa      <  1:4442683/149‑1 (MQ=255)
                                                                                                                                          ccACTGCAGTGGAGTTTCAGCTTACTCAGCAAACGCTCAACCGCTACCCGGAATGCCAGCCGAAAGCGGTGATTGAAAATTACGCGCAATATGCAGGCGTAAAACCGGAGCAGGTGCTGGTCAGCCGTGGCGCGGACGAAGGTATTGaa      <  1:3447329/149‑1 (MQ=255)
                                                                                                                                              tGCAGTGGAGTTTCAGCTTACTCAGCAAACGCTCAACCGCTACCCGGAATGCCAGCCGAAAGCGGTGATTGAAAATTACGCGCAATATGCAGGCGTAAAACCGGAGCAGGTGCTGGTCAGCCGTGGCGCGGACGAAGGTATTGAActgc  <  2:2275224/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                 
AACGCCCTTAAGGAGCAAGCATGAGCACCGTGACTATTACCGATTTAGCGCGTGAAAACGTCCGCAACCTGACGCCGTATCAGTCGGCGCGTCGTCTGGGCGGTAACGGCGATGTCTGGCTGAACGCCAACGAATACCCCACTGCCGTGGAGTTTCAGCTTACTCAGCAAACGCTCAACCGCTACCCGGAATGCCAGCCGAAAGCGGTGATTGAAAATTACGCGCAATATGCAGGCGTAAAACCGGAGCAGGTGCTGGTCAGCCGTGGCGCGGACGAAGGTATTGAACTGC  >  CP009273/2085859‑2086149

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: