Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,769,440 C→A 19.3% A166A (GCC→GCA aroD → 3‑dehydroquinate dehydratase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,769,4400CA19.3% 785.8 / 120.6 263A166A (GCC→GCAaroD3‑dehydroquinate dehydratase
Reads supporting (aligned to +/- strand):  ref base C (121/91);  new base A (25/26);  total (146/117)
Fisher's exact test for biased strand distribution p-value = 3.47e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 268 total aligned reads displayed.

TATTTACCGGTGATGATCAGGTTAAAGAAACCGTCGCCTACGCCCACGCGCATGATGTGAAAGTAGTCATGTCCAACCATGACTTCCATAAAACGCCGGAAGCCGAAGAAATCATTGCCCGTCTGCGCAAAATGCAATCCTTCGACGCCGATATTCCTAAGATTGCGCTGATGCCGCAAAGTACCAGCGATGTGCTGACGTTGCTTGCCGCGACCCTGGAGATGCAGGAGCAGTATGCCGATCGTCCAATTATCACGATGTCGATGGCAAAAACTGGCGTAATTTCTCGTCTGGCTG  >  CP009273/1769292‑1769588
                                                                                                                                                    |                                                                                                                                                    
 aTTTACCGGTGATGATCAGGTTAAAGAAACCGTCGCCTACGCCCACGCGCATGATGTGAAAGTAGTCATGTCCAACCATGACTTCCATAAAACGCCGGAAGCCGAAGAAATCATTGCCCGTCTGCGCAAAATGCAATCCTTCGACGCCGa                                                                                                                                                    <  2:2661251/150‑1 (MQ=255)
 aTTTACCGGTGATGATCAGGTTAAAGAAACCGTCGCCTACGCCCACGCGCATGATGTGAAAGTAGTCATGTCCAACCATGACTTCCATAAAACGCCGGAAGCCGAAGAAATCATTGCCCGTCTGCGCAAAATGCAATCCTTCGACGCAGa                                                                                                                                                    >  1:477820/1‑150 (MQ=255)
  tttACCGGTGATGATCAGGTTAAAGAAACCGTCGCCTACGCCCACGCGCATGATGTGAAAGTAGTCATGTCCAACCATGACTTCCATAAAACGCCGGAAGCCGAAGAAATCATTGCCCGTCTGCGCAAAATGCAATCCTTCGACGCCGat                                                                                                                                                   <  2:3351891/150‑1 (MQ=255)
  tttACCGGTGATGATCAGGTTAAAGAAACCGTCGCCTACGCCCACGCGCATGATGTGAAAGTAGTCATGTCCAACCATGACTTCCATAAAACGCCGGAAGCCGAAGAAATCATTGCCCGTCTGCGCAAAATGCAATCCTTCGACGCCGat                                                                                                                                                   <  2:3351930/150‑1 (MQ=255)
   ttACCGGTGATGATCAGGTTAAAGAAACCGTCGCCTACGCCCACGCGCATGATGTGAAAGTAGTCATGTCCAACCATGACTTCCATAAAACGCCGGAAGCCGAAGAAATCATTGCCCGTCTGCGCAAAATGCAATCCTTCGACGCCGata                                                                                                                                                  >  1:3678117/1‑150 (MQ=255)
   ttACCGGTGATGATCAGGTTAAAGAAACCGTCGCCTACGCCCACGCGCATGATGTGAAAGTAGTCATGTCCAACCATGACTTCCATAAAACGCCGGAAGCCGAAGAAATCATTGCCCGTCTGCGCAAAATGCAATCCTTCGACGCCGata                                                                                                                                                  <  1:4100383/150‑1 (MQ=255)
    tACCGGTGATGATCAGGTTAAAGAAACCGTCGCCTACGCCCACGCGCATGATGTGAAAGTAGTCATGTCCAACCATGACTTCCATAAAACGCCGGAAGCCGAAGAAATCATTGCCCGTCTGCGCAAAATGCAATCCTTCGACGCCGatat                                                                                                                                                 >  2:2862096/1‑150 (MQ=255)
          tgatgaTCAGGTTAAAGAAACCGTCGCCTACGCCCACGCGCATGATGTGAAAGTAGTCATGTCCAACCATGACTTCCATAAAACGCCGGAAGCCGAAGAAATCATTGCCCGTCTGCGCAAAATGCAATCCTTCGACGCCGATATTCCTaa                                                                                                                                           >  2:1595306/1‑150 (MQ=255)
               atCAGGTTAAAGAAACCGTCGCCTACGCCCACGCGCATGATGTGAAAGTAGTCATGTCCAACCATGACTTCCATAAAACGCCGGAAGCCGAAGAAATCATTGCCCGTCTGCGCAAAATGCAATCCTTCGACGCCGATATTCCTAAGATTg                                                                                                                                      >  1:1853918/1‑150 (MQ=255)
               atCAGGTTAAAGAAACCGTCGCCTACGCCCACGCGCATGATGTGAAAGTAGTCATGTCCAACCATGACTTCCATAAAACGCCGGAAGCCGAAGAAATCATTGCCCGTCTGCGCAAAATGCAATCCTTCGACGCCGATATTCCTAAGATTg                                                                                                                                      >  2:3393330/1‑150 (MQ=255)
               atCAGGTTAAAGAAACCGTCGCCTACGCCCACGCGCATGATGTGAAAGTAGTCATGTCCAACCATGACTTCCATAAAACGCCGGAAGCCGAAGAAATCATTGCCCGTCTGCGCAAAATGCAATCCTTCGACGCCGATATTCCTAAGATTg                                                                                                                                      >  1:1552004/1‑150 (MQ=255)
               atCAGGTTAAAGAAACCGTCGCCTACGCCCACGCGCATGATGTGAAAGTAGTCATGTCCAACCATGACTTCCATAAAACGCCGGAAGCCGAAGAAATCATTGCCCGTCTGCGCAAAATGCAATCCTTCGACGCCGATATTCCTAAGATTg                                                                                                                                      >  2:1199585/1‑150 (MQ=255)
               atCAGGTTAAAGAAACCGTCGCCTACGCCCACGCGCATGATGTGAAAGTAGTCATGTCCAACCATGACTTCCATAAAACGCCGGAAGCCGAAGAAATCATTGCCCGTCTGCGCAAAATGCAATCCTTCGACGCAGATATTCCTAAGATTg                                                                                                                                      <  2:629818/150‑1 (MQ=255)
                tCAGGTTAAAGAAACCGTCGCCTACGCCCACGCGCATGATGTGAAAGTAGTCATGTCCAACCATGACTTCCATAAAACGCCGGAAGCCGAAGAAATCATTGCCCGTCTGCGCAAAATGCAATCCTTCGACGCAGATATTCCTAAGATTgc                                                                                                                                     <  2:1643349/150‑1 (MQ=255)
                 cAGGTTAAAGAAACCGTCGCCTACGCCCACGCGCATGATGTGAAAGTAGTCATGTCCAACCATGACTTCCATAAAACGCCGGAAGCCGAAGAAATCATTGCCCGTCTGCGCAAAATGCAATCCTTCGACGCCGATATTCCTAAGATTgcg                                                                                                                                    <  1:4390267/150‑1 (MQ=255)
                     ttAAAGAAACCGTCGCCTACGCCCACGCGCATGATGTGAAAGTAGTCATGTCCAACCATGACTTCCATAAAACGCCGGAAGCCGAAGAAATCATTGCCCGTCTGCGCAAAATGCAATCCTTCGACGCCGATATTCCTAAGATTGCGCTGa                                                                                                                                >  2:692875/1‑150 (MQ=255)
                        aaGAAACCGTCGCCTACGCCCACGCGCATGATGTGAAAGTAGTCATGTCCAACCATGACTTCCATAAAACGCCGGAAGCCGAAGAAATCATTGCCCGTCTGCGCAAAATGCAATCCTTCGACGCCGATATTCCTAAGATTGCGCTGATGc                                                                                                                             >  2:3450530/1‑150 (MQ=255)
                           aaaCCGTCGCCTACGCCCACGCGCATGATGTGAAAGTAGTCATGTCCAACCATGACTTCCATAAAACGCCGGAAGCCGAAGAAATCATTGCCCGTCTGCGCAAAATGCAATCCTTCGACGCCGATATTCCTAAGATTGCGCTGATGCCGc                                                                                                                          >  2:2976878/1‑150 (MQ=255)
                           aaaCCGTCGCCTACGCCCACGCGCATGATGTGAAAGTAGTCATGTCCAACCATGACTTCCATAAAACGCCGGAAGCCGAAGAAATCATTGCCCGTCTGCGCAAAATGCAATCCTTCGACGCCGATATTCCTAAGATTGCGCTGATGCCGc                                                                                                                          <  2:531894/150‑1 (MQ=255)
                            aaCCGTCGCCTACGCCCACGCGCATGATGTGAAAGTAGTCATGTCCAACCATGACTTCCATAAAACGCCGGAAGCCGAAGAAATCATTGCCCGTCTGCGCAAAATGCAATCCTTCGACGCCGATATTCCTCAGATTGCGCTGATGACGCa                                                                                                                         >  2:2607016/1‑150 (MQ=255)
                               cGTCGCCTACGCCCACGCGCATGATGTGAAAGTAGTCATGTCCAACCATGACTTCCATAAAACGCCGGAAGCCGAAGAAATCATTGCCCGTCTGCGCAAAATGCAATCCTTCGACGCCGATATTCCTAAGATTGCGCTGATGCCGCAAAg                                                                                                                      >  1:519190/1‑150 (MQ=255)
                                   gCCTACGCCCACGCGCATGATGTGAAAGTAGTCATGTCCAACCATGACTTCCATAAAACGCCGGAAGCCGAAGAAATCATTGCCCGTCTGCGCAAAATGCAATCCTTCGACGCCGATATTCCTAAGATTGCGCTGATGCCGCAAAGTAcc                                                                                                                  >  1:107684/1‑150 (MQ=255)
                                    ccTACGCCCACGCGCATGATGTGAAAGTAGTCATGTCCAACCATGACTTCCATAAAACGCCGGAAGCCGAAGAAATCATTGCCCGTCTGCGCAAAATGCAATCCTTCGACGCCGATATTCCTAAGATTGCGCTGATGCCGCAAAGTACCa                                                                                                                 >  1:3048051/1‑150 (MQ=255)
                                    ccTACGCCCACGCGCATGATGTGAAAGTAGTCATGTCCAACCATGACTTCCATAAAACGCCGGAAGCCGAAGAAATCATTGCCCGTCTGCGCAAAATGCAATCCTTCGACGCCGATATTCCTAAGATTGCGCTGATGCCGCAAAGTACCa                                                                                                                 <  1:2214747/150‑1 (MQ=255)
                                    ccTACGCCCACGCGCATGATGTGAAAGTAGTCATGTCCAACCATGACTTCCATAAAACGCCGGAAGCCGAAGAAATCATTGCCCGTCTGCGCAAAATGCAATCCTTCGACGCCGATATTCCGAAGATTGCGCCGATGCCGCAAAGTACCa                                                                                                                 >  2:3918622/1‑150 (MQ=255)
                                      tACGCCCACGCGCATGATGTGAAAGTAGTCATGTCCAACCATGACTTCCATAAAACGCCGGAAGCCGAAGAAATCATTGCCCGTCTGCGCAAAATGCAATCCTTCGACGCCGATATTCCTAAGATTGCGCTGATGCCGCAAAGTACCAGc                                                                                                               >  1:4515127/1‑150 (MQ=255)
                                        cGCCCACGCGCATGATGTGAAAGTAGTCATGTACAACCATGACTTCCATAAAACGCCGGAAGCAGAAGAAATAATTGCCCGTCTGCGCAAAATGCAATCCTTCGAAGCCGATATTCCTAAGATTGCGCTGATGCCGCAAAGTACCAGCGa                                                                                                             <  2:1795987/150‑1 (MQ=255)
                                         gCCCACGCGCATGATGTGAAAGTAGTCATGTCCAACCATGACTTCCATAAAACGCCGGAAGCCGAAGAAATCATTGCCCGTCTGCGCAAAATGCAATCCTTCGACGCCGATATTCCTAAGATTGCGCTGATGCCGCAAAGTACCAGCGAt                                                                                                            <  1:2273475/150‑1 (MQ=255)
                                           ccACGCGCATGATGTGAAAGTAGTCATGTCCAACCATGACTTCCATAAAACGCCGGAAGCCGAAGAAATCATTGCCCGTCTGCGCAAAATGCAATCCTTCGACGCCGATATTCCTAAGATTGCGCTGATGCCGCAAAGTACAAGCGAtgt                                                                                                          >  1:316301/1‑150 (MQ=255)
                                               gcgcATGATGTGAAAGTAGTCATGTCCAACCATGACTTCCATAAAACGCCGGAAGCCGAAGAAATCATTGCCCGTCTGCGCAAAATGCAATCCTTCGACGCCGATATTCCTAAGATTGCGCTGATGCCGCAAAGTACCAGCGATGTGCTg                                                                                                      >  2:2841193/1‑150 (MQ=255)
                                                 gcATGATGTGAAAGTAGTCATGTCCAACCATGACTTCCATAAAACGCCGGAAGCCGAAGAAATCATTGCCCGTCTGCGCAAAATGCAATCCTTCGACGCAGATATTCCTAAGATTGCGCTGATGCCGCAAAGTACCAGCGATGTGCTGAc                                                                                                    <  1:149666/150‑1 (MQ=255)
                                                  cATGATGTGAAAGTAGTCATGTCCAACCATGACTTCCATAAAACGCCGGAAGCCGAAGAAATCATTGCCCGTCTGCGCAAAATGCAATCCTTCGACGCCGATATTCCTAAGATTGCGCTGATGCCGCAAAGTACCAGCGATGTGCTGACg                                                                                                   >  1:3464009/1‑150 (MQ=255)
                                                    tgatgTGAAAGTAGTCATGTCCAACCATGACTTCCATAAAACGCCGGAAGCCGAAGAAATCATTGCCCGTCTGCGCAAAATGCAATCCTTCGACGCCGATATTCCTAAGATTGCGCTGATGCCGCAAAGTACCAGCGATGTGCTGACGtt                                                                                                 >  1:2549737/1‑150 (MQ=255)
                                                    tgatgTGAAAGTAGTCATGTCCAACCATGACTTCCATAAAACGCCGGAAGCCGAAGAAATCATTGCCCGTCTGCGCAAAATGCAATCCTTCGACGCCGATATTCCTAAGATTGCGCTGATGCCGCAAAGTACCAGCGATGTGCTGACGtt                                                                                                 >  1:3732381/1‑150 (MQ=255)
                                                      atgTGAAAGTAGTCATGTCCAACCATGACTTCCATAAAACGCCGGAAGCCGAAGAAATCATTGCCCGTCTGCGCAAAATGCAATCCTTCGACGCCGATATTCCTAAGATTGCGCTGATGCCGCAAAGTACCAGCGATGTGCTGACGttgc                                                                                               <  1:1059184/150‑1 (MQ=255)
                                                       tgtgAAAGTAGTCATGTCCAACCATGACTTCCATAAAACGCCGGAAGCCGAAGAAATCATTGCCCGTCTGCGCAAAATGCAATCCTTCGACGCCGATATTCCTAAGATTGCGCTGATGCCGCAAAGTACCAGCGATGTGCTGACGttgct                                                                                              <  2:742356/150‑1 (MQ=255)
                                                        gtgAAAGTAGTCATGTCCAACCATGACTTCCATAAAACGCCGGAAGCCGAAGAAATCATTGCCCGTCTGCGCAAAATGCAATCCTTCGACGCCGATATTCCTAAGATTGCGCTGATGCCGCAAAGTACCAGCGATGTGCTGACGttgctt                                                                                             <  1:3234143/150‑1 (MQ=255)
                                                        gtgAAAGTAGTCATGTCCAACCATGACTTCCATAAAACGCCGGAAGCCGAAGAAATCATTGCCCGTCTGCGCAAAATGCAATCCTTCGACGCCGATATTCCTAAGATTGCGCTGATGCCGCAAAGTACCAGCGATGTGCTGACGttgctt                                                                                             <  1:433446/150‑1 (MQ=255)
                                                         tgAAAGTAGTCATGTCCAACCATGACTTCCATAAAACGCCGGAAGCCGAAGAAATCATTGCCCGTCTGCGCAAAATGCAATCCTTCGACGCAGATATTCCTAAGATTGCGCTGATGCCGCAAAGTACCAGCGATGTGCTGACGttgcttg                                                                                            >  2:863505/1‑150 (MQ=255)
                                                           aaaGTAGTCATGTCCAACCATGACTTCCATAAAACGCCGGAAGCCGAAGAAATCATTGCCCGTCTGCGCAAAATGCAATCCTTCGACGCCGATATTCCTAAGATTGCGCTGATGCCGCAAAGTACCAGCGATGTGCTGACGTTGCTTGcc                                                                                          >  1:234455/1‑150 (MQ=255)
                                                            aaGTAGTCATGTCCAACCATGACTTCCATAAAACGCCGGAAGCCGAAGAAATCATTGCCCGTCTGCGCAAAATGCAATCCTTCGACGCCGATATTCCTAAGATTGCGCTGATGCCGCAAAGTACCAGCGATGTGCTGACGTTGCTTGCcg                                                                                         >  1:3332294/1‑150 (MQ=255)
                                                            aaGTAGTCATGTCCAACCATGACTTCCATAAAACGCCGGAAGCCGAAGAAATCATTGCCCGTCTGCGCAAAATGCAATCCTTCGACGCCGATATTCCTAAGATTGCGCTGATGCCGCAAAGTACCAGCGATGTGCTGACGTTGCTTGCcg                                                                                         >  1:4550447/1‑150 (MQ=255)
                                                             agtagtCATGTCCAACCATGACTTCCATAAAACGCCGGAAGCCGAAGAAATCATTGCCCGTCTGCGCAAAATGCAATCCTTCGACGCCGATATTCCTAAGATTGCGCTGATGCCGCAAAGTACCAGCGATGTGCTGACGTTGCTTGCcgc                                                                                        >  2:3247407/1‑150 (MQ=255)
                                                             agtagtCATGTCCAACCATGACTTCCATAAAACGCCGGAAGCCGAAGAAATCATTGCCCGTCTGCGCAAAATGCAATCCTTCGACGCCGATATTCCTAAGATTGCGCTGATGCCGCAAAGTACCAGCGATGTGCTGACGTTGCTTGCcgc                                                                                        <  2:3592164/150‑1 (MQ=255)
                                                                 gtCATGTCCAACCATGACTTCCATAAAACGCCGGAAGCCGAAGAAATCATTGCCCGTCTGCGCAAAATGCAATCCTTCGACGCAGATATTCCTAAGATTGCGCTGATGCCGCAAAGTACCAGCGATGTGCTGACGTTGCTTGCCGCGAcc                                                                                    <  1:4004200/150‑1 (MQ=255)
                                                                      gTCCAACCATGACTTCCATAAAACGCCGGAAGCCGAAGAAATCATTGCCCGTCTGCGCAAAATGCAATCCTTCGACGCAGATATTCCTAAGATTGCGCTGATGCCGCAAAGTACCAGCGATGTGCTGACGTTGCTTGCCGCGACCCTGga                                                                               >  2:1858527/1‑150 (MQ=255)
                                                                        ccAACCATGACTTCCATAAAACGCCGGAAGCCGAAGAAATCATTGCCCGTCTGCGCAAAATGCAATCCTTCGACGCCGATATTCCTAAGATTGCGCTGATGCCGCAAAGTACCAGCGATGTGCTGACGTTGCTTGCCGCGACCCTGgaga                                                                             <  2:4170448/150‑1 (MQ=255)
                                                                         cAACCATGACTTCCATAAAACGCCGGAAGCCGAAGAAATCATTGCCCGTCTGCGCAAAATGCAATCCTTCGACGCCGATATTCCTAAGATTGCGCTGATGCCGCAAAGTACCAGCGATGTGCTGACGTTGCTTGCCGCGACCCTGGAGAt                                                                            >  2:2155965/1‑150 (MQ=255)
                                                                          aaCCATGACTTCCATAAAACGCCGGAAGCCGAAGAAATCATTGCCCGTCTGCGCAAAATGCAATCCTTCGACGCCGATATTCCTAAGATTGCGCTGATGCCGCAAAGTACCAGCGATGTGCTGACGTTGCTTGCCGCGACCCTGGAGATg                                                                           >  2:1032177/1‑150 (MQ=255)
                                                                            ccATGACTTCCATAAAACGCCGGAAGCCGAAGAAATCATTGCCCGTCTGCGCAAAATGCAATCCTTCGACGCCGATATTCCTAAGATTGCGCTGATGCCGCAAAGTACCAGCGATGTGCTGACGTTGCTTGCCGCGACCCTGGAGATGCa                                                                         <  1:3026525/150‑1 (MQ=255)
                                                                                  cTTCCATAAAACGCCGGAAGCCGAAGAAATCATTGCCCGTCTGCGCAAAATGCAATCCTTCGACGCCGATATTCCTAAGATTGCGCTGATGCCGCAAAGTACCAGCGATGTGCTGACGTTGCTTGCCGCGACCCTGGAGATGCAGGAGCa                                                                   >  2:3832151/1‑150 (MQ=255)
                                                                                            aCGCCGGAAGCCGAAGAAATCATTGCCCGTCTGCGCAAAATGCAATCCTTCGACGCAGATATTCCTAAGATTGCGCTGATGCCGCAAAGTACCAGCGATGTGCTGACGTTGCTTGCCGCGACCCTGGAGATGCAGGAGCAGTATGCCGAt                                                         <  2:477820/150‑1 (MQ=255)
                                                                                             cGCCGGAAGCCGAAGAAATCATTGCCCGTCTGCGCAAAATGCAATCCTTCGACGCCGATATTCCTAAGATTGCGCTGATGCCGCAAAGTACCAGCGATGTGCTGACGTTGCTTGCCGCGACCCTGGAGATGCAGGAGCAGTATGCCGATc                                                        <  2:1736738/150‑1 (MQ=255)
                                                                                              gCCGGAAGCCGAAGAAATCATTGCCCGTCTGCGCAAAATGCAATCCTTCGACGCCGATATTCCTAAGATTGCGCTGATGCCGCAAAGTACCAGCGATGTGCTGACGTTGCTTGCCGCGACCCTGGAGATGCAGGAGCAGTATGCCGATCg                                                       >  1:974363/1‑150 (MQ=255)
                                                                                              gCCGGAAGCCGAAGAAATCATTGCCCGTCTGCGCAAAATGCAATCCTTCGACGCCGATATTCCTAAGATTGCGCTGATGCCGCAAAGTACCAGCGATGTGCTGACGTTGCTTGCCGCGACCCTGGAGATGCAGGAGCAGTATGCCGATCg                                                       >  1:3910727/1‑150 (MQ=255)
                                                                                               ccGGAAGCCGAAGAAATCATTGCCCGTCTGCGCAAAATGCAATCCTTCGACGCCGATATTCCTAAGATTGCGCTGATGCCGCAAAGTACCAGCGATGTGCTGACGTTGCTTGCCGCGACCCTGGAGATGCAGGAGCAGTATGCCGATCGt                                                      <  2:2483853/150‑1 (MQ=255)
                                                                                                 ggAAGCCGAAGAAATCATTGCCCGTCTGCGCAAAATGCAATCCTTCGACGCCGATATTCCTAAGATTGCGCTGATGCCGCAAAGTACCAGCGATGTGCTGACGTTGCTTGCCGCGACCCTGGAGATGCAGGAGCAGTATGCCGATCGTcc                                                    <  2:4588249/150‑1 (MQ=255)
                                                                                                  gAAGCCGAAGAAATCATTGCCCGTCTGCGCAAAATGCAATCCTTCGACGCAGATATTCCTAAGATTGCGCTGATGCCGCAAAGTACCAGCGATGTGCTGACGTTGCTTGCCGCGACCCTGGAGATGCAGGAGCAGTATGCCGATCGTCCa                                                   <  1:2388368/150‑1 (MQ=255)
                                                                                                    aGCCGAAGAAATCATTGCCCGTCTGCGCAAAATGCAATCCTTCGACGCCGATATTCCTAAGATTGCGCTGATGCCGCAAAGTACCAGCGATGTGCTGACGTTGCTTGCCGCGACCCTGGAGATGCAGGAGCAGTATGCCGATCGTCCAAt                                                 >  2:1405856/1‑150 (MQ=255)
                                                                                                     gCCGAAGAAATCATTGCCCGTCTGCGCAAAATGCAATCCTTCGACGCCGATATTCCTAAGATTGCGCTGATGCCGCAAAGTACCAGCGATGTGCTGACGTTGCTTGCCGCGACCCTGGAGATGCAGGAGCAGTATGCCGATCGTCCAAtt                                                <  1:3573303/150‑1 (MQ=255)
                                                                                                       cGAAGAAATCATTGCCCGTCTGCGCAAAATGCAATCCTTCGACGCCGATATTCCTAAGATTGCGCTGATGCCGCAAAGTACCAGCGATGTGCTGACGTTGCTTGCCGCGACCCTGGAGATGCAGGAGCAGTATGCCGATCGTCCAATTAt                                              >  1:89441/1‑150 (MQ=255)
                                                                                                       cGAAGAAATCATTGCCCGTCTGCGCAAAATGCAATCCTTCGACGCAGATATTCCTAAGATTGCGCTGATGCCGCAAAGTACCAGCGATGTGCTGACGTTGCTTGCCGCGACCCTGGAGATGCAGGAGCAGTATGCCGATCGTCCAATTAt                                              >  1:3866840/1‑150 (MQ=255)
                                                                                                       cGAAGAAATCATTGCCCGTCTGCGCAAAATGCAATCCTTCGACGCAGATATTCCTAAGATTGCGCTGATGCCGCAAAGTACCAGCGATGTGCTGACGTTGCTTGCCGCGACCCTGGAGATGCAGGAGCAGTATGCCGATCGTCCAATTAt                                              >  2:4460502/1‑150 (MQ=255)
                                                                                                       cGAAGAAATCATTGCCCGTCTGCGCAAAATGCAATCCTTCGACGCAGATATTCCTAAGATTGCGCTGATGCCGCAAAGTACCAGCGATGTGCTGACGTTGCTTGCCGCGACCCTGGAGATGCAGGAGCAGTATGCCGATCGTCCAATTAt                                              >  1:3866917/1‑150 (MQ=255)
                                                                                                        gaagaaATCATTGCCCGTCTGCGCAAAATGCAATCCTTCGACGCCGATATTCCTAAGATTGCGCTGATGCCGCAAAGTACCAGCGATGTGCTGACGTTGCTTGCCGCGACCCTGGAGATGCAGGAGCAGTATGCCGATCGTCCAATTATc                                             >  2:4301739/1‑150 (MQ=255)
                                                                                                         aagaaATCATTGCCCGTCTGCGCAAAATGCAATCCTTCGACGCCGATATTCCTAAGATTGCGCTGATGCCGCAAAGTACCAGCGATGTGCTGACGTTGCTTGCCGCGACCCTGGAGATGCAGGAGCAGTATGCCGATCGTCCAATTATCa                                            <  1:1191812/150‑1 (MQ=255)
                                                                                                          agaaATCATTGCCCGTCTGCGCAAAATGCAATCCTTCGACGCAGATATTCCTAAGATTGCGCTGATGCCGCAAAGTACCAGCGATGTGCTGACGTTGCTTGCCGCGACCCTGGAGATGCAGGAGCAGTATGCCGATCGTCCAATTATCAc                                           <  1:863505/150‑1 (MQ=255)
                                                                                                           gaaATCATTGCCCGTCTGCGCAAAATGCAATCCTTCGACGCCGATATTCCTAAGATTGCGCTGATGCCGCAAAGTACCAGCGATGTGCTGACGTTGCTTGCCGCGACCCTGGAGATGCAGGAGCAGTATGCCGATCGTCCAATTATCACg                                          >  1:3238229/1‑150 (MQ=255)
                                                                                                             aaTCATTGCCCGTCTGCGCAAAATGCAATCCTTCGACGCCGATATTCCTAAGATTGCGCTGATGCCGCAAAGTACCAGCGATGTGCTGACGTTGCTTGCCGCGACCCTGGAGATGCAGGAGCAGTATGCCGATCGTCCAATTATCACGAt                                        >  2:329142/1‑150 (MQ=255)
                                                                                                              aTCATTGCCCGTCTGCGCAAAATGCAATCCTTCGACGCCGATATTCCTAAGATTGCGCTGATGCCGCAAAGTACCAGCTATGTGCTGACGTTGCTTGCCGCGACCCTGGAGATGCAGGAGCAGTATGCCGATCGTCCAATTATCACGATg                                       <  2:2421649/150‑1 (MQ=255)
                                                                                                              aTCATTGCCCGTCTGCGCAAAATGCAATCCTTCGACGCCGATATTCCTAAGATTGCGCTGATGCCGCAAAGTACCAGCGATGTGCTGACGTTGCTTGCCGCGACCCTGGAGATGCAGGAGCAGTATGCCGATCGTCCAATTATCACGATg                                       <  2:208313/150‑1 (MQ=255)
                                                                                                              aTCATTGCCCGTCTGCGCAAAATGCAATCCTTCGACGCCGATATTCCTAAGATTGCGCTGATGCCGCAAAGTACCAGCGATGTGCTGACGTTGCTTGCCGCGACCCTGGAGATGCAGGAGCAGTATGCCGATCGTCCAATTATCACGATg                                       <  2:3775709/150‑1 (MQ=255)
                                                                                                                 aTTGCCCGTCTGCGCAAAATGCAATCCTTCGACGCCGATATTCCTAAGATTGCGCTGATGCCGCAAAGTACCAGCGATGTGCTGACGTTGCTTGCCGCGACCCTGGAGATGCAGGAGCAGTATGCCGATCGTCCAATTATCACGATGTCg                                    <  1:4538496/150‑1 (MQ=255)
                                                                                                                 aTTGCCCGTCTGCGCAAAATGCAATCCTTCGACGCCGATATTCCTAAGATTGCGCTGATGCCGCAAAGTACCAGCGATGTGCTGACGTTGCTTGCCGCGACCCTGGAGATGCAGGAGCAGTATGCCGATCGTCCAATTATCACGATGTCg                                    <  1:2205409/150‑1 (MQ=255)
                                                                                                                  ttGCCCGTCTGCGCAAAATGCAATCCTTCGACGCCGATATTCCTAAGATTGCGCTGATGCCGCAAAGTACCAGCGATGTGCTGACGTTGCTTGCCGCGACCCTGGAGATGCAGGAGCAGTATGCCGATCGTCCAATTATCACGATGTCGa                                   >  2:1983642/1‑150 (MQ=255)
                                                                                                                  ttGCCCGTCTGCGCAAAATGCAATCCTTCGACGCCGATATTCCTAAGATTGCGCTGATGCCGCAAAGTACCAGCGATGTGCTGACGTTGCTTGCCGCGACCCTGGAGATGCAGGAGCAGTATGCCGATCGTCCAATTATCACGATGTCGa                                   <  2:30370/150‑1 (MQ=255)
                                                                                                                   tGCCCGTCTGCGCAAAATGCAATCCTTCGACGCCGATATTCCTAAGATTGCGCTGATGCCGCAAAGTACCAGCGATGTGCTGACGTTGCTTGCCGCGACCCTGGAGATGCAGGAGCAGTATGCCGATCGTCCAATTATCACGATGTCGAt                                  <  1:692875/150‑1 (MQ=255)
                                                                                                                          cTGCGCAAAATGCAATCCTTCGACGCCGATATTCCTAAGATTGCGCTGATGCCGCAAAGTACCAGCGATGTGCTGACGTTGCTTGCCGCGACCCTGGAGATGCAGGAGCAGTATGCCGATCGTCCAATTATCACGATGTCGATGGCaaaa                           >  1:2642165/1‑150 (MQ=255)
                                                                                                                          cTGCGCAAAATGCAATCCTTCGACGCCGATATTCCTAAGATTGCGCTGATGCCGCAAAGTACCAGCGATGTGCTGACGTTGCTTGCCGCGACCCTGGAGATGCAGGAGCAGTATGCCGATCGTCCAATTATCACGATGTCGATGGCaaaa                           >  1:3837717/1‑150 (MQ=255)
                                                                                                                           tGCGCAAAATGCAATCCTTCGACGCCGATATTCCTAAGATTGCGCTGATGCCGCAAAGTACCAGCGATGTGCTGACGTTGCTTGCCGCGACCCTGGAGATGCAGGAGCAGTATGCCGATCGTCCAATTATCACGATGTCGATGGCaaaaa                          >  2:2342841/1‑150 (MQ=255)
                                                                                                                            gcgcAAAATGCAATCCTTCGACGCCGATATTCCTAAGATTGCGCTGATGCCGCAAAGTACCAGCGATGTGCTGACGTTGCTTGCCGCGACCCTGGAGATGCAGGAGCAGTATGCCGATCGTCCAATTATCACGATGTCGATGGCAAAAAc                         <  1:3921782/150‑1 (MQ=255)
                                                                                                                             cgcAAAATGCAATCCTTCGACGCCGATATTCCTAAGATTGCGCTGATGCCGCAAAGTACCAGCGATGTGCTGACGTTGCTTGCCGCGACCCTGGAGATGCAGGAGCAGTATGCCGATCGTCCAATTATCACGATGTCGATGGCAAAAACt                        >  2:3148668/1‑150 (MQ=255)
                                                                                                                                aaaaTGCAATCCTTCGACGCCGATATTCCTAAGATTGCGCTGATGCCGCAAAGTACCAGCGATGTGCTGACGTTGCTTGCCGCGACCCTGGAGATGCAGGAGCAGTATGCCGATCGTCCAATTATCACGATGTCGATGGCAAAAACTGGc                     <  2:2172119/150‑1 (MQ=255)
                                                                                                                                    tGCAATCCTTCGACGCCGATATTCCTAAGATTGCGCTGATGCCGCAAAGTACCAGCGATGTGCTGACGTTGCTTGCCGCGACCCTGGAGATGCAGGAGCAGTATGCCGATCGTCCAATTATCACGATGTCGATGGCAAAAACTGGCGTaa                 >  2:1300553/1‑150 (MQ=255)
                                                                                                                                     gCAATCCTTCGACGCCGATATTCCTAAGATTGCGCTGATGCCGCAAAGTACCAGCGATGTGCTGACGTTGCTTGCCGCGACCCTGGAGATGCAGGAGCAGTATGCCGATCGTCCAATTATCACGATGTCGATGGCAAAAACTGGCGTAAt                <  1:1753523/150‑1 (MQ=255)
                                                                                                                                       aaTCCTTCGACGCAGATATTCCTAAGATTGCGCTGATGCCGCAAAGTACCAGCGATGTGCTGACGTTGCTTGCCGCGACCCTGGAGATGCAGGAGCAGTATGCCGATCGTCCAATTATCACGATGTCGATGGCAAAAACTGGCGTAAttt              >  2:1175213/1‑150 (MQ=255)
                                                                                                                                       aaTCCTTCGACGCAGATATTCCTAAGATTGCGCTGATGCCGCAAAGTACCAGCGATGTGCTGACGTTGCTTGCCGCGACCCTGGAGATGCAGGAGCAGTATGCCGATCGTCCAATTATCACGATGTCGATGGCAAAAACTGGCGTAAttt              >  2:716986/1‑150 (MQ=255)
                                                                                                                                         tCCTTCGACGCCGATATTCCTAAGATTGCGCTGATGCCGCAAAGTACCAGCGATGTGCTGACGTTGCTTGCCGCGACCCTGGAGATGCAGGAGCAGTATGCCGATCGTCCAATTATCACGATGTCGATGGCAAAAACTGGCGTAATTtct            >  2:36904/1‑150 (MQ=255)
                                                                                                                                            ttCGACGCCGATATTCCTAAGATTGCGCTGATGCCGCAAAGTACCAGCGATGTGCTGACGTTGCTTGCCGCGACCCTGGAGATGCAGGAGCAGTATGCCGATCGTCCAATTATCACGATGTCGATGGCAAAAACTGGCGTAATTTCTCGt         >  2:4502807/1‑150 (MQ=255)
                                                                                                                                            ttCGACGCCGATATTCCTAAGATTGCGCTGATGCCGCAAAGTACCAGCGATGTGCTGACGTTGCTTGCCGCGACCCTGGAGATGCAGGAGCAGTATGCCGATCGTCCAATTATCACGATGTCGATGGCAAAAACTGGCGTAATTTCTCGt         <  1:2909042/150‑1 (MQ=255)
                                                                                                                                              cGACGCCGATATTCCTAAGATTGCGCTGATGCCGCAAAGTACCAGCGATGTGCTGACGTTGCTTGCCGCGACCCTGGAGATGCAGGAGCAGTATGCCGATCGTCCAATTATCACGATGTCGATGGCAAAAACTGGCGTAATTTCTCGTct       <  2:1844461/150‑1 (MQ=255)
                                                                                                                                              cGACGCAGATATTCCTAAGATTGCGCTGATGCCGCAAAGTACCAGCGATGTGCTGACGTTGCTTGCCGCGACCCTGGAGATGCAGGAGCAGTATGCCGATCGTCCAATTATCACGATGTCGATGGCAAAAACTGGCGTAATTTCTCGTct       <  2:3058812/150‑1 (MQ=255)
                                                                                                                                               gACGCCGATATTCCTAAGATTGCGCTGATGCCGCAAAGTACCAGCGATGTGCTGACGTTGCTTGCCGCGACCCTGGAGATGCAGGAGCAGTATGCCGATCGTCCAATTATCACGATGTCGATGGCAAAAACTGGCGTAATTTCTCGTctg      >  1:2173330/1‑150 (MQ=255)
                                                                                                                                               gACGCCGATATTCCTAAGATTGCGCTGATGCCGCAAAGTACCAGCGATGTGCTGACGTTGCTTGCCGCGACCCTGGAGATGCAGGAGCAGTATGCCGATCGTCCAATTATCACGATGTCGATGGCAAAAACTGGCGTAATTTCTCGTctg      >  1:3364688/1‑150 (MQ=255)
                                                                                                                                               gACGCCGATATTCCTAAGATTGCGCTGATGCCGCAAAGTACCAGCGATGTGCTGACGTTGCTTGCCGCGACCCTGGAGATGCAGGAGCAGTATGCCGATCGTCCAATTATCACGATGTCGATGGCAAAAACTGGCGTAATTACTCGTctg      >  1:1990555/1‑150 (MQ=255)
                                                                                                                                                  gCCGATATTCCTAAGATTGCGCTGATGCCGCAAAGTACCAGCGATGTGCTGACGTTGCTTGCCGCGACCCTGGAGATGCAGGAGCAGTATGCCGATCGTCCAATTATCACGATGTCGATGGCAAAAACTGGCGTAATTTCTCGTctggct   <  1:1299378/150‑1 (MQ=255)
                                                                                                                                                  gCCGATATTCCTAAGATTGCGCTGATGCCGCAAAGTACCAGCGATGTGCTGACGTTGCTTGCCGCGACCCTGGAGATGCAGGAGCAGTATGCCGATCGTCCAATTATCACGATGTCGATGGCAAAAACTGGCGTAATTTCTCGTctggct   <  1:3450530/150‑1 (MQ=255)
                                                                                                                                                   ccGATATTCCTAAGATTGCGCTGATGCCGCAAAGTACCAGCGATGTGCTGACGTTGCTTGCCGCGACCCTGGAGATGCAGGAGCAGTATGCCGATCGTCCAATTATCACGATGTCGATGGCAAAAACTGGCGTAATTTCTCGTctggctg  >  2:4479772/1‑150 (MQ=255)
                                                                                                                                                   ccGATATTCCTAAGATTGCGCTGATGCCGCAAAGTACCAGCGATGTGCTGACGTTGCTTGCCGCGACCCTGGAGATGCAGGAGCAGTATGCCGATCGTCCAATTATCACGATGTCGATGGCAAAAACTGGCGTAATTTCTCGTctggctg  >  2:150829/1‑150 (MQ=255)
                                                                                                                                                   ccGATATTCCTAAGATTGCGCTGATGCCGCAAAGTACCAGCGATGTGCTGACGTTGCTTGCCGCGACCCTGGAGATGCAGGAGCAGTATGCCGATCGTCCAATTATCACGATGTCGATGGCAAAAACTGGCGTAATTTCTCGTctggctg  >  2:150339/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
TATTTACCGGTGATGATCAGGTTAAAGAAACCGTCGCCTACGCCCACGCGCATGATGTGAAAGTAGTCATGTCCAACCATGACTTCCATAAAACGCCGGAAGCCGAAGAAATCATTGCCCGTCTGCGCAAAATGCAATCCTTCGACGCCGATATTCCTAAGATTGCGCTGATGCCGCAAAGTACCAGCGATGTGCTGACGTTGCTTGCCGCGACCCTGGAGATGCAGGAGCAGTATGCCGATCGTCCAATTATCACGATGTCGATGGCAAAAACTGGCGTAATTTCTCGTCTGGCTG  >  CP009273/1769292‑1769588

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: