Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,222,642 G→T 7.8% L174I (CTC→ATC)  lysC ← lysine‑sensitive aspartokinase 3

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,222,6420GT 7.8% 1149.8 / 43.3 300L174I (CTC→ATC) lysClysine‑sensitive aspartokinase 3
Reads supporting (aligned to +/- strand):  ref base G (147/129);  new base T (8/16);  total (155/145)
Fisher's exact test for biased strand distribution p-value = 8.70e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 307 total aligned reads displayed.

GACGTCGGTCCAGATATCAACACGAGATGCGTGTAAAGCCTCCGCCAGCAAGGCTGCCGTATAATCGCTGCCTCCACGGCCAAGCGTCGTTGTACGACCTTTATTTTCGCTACCGATAAATCCCTGGGTGATCACTAAGCCTTCATTGAGACGTGGGAGCAGCTGCAGCGCGGCCAGTTCCGCCAGCGCGGCTATATCTGGCTCTGCACGACCAAATCGGTCGTTGGTACGCATCACTTTACGTACATCAAACCACTGTGCCTGAACATCGCGTTCGCGCAGGATCTCAACAAACAGCA  >  CP009273/4222493‑4222791
                                                                                                                                                     |                                                                                                                                                     
    tCGGTCCAGATATCAACACGAGATGCGTGTAAAGCCTCCGCCAGCAAGGCTGCCGTATAATCGCTGCCTCCACGGCCAAGCGTCGTTGTACGACCTTTATTTTCGCTACCGATAAATCCCTGGGTGATCACTAAGCCTTCATTGAGACGt                                                                                                                                                   <  1:2221676/150‑1 (MQ=255)
      ggTCCAGATATCAACACGAGATGCGTGTAAAGCCTCCGCCAGCAAGGCTGCCGTATAATCGCTGCCTCCACGGCCAAGCGTCGTTGTACGACCTTTATTTTCGCTACCGATAAATCCCTGGGTGATCACTAAGCCTTCATTGAGACGTgg                                                                                                                                                 >  1:3518308/1‑150 (MQ=255)
      ggTCCAGATATCAACACGAGATGCGTGTAAAGCCTCCGCCAGCAAGGCTGCCGTATAATCGCTGCCTCCACGGCCAAGCGTCGTTGTACGACCTTTATTTTCGCTACCGATAAATCCCTGGGTGATCACTAAGCCTTCATTGAGACGTgg                                                                                                                                                 <  1:374072/150‑1 (MQ=255)
         ccAGATATCAACACGAGATGCGTGTAAAGCCTCCGCCAGCAAGGCTGCCGTATAATCGCTGCCTCCACGGCCAAGCGTCGTTGTACGACCTTTATTTTCGCTACCGATAAATCCCTGGGTGATCACTAAGCCTTCATTGAGACGTGGGag                                                                                                                                              <  2:4086602/150‑1 (MQ=255)
         ccAGATATCAACACGAGATGCGTGTAAAGCCTCCGCCAGCAAGGCTGCCGTATAATCGCTGCCTCCACGGCCAAGCGTCGTTGTACGACCTTTATTTTCGCTACCGATAAATCCCTGGGTGATCACTAAGCCTTCATTGAGACGTGGGag                                                                                                                                              <  2:1236895/150‑1 (MQ=255)
          cAGATATCAACACGAGATGCGTGTAAAGCCTCCGCCAGCAAGGCTGCCGTATAATCGCTGCCTCCACGGCCAAGCGTCGTTGTACGACCTTTATTTTCGCTACCGATAAATCCCTGGGTGATCACTAAGCCTTCATTGAGACGTGGGagc                                                                                                                                             <  2:370668/150‑1 (MQ=255)
           aGATATCAACACGAGATGCGTGTAAAGCCTCCGCCAGCAAGGCTGCCGTATAATCGCTGCCTCCACGGCCAAGCGTCGTTGTACGACCTTTATTTTCGCTACCGATAAATCCCTGGGTGATCACTAAGCCTTCATTGAGACGTGGGagca                                                                                                                                            <  2:3116704/150‑1 (MQ=255)
              tatCAACACGAGAGGCGTGTAAAGCCTCCGCCAGCAAGGCTGCCGTATAATCGCTGCCTCCACGGCCAAGCGTCGTTGTACGACCTTTATTTTCGCTACCGATAAATCCCTGGGTGATCACTAAGCCTTCATTGATACGTGGGAGCAGCt                                                                                                                                         <  2:2051339/150‑1 (MQ=255)
               atCAACACGAGATGCGTGTAAAGCCTCCGCCAGCAAGGCTGCCGTATAATCGCTGCCTCCACGGCCAAGCGTCGTTGTACGACCTTTATTTTCGCTACCGATAAATCCCTGGGTGATCACTAAGCCTTCATTGAGACGTGGGAGCAGCTg                                                                                                                                        >  1:3260376/1‑150 (MQ=255)
                  aaCACGAGATGCGTGTAAAGCCTCCGCCAGCAAGGCTGCCGTATAATCGCTGCCTCCACGGCCAAGCGTCGTTGTACGACCTTTATTTTCGCTACCGATAAATCCCTGGGTGATCACTAAGCCTTCATTGAGACGTGGGAGCAGCTGCAg                                                                                                                                     <  2:2489100/150‑1 (MQ=255)
                        agaTGCGTGTAAAGCCTCCGCCAGCAAGGCTGCCGTATAATCGCTGCCTCCACGGCCAAGCGTCGTTGTACGACCTTTATTTTCGCTACCGATAAATCCCTGGGTGATCACTAAGCCTTCATTGAGACGTGGGAGCAGCTGCAGCGCGGc                                                                                                                               >  2:1023223/1‑150 (MQ=255)
                        agaTGCGTGTAAAGCCTCCGCCAGCAAGGCTGCCGTATAATCGCTGCCTCCACGGCCAAGCGTCGTTGTACGACCTTTATTTTCGCTACCGATAAATCCCTGGGTGATCACTAAGCCTTCATTGAGACGTGGGAGCAGCTGCAGCGCGGc                                                                                                                               >  2:106004/1‑150 (MQ=255)
                             cGTGTAAAGCCTCCGCCAGCAAGGCTGCCGTATAATCGCTGCCTCCACGGCCAAGCGTCGTTGTACGACCTTTATTTTCGCTACCGATAAATCCCTGGGTGATCACTAAGCCTTCATTGAGACGTGGGAGCAGCTGCAGCGCGGCCAGtt                                                                                                                          <  1:3320451/150‑1 (MQ=255)
                                gtAAAGCCTCCGCCAGCAAGGCTGCCGTATAATCGCTGCCTCCACGGCCAAGCGTCGTTGTACGACCTTTATTTTCGCTACCGATAAATCCCTGGGTGATCACTAAGCCTTCATTGAGACGTGGGAGCAGCTGCAGCGCGGCCAGTTCCg                                                                                                                       <  1:1931114/150‑1 (MQ=255)
                                 tAAAGCCTCCGCCAGCAAGGCTGCCGTATAATCGCTGCCTCCACGGCCAAGCGTCGTTGTACGACCTTTATTTTCGCTACCGATAAATCCCTGGGTGATCACTAAGCCTTCATTGAGACGTGGGAGCAGCTGCAGCGCGGCCAGTTCCGc                                                                                                                      >  1:1061988/1‑150 (MQ=255)
                                   taGCCTACTCCAGCAAGGCTGCCGTATAATCGCTGCCTCCAGGGCCAAGCGTCGTTGTACGTACTTTATTTTCGCTACAGATATATCCCTTTGTGATCACTAAGCCTTCATTGAGACTTGGGTGCAGCTGCTGCGCGGCCAGTTCAGCCa                                                                                                                    <  1:239300/149‑1 (MQ=255)
                                   aaGCCTCCGCCAGCAAGGCTGCCGTATAATCGCTGCCTCCACGGCCAAGCGTCGTTGTACGACCTTTATTTTCGCTACCGATAAATCCCTGGGTGATCACTAAGCCTTCATTGAGACGTGGGAGCAGCTGCAGCGCGGCCAGTTCCGCCa                                                                                                                    >  2:2734724/1‑150 (MQ=255)
                                      ccTCCGCCAGCAAGGCTGCCGTATAATCGCTGCCTCCACGGCCAAGCGTCGTTGTACGACCTTTATTTTCGCTACCGATAAATCCCTGGGTGATCACTAAGCCTTCATTGAGACGTGGGAGCAGCTGCAGCGCGGCCAGTTCCGCCAgcg                                                                                                                 >  2:2526319/1‑150 (MQ=255)
                                        tCCGCCAGCAAGGCTGCCGTATAATCGCTGCCTCCACGGCCAAGCGTCGTTGTACGACCTTTATTTTCGCTACCGATAAATCCCTGGGTGATCACTAAGCCTTCATTGAGACGTGGGAGCAGCTGCAGCGCGGCCAGTTCCGCCAGcgcg                                                                                                               >  1:1462734/1‑150 (MQ=255)
                                           gCCAGCAAGGCTGCCGTATAATCGCTGCCTCCACGGCCAAGCGTCGTTGTACGACCTTTATTTTCGCTACCGATAAATCCCTGGGTGATCACTAAGCCTTCATTGATACGTGGGAGCAGCTGCAGCGCGGCCAGTTCCGCCAGCGCGGCt                                                                                                            >  1:3758841/1‑150 (MQ=255)
                                            ccAGCAAGGCTGCCGTATAATCGCTGCCTCCACGGCCAAGCGTCGTTGTACGACCTTTATTTTCGCTACCGATAAATCCCTGGGTGATCACTAAGCCTTCATTGAGACGTGGGAGCAGCTGCAGCGCGGCCAGTTCCGCCAGCGCGGCta                                                                                                           <  1:252855/150‑1 (MQ=255)
                                             cAGCAAGGCTGCCGTATAATCGCTGCCTCCACGGCCAAGCGTCGTTGTACGACCTTTATTTTCGCTACCGATAAATCCCTGGGTGATCACTAAGCCTTCATTGAGACGTGGGAGCAGCTGCAGCGCGGCCAGTTCCGCCAGCGCGGCtat                                                                                                          >  1:3777182/1‑150 (MQ=255)
                                             cAGCAAGGCTGCCGTATAATCGCTGCCTCCACGGCCAAGCGTCGTTGTACGACCTTTATTTTCGCTACCGATAAATCCCTGGGTGATCACTAAGCCTTCATTGAGACGTGGGAGCAGCTGCAGCGCGGCCAGTTCCGCCAGCGCGGCtat                                                                                                          <  1:698966/150‑1 (MQ=255)
                                             cAGCAAGGCTGCCGTATAATCGCTGCCTCCACGGCCAAGCGTCGTTGTACGACCTTTATTTTCGCTACCGATAAATCCCTGGGTGATCACTAAGCCTTCATTGAGACGTGGGAGCAGCTGCAGCGCGGCCAGTTCCGCCAGCGCGGCtat                                                                                                          <  1:2617721/150‑1 (MQ=255)
                                              aGCAAGGCTGCCGTATAATCGCTGCCTCCACGGCCAAGCGTCGTTGTACGACCTTTATTTTCGCTACCGATAAATCCCTGGGTGATCACTAAGCCTTCATTGAGACGTGGGAGCAGCTGCAGCGCGGCCAGTTCCGCCAGCGCGGCtata                                                                                                         >  2:900728/1‑150 (MQ=255)
                                                 aaGGCTGCCGTATAATCGCTGCCTCCACGGCCAAGCGTCGTTGTACGACCTTTATTTTCGCTACCGATAAATCCCTGGGTGATCACTAAGCCTTCATTGAGACGTGGGAGCAGCTGCAGCGCGGCCAGTTCCGCCAGCGCGGCTATATCt                                                                                                      >  1:2436801/1‑150 (MQ=255)
                                                  aGGCTGCCGTATAATCGCTGCCTCCACGGCCAAGCGTCGTTGTACGACCTTTATTTTCGCTACCGATAAATCCCTGGGTGATCACTAAGCCTTCATTGAGACGTGGGAGCAGCTGCAGCGCGGCCAGTTCCGCCAGCGCGGCTATATCTg                                                                                                     >  2:3980038/1‑150 (MQ=255)
                                                   ggCTGCCGTATAATCGCTGCCTCCACGGCCAAGCGTCGTTGTACGACCTTTATTTCCGCTACCGATAAATCCCTGGGTGATCACTAAGCCTTCATTGAGACGTGGGAGCAGCTGCAGCGCGGCCAGTTCCGCCAGCGCGGCTATATCTgg                                                                                                    <  1:2089075/150‑1 (MQ=255)
                                                    gCTGCCGTATAATCGCTGCCTCCACGGCCAAGCGTCGTTGTACGACCTTTATTTTCGCTACCGATAAATCCCTGGGTGATCACTAAGCCTTCATTGAGACGTGGGAGCAGCTGCAGCGCGGCCAGTTCCGCCAGCGCGGCTATATCTGGc                                                                                                   <  2:215768/150‑1 (MQ=255)
                                                     cTGCCGTATAATCGCTGCCTCCACGGCCAAGCGTCGTTGTACGACCTTTATTTTCGCTACCGATAAATCCCTGGGTGATCACTAAGCCTTCATTGAGACGTGGGAGCAGCTGCAGCGCGGCCAGTTCCGCCAGCGCGGCTATATCTGGct                                                                                                  <  1:1606902/150‑1 (MQ=255)
                                                       gCCGTATAATCGCTGCCTCCACGGCCAAGCGTCGTTGTACGACCTTTATTTTCGCTACCGATAAATCCCTGGGTGATCACTAAGCCTTCATTGATACGTGGGAGCAGCTGCAGCGCGGCCAGTTCCGCCAGCGCGGCTATATCTGGctct                                                                                                <  1:1286895/150‑1 (MQ=255)
                                                        ccGTATAATCGCTGCCTCCACGGCCAAGCGTCGTTGTACGACCTTTATTTTCGCTACCGATAAATCCCTGGGTGATCACTAAGCCTTCATTGATACGTGGGAGCAGCTGCAGCGCGGCCAGTTCCGCCAGCGCGGCTATATCTGGCTCTg                                                                                               <  1:1277009/150‑1 (MQ=255)
                                                           tataATCGCTGCCTCCACGGCCAAGCGTCGTTGTACGACCTTTATTTTCGCTACCGATAAATCCCTGGGTGATCACTAAGCCTTCATTGAGACGTGGGAGCAGCTGCAGCGCGGCCAGTTCCGCCAGCGCGGCTATATCTGGCTCTGCAc                                                                                            <  1:1041500/150‑1 (MQ=255)
                                                           tataATCGCTGCCCCCACGGCCAAGCGTCGTTGTACGACCTTTATTTTCGCTACCGATAAATCCCTGGGTGATCACTAAGCCTTCATTGAGACGTGGGAGCAGCTGCAGCGCGGCCAGTTCCGCCAGCGCGGCTATATCTGGCTCTGCAc                                                                                            <  1:3928407/150‑1 (MQ=255)
                                                            ataATCGCTGCCTCCACGGCCAAGCGTCGTTGTACGACCTTTATTTTCGCTACCGATAAATCCCTGGGTGATCACTAAGCCTTCATTGAGACGTGGGAGCAGCTGCAGCGCGGCCAGTTCCGCCAGCGCGGCTATATCTGGCTCTGCACg                                                                                           >  2:2068860/1‑150 (MQ=255)
                                                              aaTCGCTGCCTCCACGGCCAAGCGTCGTTGTACGACCTTTATTTTCGCTACCGATAAATCCCTGGGTGATCACTAAGCCTTCATTGAGACGTGGGAGCAGCTGCAGCGCGGCCAGTTCCGCCAGCGCGGCTATATCTGGCTCTGCACGAc                                                                                         >  1:823328/1‑150 (MQ=255)
                                                                  gCTGCCTCCACGGCCAAGCGTCGTTGTACGACCTTTATTTTCGCTACCGATAAATCCCTGGGTGATCACTAAGCCTTCATTGATACGTGGGAGCAGCTGCAGCGCGGCCAGTTCCGCCAGCGCGGCTATATCTGGCTCTGCACGACCaaa                                                                                     >  1:4007706/1‑150 (MQ=255)
                                                                   cTGCCTCCACGGCCAAGCGTCGTTGTACGACCTTTATTTTCGCTACCGATAAATCCCTGGGTGATCACTAAGCCTTCATTGAGACGTGGGAGCAGCTGCAGCGCGGCCAGTTCCGCCAGCGCGGCTATATCTGGCTCTGCACGACCAAAt                                                                                    <  2:378549/150‑1 (MQ=255)
                                                                   cTGCCTCCACGGCCAAGCGTCGTTGTACGACCTTTATTTTCGCTACCGATAAATCCCTGGGTGATCACTAAGCCTTCATTGAGACGTGGGAGCAGCTGCAGCGCGGCCAGTTCCGCCAGCGCGGCTATATCTGGCTCTGCACGACCAAAt                                                                                    <  1:1770848/150‑1 (MQ=255)
                                                                   cTGCCTCCACGGCCAAGCGTCGTTGTACGACCTTTATTTTCGCTACCGATAAATCCCTGGGTGATCACTAAGCCTTCATTGAGACGTGGGAGCAGCTGCAGCGCGGCCAGTTCCGCCAGCGCGGCTATATCTGGCTCTGCACGACCAAAt                                                                                    <  1:2235853/150‑1 (MQ=255)
                                                                     gCCTCCACGGCCAAGCGTCGTTGTACGACCTTTATTTTCGCTACCGATAAATCCCTGGGTGATCACTAAGCCTTCATTGAGACGTGGGAGCAGCTGCAGCGCGGCCAGTTCCGCCAGCGCGGCTATATCTGGCTCTGCACGACCAAATCg                                                                                  <  2:4043967/150‑1 (MQ=255)
                                                                     gCCTCCACGGCCAAGCGTCGTTGTACGACCTTTATTTTCGCTACCGATAAATCCCTGGGTGATCACTAAGCCTTCATTGAGACGTGGGAGCAGCTGCAGCGCGGCCAGTTCCGCCAGCGCGGCTATATCTGGCTCTGCACGACCAAATCg                                                                                  <  2:3766420/150‑1 (MQ=255)
                                                                     gCCTCCACGGCCAAGCGTCGTTGTACGACCTTTATTTTCGCTACCGATAAATCCCTGGGTGATCACTAAGCCTTCATTGAGACGTGGGAGCAGCTGCAGCGCGGCCAGTTCCGCCAGCGCGGCTATATCTGGCTCTGCACGACCAAATCg                                                                                  >  2:238002/1‑150 (MQ=255)
                                                                       cTCCACGGCCAAGCGTCGTTGTACGACCTTTATTTTCGCTACCGATAAATCCCTGGGTGATCACTAAGCCTTCATTGAGACGTGGGAGCAGCTGCAGCGCGGCCAGTTCCGCCAGCGCGGCTATATCTGGCTCTGCACGACCAAATCGGt                                                                                <  2:2473453/150‑1 (MQ=255)
                                                                        tCCACGGCCAAGCGTCGTTGTACGACCTTTATTTTCGCTACCGATAAATCCCTGGGTGATCACTAAGCCTTCATTGAGACGTGGGAGCAGCTGCAGCGCGGCCAGTTCCGCCAGCGCGGCTATATCTGGCTCTGCACGACCAAATCGGTc                                                                               <  2:1490503/150‑1 (MQ=255)
                                                                         ccACGGCCAAGCGTCGTTGTACGACCTTTATTTTCGCTACCGATAAATCCCTGGGTGATCACTAAGCCTTCATTGAGACGTGGGAGCAGCTGCAGCGCGGCCAGTTCCGCCAGCGCGGCTATATCTGGCTCTGCACGACCAAATCGGTCg                                                                              >  1:2867426/1‑150 (MQ=255)
                                                                         ccACGGCCAAGCGTCGTTGTACGACCTTTATTTTCGCTACCGATAAATCCCTGGGTGATCACTAAGCCTTCATTGAGACGTGGGAGCAGCTGCAGCGCGGCCAGTTCCGCCAGCGCGGCTATATCTGGCTCTGCACGACCAAATCGGTCg                                                                              <  1:1518945/150‑1 (MQ=255)
                                                                          cACGGCCAAGCGTCGTTGTACGACCTTTATTTTCGCTACCGATAAATCCCTGGGTGATCACTAAGCCTTCATTGATACGTGGGAGCAGCTGCAGCGCGGCCAGTTCCGCCAGCGCGGCTATATCTGGCTCTGCACGACCAAATCGGTCGt                                                                             <  1:3577918/150‑1 (MQ=255)
                                                                           aCGGCCAAGCGTCGTTGTACGACCTTTATTTTCGCTACCGATAAATCCCTGGGTGATCACTAAGCCTTCATTGAGACGTGGGAGCAGCTGCAGCGCGGCCAGTTCCGCCAGCGCGGCTATATCTGGCTCTGCACGACCAAATCGGTCGtt                                                                            <  1:3123596/150‑1 (MQ=255)
                                                                              gCCAAGCGTCGTTGTACGACCTTTATTTTCGCTACCGATAAATCCCTGGGTGATCACTAAGCCTTCATTGAGACGTGGGAGCAGCTGCAGCGCGGCCAGTTCCGCCAGCGCGGCTATATCTGGCTCTGCACGACCAAATCGGTCGTTGGt                                                                         >  2:1362952/1‑150 (MQ=255)
                                                                              gCCAAGCGTCGTTGTACGACCTTTATTTTCGCTACCGATAAATCCCTGGGTGATCACTAAGCCTTCATTGAGACGTGGGAGCAGCTGCAGCGCGGCCAGTTCCGCCAGCGCGGCTATATCTGGCTCTGCACGACCAAATCGGTCGTTGGt                                                                         >  2:6791/1‑150 (MQ=255)
                                                                                cAAGCGTCGTTGTACGACCTTTATTTTCGCTACCGATAAATCCCTGGGTGATCACTAAGCCTTCATTGAGACGTGGGAGCAGCTGCAGCGCGGCCAGTTCCGCCAGCGCGGCTATATCTGGCTCTGCACGACCAAATCGGTCGTTGGTAc                                                                       >  1:1156012/1‑150 (MQ=255)
                                                                                cAAGCGTCGTTGTACGACCTTTATTTTCGCTACCGATAAATCCCTGGGTGATCACTAAGCCTTCATTGAGACGTGGGAGCAGCTGCAGCGCGGCCAGTTCCGCCAGCGCGGCTATATCTGGCTCTGCACGACCAAATCGGTCGTTGGTAc                                                                       <  1:4044875/150‑1 (MQ=255)
                                                                                 aaGCGTCGTTGTACGACCTTTATTTTCGCTACCGATAAATCCCTGGGTGATCACTAAGCCTTCATTGAGACGTGGGAGCAGCTGCAGCGCGGCCAGTTCCGCCAGCGCGGCTATATCTGGCTCTGCACGACCAAATCGGTCGTTGGTACg                                                                      >  2:3182084/1‑150 (MQ=255)
                                                                                 aaGCGTCGTTGTACGACCTTTATTTTCGCTACCGATAAATCCCTGGGTGATCACTAAGCCTTCATTGAGACGTGGGAGCAGCTGCAGCGCGGCCAGTTCCGCCAGCGCGGCTATATCTGGCTCTGCACGACCAAATCGGTCGTTGGTACg                                                                      >  2:302221/1‑150 (MQ=255)
                                                                                 aaGCGTCGTTGTACGACCTTTATTTTCGCTACCGATAAATCCCTGGGTGATCACTAAGCCTTCATTGAGACGTGGGAGCAGCTGCAGCGCGGCCAGTTCCGCCAGCGCGGCTATATCTGGCTCTGCACGACCAAATCGGTCGTTGGTACg                                                                      >  2:451672/1‑150 (MQ=255)
                                                                                  aGCTTCGTTGTACGACCATTATTTTAGCTACCGATAAATCCCGGGGTGATCAAAAAGCCTTAATTGAGACGTGTGAGCAGCTGCAGCGCGGCCAGTTCCGCCAGCGCGGCTATATCTGGCTCTGCACGAACAAATCGGTCGTTGGTACGc                                                                     >  2:2042131/1‑150 (MQ=255)
                                                                                  aGCGTCGTTGTACGACCTTTATTTTCGCTACCGATAAATCCCTGGGTGATCACTAAGCCTTCATTGAGACGTGGGAGCAGCTGCAGCGCGGCCAGTTCCGCCAGCGCGGCTATATCTGGCTCTGCACGACCAAATCGGTCGTTGGTACGc                                                                     >  2:3824519/1‑150 (MQ=255)
                                                                                   gCGTCGTTGTACGACCTTTATTTTCGCTACCGATAAATCCCTGGGTGATCACTAAGCCTTCATTGAGACGTGGGAGCAGCTGCAGCGCGGCCAGTTCCGCCAGCGCGGCTATATCTGGCTCTGCACGACCAAATCGGTCGTTGGTACGCa                                                                    >  1:2079130/1‑150 (MQ=255)
                                                                                    cgtcgtTGTACGACCTTTATTTTCGCTACCGATAAATCCCTGGGTGATCACTAAGCCTTCATTGAGACGTGGGAGCAGCTGCAGCGCGGCCAGTTCCGCCAGCGCGGCTATATCTGGCTCTGCACGACCAAATCGGTCGTTGGTACGCAt                                                                   <  2:2630347/150‑1 (MQ=255)
                                                                                       cgtTGTACGACCTTTATTTTCGCTACCGATAAATCCCTGGGTGATCACTAAGCCTTCATTGAGACGTGGGAGCAGCTGCAGCGCGGCCAGTTCCGCCAGCGCGGCTATATCTGGCTCTGCACGACCAAATCGGTCGTTGGTACGCATCAc                                                                >  1:2279528/1‑150 (MQ=255)
                                                                                        gtTGTACGACCTTTATTTTCGCTACCGATAAATCCCTGGGTGATCACTAAGCCTTCATTGAGACGTGGGAGCAGCTGCAGCGCGGCCAGTTCCGCCAGCGCGGCTATATCTGGCTCTGCACGACCAAATCGGTCGTTGGTACGCATCACt                                                               <  2:3285577/150‑1 (MQ=255)
                                                                                          tGTACGACCTTTATTTTCGCTACCGATAAATCCCTGGGTGATCACTAAGCCTTCATTGATACGTGGGAGCAGCTGCAGCGCGGCCAGTTCCGCCAGCGCGGCTATATCTGGCTCTGCACGACCAAATCGGTCGTTGGTACGCATCACttt                                                             <  2:3096432/150‑1 (MQ=255)
                                                                                            tACGACCTTTAGTTTCGCTACCGATAAATCCCTGGGTGATCACTAAGCCTTCATTGAGACGTGGGAGCAGCTGCAGCGCGGCCAGTTCCGCCAGCGCGGCTATATCTGGCTCTGCACGACCAAATCGGTCGTTGGTACGCATCACTTTAc                                                           >  2:3068398/1‑150 (MQ=255)
                                                                                             aCGACCTTTATTTTCGCTACCGATAAATCCCTGGGTGATCACTAAGCCTTCATTGATACGTGGGAGCAGCTGCAGCGCGGCCAGTTCCGCCAGCGCGGCTATATCTGGCTCTGCACGACCAAATCGGTCGTTGGTACGCATCACTTTACg                                                          <  2:3366099/150‑1 (MQ=255)
                                                                                             aCGACCTTTATTTTCGCTACCGATAAATCCCTGGGTGATCACTAAGCCTTCATTGAGACGTGGGAGCAGCTGCAGCGCGGCCAGTTCCGCCAGCGCGGCTATATCTGGCTCTGCACGACCAAATCGGTCGTTGGTACGCATCACTTTACg                                                          <  2:3609812/150‑1 (MQ=255)
                                                                                               gACCTTTATTTTCGCTACCGATAAATCCCTGGGTGATCACTAAGCCTTCATTGAGACGTGGGAGCAGCTGCAGCGCGGCCAGTTCCGCCAGCGCGGCTATATCTGGCTCTGCACGACCAAATCGGTCGTTGGTACGCATCACTTTACGTa                                                        >  2:560988/1‑150 (MQ=255)
                                                                                                aCCTTTATTTTCGCTACCGATAAATCCCTGGGTGATCACTAAGCCTTCATTGAGACGTGGGAGCAGCTGCAGCGCGGCCAGTTCCGCCAGCGCGGCTATATCTGGCTCTGCACGACCAAATCGGTCGTTGGTACGCATCACTTTACGTAc                                                       >  1:2470901/1‑150 (MQ=255)
                                                                                                 ccTTTATTTTCGCTACCGATAAATCCCTGGGTGATCACTAAGCCTTCATTGAGACGTGGGAGCAGCTGCAGCGCGGCCAGTTCCGCCAGCGCGGCTATATCTGGCTCTGCACGACCAAATCGGTCGTTGGTACGCATCACTTTACGTACa                                                      >  2:2106799/1‑150 (MQ=255)
                                                                                                  cTTTATTTTCGCTACCGATAAATCCCTGGGTGATCACTAAGCCTTCATTGAGACGTGGGAGCAGCTGCAGCGCGGCCAGTTCCGCCAGCGCGGCTATATCTGGCTCTGCACGACCAAATCGGTCGTTGGTACGCATCACTTTACGTACAt                                                     <  2:3220923/150‑1 (MQ=255)
                                                                                                      aTTTTCGCTACCGATAAATCCCTGGGTGATCACTAAGCCTTCATTGAGACGTGGGAGCAGCTGCAGCGCGGCCAGTTCCGCCAGCGCGGCTATATCTGGCTCTGAACGACCAAAACGGACGTTGGTACGCATCACTATACGTACATCaaa                                                 >  2:2585200/1‑150 (MQ=255)
                                                                                                       ttttCGCTACCGATAAATCCCTGGGTGATCACTAAGCCTTCATTGAGACGTGGGAGCAGCTGCAGCGCGGCCAGTTCCGCCAGCGCGGCTATATCTGGCTCTGCACGACCAAATCGGTCGTTGGTACGCATCACTTTACGTACATCAAAc                                                <  1:2312757/150‑1 (MQ=255)
                                                                                                       ttttCGCTACCGATAAATCCCTGGGTGATCACTAAGCCTTCATTGAGACGTGGGAGCAGCTGCAGCGCGGCCAGTTCCGCCAGCGCGGCTATATCTGGCTCTGCACGACCAAATCGGTCGTTGGTACGCATCACTTTACGTACATCAAAc                                                <  1:570951/150‑1 (MQ=255)
                                                                                                        tttCGCTACCGATAAATCCCTGGGTGATCACTAAGCCTTCATTGAGACGTGGGAGCAGCTGCAGCGCGGCCAGTTCCGCCAGCGCGGCTATATCTGGCTCTGCACGACCAAATCGGTCGTTGGTACGCATCACTTTACGTACATCAAAcc                                               <  2:1161496/150‑1 (MQ=255)
                                                                                                           cGCTACCGATAAATCCCTGGGTGATCACTAAGCCTTCATTGAGACGTGGGAGCAGCTGCAGCGCGGCCAGTTCCGCCAGCGCGGCTATATCTGGCTCTGCACGACCAAATCGGTCGTTGGTACGCATCACTTTACGTACATCAAACCACt                                            >  2:3921364/1‑150 (MQ=255)
                                                                                                               aCCGATAAATCCCTGGGTGATCACTAAGCCTTCATTGAGACGTGGGAGCAGCTGCAGCGCGGCCAGTTCCGCCAGCGCGGCTATATCTGGCTCTGCACGACCAAATCGGTCGTTGGTACGCATCACTTTACGTACATCAAACCACTGTGc                                        >  1:2754657/1‑150 (MQ=255)
                                                                                                               aCCGATAAATCCCTGGGTGATCACTAAGCCTTCATTGAGACGTGGGAGCAGCTGCAGCGCGGCCAGTTCCGCCAGCGCGGCTATATCTGGCTCTGCACGACCAAATCGGTCGTTGGTACGCATCACTTTACGTACATCAAACCACTGTGc                                        >  1:793564/1‑150 (MQ=255)
                                                                                                               aCCGATAAATCCCTGGGTGATCACTAAGCCTTCATTGAGACGTGGGAGCAGCTGCAGCGCGGCCAGTTCCGCCAGCGCGGCTATATCTGGCTCTGCACGACCAAATCGGTCGTTGGTACGCATCACTTTACGTACATCAAACCACTGTGc                                        <  1:2668391/150‑1 (MQ=255)
                                                                                                               aCCGATAAATCCCTGGGTGATCACTAAGCCTTCATTGAGACGTGGGAGCAGCTGCAGCGCGGCCAGTTCCGCCAGCGCGGCTATATCTGGCTCTGCACGACCAAATCGGTCGTTGGTACGCATCACTTTACGTACATCAAACCACTGTGc                                        >  1:2176632/1‑150 (MQ=255)
                                                                                                                ccGATAAATCCCTGGGTGATCACTAAGCCTTCATTGAGACGTGGGAGCAGCTGCAGCGCGGCCAGTTCCGCCAGCGCGGCTATATCTGGCTCTGCACGACCAAATCGGTCGTTGGTACGCATCACTTTACGTACATCAAACCACTGTGcc                                       <  1:3938946/150‑1 (MQ=255)
                                                                                                                  gATAAATCCCTGGGTGATCACTAAGCCTTCATTGAGACGTGGGAGCAGCTGCAGCGCGGCCAGTTCCGCCAGCGCGGCTATATCTGGCTCTGCACGACCAAATCGGTCGTTGGTACGCATCACTTTACGTACATCAAACCACTGTGCCTg                                     >  2:497873/1‑150 (MQ=255)
                                                                                                                   aTAAATCCCTGGGTGATCACTAAGCCTTCATTGAGACGTGGGAGCAGCTGCAGCGCGGCCAGTTCCGCCAGCGCGGCTATATCTGGCTCTGCACGACCAAATCGGTCGTTGGTACGCATCACTTTACGTACATCAAACCACTGTGCCTGa                                    >  2:3888883/1‑150 (MQ=255)
                                                                                                                       aTCCCTGGGTGATCACTAAGCCTTCATTGAGACGTGGGAGCAGCTGCAGCGCGGCCAGTTCCGCCAGCGCGGCTATATCTGGCTCTGCACGACCAAATCGGTCGTTGGTACGCATCACTTTACGTACATCAAACCACTGTGCCTGAACAt                                >  2:748426/1‑150 (MQ=255)
                                                                                                                             gggTGATCACGAAGCCTTCATTGAGACGTGGGAGCAGCTGCAGCGCGGCCAGTTCCGCCAGCGCGGCTATATCTGGCTCTGCACGACCAAATCGGTCGTTGGTACGCATCACTTTACGTACATCAAACCACTGTGCCTGAACATCGCGtt                          <  1:961316/150‑1 (MQ=255)
                                                                                                                                 gATCACTAAGCCTTCATTGAGACGTGGGAGCAGCTGCAGCGCGGCCAGTTCCGCCAGCGCGGCTATATCTGGCTCTGCACGACCAAATCGGTCGTTGGTACGCATCACTTTACGTACATCAAACCACTGTGCCTGAACATCGCGTTcgcg                      <  2:1933158/150‑1 (MQ=255)
                                                                                                                                   tCACTAAGCCTTCATTGAGACGTGGGAGCAGCTGCAGCGCGGCCAGTTCCGCCAGCGCGGCTATATCTGGCTCTGCACGACCAAATCGGTCGTTGGTACGCATCACTTTACGTACATCAAACCACTGTGCCTGAACATCGCGTTCGCGCa                    >  1:1008710/1‑150 (MQ=255)
                                                                                                                                     aCTAAGCCTTCATTGAGACGTGGGAGCAGCTGCAGCGCGGCCAGTTCCGCCAGCGCGGCTATATCTGGCTCTGCACGACCAAATCGGTCGTTGGTACGCATCACTTTACGTACATCAAACCACTGTGCCTGAACATCGCGTTCGCGCAgg                  >  1:1812744/1‑150 (MQ=255)
                                                                                                                                       tAAGCCTTCATTGAGACGTGGGAGCAGCTGCAGCGCGGCCAGTTCCGCCAGCGCGGCTATATCTGGCTCTGCACGACCAAATCGGTCGTTGGTACGCATCACTTTACGTACATCAAACCACTGTGCCTGAACATCGCGTTCGCGCAGGAt                >  2:3774037/1‑150 (MQ=255)
                                                                                                                                       tAAGCCTTCATTGAGACGTGGGAGCAGCTGCAGCGCGGCCAGTTCCGCCAGCGCGGCTATATCTGGCTCTGCACGACCAAATCGGTCGTTGGTACGCATCACTTTACGTACATCAAACCACTGTGCCTGAACATCGCGTTCGCGCAGGAt                >  1:3251094/1‑150 (MQ=255)
                                                                                                                                       tAAGCCTTCATTGAGACGTGGGAGCAGCTGCAGCGCGGCCAGTTCCGCCAGCGCGGCTATATCTGGCTCTGCACGACCAAATCGGTCGTTGGTACGCATCACTTTACGTACATCAAACCACTGTGCCTGAACATCGCGTTCGCGCAGGAt                >  1:411163/1‑150 (MQ=255)
                                                                                                                                       tAAGCCTTCATTGAGACGTGGGAGCAGCTGCAGCGCGGCCAGTTCCGCCAGCGCGGCTATATCTGGCTCTGCACGACCAAAACGGTCGTTGGTACGCATCACTTTACGTACATCAAACCACTGTGCCTGAACATCGCGATCGCGCAGGAt                >  2:2884761/1‑150 (MQ=255)
                                                                                                                                         aGCCTTCATTGAGACGTGGGAGCAGCTGCAGCGCGGCCAGTTCCGCCAGCGCGGCTATATCTGGCTCTGCACGACCAAATCGGTCGTTGGTACGCATCACTTTACGTACATCAAACCACTGTGCCTGAACATCGCGTTCGCGCAGGAtct              <  1:3866770/150‑1 (MQ=255)
                                                                                                                                          gCCTTCATTGATACGTGGGAGCAGCTGCAGCGCNGCCAGTTCCNCCAGCGCGGCTATATCTGGCTCTGCACGACCAAATCGGTCGTTGGTACGCATCACTTTACGTACATCAAACCACTGTGCCTGAACATCGCGTTCGCGCAGGAtctc             <  1:2937874/150‑1 (MQ=255)
                                                                                                                                                aTTGAGACGTGGGAGCAGCTGCAGCGCGGCCAGTTCCGCCAGCGCGGCTATATCTGGCTCTGCACGACCAAATCGGTCGTTGGTACGCATCACTTTACGTACATCAAACCACTGTGCCTGAACATCGCGTTCGCGCAGGATCTCaacaaa       <  2:2290424/150‑1 (MQ=255)
                                                                                                                                                aTTGAGAAGTGGGAGCAGCTGCAGCGCGGCCAGTTCCGCCAGCGCGGCTATATCTGGCTCTGCACGACCAAATCGGTCGTTGGTACGCATCACTTTACGTACATCAAACCACTGTGCCTGAACATCGCGTTCGCGCAGGATCTCaacaaa       >  2:3817246/1‑150 (MQ=255)
                                                                                                                                                 ttGAGACGTGGGAGCAGCTGCAGCGCGGCCAGTTCCGCCAGCGCGGCTATATCTGGCTCTGCACGACCAAATCGGTCGTTGGTACGCATCACTTTACGTACATCAAACCACTGTGCCTGAACATCGCGTTCGCGCAGGATCTCaacaaac      >  2:1610897/1‑150 (MQ=255)
                                                                                                                                                  tGAGACGTGGGAGCAGCTGCAGCGCGGCCAGTTCCGCCAGCGCGGCTATATCTGGCTCTGCACGACCAAATCGGTCGTTGGTACGCATCACTTTACGTACATCAAACCACTGTGCCTGAACATCGCGTTCGCGCAGGATCTCaacaaaca     <  2:42735/150‑1 (MQ=255)
                                                                                                                                                    agaCGTGGGAGCAGCTGCAGCGCGGCCAGTTCCGCCAGCGCGGCTATATCTGGCTCTGCACGACCAAATCGGTCGTTGGTACGCATCACTTTACGTACATCAAACCACTGTGCCTGAACATCGCGTTCGCGCAGGATCTCAACAAAcagc   >  1:3023218/1‑150 (MQ=255)
                                                                                                                                                    agaCGTGGGAGCAGCTGCAGCGCGGCCAGTTCCGCCAGCGCGGCTATATCTGGCTCTGCACGACCAAATCGGTCGTTGGTACGCATCACTTTACGTACATCAAACCACTGTGCCTGAACATCGCGTTCGCGCAGGATCTCAACAAAcagc   >  1:136738/1‑150 (MQ=255)
                                                                                                                                                     gaCGTGGGAGCAGCTGCAGCGCGGCCAGTTCCGCCAGCGCGGCTATATCTGGCTCTGCACGACCAAATCGGTCGTTGGTACGCATCACTTTACGTACATCAAACCACTGTGCCTGAACATCGCGTTCGCGCAGGATCTCAACAAAcagca  >  1:51753/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GACGTCGGTCCAGATATCAACACGAGATGCGTGTAAAGCCTCCGCCAGCAAGGCTGCCGTATAATCGCTGCCTCCACGGCCAAGCGTCGTTGTACGACCTTTATTTTCGCTACCGATAAATCCCTGGGTGATCACTAAGCCTTCATTGAGACGTGGGAGCAGCTGCAGCGCGGCCAGTTCCGCCAGCGCGGCTATATCTGGCTCTGCACGACCAAATCGGTCGTTGGTACGCATCACTTTACGTACATCAAACCACTGTGCCTGAACATCGCGTTCGCGCAGGATCTCAACAAACAGCA  >  CP009273/4222493‑4222791

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: