Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,208,251 A→T 9.9% 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 9.9% 639.4 / 50.5 213T405T (ACA→ACTaceAisocitrate lyase
Reads supporting (aligned to +/- strand):  ref base A (107/85);  new base T (11/10);  total (118/95)
Fisher's exact test for biased strand distribution p-value = 8.20e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.64e-01

Only 100 of 214 total aligned reads displayed.

CAGCATGTGGTTCAACATGTTTGACCTGGCAAACGCCTATGCCCAGGGCGAGGGTATGAAGCACTACGTTGAGAAAGTGCAGCAGCCGGAATTTGCCGCCGCGAAAGATGGCTATACCTTCGTATCTCACCAGCAGGAAGTGGGTACAGGTTACTTCGATAAAGTGACGACTATTATTCAGGGCGGCACGTCTTCAGTCACCGCGCTGACCGGCTCCACTGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCGTTCACGCCGCATC  >  CP009273/4208104‑4208393
                                                                                                                                                   |                                                                                                                                              
caGCATGTGGTTCAACATGTTTGACCTGGCAAACGCCTATGCCCAGGGCGAGGGTATGAAGCACTACGTTGAGAAAGTGCAGCAGCCGGAATTTGCCGCCGCGAAAGATGGCTATACCTTCGTATCTCACCAGCAGGAAGTGGGTACa                                                                                                                                                >  2:121084/1‑148 (MQ=255)
   cATGTGGTTCAACATGTTTGACCTGGCAAACGCCTATGCCCAGGGCGAGGGTATGAAGCACTACGTTGAGAAAGTGCAGCAGCCGGAATTTGCCGCCGCGAAAGATGGCTATACCTTCGTATCTCACCAGCAGGAAGTGGGTACAGGt                                                                                                                                             <  1:2874743/148‑1 (MQ=255)
   cATGTGGTTCAACATGTTTGACCTGGCAAACGCCTATGCCCAGGGCGAGGGTATGAAGCACTACGTTGAGAAAGTGCAGCAGCCGGAATTTGCCGCCGCGAAAGATGGCTATACCTTCGTATCTCACCAGCAGGAAGTGGGTACAGGt                                                                                                                                             >  2:521035/1‑148 (MQ=255)
     tgtgGTTCAACATGTTTGACCTGGCAAACGCCTATGCCCAGGGCGAGGGTATGAAGCACTACGTTGAGAAAGTGCAGCAGCCGGAATTTGCCGCCGCGAAAGATGGCTATACCTTCGTATCTCACCAGCAGGAAGTGGGTACAGGTTa                                                                                                                                           <  1:881353/148‑1 (MQ=255)
     tgtgGTTCAACATGTTTGACCTGGCAAACGCCTATGCCCAGGGCGAGGGTATGAAGCACTACGTTGAGAAAGTGCAGCAGCCGGAATTTGCCGCCGCGAAAGATGGCTATACCTTCGTATCTCACCAGCAGGAAGTGGGTACAGGTTa                                                                                                                                           >  2:365242/1‑148 (MQ=255)
     tgtgGTTCAACATGTTTGACCTGGCAAACGCCTATGCCCAGGGCGAGGGTATGAAGCACTACGTTGAGAAAGTGCAGCAGCCGGAATTTGCCGCCGCGAAAGATGGCTATACCTTCGTATCTCACCAGCAGGAAGTGGGTACAGGTTa                                                                                                                                           >  1:1179547/1‑148 (MQ=255)
      gtgGTTCAACATGTTTGACCTGGCAAACGCCTATGCCCAGGGCGAGGGTATGAAGCACTACGTTGAGAAAGTGCAGCAGCCGGAATTTGCCGCCGCGAAAGATGGCTATACCTTCGTATCTCACCAGCAGGAAGTGGGTCCAGGTTAc                                                                                                                                          >  1:3120874/1‑148 (MQ=255)
      gtgGTTCAACATGTTTGACCTGGCAAACGCCTATGCCCAGGGCGAGGGTATGAAGCACTACGTTGAGAAAGTGCAGCAGCCGGAATTTGCCGCCGCGAAAGATGGCTATACCTTCGTATCTCACCAGCAGGAAGTGGGTACAGGTTAc                                                                                                                                          <  2:2782636/148‑1 (MQ=255)
      gtgGTTCAACATGTTTGACCTGGCAAACGCCTATGCCCAGGGCGAGGGTATGAAGCACTACGTTGAGAAAGTGCAGCAGCCGGAATTTGCCGCCGCGAAAGATGGCTATACCTTCGTATCTCACCAGCAGGAAGTGGGTACAGGTTAc                                                                                                                                          >  2:1299078/1‑148 (MQ=255)
      gtgGTTCAACATGTTTGACCTGGCAAACGCCTATGCCCAGGGCGAGGGTATGAAGCACTACGTTGAGAAAGTGCAGCAGCCGGAATTTGCCGCCGCGAAAGATGGCTATACCTTCGTATCTCACCAGCAGGAAGTGGGTACAGGTTAc                                                                                                                                          >  1:3030317/1‑148 (MQ=255)
      gtgGTTCAACATGTTTGACCTGGCAAACGCCTATGCCCAGGGCGAGGGTATGAAGCACTACGTTGAGAAAGTGCAGCAGCCGGAATTTGCCGCCGCGAAAGATGGCTATACCTTCGTATCTCACCAGCAGGAAGTGGGTACAGGTTAc                                                                                                                                          >  2:1588963/1‑148 (MQ=255)
      gtgGTTCAACATGTTTGACCTGGCAAACGCCTATGCCCAGGGCGAGGGTATGAAGCACTACGTTGAGAAAGTGCAGCAGCCGGAATTTGCCGCCGCGAAAGATGGCTATACCTTCGTATCTCACCAGCAGGAAGTGGGTACAGGTTAc                                                                                                                                          <  2:231204/148‑1 (MQ=255)
           tCAACATGTTTGACCTGGCAAACGCCTATGCCCAGGGCGAGGGTATGAAGCACTACGTTGAGAAAGTGCAGCAGCCGGAATTTGCCGCCGCGAAAGATGGCTATACCTTCGTATCTCACCAGCAGGAAGTGGGTACAGGTTACTTCGa                                                                                                                                     <  1:1654013/148‑1 (MQ=255)
           tCAACATGTTTGACCTGGCAAACGCCTATGCCCAGGGCGAGGGTATGAAGCACTACGTTGAGAAAGTGCAGCAGCCGGAATTTGCCGCCGCGAAAGATGGCTATACCTTCGTATCTCACCAGCAGGAAGTGGGTACAGGTTACTTCGa                                                                                                                                     <  2:2520500/148‑1 (MQ=255)
           tCAACATGTTTGACCTGGCAAACGCCTATGCCCAGGGCGAGGGTATGAAGCACTACGTTGAGAAAGTGCAGCAGCCGGAATTTGCCGCCGCGAAAGATGGCTATACCTTCGTATCTCACCAGCAGGAAGTGGGTACAGGTTACTTCGa                                                                                                                                     >  2:2645586/1‑148 (MQ=255)
             aaCATGTTTGACCTGGCAAACGCCTATGCCCAGGGCGAGGGTATGAAGCACTACGTTGAGAAAGTGCAGCAGCCGGAATTTGCCGCCGCGAAAGATGGCTATACCTGCGTATCTCACCAGCAGGACGTGGGTACAGGTTACTTCGATa                                                                                                                                   >  1:1416942/1‑148 (MQ=255)
              aCATGTTTGACCTGGCAAACGCCTATGCCCAGGGCGAGGGTATGAAGCACTACGTTGAGAAAGTGCAGCAGCCGGAATTTGCCGCCGCGAAAGATGGCTATACCTTCGTATCTCACCAGCAGGAAGTGGGTACAGGTTACTTCGATaa                                                                                                                                  <  2:1932648/148‑1 (MQ=255)
               cATGTTTGACCTGGCAAACGCCTATGCCCAGGGCGAGGGTATGAAGCACTACGTTGAGAAAGTGCAGCAGCCGGAATTTGCCGCCGCGAAAGATGGCTATACCTTCGTATCTCACCAGCAGGAAGTGGGTACAGGTTACTTCGATaaa                                                                                                                                 >  2:1132872/1‑148 (MQ=255)
                aTGTTTGACCTGGCAAACGCCTATGCCCAGGGCGAGGGTATGAAGCACTACGTTGAGAAAGTGCAGCAGCCGGAATTTGCCGCCGCGAAAGATGGCTATACCTTCGTATCTCACCAGCAGGAAGTGGGTACAGGTTACTTCGATAAAg                                                                                                                                >  2:692598/1‑148 (MQ=255)
                aTGTTTGACCTGGCAAACGCCTATGCCCAGGGCGAGGGTATGAAGCACTACGTTGAGAAAGTGCAGCAGCCGGAATTTGCCGCCGCGAAAGATGGCTATACCTTCGTATCTCACCAGCAGGAAGTGGGTACAGGTTACTTCGATAAAg                                                                                                                                <  2:3255175/148‑1 (MQ=255)
                aTGTTTGACCTGGCAAACGCCTATGCCCAGGGCGAGGGTATGAAGCACTACGTTGAGAAAGTGCAGCAGCCGGAATTTGCCGCCGCGAAAGATGGCTATACCTTCGTATCTCACCAGCAGGAAGTGGGTACAGGTTACTTCGAGAAAg                                                                                                                                >  2:352109/1‑148 (MQ=255)
                 tGTTTGACCTGGCAAACGCCTATGCCCAGGGCGAGGGTATGAAGCACTACGTTGAGAAAGTGCAGCAGCCGGAATTTGCCGCCGCGAAAGATGGCTATACCTTCGTATCTCACCAGCAGGAAGTGGGTACAGGTTACTTCGATAAAGt                                                                                                                               >  1:508031/1‑148 (MQ=255)
                 tGTTTGACCTGGCAAACGCCTATGCCCAGGGCGAGGGTATGAAGCACTACGTTGAGAAAGTGCAGCAGCCGGAATTTGCCGCCGCGAAAGATGGCTATACCTTCGTATCTCACCAGCAGGAAGTGGGTACAGGTTACTTCGATAAAGt                                                                                                                               >  1:681702/1‑148 (MQ=255)
                  gTTTGACCTGGCAAACGCCTATGCCCAGGGCGAGGGTATGAAGCACTACGTTGAGAAAGTGCAGCAGCCGGAATTTGCCGCCGCGAAAGATGGCTATACCTTCGTATCTCACCAGCAGGAAGTGGGTACAGGTTACTTCGATAAAGTg                                                                                                                              <  2:1948503/148‑1 (MQ=255)
                  gTTTGACCTGGCAAACGCCTATGCCCAGGGCGAGGGTATGAAGCACTACGTTGAGAAAGTGCAGCAGCCGGAATTTGCCGCCGCGAAAGATGGCTATACCTTCGTATCTCACCAGCAGGAAGTGGGTACAGGTTACTTCGATAAAGTg                                                                                                                              >  1:1961813/1‑148 (MQ=255)
                   tttGACCTGGCAAACGCCTATGCCCAGGGCGAGGGTATGAAGCACTACGTTGAGAAAGTGCAGCAGCCGGAATTTGCCGCCGCGAAAGATGGCTATACCTTCGTATCTCACCAGCAGGAAGTGGGTACTGGTTACTTCGATAAAGTga                                                                                                                             <  1:529327/148‑1 (MQ=255)
                        ccTGGCAAACGCCTATGCCCAGGGCGAGGGTATGAAGCACTACGTTGAGAAAGTGCAGCAGCCGGAATTTGCCGCCGCGAAAGATGGCTATACCTTCGTATCTCACCAGCAGGAAGTGGGTACAGGTTACTTCGATAAAGTGACGACt                                                                                                                        <  1:63300/148‑1 (MQ=255)
                          tGGCAAACGCCTATGCCCAGGGCGAGGGTATGAAGCACTACGTTGAGAAAGTGCAGCAGCCGGAATTTGCCGCCGCGAAAGATGGCTATACCTTCGTATCTCACCAGCAGGAAGTGGGTACTGGTTACTTCGATAAAGTGACGACtat                                                                                                                      >  2:1334361/1‑148 (MQ=255)
                          tGGCAAACGCCTATGCCCAGGGCGAGGGTATGAAGCACTACGTTGAGAAAGTGCAGCAGCCGGAATTTGCCGCCGCGAAAGATGGCTATACCTTCGTATCTCACCAGCAGGAAGTGGGTACTGGTTACTTCGATAAAGTGACGACtat                                                                                                                      <  1:2211930/148‑1 (MQ=255)
                           ggCAAACGCCTATGCCCAGGGCGAGGGTATGAAGCACTACGTTGAGAAAGTGCAGCAGCCGGAATTTGCCGCCGCGAAAGATGGCTATACCTTCGTATCTCACCAGCAGGAAGTGGGTACAGGTTACTTCGATAAAGTGACGACtatt                                                                                                                     <  2:1840399/148‑1 (MQ=255)
                            gCAAACGCCTATGCCCAGGGCGAGGGTATGAAGCACTACGTTGAGAAAGTGCAGCAGCCGGAATTTGCCGCCGCGAAAGATGGCTATACCTTCGTATCTCACCAGCAGGAAGTGGGTACTGGTTACTTCGATAAAGTGACGACtatta                                                                                                                    >  1:2647490/1‑148 (MQ=255)
                            gCAAACGCCTATGCCCAGGGCGAGGGTATGAAGCACTACGTTGAGAAAGTGCAGCAGCCGGAATTTGCCGCCGCGAAAGATGGCTATACCTTCGTATCTCACCAGCAGGAAGTGGGTACAGGTTACTTCGATAAAGTGACGACtatta                                                                                                                    <  2:751779/148‑1 (MQ=255)
                            gCAAACGCCTATGCCCAGGGCGAGGGTATGAAGCACTACGTTGAGAAAGTGCAGCAGCCGGAATTTGCCGCCGCGAAAGATGGCTATACCTTCGTATCTCACCAGCAGGAAGTGGGTACAGGTTACTTCGATAAAGTGACGACtatta                                                                                                                    <  1:3280735/148‑1 (MQ=255)
                                aCGCCTATGCCCAGGGCGAGGGTATGAAGCACTACGTTGAGAAAGTGCAGCAGCCGGAATTTGCCGCCGCGAAAGATGGCTATACCTTCGTATCTCACCAGCAGGAAGTGGGTACTGGTTACTTCGATAAAGTGACGACTATTATTCa                                                                                                                >  1:2620799/1‑148 (MQ=255)
                                 cGCCTATGCCCAGGGCGAGGGTATGAAGCACTACGTTGAGAAAGTGCAGCAGCCGGAATTTGCCGCCGCGAAAGATGGCTATACCTTCGTATCTCACCAGCAGGAAGTGGGTACAGGTTACTTCGATAAAGTGACGACTATTATTCAg                                                                                                               <  1:1548740/148‑1 (MQ=255)
                                 cGCCTATGCCCAGGGCGAGGGTATGAAGCACTACGTTGAGAAAGTGCAGCAGCCGGAATTTGCCGCCGCGAAAGATGGCTATACCTTCGTATCTCACCAGCAGGAAGTGGGTACAGGTTACTTCGATAAAGTGACGACTATTATTCAg                                                                                                               <  1:1548465/148‑1 (MQ=255)
                                 cGCCTATGCCCAGGGCGAGGGTATGAAGCACTACGTTGAGAAAGTGCAGCAGCCGGAATTTGCCGCCGCGAAAGATGGCTATACCTTCGTATCTCACCAGCAGGAAGTGGGTACAGGTTACTTCGATAAAGTGACGACTATTATTCAg                                                                                                               <  1:1138147/148‑1 (MQ=255)
                                 cGCCTATGCCCAGGGCGAGGGTATGAAGCACTACGTTGAGAAAGTGCAGCAGCCGGAATTTGCCGCCGCGAAAGATGGCTATACCTTCGTATCTCACCAGCAGGAAGTGGGTACAGGTTACTTCGATAAAGTGACGACTATTATTCAg                                                                                                               >  2:2455146/1‑148 (MQ=255)
                                              ggCGAGGGTATGAAGCACTACGTTGAGAAAGTGCAGCAGCCGGAATTTGCCGCCGCGAAAGATGGCTATACCTTCGTATCTCACCAGCAGGAAGTGGGTACAGGTTACTTCGATAAAGTGACGACTATTATTCAGGGCGGCACGTCtt                                                                                                  >  1:1916214/1‑148 (MQ=255)
                                                       aTGAAGCACTACGTTGAGAAAGTGCAGCAGCCGGAATTTGCCGCCGCGAAAGATGGCTATACCTTCGTATCTCACCAGCAGGAAGTGGGTACTGGTTACTTCGATAAAGTGACGACTATTATTCAGGGCGGCACGTCTTCAGTCACcg                                                                                         >  2:1629663/1‑148 (MQ=255)
                                                       aTGAAGCACTACGTTGAGAAAGTGCAGCAGCCGGAATTTGCCGCCGCGAAAGATGGCTATACCTTCGTATCTCACCAGCAGGAAGTGGGTACTGGTTACTTCGATAAAGTGACGACTATTATTCAGGGCGGCACGTCTTCAGTCACcg                                                                                         >  2:1983156/1‑148 (MQ=255)
                                                       aTGAAGCACTACGTTGAGAAAGTGCAGCAGCCGGAATTTGCCGCCGCGAAAGATGGCTATACCTTCGTATCTCACCAGCAGGAAGTGGGTACAGGTTACTTCGATAAAGTGACGACTATTATTCAGGGCGGCACGTCTTCAGTCACcg                                                                                         >  1:3198317/1‑148 (MQ=255)
                                                          aaGCACTACGTTGAGAAAGTGCAGCAGCCGGAATTTGCCGCCGCGAAAGATGGCTATACCTTCGTATCTCACCAGCAGGAAGTGGGTACAGGTTACTTCGATAAAGTGACGACTATTATTCAGGGCGGCACGTCTTCAGTCACCgcgc                                                                                      >  2:779896/1‑148 (MQ=255)
                                                            gCACTACGTTGAGAAAGTGCAGCAGCCGGAATTTGCCGCCGCGAAAGATGGCTATACCTTCGTATCTCACCAGCAGGAAGTGGGTACTGGTTACTTCGATAAAGTGACGACTATTATTCAGGGCGGCACGTCTTCAGTCACCGCGCTg                                                                                    <  2:2789441/148‑1 (MQ=255)
                                                             cACTACGTTGAGAAAGTGCAGCAGCCGGAATTTGCCGCCGCGAAAGATGGCTATACCTTCGTATCTCACCAGCAGGAAGTGGGTACAGGTTACTTCGATAAAGTGACGACTATTATTCAGGGCGGCACGTCTTCAGTCACCGCGCTGa                                                                                   >  1:1658744/1‑148 (MQ=255)
                                                              aCTACGTTGAGAAAGTGCAGCAGCCGGAATTTGCCGCCGCGAAAGATGGCTATACCTTCGTATCTCACCAGCAGGAAGTGGGTACTGGTTACTTCGATAAAGTGACGACTATTATTCAGGGCGGCACGTCTTCAGTCACCGCGCTGAc                                                                                  <  2:310626/148‑1 (MQ=255)
                                                                    ttGAGAAAGTGCAGCAGCCGGAATTTGCCGCCGCGAAAGATGGCTATACCTTCGTATCTCACCAGCAGGAAGTGGGTACAGGTTACTTCGATAAAGTGACGACTATTATTCAGGGCGGCACGTCTTCAGTCACCGCGCTGACCGGCTc                                                                            <  1:1329036/148‑1 (MQ=255)
                                                                     tGAGAAAGTGCAGCAGCCGGAATTTGCCGCCGCGAAAGATGGCTATACCTTCGTATCTCACCAGCAGGAAGTGGGTACAGGTTACTTCGATAAAGTGACGACTATTATTCAGGGCGGCACGTCTTCAGTCACCGCGCTGACCGGCTcc                                                                           >  2:1218130/1‑148 (MQ=255)
                                                                       agaAAGTGCAGCAGCCGGAATTTGCCGCCGCGAAAGATGGCTATACCTTCGTATCTCACCAGCAGGAAGTGGGTACAGGTTACTTCGATAAAGTGACGACTATTATTCAGGGCGGCACGTCTTCAGTCACCGCGCTGACCGGCTCCAc                                                                         <  1:3237522/148‑1 (MQ=255)
                                                                         aaaGTGCAGCAGCCGGAATTTGCCGCCGCGAAAGATGGCTATACCTTCGTATCTCACCAGCAGGAAGTGGGTACAGGTTACTTCGATAAAGTGACGACTATTATTCAGGGCGGCACGTCTTCAGTCACCGCGCTGACCGGCTCCACTg                                                                       >  2:1666846/1‑148 (MQ=255)
                                                                         aaaGTGCAGCAGCCGGAATTTGCCGCCGCGAAAGATGGCTATACCTTCGTATCTCACCAGCAGGAAGTGGGTACAGGTTACTTCGATAAAGTGACGACTATTATTCAGGGCGGCACGTATTCAGTCACCGCGCTGACCGGCTCCACTg                                                                       >  2:1666183/1‑148 (MQ=255)
                                                                           aGTGCAGCAGCCGGAATTTGCCGCCGCGAAAGATGGCTATACCTTCGTATCTCACCAGCAGGAAGTGGGTACAGGTTACTTCGATAAAGTGACGACTATTATTCAGGGCGGCACGTCTTCAGTCACCGCGCTGACCGGCTCCACTgaa                                                                     >  2:1308130/1‑148 (MQ=255)
                                                                            gTGCAGCAGCCGGAATTTGCCGCCGCGAAAGATGGCTATACCTTCGTATCTCACCAGCAGGAAGTGGGTACAGGTTACTTCGATAAAGTGACGACTATTATTCAGGGCGGCACGTCTTCAGTCACCGCGCTGACCGGCTCCACTgaag                                                                    >  2:2175194/1‑148 (MQ=255)
                                                                            gTGCAGCAGCCGGAATTTGCCGCCGCGAAAGATGGCTATACCTTCGTATCTCACCAGCAGGAAGTGGGTACAGGTTACTTCGATAAAGTGACGACTATTATTCAGGGCGGCACGTCTTCAGTCACCGCGCTGACCGGCTCCACTgaag                                                                    >  2:1514136/1‑148 (MQ=255)
                                                                             tGCAGCAGCCGGAATTTGCCGCCGCGAAAGATGGCTATACCTTCGTATCTCACCAGCAGGAAGTGGGTACTGGTTACTTCGATAAAGTGACGACTATTATTCAGGGCGGCACGTCTTCAGTCACCGCGCTGACCGGCTCCACTgaaga                                                                   >  2:1941655/1‑148 (MQ=255)
                                                                                 gcagcCGGAATTTGCCGCCGCGAAAGATGGCTATACCTTCGTATCTCACCAGCAGGAAGTGGGTACAGGTTACTTCGATAAAGTGACGACTATTATTCAGGGCGGCACGTCTTCAGTCACCGCGCTGACCGGCTCCACTGAAGAATCg                                                               <  2:629970/148‑1 (MQ=255)
                                                                                  cagcCGGAATTTGCCGCCGCGAAAGATGGCTATACCTTCGTATCTCACCAGCAGGAAGTGGGTACAGGTTACTTCGATAAAGTGACGACTATTATTCAGGGCGGCACGTCTTCAGTCACCGCGCTGACCGGCTCCACTGAAGAATCGc                                                              >  1:1367592/1‑148 (MQ=255)
                                                                                     ccGGAATTTGCCGCCGCGAAAGATGGCTATACCTTCGTATCTCACCAGCAGGAAGTGGGTACAGGTTACTTCGATAAAGTGACGACTATTATTCAGGGCGGCACGTCTTCAGTCACCGCGCTGACCGGCTCCACTGAAGAATCGCAGt                                                           >  1:476943/1‑148 (MQ=255)
                                                                                      cGGAATTTGCCGCCGCGAAAGATGGCTATACCTTCGTATCTCACCAGCAGGAAGTGGGTACAGGTTACTTCGATAAAGTGACGACTATTATTCAGGGCGGCACGTCTTCAGTCACCGCGCTGACCGGCTCCACTGAAGAATCGCAGtt                                                          <  2:411627/148‑1 (MQ=255)
                                                                                        gAATTTGCCGCCGCGAAAGATGGCTATACCTTCGTATCTCACCAGCAGGAAGTGGGTACATGTTACTTCGATAAAGTGACGACTATTATTCAGGGCGGCACGTCTTCAGTCACCGCGCTGACCGGCTCCACTGAAGAATCGCAGTTCt                                                        >  1:2655557/1‑148 (MQ=255)
                                                                                           tttGCCGCCGCGAAAGATGGCTATACCTTCGTATCTCACCAGCAGGAAGTGGGTACAGGTTACTTCGATAAAGTGACGACTATTATTCAGGGCGGCACGTCTTCAGTCACCGCGCTGACCGGCTCCACTGAAGAATCGCAGTTCTAAg                                                     >  1:892258/1‑148 (MQ=255)
                                                                                             tGCCGCCGCGAAAGATGGCTATACCTTCGTATCTCACCAGCAGGAAGTGGGTACAGGTTACTTCGATAAAGTGACGACTATTATTCAGGGCGGCACGTCTTCAGTCACCGCGCTGACCGGCTCCACTGAAGAATCGCAGTTCTAAGca                                                   <  2:2072544/148‑1 (MQ=255)
                                                                                             tGCCGCCGCGAAAGATGGCTATACCTTCGTATCTCACCAGCAGGAAGTGGGTACAGGTTACTTCGATAAAGTGACGACTATTATTCAGGGCGGCACGTCTTCAGTCACCGCGCTGACCGGCTCCACTGAAGAATCGCAGTTCTAAGca                                                   <  1:179237/148‑1 (MQ=255)
                                                                                              gccgccGCGAAAGATGGCTATACCTTCGTATCTCACCAGCAGGAAGTGGGTACAGGTTACTTCGATAAAGTGACGACTATTATTCAGGGCGGCACGTCTTCAGTCACCGCGCTGACCGGCTCCACTGAAGAATCGCAGTTCTAAGcaa                                                  <  1:851484/148‑1 (MQ=255)
                                                                                                 gccgcGAAAGATGGCTATACCTTCGTATCTCACCAGCAGGAAGTGGGTACAGGTTACTTCGATAAAGTGACGACTATTATTCAGGGCGGCACGTCTTCAGTCACCGCGCTGACCGGCTCCACTGAAGAATCGCAGTTCTAAGcaacaa                                               <  1:1583385/148‑1 (MQ=255)
                                                                                                       aaaGGTGGCTATACCTTCGTATCTCACCAGCAGGAAGTGGGTACAGGTTACTTCGATAAAGTGACGACTATTATTCAGGGCGGCACGTCTTCAGTCACCGCGCTGACCGGCTCCACTGAAGAATCGCAGTTCTAAGCAACAACAACCg                                         >  1:3161130/1‑148 (MQ=255)
                                                                                                       aaaGATGGCTATACCTTCGTATCTCACCAGCAGGAAGTGGGTACAGGTTACTTCGATAAAGTGACGACTATTATTCAGGGCGGCACGTCTTCAGTCACCGCGCTGACCGGCTCCACTGAAGAATCGCAGTTCTAAGCAACAACAACCg                                         >  1:2021153/1‑148 (MQ=255)
                                                                                                       aaaGATGGCTATACCTTCGTATCTCACCAGCAGGAAGTGGGTACAGGTTACTTCGATAAAGTGACGACTATTATTCAGGGCGGCACGTCTTCAGTCACCGCGCTGACCGGCTCCACTGAAGAATCGCAGTTCTAAGCAACAACAACCg                                         >  2:900744/1‑148 (MQ=255)
                                                                                                        aaGATGGCTATACCTTCGTATCTCACCAGCAGGAAGTGGGTACAGGTTACTTCGATAAAGTGACGACTATTATTCAGGGCGGCACGTCTTCAGTCACCGCGCTGACCGGCTCCACTGAAGAATCGCAGTTCTAAGCAACAACAACCGt                                        <  2:955759/148‑1 (MQ=255)
                                                                                                        aaGATGGCTATACCTTCGTATCTCACCAGCAGGAAGTGGGTACAGGTTACTTCGATAAAGTGACGACTATTATTCAGGGCGGCACGTCTTCAGTCACCGCGCTGACCGGCTCCACTGAAGAATCGCAGTTCTAAGCAACAACAACCGt                                        >  2:2600594/1‑148 (MQ=255)
                                                                                                          gATGGCTATACCTTCGTATCTCACCAGCAGGAAGGGGGTACAGGTTACTTCGATAAAGTGACGACTATTATTCAGGGCGGCACGTCTTCAGTCACCGCGCTGACCGGCTCCACTGAAGAATCGCAGTTCTAAGCAACAACAACCGTTg                                      >  2:2130939/1‑148 (MQ=255)
                                                                                                               cTATACCTTCGTATCTCACCAGCAGGAAGTGGGTACAGGTTACTTCGATAAAGTGACGACTATTATTCAGGGCGGCACGTCTTCAGTCACCGCGCTGACCGGCTCCACTGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCTGAc                                 >  1:241199/1‑148 (MQ=255)
                                                                                                                 ataCCTTCGTATCTCACCAGCAGGAAGTGGGTACAGGTTACTTCGATAAAGTGACGACTATTATTCAGGGCGGCACGTCTTCAGTCACCGCGCTGACCGGCTCCACTGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCTGACTg                               >  2:2590555/1‑148 (MQ=255)
                                                                                                                  taCCTTCGTATCTCACCAGCAGGAAGTGGGTACAGGTTACTTCGATAAAGTGACGACTATTATTCAGGGCGGCACGTCTTCAGTCACCGCGCTGACCGGCTCCACTGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGt                              <  1:1183271/148‑1 (MQ=255)
                                                                                                                  taCCTTCGTATCTCACCAGCAGGAAGTGGGTACAGGTTACTTCGATAAAGTGACGACTATTATTCAGGGCGGCACGTCTTCAGTCACCGCGCTGACCGGCTCCACTGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGt                              >  1:2391099/1‑148 (MQ=255)
                                                                                                                  taCCTTCGTATCTCACCAGCAGGAAGTGGGTACAGGTTACTTCGATAAAGTGACGACTATTATTCAGGGCGGCACGTCTTCAGTCACCGCGCTGACCGGCTCCACTGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGt                              >  1:729741/1‑148 (MQ=255)
                                                                                                                   aCCTTCGTATCTCACCAGCAGGAAGTGGGTACAGGTTACTTCGATAAAGTGACGACTATTATTCAGGGCGGCACGTCTTCAGTCACCGCGCTGACCGGCTCCACTGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTa                             <  1:3070777/148‑1 (MQ=255)
                                                                                                                   aCCTTCGTATCTCACCAGCAGGAAGTGGGTACAGGTTACTTCGATAAAGTGACGACTATTATTCAGGGCGGCACGTCTTCAGTCACCGCGCTGACCGGCTCCACTGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTa                             <  2:3192327/148‑1 (MQ=255)
                                                                                                                     cTTCGTATCTCACCAGCAGGAAGTGGGTACAGGTTACTTCGATAAAGTGACGACTATTATTCAGGGCGGCACGTCTTCAGTCACCGCGCTGACCGGCTCCACTGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAgg                           >  2:1172720/1‑148 (MQ=255)
                                                                                                                     cTTCGTATCTCACCAGCAGGAAGTGGGTACAGGTTACTTCGATAAAGTGACGACTATTATTCAGGGCGGCACGTCTTCAGTCACCGCGCTGACCGGCTCCACTGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAgg                           >  1:3144828/1‑148 (MQ=255)
                                                                                                                     cTTCGTATCTCACCAGCAGGAAGTGGGTACAGGTTACTTCGATAAAGTGACGACTATTATTCAGGGCGGCACGTCTTCAGTCACCGCGCTGACCGGCTCCACTGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAgg                           >  1:428540/1‑148 (MQ=255)
                                                                                                                     cTTCGTATCTCACCAGCAGGAAGTGGGTACAGGTTACTTCGATAAAGTGACGACTATTATTCAGGGCGGCACGTCTTCAGTCACCGCGCTGACCGGCTCCACTGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAgg                           >  1:428479/1‑148 (MQ=255)
                                                                                                                      ttCGTATCTCACCAGCAGGAAGTGGGTACAGGTTACTTCGATAAAGTGACGACTATTATTCAGGGCGGCACGTCTTCAGTCACCGCGCTGACCGGCTCCACTGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGc                          <  2:1499737/148‑1 (MQ=255)
                                                                                                                       tCGTATCTCACCAGCAGGAAGTGGGTACAGGTTACTTCGATAAAGTGACGACTATTATTCAGGGCGGCACGTCTTCAGTCACCGCGCTGACCGGCTCCACTGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGcc                         >  2:2941976/1‑148 (MQ=255)
                                                                                                                         gTATCTCACCAGCCGGAAGTGGGTACAGGTTACTTCGATAAAGTGACGACTATTATTCCGGGCGGCACCTCTTCAGTCACCGCGCGGACCGGCTCCACTGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCgg                       <  1:786451/148‑1 (MQ=255)
                                                                                                                          tATCTCACCAGCAGGAAGTGGGTACAGGTTACTTCGATAAAGTGACGACTATTATTCAGGGCGGCACGTCTTCAGTCACCGCGCTGACCGGCTCCACTGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGa                      >  1:1648140/1‑148 (MQ=255)
                                                                                                                          tATCTCACCAGCAGGAAGTGGGTACAGGTTACTTCGATAAAGTGACGACTATTATTCAGGGCGGCACGTCTTCAGTCACCGCGCTGACCGGCTCCACTGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGa                      >  1:1648128/1‑148 (MQ=255)
                                                                                                                           aTCTCACCAGCAGGAAGTGGGTACAGGTTACTTCGATAAAGTGACGACTATTATTCAGGGCGGCACGTCTTCAGTCACCGCGCTGACCGGCTCCACTGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGAt                     >  1:1024123/1‑148 (MQ=255)
                                                                                                                            tctcACCAGCAGGAAGTGGGTACTGGTTACTTCGATAAAGTGACGACTATTATTCAGGGCGGCACGTCTTCAGTCACCGCGCTGACCGGCTCCACTGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATa                    >  2:2447364/1‑148 (MQ=255)
                                                                                                                            tctcACCAGCAGGAAGTGGGTACAGGTTACTTCGATAAAGTGACGACTATTATTCAGGGCGGCACGTCTTCAGTCACCGCGCTGACCGGCTCCACTGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATa                    <  1:1062097/148‑1 (MQ=255)
                                                                                                                             ctcACCAGCAGGAAGTGGGTACAGGTTACTTCGATAAAGTGACGACTATTATTCAGGGCGGCACGTCTTCAGTCACCGCGCTGACCGGCTCCACTGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATaa                   <  1:2919995/148‑1 (MQ=255)
                                                                                                                              tcACCAGCAGGAAGTGGGTACAGGTTACTTCGATAAAGTGACGACTATTATTCAGGGCGGCACGTCTTCAGTCACCGCGCTGACCGGCTCCACTGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATAAg                  >  2:27318/1‑148 (MQ=255)
                                                                                                                                aCCAGCAGGAAGTGGGTACAGGTTACTTCGATAAAGTGACGACTATTATTCAGGGCGGCACGTCTTCAGTCACCGCGCTGACCGGCTCCACTGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGc                >  2:2526764/1‑148 (MQ=255)
                                                                                                                                  cagcagGAAGTGGGTACAGGTTACTTCGATAAAGTGACGACTATTATTCAGGGCGGCACGTCTTCAGTCACCGCGCTGACCGGCTCCACTGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCGt              <  1:2093089/148‑1 (MQ=255)
                                                                                                                                    gcagGAAGTGGGTACAGGTTACTTCGATAAAGTGACGACTATTATTCAGGGCGGCACGTCTTCAGTCACCGCGCTGACCGGCTCCACTGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCGTTc            >  1:1464732/1‑148 (MQ=255)
                                                                                                                                     cagGAAGTGGGTACAGGTTACTTCGATAAAGTGACGACTATTATTCAGGGCGGCACGTCTTCAGTCACCGCGCTGACCGGCTCCACTGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCGTTCa           >  2:2368223/1‑148 (MQ=255)
                                                                                                                                       ggAAGTGGGTACAGGTTACTTCGATAAAGTGACGACTATTATTCAGGGCGGCACGTCTTCAGTCACCGCGCTGACCGGCTCCACTGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCGTTCAcg         >  2:1401087/1‑148 (MQ=255)
                                                                                                                                          aGTGGGTACAGGTTACTTCGATAAAGTGACGACTATTATTCAGGGCGGCACGTCTTCAGTCACCGCGCTGACCGGCTCCACTGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCGTTCAcgccg      >  1:751936/1‑148 (MQ=255)
                                                                                                                                              ggTACAGGTTACTTCGATAAAGTGACGACTATTATTCAGGGCGGCACGTCTTCAGTCACCGCGCTGACCGGCTCCACTGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCGTTCACGCCGCATc  <  1:365033/148‑1 (MQ=255)
                                                                                                                                              ggTACAGGTTACTTCGATAAAGTGACGACTATTATTCAGGGCGGCACGTCTTCAGTCACCGCGCTGACCGGCTCCACTGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCGTTCACGCCGCATc  <  2:868320/148‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                              
CAGCATGTGGTTCAACATGTTTGACCTGGCAAACGCCTATGCCCAGGGCGAGGGTATGAAGCACTACGTTGAGAAAGTGCAGCAGCCGGAATTTGCCGCCGCGAAAGATGGCTATACCTTCGTATCTCACCAGCAGGAAGTGGGTACAGGTTACTTCGATAAAGTGACGACTATTATTCAGGGCGGCACGTCTTCAGTCACCGCGCTGACCGGCTCCACTGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCGTTCACGCCGCATC  >  CP009273/4208104‑4208393

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: