Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,512,114 C→A 9.3% E104* (GAA→TAA)  aroK ← shikimate kinase I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,512,1140CA 9.3% 519.1 / 4.5 219E104* (GAA→TAA) aroKshikimate kinase I
Reads supporting (aligned to +/- strand):  ref base C (76/109);  new base A (16/12);  total (98/121)
Fisher's exact test for biased strand distribution p-value = 1.51e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 224 total aligned reads displayed.

ACGAATGGTCACGTCGGCAATCTCTTCATACAGCGGATTGCGTTCATTGGCCAACGCTTCCAGAACTTCACGCGGCGGTGTTTCAACGTGCAGCAACGGGCGTTTTTTATCACGCTGCGTGCGTGCAAGTTGCTTTTCGATGGTCGTTTCAAGATAAACGACAACGCCACGAGCGGAAAGACGGTTACGCGTTTCACGGGATTTCACAGAGCCGCCGCCAGTAGCCAGCACAATACCCTGTTTCTCGGTCAACTCATTGATGACCTTTTCTTCGCGATCGCGGAAGCCTTCTTCGCCTT  >  CP009273/3511965‑3512263
                                                                                                                                                     |                                                                                                                                                     
aCGAATGGTCACGTCGGCAATCTCTTCATACAGCGGATTGCGTTCATTGGCCAACGCTTCCAGAACTTCACGCGGCGGTGTTTCAACGTGCAGCAACGGGCGTTTTTTATCACGCTGCGTGCGTGCAAGTTGCTTTTCGATGGTCGTTTc                                                                                                                                                       >  2:403754/1‑150 (MQ=255)
   aaTGGTCTCGTCGGCAATCTCTTCATACAGCGGATTGCGTTCATTGGCCAACGCTTCCAGAACTTCACGCGGCGGTGTTTCATCGTGCAGCAACGGGCGTTTTTTATCACGCTGCGTGCGTGCAAGTTGCTTTTCGATGGTCGTTTCAAg                                                                                                                                                    <  1:1554245/150‑1 (MQ=255)
    aTGGTCACGTCGGCAATCTCTTCATACAGCGGATTGCGTTCATTGGCCAACGCTTCCAGAACTTCACGCGGCGGTGTTTCAACGTGCAGCAACGGGCGTTTTTTATAACGCTGCGTGCGTGCAAGTTGCTTTTCGATGGTCGTTTCAAGa                                                                                                                                                   >  1:302632/1‑150 (MQ=255)
      ggTCACGTCGGCAATCTCTTCATACAGCGGATTGCGTTCATTGGCCAACGCTTCCAGAACTTCACGCGGCGGTGTTTCAACGTGCAGCAACGGGCGTTTTTTATCACGCTGCGTGCGTGCAAGTTGCTTTTCGATGGTCGTTTCAAGATa                                                                                                                                                 <  1:136063/150‑1 (MQ=255)
      ggTCACGTCGGCAATCTCTTCATACAGCGGATTGCGTTCATTGGCCAACGCTTCCAGAACTTCACGCGGCGGTGTTTCAACGTGCAGCAACGGGCGTTTTTTATCACGCTGCGTGCGTGCAAGTTGCTTTTCGATGGTCGTTTCAAGATa                                                                                                                                                 >  2:2960888/1‑150 (MQ=255)
      ggTCACGTCGGCAATCTCTTCATACAGCGGATTGCGTTCATTGGCCAACGCTTCCAGAACTTCACGCGGCGGTGTTTCAACGTGCAGCAACGGGCGTTTTTTATAACGCTGCGTGCGTGCAAGTTGCTTTTAGTAGGTCGTTTAAAtata                                                                                                                                                 >  1:474046/1‑150 (MQ=255)
         cACGTCGGCAATCTCTTCATACAGCGGATTGCGTTCATTGGCCAACGCTTCCAGAACTTCACGCGGCGGTGTTTCAACGTGCAGCAACGGGCGTTTTTTATCACGCTGCGTGCGTGCAAGTTGCTTTTCGATGGTCGTTTCAAGATAAAc                                                                                                                                              <  2:4240397/150‑1 (MQ=255)
          aCGTCGGCAATCTCTTCATACAGCGGATTGCGTTCATTGGCCAACGCTTCCAGAACTTCACGCGGCGGTGTTTCAACGTGCAGCAACGGGCGTTTTTTATCACGCTGCGTGCGTGCAAGTTGCTTTTCGATGGTCGTTTCAAGATAAACg                                                                                                                                             <  2:1310858/150‑1 (MQ=255)
           cGTCGGCAATCTCTTCATACAGCGGATTGCGTTCATTGGCCAACGCTTCCAGAACTTCACGCGGCGGTGTTTCAACGTGCAGCAACGGGCGTTTTTTATCACGCTGCGTGCGTGCAAGTTGCTTTTCGATGGTCGTTTCAAGATAAAAGa                                                                                                                                            >  2:3438604/1‑150 (MQ=255)
            gTCGGCAATCTCTTCATACAGCGGATTGCGTTCATTGGCCAACGCTTCCAGAACTTCACGCGGCGGTGTTTCAACGTGCAGCAACGGGCGTTTTTTATCACGCTGCGTGCGTGCAAGTTGCTTTTCGATGGTCGTTTCAAGATAAACGAc                                                                                                                                           >  2:1079486/1‑150 (MQ=255)
            gTCGGCAATCTCTTCATACAGCGGATTGCGTTCATTGGCCAACGCTTCCAGAACTTCACGCGGCGGTGTTTCAACGTGCAGCAACGGGCGTTTTTTATCACGCTGCGTGCGTGCAAGTTGCTTTTCGATGGTCGTTTCAAGATAAACGAc                                                                                                                                           <  2:4265000/150‑1 (MQ=255)
               ggCAATCTCTTCATACAGCGGATTGCGTTCATTGGCCAACGCTTCCAGAACTTCACGCGGCGGTGTTTCAACGTGCAGCACCGGGCGTTTTTTATCCCGCTGCGTGCGTGCAAGTTGCTTTTCGATGGTCGTTTCAAGATAAACGACAAc                                                                                                                                        >  2:1427633/1‑150 (MQ=255)
                  aaTCTCTTCATACAGCGGATTGCGTTCATTGGCCAACGCTTCCAGAACTTCACGCGGCGGTGTTTCAACGTGCAGCAACGGGCGTTTTTTATCACGCTGCGTGCGTGCAAGTTGCTTTTCGATGGTCGTTTCAAGATAAACGACAACGcc                                                                                                                                     >  1:4799838/1‑150 (MQ=255)
                  aaTCTCTTCATACAGCGGATTGCGTTCATTGGCCAACGCTTCCAGAACTTCACGCGGCGGTGTTTCAACGTGCAGCAACGGGCGTTTTTTATAACGCTGCGTGCGTGCAAGTTGCTTTTCGATGGTCGTTTCAAGATAAACGACAAcccc                                                                                                                                     >  1:1109378/1‑150 (MQ=255)
                  aaTCTATTCATACAGCGTATTGCGTTCATTGGCAAACGCTTCCATAACTTCACGCGGCGGTGTTTCAACGTGCAGCAACGGGCGTTTTTTATCACGCTGCGTGCGTGCAAGTTGCTTTTCGATGGTCGTTTCAAGATAAACGACAACGcc                                                                                                                                     <  2:3087399/150‑1 (MQ=255)
                    tctcTTCATACAGCGGATTGCGTTCATTGGCCAACGCTTCCAGAACTTCACGCGGCGGTGTTTCAACGTGCAGCAACGGGCGTTTTTTATAAAGCTGCGTGCGTGCAAGTTGCTTTTAGATGGTCGTTTAAATTTTAAAGACAAAGCCAc                                                                                                                                   >  1:4433177/1‑150 (MQ=255)
                       ctTCATACAGCGGATTGCGTTCATTGGCCAACGCTTCCAGAACTTCACGCGGCGGTGTTTCAACGTGCAGCAACGGGCGTTTTTTATAACGCTGCGTGCGTGCAAGTTGCTTTTCGATTGTCGTTTCAAGATAAACGACAACGCCACGAg                                                                                                                                >  1:4814769/1‑150 (MQ=255)
                        ttCATACAGCGGATTGCGTTCATTGGCCAACGCTTCCAGAACTTCACGCGGCGGTGTTTCAACGTGCAGCAACGGGCGTTTTTTATCACGCTGCGTGCGTGCAAGTTGCTTTTCGATGGTCGTTTAAAGATAAACGACAACGCCACGAGc                                                                                                                               >  2:140860/1‑150 (MQ=255)
                          cATACAGCGGATTGCGTTCATTGGCCAACGCTTCCAGAACTTCACGCGGCGGTGTTTCAACGTGCAGCAACGGGCGTTTTTTATCACGCTGCGTGCGTGCAAGTTGCTTTTAGATTGTCGTTTCAAGATAAACGACAACCCCCCCAGCgg                                                                                                                             >  1:239716/1‑150 (MQ=255)
                           aTACAGCGGATTGCGTTCATTGGCCAACGCTTCCAGAACTTCACGCGGCGGTGTTTCAACGTGCAGCAACGGGCGTTTTTTATCACGCTGCGTGCGTGCAAGTTGCTTTTCGATGGTCGTTTCAAGATAAACGACAACGCCACGAGCGGa                                                                                                                            >  2:634083/1‑150 (MQ=255)
                            tACAGCGGATTGCGTTCATTGGCCAACGCTTCCAGAACTTCACGCGGCGGTGTTTCAACGTGCAGCAACGGGCGTTTTTTATCACGCTGCGTGCGTTCAAATTTCTTTTCGATGGTCGTTTCAATAAAAAAGACAAACCCACCCGCggca                                                                                                                           >  1:3766301/1‑148 (MQ=255)
                             aCAGCGGATTGCGTTCATTGGCCAACGCTTCCAGAACTTCACGCGGCGGTGTTTCAACGTGCAGCAACGGGCGTTTTTTATCACGCTGCGTGCGTGCAAGTTGCTTTTCGATGGTCGTTTCAAGATAAACGACAACGCCACGAGCGGaaa                                                                                                                          <  2:1734951/150‑1 (MQ=255)
                              cAGCGGATTGCGTTCATTGGCCAACGCTTCCAGAACTTCACGCGGCGGTGTTTCAACGTGCAGCAACGGGCGTTTTTTATCACGCTGCGTGCGTGCAAGTTGCTTTTCGAGGGTCGTTTCAAGATAAACGACAACGCCACGAGCGGAAAg                                                                                                                         <  2:274927/150‑1 (MQ=255)
                               aGCGGATTGCGTTCATTGGCCAACGCTTCCAGAACTTCACGCGGCGGTGTTTCAACGTGCAGCAACGGGCGTGTTTTATAACGCTGCGTGCGTGCAAGTTGCTTTTAGAAGGTCGGTTAAACATAAAAGACAACCCCACCAGCGGaaaaa                                                                                                                        >  1:2935636/1‑148 (MQ=255)
                                gCGGATTGCGTTCATTGGCCAACGCTTCCAGAACTTCACGCGGCGGTGTTTCAACGTGCAGCAACGGGCGTTTTTTATCACGCTGCGTGCGTGCAAGTTGCTTTTCGATGGTCGTTTCAAGAAAAACGACAACGCCACGAGCGGAAAGAc                                                                                                                       <  1:1965658/150‑1 (MQ=255)
                                    aTTGCGTTCATTGGCCAACGCTTCCAGAACTTCACGCGGCGGTGTTTCAACGTGCAGCAACGGGCGTTTTTTATAACGCTGCGTGCGTGCAAGTTGCTTTTCGATGGTCGTTTCAATATAAAAGAAAAAGGCAACAGAGGAAAGaacgtt                                                                                                                   >  2:1055661/1‑145 (MQ=255)
                                     ttGCGTTCATTGGCCAACGCTTCCAGAACTTCACGCGGCGGTGTTTCAACGTGCAGCAACGGGCGTTTTTTATCACGCTGCGTGCGTGCAAGTTGCTTTTCGATGGTCGTTTCAAGATAAACGACAACGCCACGAGCGGAAAGACGGTTa                                                                                                                  <  2:2548605/150‑1 (MQ=255)
                                      tGCGTTCATTGGCCAACGCTTCCAGAACTTCACGCGGCGGTGTTTCAACGTGCAGCAACGGGCGTTTTTTATCACGCTGCGTGCGTGCCAGTTGCTTTTCGATTGTCGTTTCAATATAAAAGAAAACGCCAAGAGCGGAAAGACGGTTAc                                                                                                                 >  2:1688367/1‑150 (MQ=255)
                                         gttatttGGCCAACGCTTCCATAACTTCCCGGGGCGGTTTTTAAACGTGCAGCAAGGGGCTTTTTTTATCCCGCTGCGTGCGTGAATTTTGCTTTTCGTTGGTCGTTTAAAGATAAACGACAACGCCACGAGCGGAAAGACGGTTAcgcg                                                                                                              <  1:4331390/145‑1 (MQ=255)
                                           tCATTGGCCAACGCTTCCTGAACTGCACGCGGGGGTGTTTCAACGTGCAGCAACGGGCGTTTTTTATCACGCTGCGTGCGTGCAAGTTGCTTTTCGATGGTCGTTTCATGATAAACGACAACGCCACGAGCGGAAAGACGGTTACGCGtt                                                                                                            <  1:4011821/150‑1 (MQ=255)
                                           tCATTGGCCAACGCTTCCAGAACTTCACGCGGCGGTTTTTCAACGTGCAGCAACGGGCTTTTTTTATCACGCTGCGTGCGTGAAAGTTGCTTTTCGAGGGTCGTTTCAAGATAAACGACAACGCCACGAGCGGAAAGACGGTTACGCGtt                                                                                                            <  1:3186932/150‑1 (MQ=255)
                                           tCATTGGCCAACGCTTCCAGAACTTCACGCGGCGGTGTTTCAACGTGCAGCAACGGGCGTTTTTTATCACGCTGCGTGCGTGCAAGTTGCTTTTCGATGGTCGTTTCAAGATAAAACACAACCCCACGAGCGGAAAGACGGTTACGCGtt                                                                                                            >  1:3853475/1‑150 (MQ=255)
                                            cATTGGCCAACGCTTCCAGAACTTCACGCGGCGGTGTTTCAACGTGCAGCAACGGGCGTTTTTTATCACGCTGCGTGCGTGCAAGTTGCTTTTCGATGGTCGTTTCAAGAAAAACGACAACGCCACGAGCGGAAAGACGGTTACGCGttt                                                                                                           <  2:1823480/150‑1 (MQ=255)
                                             aTTGGCCAACGCTTCCAGAACTTCACGCGGCGGTGTTTCAACGTGCAGCAACGGGCGTTTTTTATCACGCTGCGTGCGTGCAAGTTGCTTTTCGATGGTCGTTTCAAGATAAACGACAACGCCACGAGCGGAAAGACGGTTACGCGTTTc                                                                                                          <  2:256057/150‑1 (MQ=255)
                                              ttGGCCACCGCTTCCAGAACTTCACGCGGCGGTGTTTCACCGTGCAGCAACGGGCGTTTTTTATCACGCTGCGTGCGTGCATGTTGCTTTTCGATGGTCGTTTCAAGATACACGACAACGCCACGAGCGGAAAGACGGTTACGCGTTTCa                                                                                                         <  2:3286899/150‑1 (MQ=255)
                                              ttGGCCAACGCTTCCAGAACTTCACGCGGCGGTGTTTCAACGTGCAGCAACGGGCGTTTTTTATCACGCTGCGTGCGTGCAAGTTGCGTTTCGATGGTCGTTTCAATATAAACGACAACGCCACGAGCGGAAAGACGGTTACGCGTTTCa                                                                                                         >  1:2807740/1‑150 (MQ=255)
                                                  ccAACGCTTCCAGAACTTCACGCGGCGGTGTTTCAACGTGCAGCAACGGGCGTTTTTTATCACGCTGCGTGCGTGCAAGTTGCTTTTCGATGGTCGTTTCAAGATAAACGAAAACCCCACCAGCGGAAAGACGGTTAAGCGTTTAACggg                                                                                                     >  2:3918896/1‑150 (MQ=255)
                                                  ccAACGCTTCCAGAACTTCACGCGGCGGTGTTTCAACGTGCAGCAACGGGCGTTTTTTATCACGCTGCGTGCGTGCAAGTTGCTTTTCGATGGTCGTTTCAAGATAAAAGACAACGCCACGAGCGGAAAGACGGTTACGCGTTTCACggg                                                                                                     >  2:3921578/1‑150 (MQ=255)
                                                   cAACGCTTCCAGAACTTCACGCGGCGGTGTTTCAACGTGCAGCAACGGGCGTTTTTTTTAACGCTGCGTGGGTGGAAGTTGCTTTTCGATTTTAGTTTTAAAATAAAAAAAAAACCCACCAGCGGAAAAAAAGTTTAGCGTTTTAAGGGa                                                                                                    >  1:3733919/1‑150 (MQ=255)
                                                    aaCGCTTCCAGAACTTCACGCGGCGGTGTTTCAACGTGCAGCAACGGGCGTTTTTTATCACGCTGCGTGCGTGCAAGTTGCTTTTCGATGGTCGTTTCAAGATAAACGACAACGCCACGAGCGGAAAGACGGTTACGCGTTTAAAGGGAt                                                                                                   >  2:2076508/1‑150 (MQ=255)
                                                    aaCGCTTCCAGAACTTCACGCGGCGGTGTTTCAACGTGCAGCAACGGGCGTTTTTTATCACGCTGCGTGCGTGCAAGTTGCTTTTAGATGGTCGTTTAAATATAAACGACAACGCCACGAGCGGAAAGACGGTTACGCGTTTCAAGGGAt                                                                                                   >  1:1853397/1‑150 (MQ=255)
                                                     aCGCTTCCAGAACTTCACGCGGCGGTGTTTCAACGTGCAGCAACGGGCGTTTTTTATCACGCTGCGTGCGTGCAAGTTGCTTTTCGATGGTCGTTTCACGATAAACGACAACGCCACGAGCGGAAAGACGGTTACGCGTTTCACGGGAtt                                                                                                  <  2:2499046/150‑1 (MQ=255)
                                                      cGCTTCCAGAACTTCCCGCGGCGGTGTTTCAACGTGCAGCAACGGGCGTTTTTTATCACGCTGCGTGCGTGCAAGTTGCTTTTCGATGGTCGTTTCAAGATAAACGACAACGCCACGAGCGGAAAGACGGTTACGCGTTTCACGGGAttt                                                                                                 <  2:4729439/150‑1 (MQ=255)
                                                       gCTTCCAGAACTTCACGCGGCGGTGTTTCAACGTGCAGCAACGGGCGTTTTTTATCACGCTGCGTGCGTGCATGTTGCTTTTCGATGGTCGTTTCAAGATAAACGACAACGCCACGAGCGGAAAGACGGTTACGCGTTTCACGGGATTTc                                                                                                <  1:1902995/150‑1 (MQ=255)
                                                        cTTCCAGAACTTCACGCGGCGGTGTTTCAACGTGCAGCAACGGGCGTTTTTTATCACGCTGCGTGCGTGCAAGTTGCTTTTCTATGGTCGTTTAAATATAAAAAACAAACCCACCAGCCGCAAGACGGTTTAGCGGTTCAAGGGAtttaa                                                                                               >  1:4242637/1‑148 (MQ=255)
                                                          tCCAGAACTTCCGGGGGGGGTGTTACAACGTGCAGCACGGGGGTTTTTTTATCCCGCGGGGTGCGTGCATGTTGCTTTTCGTTGGTGGTTTAAATAAAAACGACAACGCCACGAGCGGAAAGACGGTTACGCGTTTCACGGGATTTcaca                                                                                             <  1:3781389/150‑1 (MQ=255)
                                                            cAGAACTTCACGCGGCGGTGTTTCAACGTGCAGCAACGGGCGTTTTTTATCACGCTGCGTGCGTGCAAGTTTCTTTTCGATGGTCGTTTAAATATAAAAGAAAAACCCACCAGCGGAAAGAAGGTTAAGCGTTTAAACGGAtttaaaaca                                                                                           >  1:3096549/1‑144 (MQ=255)
                                                            cAGAACTTCACGCGGCGGTGTTTCAACGTGCAGCAACGGGCGTTTTTTATCACGCTGCGTGCGTGCAAGTTGCTTTTCGATGGTCGTTTCAAGATAAACGACAACGCCACGAGCGGAAAGACGGTTACGCGTTTCACGGGATTTCACaga                                                                                           <  2:3413868/150‑1 (MQ=255)
                                                             aGAACTTCCCGCGGCTGTGTTTCAACGTGCAGCAACGGGCGTTTTTTATCACGCTGCGTGCGTGCAAGGTGCTTTTCGATGGTCGTTTCAAGATAAACGACAACGCCACGAGCGGAAAGACGGTTACGCGTTTCACGGGATTTCACagag                                                                                          <  2:4571612/150‑1 (MQ=255)
                                                                 cTTCACGCGGCGGTGTTTCAACGTGCAGCAACGGGCGTTTTTTATCACGCTGCGTGCGTGCAAGTTGCTTTTCGATGGGCGTTTCAAGATAAACGACAACGCCACGAGCGGAAAGACGGTTCAGCGGTTCACGGGATTTCACAGAgccgc                                                                                      >  2:4083588/1‑150 (MQ=255)
                                                                        ggggggTGTTTAACGGTGCAGCAAGGGTGTTTTTTTATCTCGCTGGGTGGGTGAATGTTGTTTTTCTTTGGTATTTAAAACAAAAACGACAACGCCACGAGCGGAAAGACGGTTACGCGTTTCACGGGATTTCACAGAGCCGCCGCCAGt                                                                               <  1:4106547/146‑1 (MQ=255)
                                                                        cggcggTGTTTCAACGTGCAGCAACGGGCGTTTTTTATCACGCTGCGTGCGTGCAAGTTGCTTTTCGATGGTCGTTTCAAGATAAACGACAACGCCACGAGCGGAAAGACGGTTACGCGTTTCACGGGATTTCACAGAGCCGCCGCCAGt                                                                               <  1:2829147/150‑1 (MQ=255)
                                                                         ggcggTGTTTCAACGTGCAGCAACGGGCGTTTTTTATAACGCTGCGTGCGTGCAAGTTGCTTTTCGATGGTCGTTTCAAGATAAAACACAAACCCACCAGCGGAAAGACGGTTACGCGGTTCACCGGAGTTGACAGAgccccccccccga                                                                              >  1:2672538/1‑146 (MQ=255)
                                                                          gcggTGTTTCAACGTGCAGCAAGGGGCGTTTTTTATCACGCTGCTTGCGTGCAAGTTGCTTTTCGATGGTCGTTTCAAGATAAACGACAACGCCACGAGCGGAAAGACGGTTACGCGTTTCACGGGATTTCACAGAGCCGCCGCCAGTAg                                                                             <  2:1790268/150‑1 (MQ=255)
                                                                             gtttttCAACGTGCAGCAACGGGCGTTTTTTATCACGCTGCGTGGGTGCAAGTTGCTTTTCGATGGTGTTTTCAATATAAACGACAACGCCACGAGCGGAAAGACGGTTACGCGTTTCACGGGATTTCACAGAGCCGCCGCCAGTAGCCa                                                                          <  2:904886/147‑1 (MQ=255)
                                                                             gtgtTTCAACGTGCAGCACCGGGCGTTTTTTATCACGCGGCGGGCGTGCAAGTTGCTTTTCGATGGTCTTTTCAAGATAAACGACAACGCCACGAGCGGAAAGACGGTTACGCGTTTCACGGGATTTCACAGAGCCGCCGCCAGTAGCCa                                                                          <  2:3252135/150‑1 (MQ=255)
                                                                                  tCAACGTGCAGCAACGGGCGTTTTTTATCACGCTGCGTGCGTGCAAGTTGCTTTTCGATGGTCGTTTCAAGATAAACGACAACGCCACGAGCGGAAAGACGGTTACGCGTTTCACGGGATTTCACAGAGCCGCCGCCAGTAGCCAGcaca                                                                     <  2:36480/150‑1 (MQ=255)
                                                                                    aaggTGCTGAAACGGGCTTTTTTTATCACGCTGGGTGGTTGATTGTTGGTTTTCTTTTGTCGTTTAAATATAAACGACAACGCCACGAGCGGAAAGACGGTTACGCGTTTCACGGGATTTCACAGAGCCGCCGCCAGTAGCCAGCACAAt                                                                   <  1:3109904/147‑1 (MQ=255)
                                                                                      ctggcagcaACGGGCGTTTTTTATCACGCTGCGTGCGTGCAAGTTGCTTTTCGATGGTCGTTTCAAGATAAACGACAACGCCACGAGCGGAAAGACGGTTACGCGTTTCACGGGATTTCACAGAGCCGCCGCCAGTAGCCAGCACAATAc                                                                 <  2:3683907/147‑1 (MQ=255)
                                                                                      cGTGCAGCAACGGGCGTTTTTTATCACGCTGCGTGCGTGCAAGTTGCTTTTCGATGGTCGTTTCAAGATAAACGACAACGCCACGAGCGGAAAGACGGTTACGCGTTTCACGGGATTTCACAGAGCCGCCGCCAGTAGCCAGCACAATAc                                                                 <  1:1079486/150‑1 (MQ=255)
                                                                                        tGCAGCAACGGGCGTTTTTTATCACGCTGCGTGCGTGCAAGTTGCTTTTCGATGGTCGTTTCAAGATAAACGACAACGCCACGAGCGGAAAGACGGTTACGCGTTTCACGGGATTTCACAGAGCCGCCGCCAGTAGCCAGCACAATAccc                                                               <  2:3348634/150‑1 (MQ=255)
                                                                                        tGCAGCAACGGGCGTTTTTTATCACGCTGCGTGCGTGCAAGTTGCTTTTCGATGGTCGTTTAAACAAAAACGACAAACCCACCCGCGGCAACACCGGTTCGCGTTTAACGGGATTTAAAAGAGCCCCCGCCACGAGCCACCACaaaaacc                                                               >  1:2597424/1‑145 (MQ=255)
                                                                                          cagcaACGGGCGTTTTTTATCACGCTGCGTGCGTGCAAGTTGCTTTTCGATGGTCGTTTCAAGATAAACGACAACGCCACGAGCGGAAAGACGGTTACGCGTTTCACGGGATTTCACAGAGCCGCCGCCAGTAGCCAGCACAATACCCTg                                                             <  1:3149457/150‑1 (MQ=255)
                                                                                           agcaACGGGCGTTTTTTATCACGCTGCGTGCGTGCAAGTTGCTTTTCGAGGGTCGTTTCAAGATAAACGACAACGCCACGAGCGGAAAGACGGTTACGCGTTTCACGGGATTTCACAGAGCCGCCGCCAGTAGCCAGCACAATACCCTGt                                                            <  1:3572776/150‑1 (MQ=255)
                                                                                           agcaACGGGCGTTTTTTAACACGCTGCGTGCGTGCAAGTTGCTTTTCGATGGTCGTTTCAAGATAAACGACAACGCCACGAGCGGAAAGACGGTTACGCGTTTCACGGGATTTCACAGAGCCGCCGCCAGTAGCCAGCACAATACCCTGt                                                            <  1:3570922/150‑1 (MQ=255)
                                                                                            gcaTCGGGGGTTTTTTATCACGCTGCGTGCGTGAATGTTGTTTTTCGTTGGTCGTTTCAAGATAAACGACAACGCCACGAGCGGAAAGACGGTTACGCGTTTCACGGGATTTCACAGAGCCGCCGCCAGTAGCCAGCACAATACCCTGtt                                                           <  1:4545173/150‑1 (MQ=255)
                                                                                               aCGGGCGTTTTTTATCACGCTGCGTGCGTGCAAGTTGCTTTTCGATGGTCGTTTCAAGATAAACGACAACGCCACGAGCGGAAAGACGGTTACGCGTTTCACGGGATTTCACAGAGCCGCCACCAGTAGCCAGCACAATACCCTGTTtct                                                        <  2:4481153/150‑1 (MQ=255)
                                                                                                cGGGCTTTTTTTATCACGCTGCGTGCGTGAAAGTTGCTTTTCGTTGGTCTTTTAAATATAAACGACAACGCCACGAGCGGAAAGACGGTTACGCGTTTCACGGGATTTCACAGAGCCGCCGCCAGTAGCCAGCACAATACCCTGTTtctc                                                       <  1:3469431/150‑1 (MQ=255)
                                                                                                cGGGCGTTTTTTATCACGCTGCGTGCGTGCAAGTTGCTTTTCGATGGTCGTTTCAAGATAAACGACAACGCCACGAGCGGAAAGACGGTTACGCGTTTCACGGGATTTCACAGAGCCGCCGCCAGTAGCCAGCACAATACCCTGTTtctc                                                       <  1:793029/150‑1 (MQ=255)
                                                                                                 gggCTTTTTTTATCACGCTGCGTGCGTGAAATTTGCTTTTCGAGGGTCGTTACAAGATAAACGACAACGCCACGAGCGGAAAGACGGTTACGCGTTTCACGGGATTTCACAGAGCCGCCGCCAGTAGCCAGCACAATACCCTGTTTCTCg                                                      <  2:831567/150‑1 (MQ=255)
                                                                                                 gggCGTTTTTTATCACGCTGCGTGCGTGCAAGTTGCTTTTCGATGGTCGTTTCAAGATAAACGACAACGCCACGAGCGGAAAGACGGTTACGCGTTTCACGGGATTTCACAGAGCCGCCGCCAGTAGCCAGCACAATACCCTGTTTCTCg                                                      <  1:2422972/150‑1 (MQ=255)
                                                                                                   gCGTTTTTTATCACGCTGCGTGCGTGCAAGTTGCTTTTCGATGGTCGTTTCAAGATAAACGACAACGCCACGAGCGGAAAGACGGTTACGCGTTTCACGGGATTTCACAGAGCCGCCGCCAGTAGCCAGCACAATACCCTGTTTCTCGGt                                                    >  1:367328/1‑150 (MQ=255)
                                                                                                    ctttttttATCACGCTGCGTGCGTGCATGTTGATTTTCGATGTTCGTTTCAAGATAAACGACAACGCCACGAGCGGAAAGACGGTTACGCGTTTCACGGGATTTCACAGAGCCGCCGCCAGTAGCCAGCACAATACCCTGTTTCTCGGTc                                                   <  1:4302287/148‑1 (MQ=255)
                                                                                                    cGTTTTTTATCACGCTGCGTGCGTGCAAGTTGCTTTTCGATGGTCGTTTCAAGATAAACGACAACGCCACGAGCGGAAAGACGGTTACGCGTTTCACGGGATTTCACAGAGCCGCCGCCAGTAGCCAGCACAATACCCTGTTTCTCGGTc                                                   <  2:1369869/150‑1 (MQ=255)
                                                                                                          ttATCACGCTGGGGGCGTGAATGTTGTTTTTGGTTGGTCGTTTAAAGATAAACGACAACGCCACGAGCGGAAAGACGGTTACGCGTTTCACGGGATTTCACAGAGCCGCCGCCAGTAGCCAGCACAATACCCTGTTTCTCGGTCAACTCa                                             <  1:2076555/150‑1 (MQ=255)
                                                                                                            aTCACGCTGCGTGCGTGAAAGTTGATTTTCGATGGTCGTTTCAAGATAAACGACAACGCCACGAGCGGAAAGACGGTTACGCGTTTCACGGGATTTCACAGAGCCGCCGCCAGTAGCCAGCACAATACCCTGTTTCTCGGTCAACTCAtt                                           <  1:1476911/150‑1 (MQ=255)
                                                                                                             tCACGCTGCGTGCGTGCAAGTTGCTTTTCGATGGTCGTTTCAAGATAAACGACAACGCCACGAGCGGAAAGACGGTTACGCGTTTCACGGGATTTCACAGAGCCGCCGCCAGTAGCCAGCACAATACCCTGTTTCTCGGTCAACTCATtg                                          <  2:1566903/150‑1 (MQ=255)
                                                                                                             tCACGCTGCGTGCGTGCAAGTTGCTTTTCGATGGTCGTTTCAAGATAAACGACAACGCCACGAGCGGAAAGACGGTTACGCGTTTCACGGGATTTCACAGAGCCGCCGCCAGTAGCCAGCACAATACCCTGTTTCTCGGTCAACTCATtg                                          >  2:2098702/1‑150 (MQ=255)
                                                                                                             tCACGCTGCGTGCGTGCAAGTTGCTTTTCGATGGTCGTTTCAAGATAAAAGACAACGCCACGAGCGGAAAGACGGTTACGCGTTTCACGGGATTTAACAGAGCCGCCGCCAGTAGCCAGCACAATACCCTGTTTCTCGGTCAACTCATtg                                          >  1:2566715/1‑150 (MQ=255)
                                                                                                              cACGCTGCGTGCGTGCAAGTTGCTTTTCGATTGTCGTTTCAAGATAAAACACAACCCCACGAGCCGAAAGAAGGTTACGCGTTTCCCGGGATTTAACAGAGCCGCCCCCCGTAGCCAGCACAAAAACCTGTTTCTCGGTCAACTCagtta                                         >  1:33446/1‑146 (MQ=255)
                                                                                                               aCGCTGCGTGCGTGCAAGTTGCTTTTCGATGGTCGTGTCAAGATAAACGACACAGCCACGAGCGGAAAGACGGTTACGCGTTTCACGGGATTTCACAGAGCCGCCGCCAGTAGCCAGCACAATACCCGGTTTCTCGGTCAACTCATtgat                                        >  2:3371526/1‑150 (MQ=255)
                                                                                                                       tgcgtgcAAGTTGCTTTTCGATGGTCGTTTCAAGATAAACGACAACGCCACGAGCGGAAAGACGGTTACGCGTTTCACGGGATTTCACAGAGCCGCCGCCAGTAGCCAGCACAATACCCTGTTTCTCGGTCAACTCATTGATGACCtttt                                <  1:3532688/150‑1 (MQ=255)
                                                                                                                       tgcgtgcAAGTTGCTTTTCGATGGTCGTTTCAAGATAAACGACAACGCCACGAGCGGAAAGACGGTTACGCGTTTCACGGGATTTCACAGAGCCGCCGCCAGTAGCCAGCACAATACCCTGTTTCTCGGTCAACTCATTGATGACCtttt                                <  2:4592809/150‑1 (MQ=255)
                                                                                                                              aaGTTGCTTTTCGATGGTCGTTTCAAGATAAACGACAACGCCACGAGCGGAAAGACGGTTACGCGTTTCACGGGATTTCACAGAGCCGCCGCCAGTAGCCAGCACAATACCCTGTTTCTCGGTCAACTCATTGATGACCTTTTCTTcgcg                         >  2:781392/1‑150 (MQ=255)
                                                                                                                              aaGTTGCTTTTCGATGGTCGTTTCAAGATAAACGACAACGCCACGAGCGGAAAGACGGTTACGCGTTTCACGGGATTTCACAGAGCCGCCGCCAGTAGCCAGCACAATACCCTGTTTCTCGGTCAACTCATTGATGACCTTTTCTTcgcg                         >  1:2404386/1‑150 (MQ=255)
                                                                                                                                gTTGCTTTTCGATGGTCGTTTCAAGATAAACGACAACGCCACGAGCGGAAAGACGGTTACGCGTTTCACGGGATTTCACAGAGCCGCCGCCAGTAGCCAGCACAATACCCTGTTTCTCGGTCAACTCATTGATGACCTTTTCTTCGCGAt                       <  2:2445329/150‑1 (MQ=255)
                                                                                                                                  tGCTTTTCGATGGTCGTTTCAAGATAAACGACAACGCCACGAGCGGAAAGACGGTTACGCGTTTCACGGGATTTCACAGAGCCGCCGCCAGTAGCCAGCACAATACCCTGTTTCTCGGTCAACTCATTGATGACCTTTTCTTCGCGATcg                     >  2:4414415/1‑150 (MQ=255)
                                                                                                                                    cTTTTCGATGGTCGTTTCAAGATAAACGACAACGCCACGAGCGGAAAGACGGTTACGCGTTTCACGGGATTTCACAGAGCCGCCGCCAGTAGCCAGCACAAGACCCTGTTTCTCGGTCAACTCATTGATGACCTTTTCTTCGCGATcgcg                   <  2:3706481/150‑1 (MQ=255)
                                                                                                                                     ttttCGATGGTCGTTTCAAGATAAAAGACAAAGCCACGAGCGGAAAGACGGTTACGCGTTTCACGGGATTTAACAGAGCCGCCGCCCCTAGCCAGCCCAAAACCCTGTTTCTCGGTCAACTCATTGATGACCGTTTCTTCGCGATCGCgg                  >  1:3294184/1‑150 (MQ=255)
                                                                                                                                      tttCGATGGTCGTTTCAAGATAAACGACAACGCCACGAGCGGAAAGACGGTTACGCGTTTCACGGGATTTCACAGAGCCGCCGCCAGTAGCCAGCACAATACCCTGTTTCTCGGTCAACTCATTGATGACCTTTTCTTCGCGATCGCGGa                 >  2:4188266/1‑150 (MQ=255)
                                                                                                                                      ggtCGAGGGGCGTGTCAAGATAAACGACAACGCCACGAGCGGAAAGACGGTTACGCGTTTCACGGGATTTCACAGAGCCGCCGCCAGTAGCCAGCACAAGACCCTGTTTCTCGGTCAACTCATTGATGACCTTTTCTTCGCGATCGCGGa                 <  2:4244006/148‑1 (MQ=255)
                                                                                                                                       ttCGATGGTCGTTTCAAGATAAACGACAACGCCACGAGCGGAAAGACGGTTACGCGTTTCACGGGATTTCACAGAGCCGCCGCCAGTAGCCAGCACAATACCCTGTTTCTCGGTCAACTCATTGATGACCTTTTCTTCGCGATCGCGGaa                <  2:4452995/150‑1 (MQ=255)
                                                                                                                                           aTGGTCGTTTCAAGATAAACGACAACGCCACGAGCGGAAAGACGGTTACGCGTTTCACGGGATTTCACAGAGCCGCCGCCAGTAGCCAGCACAATACCCTGTTTCTCGGTCAACTCATTGATGACCTTTTCTTCGCGATCGCGGAAGCct            <  1:2745874/150‑1 (MQ=255)
                                                                                                                                            tgttCGTTTAAAGAAAAACGACAACGCCACGAGCGGAAAGACGGTTACGCGTTTCACGGGATTTCACAGAGCCGCCGCCAGTAGCCAGCACAATACCCTGTTTCTCGGTCAACTCATTGATGACCTTTTCTTCGCGATCGCGGAAGCctt           <  1:1808355/147‑1 (MQ=255)
                                                                                                                                            tGGTCGTTTCAAGATAAACGACAACGCCACGAGCGGAAAGACGGTTACGCGTTTCACGGGATTTCACAGAGCCGCCGCCAGTAGCCAGCACAATACCCTGTTTCTCGGTCAACTCATTGATGACCTTTTCTTCGCGATCGCGGAAGCctt           >  1:3574658/1‑150 (MQ=255)
                                                                                                                                              gTCGTTTCAAGATACACGACAACGCCACGAGCGGAAAGACGGTTACGCGTTTCACGGGATTTCACAGAGCCGCCGCCAGTAGCCAGCACAATACCCTGTTTCTCGGTCAACTCATTGATGACCTTTTCTTCGCGATCGCGGAAGCcttct         <  1:471418/150‑1 (MQ=255)
                                                                                                                                                cGTTTCAAGATAAACGACAACGCCACGAGCGGAAAGACGGTTACGCGTTTCACGGGATTTCACAGAGCCGCCGCCAGTAGCCAGCACAATACCCTGTTTCTCGGTCAACTCATTGATGACCTTTTCTTCGCGATCGCGGAAGCCttcttc       <  2:1813394/150‑1 (MQ=255)
                                                                                                                                                   ttCAAGATAAACGACAACGCCACGAGCGGAAAGACGGTTACGCGTTTCACGGGATTTCACAGAGCCGCCGCCAGTAGCCAGCACAATACCCTGTTTCTCGGTCAACTCATTGATGACCTTTTCTTCGCGATCGCGGAAGCCTTCTTCGcc    <  2:4538990/150‑1 (MQ=255)
                                                                                                                                                    tCAAGATAAACGACAACGCCACGAGCGGAAAGACGGTTACGCGTTTCACGGGATTTCACAGAGCCGCCGCCAGTAGCCAGCACAATACCCTGTTTCTCGGTCAACTCATTGATGACCTTTTCTTCGCGATCGCGGAAGCCTTCTTCGCCt   <  1:291022/150‑1 (MQ=255)
                                                                                                                                                     cAAGATAAACGACAACGCCACGAGCGGAAAGACGGTTACGCGTTTCACGGGATTTCACAGAGCCGCCGCCAGTAGCCAGCACAATACCCTGTTTCTCGGTCAACTCATTGATGACCTTTTCTTCGCGATCGCGGAAGCCTTCTTCGCCtt  <  1:4383162/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
ACGAATGGTCACGTCGGCAATCTCTTCATACAGCGGATTGCGTTCATTGGCCAACGCTTCCAGAACTTCACGCGGCGGTGTTTCAACGTGCAGCAACGGGCGTTTTTTATCACGCTGCGTGCGTGCAAGTTGCTTTTCGATGGTCGTTTCAAGATAAACGACAACGCCACGAGCGGAAAGACGGTTACGCGTTTCACGGGATTTCACAGAGCCGCCGCCAGTAGCCAGCACAATACCCTGTTTCTCGGTCAACTCATTGATGACCTTTTCTTCGCGATCGCGGAAGCCTTCTTCGCCTT  >  CP009273/3511965‑3512263

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑