Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,208,251 A→T 7.2% T405T (ACA→ACT)  aceA → isocitrate lyase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,208,2510AT 7.2% 641.6 / 31.8 209T405T (ACA→ACT) aceAisocitrate lyase
Reads supporting (aligned to +/- strand):  ref base A (95/99);  new base T (5/10);  total (100/109)
Fisher's exact test for biased strand distribution p-value = 2.91e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.42e-01

Only 100 of 209 total aligned reads displayed.

GCATGTGGTTCAACATGTTTGACCTGGCAAACGCCTATGCCCAGGGCGAGGGTATGAAGCACTACGTTGAGAAAGTGCAGCAGCCGGAATTTGCCGCCGCGAAAGATGGCTATACCTTCGTATCTCACCAGCAGGAAGTGGGTACAGGTTACTTCGATAAAGTGACGACTATTATTCAGGGCGGCACGTCTTCAGTCACCGCGCTGACCGGCTCCACTGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGC  >  CP009273/4208106‑4208397
                                                                                                                                                 |                                                                                                                                                  
gCATGTGGTTCAACATGTTTGACCTGGCAAACGCCTATGCCCAGGGCGAGGGTATGAAGCACTACGTTGAGAAAGTGCAGCAGCCGGAATTTGCCGCCGCGAAAGATGGCTATACCTTCGTATCTCACCAGCAGGAAGTGGGTACAgg                                                                                                                                                  >  2:2730713/1‑148 (MQ=255)
gCATGTGGTTCAACATGTTTGACCTGGCAAACGCCTATGCCCAGGGCGAGGGTATGAAGCACTACGTTGAGAAAGTGCAGCAGCCGGAATTTGCCGCCGCGAAAGATGGCTATACCTTCGTATCTCACCAGCAGGAAGTGGGTACAgg                                                                                                                                                  <  1:424741/148‑1 (MQ=255)
   tgtgGTTCAACATGTTTGACCTGGCAAACGCCTATGCCCAGGGCGAGGGTATGAAGCACTACGTTGAGAAAGTGCAGCAGCCGGAATTTGCCGCCGCGAAAGATGGCTATACCTTCGTATCTCACCAGCAGGAAGTGGGTACAGGTTa                                                                                                                                               >  1:2721110/1‑148 (MQ=255)
    gtgGTTCAACATGTTTGACCTGGCAAACGCCTATGCCCAGGGCGAGGGTATGAAGCACTACGTTGAGAAAGTGCAGCAGCCGGAATTTGCCGCCGCGAAAGATGGCTATACCTTCGTATCTCACCAGCAGGAAGTGGGTACAGGTTAc                                                                                                                                              >  2:2702859/1‑148 (MQ=255)
    gtgGTTCAACATGTTTGACCTGGCAAACGCCTATGCCCAGGGCGAGGGTATGAAGCACTACGTTGAGAAAGTGCAGCAGCCGGAATTTGCCGCCGCGAAAGATGGCTATACCTTCGTATCTCACCAGCAGGAAGTGGGTACAGGTTAc                                                                                                                                              >  1:3507533/1‑148 (MQ=255)
    gtgGTTCAACATGTTTGACCTGGCAAACGCCTATGCCCAGGGCGAGGGTATGAAGCACTACGTTGAGAAAGTGCAGCAGCCGGAATTTGCCGCCGCGAAAGATGGCTATACCTTCGTATCTCACCAGCAGGAAGTGGGTACAGGTTAc                                                                                                                                              <  2:682217/148‑1 (MQ=255)
     tgGTTCAACATGTTTGACCTGGCAAACGCCTATGCCCAGGGCGAGGGTATGAAGCACTACGTTGAGAAAGTGCAGCAGCCGGAATTTGCCGCCGCGAAAGATGGCTATACCTTCGTATCTCACCAGCAGGAAGTGGGTACAGGTTACt                                                                                                                                             >  2:2952965/1‑148 (MQ=255)
      ggTTCAACATGTTTGACCTGGCAAACGCCTATGCCCAGGGCGAGGGTATGAAGCACTACGTTGAGAAAGTGCAGCAGCCGGAATTTGCCGCCGCGAAAGATGGCTATACCTTCGTATCTCACCAGCAGGAAGTGGGTACAGGTTACtt                                                                                                                                            <  2:3004600/148‑1 (MQ=255)
       gTTCAACATGTTTGACCTGGCAAACGCCTATGCCCAGGGCGAGGGTATGAAGCACTACGTTGAGAAAGTGCAGCAGCCGGAATTTGCCGCCGCGAAAGATGGCTATACCTTCGTATCTCACCAGCAGGAAGTGGGTACAGGTTACTTc                                                                                                                                           >  2:3382776/1‑148 (MQ=255)
         tCAACATGTTTGACCTGGCAAACGCCTATGCCCAGGGCGAGGGTATGAAGCACTACGTTGAGAAAGTGCAGCAGCCGGAATTTGCCGCCGCGAAAGATGGCTATACCTTCGTATCTCACCAGCAGGAAGTGGGTACAGGTTACTTCGa                                                                                                                                         <  1:119447/148‑1 (MQ=255)
         tCAACATGTTTGACCTGGCAAACGCCTATGCCCAGGGCGAGGGTATGAAGCACTACGTTGAGAAAGTGCAGCAGCCGGAATTTGCCGCCGCGAAAGATGGCTATACCTTCGTATCTCACCAGCAGGAAGTGGGTACAGGTTACTTCGa                                                                                                                                         <  1:3009377/148‑1 (MQ=255)
            aCATGTTTGACCTGGCAAACGCCTATGCCCAGGGCGAGGGTATGAAGCACTACGTTGAGAAAGTGCAGCAGCCGGAATTTGCCGCCGCGAAAGATGGCTATACCTTCGTATCTCACCAGCAGGAAGTGGGTACAGGTTACTTCGATaa                                                                                                                                      <  2:1991754/148‑1 (MQ=255)
               tGTTTGACCTGGCAAACGCCTATGCCCAGGGCGAGGGTATGAAGCACTACGTTGAGAAAGTGCAGCAGCCGGAATTTGCCGCCGCGAAAGATGGCTATACCTTCGTATCTCACCAGCAGGAAGTGGGTACTGGTTACTTCGATAAAGt                                                                                                                                   <  2:744773/148‑1 (MQ=255)
                 tttGACCTGGCAAACGCCTATGCCCAGGGCGAGGGTATGAAGCACTACGTTGAGAAAGTGCAGCAGCCGGAATTTGCCGCCGCGAAAGATGGCTATACCTTCGTATCTCACCAGCAGGAAGTGGGTACAGGTTACTTCGATAAAGTga                                                                                                                                 >  1:2587798/1‑148 (MQ=255)
                 tttGACCTGGCAAACGCCTATGCCCAGGGCGAGGGTATGAAGCACTACGTTGAGAAAGTGCAGCAGCCGGAATTTGCCGCCGCGAAAGATGGCTATACCTTCGTATCTCACCAGCAGGAAGTGGGTACAGGTTACTTCGATAAAGTga                                                                                                                                 <  2:611839/148‑1 (MQ=255)
                       cTGGCAAACGCCTATGCCCAGGGCGAGGGTATGAAGCACTACGTTGAGAAAGTGCAGCAGCCGGAATTTGCCGCCGCGAAAGATGGCTATACCTTCGTATCTCACCAGCAGGAAGTGGGTACAGGTTACTTCGATAAAGTGACGACta                                                                                                                           <  1:2812923/148‑1 (MQ=255)
                        tGGCAAACGCCTATGCCCAGGGCGAGGGTATGAAGCACTACGTTGAGAAAGTGCAGCAGCCGGAATTTGCCGCCGCGAAAGATGGCTATACCTTCGTATCTCACCAGCAGGAAGTGGGTACAGGTTACTTCGATAAAGTGACGACtat                                                                                                                          <  1:17631/148‑1 (MQ=255)
                         ggCAAACGCCTATGCCCAGGGCGAGGGTATGAAGCACTACGTTGAGAAAGTGCAGCAGCCGGAATTTGCCGCCGCGAAAGATGGCTATACCTTCGTATCTCACCAGCAGGAAGTGGGTACAGGTTACTTCGATAAAGTGACGACtatt                                                                                                                         <  2:3533720/148‑1 (MQ=255)
                          gCAAACGCCTATGCCCAGGGCGAGGGTATGAAGCACTACGTTGAGAAAGTGCAGCAGCCGGAATTTGCCGCCGCGAAAGATGGCTATACCTTCGTATCTCACCAGCAGGAAGTGGGTACAGGTTACTTCGATAAAGTGACGACtatta                                                                                                                        >  2:3131721/1‑148 (MQ=255)
                            aaaCGGCTATGCCCTGGGCGAGGGTATGAAGCACTACGTTGAGAATGTGCAGCAGCCTGAATTTGCCGCCGCGAAAGATGGCTATACCTTCGTATCTCACCAGCAGGAAGTGGGTACAGGTTACTTCGATAAAGTGACGACTattatt                                                                                                                      <  2:2440848/148‑1 (MQ=255)
                            aaaCGCCTATGCCCAGGGCGAGGGTATGAAGCACTACGTTGAGAAAGTGCAGCAGCCGGAATTTGCCGCCGCGAAAGATGGCTATACCTTCGTATCTCACCAGCAGGAAGTGGGTACAGGTTACTTCGATAAAGTGACGACTattatt                                                                                                                      <  1:3248442/148‑1 (MQ=255)
                             aaCGCCTATGCCCAGGGCGAGGGTATGAAGCACTACGTTGAGAAAGTGCAGCAGCCGGAATTTGCCGCCGCGAAAGATGGCTATACCTTCGTATCTCACCAGCAGGAAGTGGGTACAGGTTACTTCGATAAAGTGACGACTATTATTc                                                                                                                     <  1:2108802/148‑1 (MQ=255)
                             aaCGCCTATGCCCAGGGCGAGGGTATGAAGCACTACGTTGAGAAAGTGCAGCAGCCGGAATTTGCCGCCGCGAAAGATGGCTATACCTTCGTATCTCACCAGCAGGAAGTGGGTACAGGTTACTTCGATAAAGTGACGACTATTATTc                                                                                                                     <  2:873948/148‑1 (MQ=255)
                              aCGCCTATGCCCAGGGCGAGGGTATGAAGCACTACGTTGAGAAAGTGCAGCAGCCGGAATTTGCCGCCGCGAAAGATGGCTATACCTTCGTATCTCACCAGCAGGAAGTGGGTACAGGTTACTTCGATAAAGTGACGACTATTATTCa                                                                                                                    <  1:557568/148‑1 (MQ=255)
                                gCCTATGCCCAGGGCGAGGGTATGAAGCACTACGTTGAGAAAGTGCAGCAGCCGGAATTTGCCGCCGCGAAAGATGGCTATACCTTCGTATCTCACCAGCAGGAAGTGGGTACAGGTTACTTCGATAAAGTGACGACTATTATTCAgg                                                                                                                  >  1:743028/1‑148 (MQ=255)
                                 ccTATGCCCAGGGCGAGGGTATGAAGCACTACGTTGAGAAAGTGCAGCAGCCGGAATTTGCCGCCGCGAAAGATGGCTATACCTTCGTATCTCACCAGCAGGAAGTGGGTACAGGTTACTTCGATAAAGTGACGACTATTATTCAggg                                                                                                                 <  2:1636767/148‑1 (MQ=255)
                                 ccTATGCCCAGGGCGAGGGTATGAAGCACTACGTTGAGAAAGTGCAGCAGCCGGAATTTGCCGCCGCGAAAGATGGCTATACCTTCGTATCTCACCAGCAGGAAGTGGGTACAGGTTACTTCGATAAAGTGACGACTATTATTCAggg                                                                                                                 >  1:111114/1‑148 (MQ=255)
                                   tATGCCCAGGGCGAGGGTATGAAGCACTACGTTGAGAAAGTGCAGCAGCCGGAATTTGCCGCCGCGAAAGATGGCTATACCTTCGTATCTCACCAGCAGGAAGTGGGTACAGGTTACTTCGATAAAGTGACGACTATTATTCAGggcg                                                                                                               >  2:2250038/1‑148 (MQ=255)
                                    aTGCCCAGGGCGAGGGTATGAAGCACTACGTTGAGAAAGTGCAGCAGCCGGAATTTGCCGCCGCGAAAGATGGCTATACCTTCGTATCTCACCAGCAGGAAGTGGGTACAGGTTACTTCGATAAAGTGACGACTATTATTCAGggcgg                                                                                                              >  2:1175829/1‑148 (MQ=255)
                                      gCCCAGGGCGAGGGTATGAAGCACTACGTTGAGAAAGTGCAGCAGCCGGAATTTGCCGCCGCGAAAGATGGCTATACCTTCGTATCTCACCAGCAGGAAGTGGGTACAGGTTACTTCGATAAAGTGACGACTATTATTCAGGGCGGCa                                                                                                            >  1:1954855/1‑148 (MQ=255)
                                         cAGGGCGAGGGTATGAAGCACTACGTTGAGAAAGTGCAGCAGCCGGAATTTGCCGCCGCGAAAGATGGCTATACCTTCGTATCTCACCAGCAGGAAGTGGGTACAGGTTACTTCGATAAAGTGACGACTATTATTCAGGGCGGCACGt                                                                                                         <  1:1015682/148‑1 (MQ=255)
                                          aGGGCGAGGGTATGAAGCACTACGTTGAGAAAGTGCAGCAGCCGGAATTTGCCGCCGCGAAAGATGGCTATACCTTCGTATCTCACCAGCAGGAAGTGGGTACAGGTTACTTCGATAAAGTGACGACTATTATTCAGGGCGGCACGTc                                                                                                        <  1:1867652/148‑1 (MQ=255)
                                             gCGAGGGTATGAAGCACTACGTTGAGAAAGTGCAGCAGCCGGAATTTGCCGCCGCGAAAGATGGCTAGACCTTCGTATCTCACCAGCAGGAAGTGGGTACAGGTTACTTCGATAAAGTGACGACTATTATTCAGGGCGGCACGTCTTc                                                                                                     <  2:2978682/148‑1 (MQ=255)
                                               gAGGGTATGAAGCACTACGTTGAGAAAGTGCAGCAGCCGGAATTTGCCGCCGCGAAAGATGGCTATACCTTCGTATCTCACCAGCAGGAAGTGGGTACTGGTTACTTCGATAAAGTGACGACTATTATTCAGGGCGGCACGTCTTCAg                                                                                                   <  2:3372917/148‑1 (MQ=255)
                                                  ggTATGAAGCACTACGTTGAGAAAGTGCAGCAGCCGGAATTTGCCGCCGCGAAAGATGGCTATACCTTCGTATCTCACCAGCAGGAAGTGGGTACAGGTTACTTCGATAAAGTGACGACTATTATTCAGGGCGGCACGTCTTCAGTCa                                                                                                <  2:2994289/148‑1 (MQ=255)
                                                     aTGAAGCACTACGTTGAGAAAGTGCAGCAGCCGGAATTTGCCGCCGCGAAAGATGGCTATACCTTCGTATCTCACCAGCAGGAAGTGGGTACAGGTTACTTCGATAAAGTGACGACTATTATTAAGGGCGGCACGTCTTCAGTCACcg                                                                                             <  1:1145428/148‑1 (MQ=255)
                                                       gAAGCACTACGTTGAGAAAGTGCAGCAGCCGGAATTTGCCGCCGCGAAAGATGGCTATACCTTCGTATCTCACCAGCAGGAAGTGGGTACAGGTTACTTCGATAAAGTGACGACTATTATTCAGGGCGGCACTTCTTCAGTCACcgcg                                                                                           >  1:1508362/1‑148 (MQ=255)
                                                        aaGCACTACGTTGAGAAAGTGCAGCAGCCGGAATTTGCCGCCGCGAAAGATGGCTATACCTTCGTATCTCACCAGCAGGAAGTGGGTACAGGTTACTTCGATAAAGTGACGACTATTATTCAGGGCGGCACGTCTTCAGTCACCgcgc                                                                                          >  1:1165653/1‑148 (MQ=255)
                                                          gCACTACGTTGAGAAAGTGCAGCAGCCGGAATTTGCCGCCGCGAAAGATGGCTATACCTTCGTATCTCACCAGCAGGAAGTGGGTACAGGTTACTTCGATAAAGTGACGACTATTATTCAGGGCGGCACGTCTTCAGTCACCGCGCTg                                                                                        <  1:3341269/148‑1 (MQ=255)
                                                          gCACTACGTTGAGAAAGTGCAGCAGCCGGAATTTGCCGCCGCGAAAGATGGCTATACCTTCGTATCTCACCAGCAGGAAGTGGGTACAGGTTACTTCGATAAAGTGACGACTATTATTCAGGGCGGCACGTCTTCAGTCACCGCGCTg                                                                                        >  1:746252/1‑148 (MQ=255)
                                                          gCACTACGTTGAGAAAGTGCAGCAGCCGGAATTTGCCGCCGCGAAAGATGGCTATACCTTCGTATCTCACCAGCAGGAAGTGGGTACAGGTTACTTCGATAAAGTGACGACTATTATTCAGGGCGGCACGTCTTCAGTCACCGCGCTg                                                                                        <  1:752252/148‑1 (MQ=255)
                                                          gCACTACGTTGAGAAAGTGCAGCAGCCGGAATTTGCCGCCGCGAAAGATGGCTATACCTTCGTATCTCACCAGCAGGAAGTGGGTACAGGTTACTTCGATAAAGTGACGACTATTATTCAGGGCGGCACGTCTTCAGTCACCGCGCTg                                                                                        <  1:2534134/148‑1 (MQ=255)
                                                          gCACTACGTTGAGAAAGTGCAGCAGCCGGAATTTGCCGCCGCGAAAGATGGCTATACCTTCGTATCTCACCAGCAGGAAGTGGGTACAGGTTACTTCGATAAAGTGACGACTATTATTCAGGGCGGCACGTCTTCAGTCACCGCGCTg                                                                                        >  1:3167206/1‑148 (MQ=255)
                                                           cACTACGTTGAGAAAGTGCAGCAGCCGGAATTTGCCGCCGCGAAAGATGGCTATACCTTCGTATCTCACCAGCAGGAAGTGGGTACAGGTTACTTCGATAAAGTGACGACTATTATTCAGGGCGGCACGTCTTCAGTCACCGCGCTGa                                                                                       <  1:1128590/148‑1 (MQ=255)
                                                              tACGTTGAGAAAGTGCAGCAGCCGGAATTTGCCGCCGCGAAAGATGGCTATACCTTCGTATCTCACCAGCAGGAAGTGGGTACAGGTTACTTCGATAAAGTGACGACTATTATTCAGGGCGGCACGTCTTCAGTCACCGCGCTGACCg                                                                                    >  1:1841541/1‑148 (MQ=255)
                                                              tACGTTGAGAAAGTGCAGCAGCCGGAATTTGCCGCCGCGAAAGATGGCTATACCTTCGTATCTCACCAGCAGGAAGTGGGTACAGGTTACTTCGATAAAGTGACGACTATTATTCAGGGCGGCACGTCTTCAGTCACCGCGCTGACCg                                                                                    >  1:1460237/1‑148 (MQ=255)
                                                                  ttGAGAAAGTGCAGCAGCCGGAATTTGCCGCCGCGAAAGATGGCTATACCTTCGTATCTCACCAGCAGGAAGTGGGTACAGGTTACTTCGATAAAGTGACGACTATTATTCAGGGCGGCACGTCTTCAGTCACCGCGCTGACCGGCTc                                                                                <  2:794579/148‑1 (MQ=255)
                                                                  ttGAGAAAGTGCAGCAGCCGGAATTTGCCGCCGCGAAAGATGGCTATACCTTCGTATCTCACCAGCAGGAAGTGGGTACAGGTTACTTCGATAAAGTGACGACTATTATTCAGGGCGGCACGTCTTCAGTCACCGCGCTGACCGGCTc                                                                                >  1:2360155/1‑148 (MQ=255)
                                                                   tGAGAAAGTGCAGCAGCCGGAATTTGCCGCCGCGAAAGATGGCTATACCTTCGTATCTCACCAGCAGGAAGTGGGTACAGGTTACTTCGATAAAGTGACGACTATTATTCAGGGCGGCACGTCTTCAGTCACCGCGCTGACCGGCTcc                                                                               <  2:2539085/148‑1 (MQ=255)
                                                                    gagaAAGTGCAGCAGCCGGAATTTGCCGCCGCGAAAGATGGCTATACCTTCGTATCTCACCAGCAGGAAGTGGGTACAGGTTACTTCGATAAAGTGACGACTATTATTCAGGGCGGCACGTCTTCAGTCACCGCGCTGACCGGCTCCa                                                                              <  2:453961/148‑1 (MQ=255)
                                                                    gagaAAGTGCAGCAGCCGGAATTTGCCGCCGCGAAAGATGGCTATACCTTCGTATCTCACCAGCAGGAAGTGGGTACAGGTTACTTCGATAAAGTGACGACTATTATTCAGGGCGGCACGTCTTCAGTCACCGCGCTGACCGGCTCCa                                                                              <  1:1075806/148‑1 (MQ=255)
                                                                    gagaAAGTGCAGCAGCCGGAATTTGCCGCCGCGAAAGATGGCTATACCTTCGTATCTCACCAGCAGGAAGTGGGTACAGGTTACTTCGATAAAGTGACGACTATTATTCAGGGCGGCACGTCTTCAGTCACCGCGCTGACCGGCTCCa                                                                              <  1:3547816/148‑1 (MQ=255)
                                                                    gagaAAGTGCAGCAGCCGGAATTTGCCGCCGCGAAAGATGGCTATACCTTCGTATCTCACCAGCAGGAAGTGGGTACAGGTTACTTCGATAAAGTGACGACTATTATTCAGGGCGGCACGTCTTCAGTCACCGCGCTGACCGGCTCCa                                                                              <  1:264902/148‑1 (MQ=255)
                                                                      gaAAGTGCAGCAGCCGGAATTTGCCGCCGCGAAAGATGGCTATACCTTCGTATCTCACCAGCAGGAAGTGGGTACAGGTTACTTCGATAAAGTGACGACTATTATTCAGGGCGGCACGTCTTCAGTCACCGCGCTGACCGGCTCCACt                                                                            >  1:1883720/1‑148 (MQ=255)
                                                                      gaAAGTGCAGCAGCCGGAATTTGCCGCCGCGAAAGATGGCTATACCTTCGTATCTCACCAGCAGGAAGTGGGTACAGGTTACTTCGATAAAGTGACGACTATTATTCAGGGCGGCACGTCTTCAGTCACCGCGCTGACCGGCTCCACt                                                                            >  1:1884235/1‑148 (MQ=255)
                                                                        aaGTGCAGCAGCCGGAATTTGCCGCCGCGAAAGATGGCTATACCTTCGTATCTCACCAGCAGGAAGTGGGTACAGGTTACTTCGATAAAGTGACGAATATTATTCAGGGCGGCACGTCTTCAGTCACCGCGCTGACCGGCTCCACTga                                                                          <  1:482491/148‑1 (MQ=255)
                                                                          gTGCAGCAGCCGGAATTTGCCGCCGCGAAAGATGGCTATACCTTCGTATCTCACCAGCAGGAAGTGGGTACAGGTTACTTCGATAAAGTGACGACTATTATTCAGGGCGGCACGTCTTCAGTCACCGCGCTGACCGGCTCCACTgaag                                                                        <  2:3213343/148‑1 (MQ=255)
                                                                           tGCAGCAGCCGGAATTTGCCGCCGCGAAAGATGGCTATACCTTCGTATCTCACCAGCAGGAAGTGGGTACAGGTTACTTCGATAAAGTGACGACTATTATTCAGGGCGGCACGTCTTCAGTCACCGCGCTGACCGGCTCCACTgaaga                                                                       <  2:2590287/148‑1 (MQ=255)
                                                                            gcagcaGCCGGAATTTGCCGCCGCGAAAGATGGCTATACCTTCGTATCTCACCAGCAGGAAGTGGGTACTGGTTACTTCGATAAAGTGACGACTATTATTCAGGGCGGCACGTCTTCAGTCACCGCGCTGACCGGCTCCACTgaagaa                                                                      <  2:1590380/148‑1 (MQ=255)
                                                                            gcagcaGCCGGAATTTGCCGCCGCGAAAGATGGCTATACCTTCGTATCTCACCAGCAGGAAGTGGGTACTGGTTACTTCGATAAAGTGACGACTATTATTCAGGGCGGCACGTCTTCAGTCACCGCGCTGACCGGCTCCACTgaagaa                                                                      <  2:1588230/148‑1 (MQ=255)
                                                                               gcagcCGGAATTTGCCGCCGCGAAAGATGGCTATACCTTCGTATCTCACCAGCAGGAAGTGGGTACAGGTTACTTCGATAAAGTGACGACTATTATTCAGGGCGGCACGTCTTCAGTCACCGCGCTGACCGGCTCCACTGAAGAATCg                                                                   <  2:2324455/148‑1 (MQ=255)
                                                                                cagcCGGAATTTGCCGCCGCGAAAGATGGCTATACCTTCGTATCTCACCAGCAGGAAGTGGGTACAGGTTACTTCGATAAAGTGACGACTATTATTCAGGGCGGCACGTCTTCAGTCACCGCGCTGACCGGCTCCACTGAAGAATCGc                                                                  >  1:133269/1‑148 (MQ=255)
                                                                                 agcCGGAATTTGCCGCCGCGAAAGATGGCTATACCTTCGTATCTCACCAGCAGGAAGTGGGTACAGGTTACTTCGATAAAGTGACGACTATTATTCAGGGCGGCACGTCTTCAGTCACCGCGCTGACCGGCTCCACTgaaga                                                                       >  1:399063/1‑142 (MQ=255)
                                                                                 agcCGGAATTTGCCGCCGCGAAAGATGGCTATACCTTCGTATCTCACCAGCAGGAAGTGGGTACAGGTTACTTCGATAAAGTGACGACTATTATTCAGGGCGGCACGTCTTCAGTCACCGCGCTGACCGGCTCCACTgaaga                                                                       <  2:399063/142‑1 (MQ=255)
                                                                                   ccGGAATTTGCCGCCGCGAAAGATGGCTATACCTTCGTATCTCACCAGCAGGAAGTGGGTACAGGTTACTTCGATAAAGTGACGACTATTATTCAGGGCGGCACGTCTTCAGTCACCGCGCTGACCGGCTCCACTGAAGAATCGCAGt                                                               >  1:3322686/1‑148 (MQ=255)
                                                                                    cGGAATTTGCCGCCGCGAAAGATGGCTATACCTTCGTATCTCACCAGCAGGAAGTGGGTACAGGTTACTTCGATAAAGTGACGACTATTATTCAGGGCGGCACGTCTTCAGTCACCGCGCTGACCGGCTCCACTTAAGAATCGCAGtt                                                              <  2:1417367/148‑1 (MQ=255)
                                                                                    cGGAATTTGCCGCCGCGAAAGATGGCTATACCTTCGTATCTCACCAGCAGGAAGTGGGTACAGGTTACTTCGATAAAGTGACGACTATTATTCAGGGCGGCACGTCTTCAGTCACCGCGCTGACCGGCTCCACTGAAGAATCGCAGtt                                                              <  2:214121/148‑1 (MQ=255)
                                                                                       aaTTTGCCGCCGCGAAAGATGGCTATACCTTCGTATCTCACCAGCAGGAAGTGGGTACAGGTTACTTCGATAAAGTGACGACTATTATTCAGGGCGGCACGTCTTCAGTCACCGCGCTGACCGGCTCCACTGAAGAATCGCAGTTCTa                                                           <  2:1507267/148‑1 (MQ=255)
                                                                                         tttGCCGCCGCGAAAGATGGCTATACCTTCGTATCTCACCAGCAGGAAGTGGGTACTGGTTACTTCGATAAAGTGACGACTATTATTCAGGGCGGCACGTCTTCAGTCACCGCGCTGACCGGCTCCACTGAAGAATCGCAGTTCTAAg                                                         >  1:2153223/1‑148 (MQ=255)
                                                                                             ccgccgCGAAAGATGGCTATACCTTCGTATCTCACCAGCAGGAAGTGGGTACAGGTTACTTCGATAAAGTGACGACTATTATTCAGGGCGGCACGTCTTCAGTCACCGCGCTGACCGGCTCCACTGAAGAATCGCAGTTCTAAGcaac                                                     <  2:3298013/148‑1 (MQ=255)
                                                                                                 cgcgAAAGATGGCTATACCTTCGTATCTCACCAGCAGGAAGTGGGTACAGGTTACTTCGATAAAGTGACGACTATTATTCAGGGCGGCACGTCTTCAGTCACCGCGCTGACCGGCTCCACTGAAGAATCGCAGTTCTAAGcaacaaca                                                 <  2:474643/148‑1 (MQ=255)
                                                                                                     aaaGATGGCTATACCTTCGTATCTCACCAGCAGGAAGTGGGTACAGGTTACTTCGATAAAGTGACGACTATTATTCAGGGCGGCACGTCTTCAGTCACCGCGCTGACCGGCTCCACTGAAGAATCGCAGTTCTAAGCAACAACAACCg                                             <  2:844000/148‑1 (MQ=255)
                                                                                                        gATGGCTATACCTTCGTATCTCACCAGCAGGAAGTGGGTACAGGTTACTTCGATAAAGTGACGACTATTATTCAGGGCGGCACGTCTTCAGTCACCGCGCTGACCGGCTCCACTGAAGAATCGCAGTTCTAAGCAACAACAACCGTTg                                          >  2:1402106/1‑148 (MQ=255)
                                                                                                          tGGCTATACCTTCGTATCTCACCAGCAGGAAGTGGGTACTGGTTACTTCGATAAAGTGACGACTATTATTCAGGGCGGCACGTCTTCAGTCACCGCGCTGACCGGCTCCACTGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCt                                        <  2:1711741/148‑1 (MQ=255)
                                                                                                           ggCTATACCTTCGTATCTCACCAGCAGGAAGTGGGTACAGGTTACTTCGATAAAGTGACGACTATTATTCAGGGCGGCACGTCTTCAGTCACCGCGCTGACCGGCTCCACTGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCTg                                       <  2:1801480/148‑1 (MQ=255)
                                                                                                            gCTATACCTTCGTATCTCACCAGCAGGAAGTGGGTACTGGTTACTTCGATAAAGTGACGACTATTATTCAGGGCGGCACGTCTTCAGTCACCGCGCTGACCGGCTCCACTGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCTGa                                      >  1:585959/1‑148 (MQ=255)
                                                                                                            gCTATACCTTCGTATCTCACCAGCAGGAAGTGGGTACAGGTTACTTCGATAAAGTGACGACTATTATTCAGGGCGGCACGTCTTCAGTCACCGCGCTGACCGGCTCCACTGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCTGa                                      <  1:2212613/148‑1 (MQ=255)
                                                                                                            gCTATACCTTCGTATCTCACCAGCAGGAAGTGGGTACAGGTTACTTCGATAAAGTGACGACTATTATTCAGGGCGGCACGTCTTCAGTCACCGCGCTGACCGGCTCCACTGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCTGa                                      <  1:2381326/148‑1 (MQ=255)
                                                                                                              tataCCTTCGTATCTCACCAGCAGGAAGTGGGTACAGGTTACTTCGATAAAGTGACGACTATTATTCAGGGCGGCACGTCTTCAGTCACCGCGCTGACCGGCTCCACTGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCTGACt                                    >  1:96769/1‑148 (MQ=255)
                                                                                                                 aCCTTCGTATCTCACCAGCAGGAAGTGGGTACAGGTTACTTCGATAAAGTGACGACTATTATTCAGGGCGGCACGTCTTCAGTCACCGCGCTGACCGGCTCCACTGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTa                                 >  1:2217575/1‑148 (MQ=255)
                                                                                                                 aCCTTCGTATCTCACCAGCAGGAAGTGGGTACAGGTTACTTCGATAAAGTGACGACTATTATTCAGGGCGGCACGTCTTCAGTCACCGCGCTGACCGGCTCCACTGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTa                                 >  1:71069/1‑148 (MQ=255)
                                                                                                                 aCCTTCGTATCTCACCAGCAGGAAGTGGGTACAGGTTACTTCGATAAAGTGACGACTATTATTCAGGGCGGCACGTCTTCAGTCACCGCGCTGACCGGCTCCACTGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTa                                 >  1:2536159/1‑148 (MQ=255)
                                                                                                                 aCCTTCGTATCTCACCAGCAGGAAGTGGGTACAGGTTACTTCGATAAAGTGACGACTATTATTCAGGGCGGAACGTCTTCAGTCACCGCGCTGACCGGCTCCACTGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTa                                 >  2:690934/1‑148 (MQ=255)
                                                                                                                          tctcACCAGCAGGAAGTGGGTACAGGTTACTTCGATAAAGTGACGACTATTATTCAGGGCGGCACGTCTTCAGTCACCGCGCTGACCGGCTCCACTGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATa                        >  2:1777411/1‑148 (MQ=255)
                                                                                                                           ctcACCAGCAGGAAGTGGGTACAGGTTACTTCGATAAAGTGACGACTATTATTCAGGGCGGCACGTCTTCAGTCACCGCGCTGACCGGCTCCACTGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATaa                       >  2:2277220/1‑148 (MQ=255)
                                                                                                                           ctcACCAGCAGGAAGTGGGTACAGGTTACTTCGATAAAGTGACGACTATTATTCAGGGCGGCACGTCTTCAGTCACCGCGCTGACCGGCTCCACTGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATaa                       >  1:3344847/1‑148 (MQ=255)
                                                                                                                               ccAGCAGGAAGTGGGTACAGGTTACTTCGATAAAGTGACGACTATTATTCAGGGCGGCACGTCTTCAGTCACCGCGCTGACCGGCTCCACTGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCg                   <  1:233224/148‑1 (MQ=255)
                                                                                                                               ccAGCAGGAAGTGGGTACAGGTTACTTCGATAAAGTGACGACTATTATTCAGGGCGGCACGTCTTCAGTCACCGCGCTGACCGGCTCCACTGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCg                   >  1:2095513/1‑148 (MQ=255)
                                                                                                                                cagcagGAAGTGGGTACAGGTTACTTCGATAAAGTGACGACTATTATTCAGGGCGGCACGTCTTCAGTCACCGCGCTGACCGGCTCCACTGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCGt                  <  1:3056545/148‑1 (MQ=255)
                                                                                                                                 agcagGAAGTGGGTACAGGTTACTTCGATAAAGTGACGACTATTATTCAGGGCGGCACGTCTTCAGTCACCGCGCTGACCGGCTCCACTGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCGtt                 >  1:297092/1‑148 (MQ=255)
                                                                                                                                    agGAAGTGGGTACAGGTTACTTCGATAAAGTGACGACTATTATTCAGGGCGGCACGTCTTCAGTCACCGCGCTGACCGGCTCCACTGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCGTTCAc              <  1:1372147/148‑1 (MQ=255)
                                                                                                                                     ggAAGTGGGTACAGGTTACTTCGATAAAGTGACGACTATTATTCAGGGCGGCACGTCTTCAGTCACCGCGCTGACCGGCTCCACTGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCGTTCAcg             <  1:2322312/148‑1 (MQ=255)
                                                                                                                                     ggAAGTGGGTACAGGTTACTTCGATAAAGTGACGACTATTATTCAGGGCGGCACGTCTTCAGTCACCGCGCTGACCGGCTCCACTGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCGTTCAcg             >  2:2737261/1‑148 (MQ=255)
                                                                                                                                      gAAGTGGGTACAGGTTACTTCGATAAAGTGACGACTATTATTCAGGGCGGCACGTCTTCAGTCACCGCGCTGACCGGCTCCACTGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCGTTCAcgc            >  2:849072/1‑148 (MQ=255)
                                                                                                                                      gAAGTGGGTACAGGTTACTTCGATAAAGTGACGACTATTATTCAGGGCGGCACGTCTTCAGTCACCGCGCTGACCGGCTCCACTGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCGTTCAcgc            >  2:1919626/1‑148 (MQ=255)
                                                                                                                                             gTACAGGTTACTTCGATAAAGTGACGACTATTATTCAGGGCGGCACGTCTTCAGTCACCGCGCTGACCGGCTCCACTGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCGTTCACGCCGCATcc     <  1:2361769/148‑1 (MQ=255)
                                                                                                                                               aCAGGTTACTTCGATAAAGTGACGACTATTATTCAGGGCGGCACGTCTTCAGTCACCGCGCTGACCGGCTCCACTGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCGTTCACGCCGCATCCgg   >  1:1501644/1‑148 (MQ=255)
                                                                                                                                               aCAGGTTACTTCGATAAAGTGACGACTATTATTCAGGGCGGCACGTCTTCAGTCACCGCGCTGACCGGCTCCACTGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCGTTCACGCCGCATCCgg   >  2:1243986/1‑148 (MQ=255)
                                                                                                                                               aCAGGTTACTTCGATAAAGTGACGACTATTATTCAGGGCGGCACGTCTTCAGTCACCGCGCTGACCGGCTCCACTGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCGTTCACGCCGCATCCgg   >  1:2961695/1‑148 (MQ=255)
                                                                                                                                                cAGGTTACTTCGATAAAGTGACGACTATTATTCAGGGCGGCACGTCTTCAGTCACCGCGCTGACCGGCTCCACTGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGc  >  1:3442093/1‑148 (MQ=255)
                                                                                                                                                 |                                                                                                                                                  
GCATGTGGTTCAACATGTTTGACCTGGCAAACGCCTATGCCCAGGGCGAGGGTATGAAGCACTACGTTGAGAAAGTGCAGCAGCCGGAATTTGCCGCCGCGAAAGATGGCTATACCTTCGTATCTCACCAGCAGGAAGTGGGTACAGGTTACTTCGATAAAGTGACGACTATTATTCAGGGCGGCACGTCTTCAGTCACCGCGCTGACCGGCTCCACTGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGC  >  CP009273/4208106‑4208397

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