Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,207,552 G→T 100% M172I (ATG→ATT aceA → isocitrate lyase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,207,5520GT98.8% 575.5 / ‑1.7 180M172I (ATG→ATTaceAisocitrate lyase
Reads supporting (aligned to +/- strand):  ref base G (1/2);  new base T (84/93);  total (85/95)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

Only 100 of 183 total aligned reads displayed.

TCAACAACACCTTCCGTCGTGCCGATCAGATCCAATGGTCCGCGGGCATTGAGCCGGGCGATCCGCGCTATGTCGATTACTTCCTGCCGATCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATGAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGTGAAGAAATGCGGTCACATGGGCGGCAAAGTTTTAGTGCCAACTCAGGAAGCTATTCAGAAACTGGTCGCGGCGCGTCTGGCAGCT  >  CP009273/4207407‑4207696
                                                                                                                                                 |                                                                                                                                                
 caacaaCACCTTCCGTCGTGCCGATCAGATCCAATGGTCCGCGGGCATTGAGCCGGGCGATCCGCGCTATGTCGATTACTTCCTGCCGATCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATTaaa                                                                                                                                               >  2:2227644/1‑148 (MQ=255)
  aacaacACCTTCCGTCGTGCCGATCAGATCCAATGGTCCGCGGGCATTGAGCCGGGCGATCCGCGCTATGTCGATTACTTCCTGCCGATCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATTAAAg                                                                                                                                              <  1:1065107/148‑1 (MQ=255)
    caacaCCTTCCGTCGTGCCGATCAGATCCAATGGTCCGCGGGCATTGAGCCGGGCGATCCGCGCTATGTCGATTACTTCCTGCCGATCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATTAAAGCg                                                                                                                                            <  2:404028/148‑1 (MQ=255)
     aacaCCTTCCGTCGTGCCGATCAGATCCAATGGTCCGCGGGCATTGAGCCGGGCGATCCGCGCTATGTCGATTACTTCCTGCCGATCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATTAAAGCga                                                                                                                                           >  1:1494625/1‑148 (MQ=255)
       cacCTTCCGTCGTGCCGATCAGATCCAATGGTCCGCGGGCATTGAGCCGGGCGATCCGCGCTATGTCGATTACTTCCTGCCGATCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATTAAAGCgatg                                                                                                                                         >  1:899656/1‑148 (MQ=255)
          cTTCCGTCGTGCCGATCAGATCCAATGGTCCGCGGGCATTGAGCCGGGCGATCCGCGCTATGTCGATTACTTCCTGCCGATCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATTAAAGCGATGAtt                                                                                                                                      >  1:1102444/1‑148 (MQ=255)
          cTTCCGTCGTGCCGATCAGATCCAATGGTCCGCGGGCATTGAGCCGGGCGATCCGCGCTATGTCGATTACTTCCTGCCGATCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATTAAAGCGATGAtt                                                                                                                                      <  1:2563382/148‑1 (MQ=255)
          cTTCCGTCGTGCCGATCAGATCCAATGGTCCGCGGGCATTGAGCCGGGCGATCCGCGCTATGTCGATTACTTCCTGCCGATCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATTAAAGCGATGAtt                                                                                                                                      <  1:2117214/148‑1 (MQ=255)
           ttCCGTCGTGCCGATCAGATCCAATGGTCCGCGGGCATTGAGCCGGGCGATCCGCGCTATGTCGATTACTTCCTGCCGATCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATTAAAGCGATGATTg                                                                                                                                     <  1:342065/148‑1 (MQ=255)
             ccGTCGTGCCGATCAGATCCAATGGTCCGCGGGCATTGAGCCGGGCGATCCGCGCTATGTCGATTACTTCCTGCCGATCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATTAAAGCGATGATTGaa                                                                                                                                   >  1:425012/1‑148 (MQ=255)
             ccGTCGTGCCGATCAGATCCAATGGTCCGCGGGCATTGAGCCGGGCGATCCGCGCTATGTCGATTACTTCCTGCCGATCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATTAAAGCGATGATTGaa                                                                                                                                   <  1:2734728/148‑1 (MQ=255)
             ccGTCGTGCCGATCAGATCCAATGGTCCGCGGGCATTGAGCCGGGCGATCCGCGCTATGTCGATTACTTCCTGCCGATCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATTAAAGCGATGATTGaa                                                                                                                                   >  1:1365011/1‑148 (MQ=255)
             ccGTCGTGCCGATCAGATCCAATGGTCCGCGGGCATTGAGCCGGGCGATCCGCGCTATGTCGATTACTTCCTGCCGATCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATTAAAGCGATGATTGaa                                                                                                                                   <  2:3048650/148‑1 (MQ=255)
                tcgtGCCGATCAGATCCAATGGTCCGCGGGCATTGAGCCGGGCGATCCGCGCTATGTCGATTACTTCCTGCCGATCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATTAAAGCGATGATTGAAGcc                                                                                                                                >  2:1480999/1‑148 (MQ=255)
                tcgtGCCGATCAGATCCAATGGTCCGCGGGCATTGAGCCGGGCGATCCGCGCTATGTCGATTACTTCCTGCCGATCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATTAAAGCGATGATTGAAGcc                                                                                                                                <  2:384987/148‑1 (MQ=255)
                     ccGATCAGATCCAATGGTCCGCGGGCATTGAGCCGGGCGATCCGCGCTATGTCGATTACTTCCTGCCGATCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATTAAAGCGATGATTGAAGCCGGTGc                                                                                                                           >  1:3464949/1‑148 (MQ=255)
                       gATCAGATCCAATGGTCCGCGGGCATTGAGCCGGGCGATCCGCGCTATGTCGATTACTTCCTGCCGATCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATTAAAGCGATGATTGAAGCCGGTGCAg                                                                                                                         <  2:2052550/148‑1 (MQ=255)
                         tCAGATCCAATGGTCCGCGGGCATTGAGCCGGGCGATCCGCGCTATGTCGATTACTTCCTGCCGATCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATTAAAGCGATGATTGAAGCCGGTGCAGCg                                                                                                                       >  2:2565940/1‑148 (MQ=255)
                            gATCCAATGGTCCGCGGGCATTGAGCCGGGCGATCCGCGCTATGTCGATTACTTCCTGCCGATCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATTAAAGCGATGATTGAAGCCGGTGCAGCGGCa                                                                                                                    <  2:2687268/148‑1 (MQ=255)
                               ccAATGGTCCGCGGGCATTGAGCCGGGCGATCCGCGCTATGTCGATTACTTCCTGCCGATCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATTAAAGCGATGATTGAAGCCGGTGCAGCGGCAGtt                                                                                                                 <  1:1255158/148‑1 (MQ=255)
                               ccAATGGTCCGCGGGCATTGAGCCGGGCGATCCGCGCTATGTCGATTACTTCCTGCCGATCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATTAAAGCGATGATTGAAGCCGGTGCAGCGGCAGtt                                                                                                                 <  2:2605640/148‑1 (MQ=255)
                                 aaTGGGCCGCGGGCATTGAGCCGGGCGATCCGCGCTATGTCGATTACTTCCTGCCGATCGTTGCCGATGCGGTAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATTAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCa                                                                                                               <  2:3156050/148‑1 (MQ=255)
                                  aTGGTCCGCGGGCATTGAGCCGGGCGATCCGCGCTATGTCGATTACTTCCTGCCGATCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATTAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCAc                                                                                                              >  1:1350939/1‑148 (MQ=255)
                                   tGGTCCGCGGGCATTGAGCCGGGCGATCCGCGCTATGTCGATTACTTCCTGCCGATCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATTAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACt                                                                                                             <  2:1640379/148‑1 (MQ=255)
                                     gTCCGCGGGCATTGAGCCGGGCGATCCGCGCTATGTCGATTACTTCCTGCCGATCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATTAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTc                                                                                                           >  2:2296770/1‑148 (MQ=255)
                                          cgGGCATTGAGCCGGGCGATCCGCGCTATGTCGATTACTTCCTGCCGATCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATTAAAGCGATGATTGAAGCCGGTGCAGCCCCAGTTCACTTCGAAGa                                                                                                      >  1:653119/1‑148 (MQ=255)
                                           gggCATTGAGCCGGGCGATCCGCGCTATGTCGATTACTTCCTGCCGATCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATTAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGAt                                                                                                     <  1:1007662/148‑1 (MQ=255)
                                            ggCATTGAGCCGGGCGATCCGCGCTATGTCGATTACTTCCTGCCGATCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATTAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATc                                                                                                    <  1:2184561/148‑1 (MQ=255)
                                            ggCATTGAGCCGGGCGATCCGCGCTATGTCGATTACTTCCTGCCGATCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATTAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATc                                                                                                    <  2:634507/148‑1 (MQ=255)
                                              cATTGAGCCGGGCGATCCGCGCTATGTCGATTACTTCCTGCCGATCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATTAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAg                                                                                                  <  2:1103965/148‑1 (MQ=255)
                                                ttGAGCCGGGCGATCCGCGCTATGTCGATTACTTCCTGCCGATCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATTAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCt                                                                                                <  1:855082/148‑1 (MQ=255)
                                                 tGAGCCGGGCGATCCGCGCTATGTCGATTACTTCCTGCCGATCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATTAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTg                                                                                               >  2:196148/1‑148 (MQ=255)
                                                            aTCCGCGCTATGTCGATTACTTCCTGCCGATCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATTAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGTgaa                                                                                    >  1:1965513/1‑148 (MQ=255)
                                                             tCCGCGCTATGTCGATTACTTCCTGCCGATCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATTAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGTgaag                                                                                   >  1:2397556/1‑148 (MQ=255)
                                                              ccGCGCTATGTCGATTACTTCCTGCCGATCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATTAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGTgaaga                                                                                  >  1:3106358/1‑148 (MQ=255)
                                                                    tATGTCGATTACTTACTGCCGATCGTTTCCCATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAAATGATTAAAGCGATGATTGAATCCGGTGCAGCGGCAGTTCACTTCGAAGAACAGCTGGCGTCAGTGAAGAAATGCg                                                                            >  1:2839898/1‑148 (MQ=255)
                                                                     aTGTCGATTACTTCCTGCCGATCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATTAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGTGAAGAAATGCgg                                                                           <  1:2103279/148‑1 (MQ=255)
                                                                     aTGTCGATTACTTCCTGCCGATCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATTAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGTGAAGAAATGCgg                                                                           <  2:3106358/148‑1 (MQ=255)
                                                                     aTGTCGATTACTTCCTGCCGATCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATTAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGTGAAGAAATGCgg                                                                           >  1:2561570/1‑148 (MQ=255)
                                                                       gTCGATTACTTCCCGCCGAGCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATTAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGTGAAGAAATGCGGTc                                                                         <  2:2902323/148‑1 (MQ=255)
                                                                          gATTACTTCCTGCCGATCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATTAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGTGAAGAAATGCGGTcaca                                                                      <  2:645003/148‑1 (MQ=255)
                                                                          gATTACTTCCTGCCGATCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATTAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGTGAAGAAATGCGGTcaca                                                                      >  2:2987172/1‑148 (MQ=255)
                                                                           aTTACTTCCTGCCGATCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATTAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGTGAAGAAATGCGGTCACAt                                                                     <  2:2057762/148‑1 (MQ=255)
                                                                              aCTTCCTGCCGATCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATTAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGTGAAGAAATGCGGTCACATggg                                                                  <  2:2753764/148‑1 (MQ=255)
                                                                                    tGCCGATCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATTAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGTGAAGAAATGCGGTCACATGGGCGGCaa                                                            >  2:2371412/1‑148 (MQ=255)
                                                                                      ccGATCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATTAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGTGAAGAAATGCGGTCACATGGGCGGCAAAg                                                          <  1:945563/148‑1 (MQ=255)
                                                                                      ccGATCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATTAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGTGAAGAAATGCGGTCACATGGGCGGCAAAg                                                          <  2:1664809/148‑1 (MQ=255)
                                                                                        gATCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATTAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGTGAAGAAATGCGGTCACATGGGCGGCAAAGtt                                                        >  2:3006046/1‑148 (MQ=255)
                                                                                         aTCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATTAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGTGAAGAAATGCGGTCACATGGGCGGCAAAGttt                                                       >  2:673178/1‑148 (MQ=255)
                                                                                          tCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATTAAAGCGATGATTGAAGCCGGTTCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGTGAAGAAATGCGGTCACATGGGCGGCAAAGtttt                                                      >  1:1544071/1‑148 (MQ=255)
                                                                                          tCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATTAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGTGAAGAAATGCGGTCACATGGGCGGCAAAGtttt                                                      >  1:712195/1‑148 (MQ=255)
                                                                                           cGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATTAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGTGAAGAAATGCGGTCACATGGGCGGCAAAGTTTTa                                                     <  2:1835299/148‑1 (MQ=255)
                                                                                           cGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATTAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGTGAAGAAATGCGGTCACATGGGCGGCAAAGTTTTa                                                     >  2:979472/1‑148 (MQ=255)
                                                                                           cGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATTAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGTGAAGAAATGCGGTCACATGGGCGGCAAAGTTTTa                                                     <  2:1168066/148‑1 (MQ=255)
                                                                                            gTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATTAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGTGAAGAAATGCGGTCACATGGGCGGCAAAGTTTTAg                                                    <  2:1401405/148‑1 (MQ=255)
                                                                                              tGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATTAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGTGAAGAAATGCGGTcaca                                                                      >  2:569784/1‑128 (MQ=255)
                                                                                              tGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATTAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGTGAAGAAATGCGGTcaca                                                                      <  1:569784/128‑1 (MQ=255)
                                                                                              tGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATTAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGTGAAGAAATGCGGTCACATGGGCGGCAAAGTTTTAGTg                                                  >  1:3170495/1‑148 (MQ=255)
                                                                                              tGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATTAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGTGAAGAAATGCGGTCACATGGGCGGCAAAGTTTTAGTg                                                  >  1:2708656/1‑148 (MQ=255)
                                                                                              tGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATTAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGTGAAGAAATGCGGTCACATGGGCGGCAAAGTTTTAGTg                                                  >  1:2708543/1‑148 (MQ=255)
                                                                                              tGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATTAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGTGAAGAAATGCGGTCACATGGGCGGCAAAGTTTTAGTg                                                  >  1:1845927/1‑148 (MQ=255)
                                                                                                   aTGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATTAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGTGAAGAAATGCGGTCACATGGGCGGCAAAGTTTTAGTGCCAAc                                             >  2:1090600/1‑148 (MQ=255)
                                                                                                   aTGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATTAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGTGAAGAAATGCGGTCACATGGGCGGCAAAGTTTTAGTGCCAAc                                             >  2:469447/1‑148 (MQ=255)
                                                                                                   aTGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATTAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGTGAAGAAATGCGGTCACATGGGCGGCAAAGTTTTAGTGCCAAc                                             >  2:1090744/1‑148 (MQ=255)
                                                                                                     gCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATTAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGTGAAGAAATGCGGTCACATGGGCGGCAAAGTTTTAGTGCCAACTc                                           <  2:757708/148‑1 (MQ=255)
                                                                                                      cGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATTAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGTGAAGAAATGCGGTCACATGGGCGGCAAAGTTTTAGTGCCAACTCa                                          >  1:732356/1‑148 (MQ=255)
                                                                                                      cGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATTAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGTGAAGAAATGCGGTCACATGGGCGGCAAAGTTTTAGTGCCAACTCa                                          <  2:425012/148‑1 (MQ=255)
                                                                                                      cGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATTAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGTGAAGAAATGCGGTCACATGGGCGGCAAAGTTTTAGTGCCAACTCa                                          <  1:3357506/148‑1 (MQ=255)
                                                                                                      cGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATTAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGTGAAGAAATGCGGTCACATGGGCGGCAAAGTTTTAGTGCCAACTCa                                          <  2:1365011/148‑1 (MQ=255)
                                                                                                      cGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATTAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGTGAAGAAATGCGGTCACATGGGCGGCAAAGTTTTAGTGCCAACTCa                                          <  2:1365010/148‑1 (MQ=255)
                                                                                                        gAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATTAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGTGAAGAAATGCGGTCACATGGGCGGCAAAGTTTTAGTGCCAACTCAgg                                        >  2:3395961/1‑148 (MQ=255)
                                                                                                         aaGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATTAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGTGAAGAAATGCGGTCACATGGGCGGCAAAGTTTTAGTGCCAACTCAGGa                                       <  2:438994/148‑1 (MQ=255)
                                                                                                         aaGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATTAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGTGAAGAAATGCGGTCACATGGGCGGCAAAGTTTTAGTGCCAACTCAGGa                                       <  2:438986/148‑1 (MQ=255)
                                                                                                         aaGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATGAAAGCTATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGTGAAGAAATGCGGTCACATGGGCGGCAAAGTTTTAGTGCCAACTCAGGa                                       <  2:2160785/148‑1 (MQ=255)
                                                                                                           gCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATTAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGTGAAGAAATGCGGTCACATGGGCGGCAAAGTTTTAGTGCCAACTCAGGAAg                                     >  1:1114122/1‑148 (MQ=255)
                                                                                                              ggTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATTAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGTGAAGAAATGCGGTCACATGGGCGGCAAAGTTTTAGTGCCAACTCAGGAAGCTa                                  <  1:979472/148‑1 (MQ=255)
                                                                                                                ttttGGCGGTGTCCTGAATGCCTTTGAACTGATTAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACGTCGAAGATCAGCTGGCGTCAGTGAAGAAATGCGGTCACATGGGCGGCAAAGTTTTAGTGCCAACTCAGGAAGCtctg                                >  2:3396950/1‑145 (MQ=255)
                                                                                                                 tttGGCGGTGTCCTGAATGCCTTTGAACTGATTAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGTGAAGAAATGCGGTCACATGGGCGGCAAAGTTTTAGTGCCAACTCAGGAAGCTATTc                               <  2:93941/148‑1 (MQ=255)
                                                                                                                 tttGGCGGTGTCCTGAATGCCTTTGAACTGATTAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGTGAAGAAATGCGGTCACATGGGCGGCAAAGTTTTAGTGCCAACTCAGGAAGCTATTc                               <  2:2704271/148‑1 (MQ=255)
                                                                                                                      cGGTGTCCTGAATGCCTTTGAACTGATTAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGTGAAGAAATGCGGTCACATGGGCGGCAAAGTTTTAGTGCCAACTCAGGAAGCTATTCAGaaa                          >  2:1799594/1‑148 (MQ=255)
                                                                                                                      cGGTGTCCTGAATGACTTTGAACTGATTAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGTGAAGAAATGCGGTCACATGGGCGGCAAAGTTTTAGTGCCAACTCAGGAAGCTATTCAGaaa                          <  2:3301489/148‑1 (MQ=255)
                                                                                                                        gtgtCCTGAATGCCTTTGAACTGATTAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGTGAAGAAATGCGGTCACATGGGCGGCAAAGTTTTAGTGCCAACTCAGGAAGCTATTCAGAAACt                        >  1:3283517/1‑148 (MQ=255)
                                                                                                                         tgtCCTGAATGCCTTTGAACTGATTAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGTGAAGAAATGCGGTCACATGGGCGGCAAAGTTTTAGTGCCAACTCAGGAAGCTATTCAGAAACTg                       <  2:2887872/148‑1 (MQ=255)
                                                                                                                         tgtCCTGAATGCCTTTGAACTGATTAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGTGAAGAAATGCGGTCACATGGGCGGCAAAGTTTTAGTGCCAACTCAGGAAGCTATTCAGAAACTg                       >  1:2703620/1‑148 (MQ=255)
                                                                                                                         tgtCCTGAATGCCTTTGAACTGATTAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGTGAAGAAATGCGGTCACATGGGCGGCAAAGTTTTAGTGCCAACTCAGGAAGCTATTCAGAAACTg                       <  2:1015070/148‑1 (MQ=255)
                                                                                                                         tgtCCTGAATGCCTTTGAACTGATTAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGTGAAGAAATGCGGTCACATGGGCGGCAAAGTTTTAGTGCCAACTCAGGAAGCTATTCAGAAAATg                       >  2:90659/1‑148 (MQ=255)
                                                                                                                           tCCTGAATGCCTTTGAACTGATTAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGTGAAGAAATGCGGTCACATGGGCGGCAAAGTTTTAGTGCCAACTCAGGAAGCTATTCAGAAACTGGt                     <  1:2857369/148‑1 (MQ=255)
                                                                                                                           tCCTGAATGCCTTTGAACTGATTAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGTGAAGAAATGCGGTCACATGGGCGGCAAAGTTTTAGTGCCAACTCAGGAAGCTATTCAGAAACTGGt                     <  2:13077/148‑1 (MQ=255)
                                                                                                                              tGAATGCCTTTGAACTGATTAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGTGAAGAAATGCGGTCACATGGGCGGCAAAGTTTTAGTGCCAACTCAGGAAGCTATTCAGAAACTGGTcgc                  >  1:1257784/1‑148 (MQ=255)
                                                                                                                                     cTTTGAACTGATTAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGTGAAGAAATGCGGTCACATGGGCGGCAAAGTTTTAGTGCCAACTCAGGAAGCTATTCAGAAACTGGTCGCGGCGCGt           >  1:2239182/1‑148 (MQ=255)
                                                                                                                                     cTTTGAACTGATTAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGTGAAGAAATGCGGTCACATGGGCGGCAAAGTTTTAGTGCCAACTCAGGAAGCTATTCAGAAACTGGTCGCGGCGCGt           <  1:1702518/148‑1 (MQ=255)
                                                                                                                                     cTTTGAACTGATTAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGTGAAGAAATGCGGTCACATGGGCGGCAAAGTTTTAGTGCCAACTCAGGAAGCTATTCAGAAACTGGTCGCGGCGCGt           >  1:3239785/1‑148 (MQ=255)
                                                                                                                                     cTTTGAACTGATTAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGTGAAGAAATGCGGTCACATGGGCGGCAAAGTTTTAGTGCCAACTCAGGAAGCTATTCAGAAACTGGTCGCGGCGCGt           >  2:2808358/1‑148 (MQ=255)
                                                                                                                                       ttGAACTGATTAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGTGAAGAAATGCGGTCACATGGGCGGCAAAGTTTTAGTGCCAACTCAGGAAGCTATTCAGAAACTGGTCGCGGCGCGTCt         >  1:621991/1‑148 (MQ=255)
                                                                                                                                       ttGAACTGATTAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGTGAAGAAATGCGGTCACATGGGCGGCAAAGTTTTAGTGCCAACTCAGGAAGCTATTCAGAAACTGGTCGCGGCGCGTCt         >  1:2266127/1‑148 (MQ=255)
                                                                                                                                         gAACTGATTAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGTGAAGAAATGCGGTCACATGGGCGGCAAAGTTTTAGTGCCAACTCAGGAAGCTATTCAGAAACTGGTCGCGGCGCGTCTgg       <  2:608075/148‑1 (MQ=255)
                                                                                                                                         gAACTGATTAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGTGAAGAAATGCGGTCACATGGGCGGCAAAGTTTTAGTGCCAACTCAGGAAGCTATTCAGAAACTGGTCGCGGCGCGTCTgg       <  2:1201263/148‑1 (MQ=255)
                                                                                                                                            cTGATTAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGTGAAGAAATGCGGTCACATGGGCGGCAAAGTTTTAGTGCCAACTCAGGAAGCTATTCAGAAACTGGTCGCGGCGCGTCTGGCAg    <  2:3464949/148‑1 (MQ=255)
                                                                                                                                            cTGATTAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGTGAAGAAATGCGGTCACATGGGCGGCAAAGTTTTAGTGCCAACTCAGGAAGCTATTCAGAAACTGGTCGCGGCGCGTCTGGCAg    >  1:1220848/1‑148 (MQ=255)
                                                                                                                                              gattaAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGTGAAGAAATGCGGTCACATGGGCGGCAAAGTTTTAGTGCCAACTCAGGAAGCTATTCAGAAACTGGTCGCGGCGCGTCTGGCAGCt  >  2:222322/1‑148 (MQ=255)
                                                                                                                                                 |                                                                                                                                                
TCAACAACACCTTCCGTCGTGCCGATCAGATCCAATGGTCCGCGGGCATTGAGCCGGGCGATCCGCGCTATGTCGATTACTTCCTGCCGATCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATGAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGTGAAGAAATGCGGTCACATGGGCGGCAAAGTTTTAGTGCCAACTCAGGAAGCTATTCAGAAACTGGTCGCGGCGCGTCTGGCAGCT  >  CP009273/4207407‑4207696

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: