Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,258,871 G→T 100% L173I (CTC→ATC)  tdcB ← L‑threonine dehydratase, catabolic

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,258,8710GT98.4% 805.5 / 2.4 251L173I (CTC→ATC) tdcBL‑threonine dehydratase, catabolic
Reads supporting (aligned to +/- strand):  ref base G (3/1);  new base T (130/116);  total (133/118)
Fisher's exact test for biased strand distribution p-value = 6.25e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.59e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

Only 100 of 255 total aligned reads displayed.

AGTGGAAAGAAGCCGCCATGCCGTGAACGTTTTCAGACTGTACGCCAATAACACGAATGGTCGGGTTAATAGATTTAATTGCCACCGCAATACCAGCAATTAAACCGCCACCACCAATTGGCACAATCACGTTATCGACATCATAGAGATCTTCCATAATTTCCAGACCAATCGTTCCCTGGCCAGCAATCACTTTCGGATCATCGTAAGGTGGGATAAAAATACGGCCTTCCATTTCGACAATTTCGCTCACTTTAGCGATAGTGTCGTTGAAGTTATCACCATGCAGAACG  >  CP009273/3258724‑3259016
                                                                                                                                                   |                                                                                                                                                 
aGTGGAAAGAAGCCGCCATGCCGTGAACGTTTTCAGACTGTACGCCAATAACACGAATGGTCGGGTTAATAGATTTAATTGCCACCGCAATACCAGCAATTAAACCGCCACCACCAATTGGCACAATCACGTTATCGACATCATAGatat                                                                                                                                                 <  2:4061498/150‑1 (MQ=255)
     aaaGAAGCCGCCATGCCGTGAACGTTTTCAGACTGTACGCCAATAACACGAATGGTCGGGTTAATAGATTTAATTGCCACCGCAATACCAGCAATTAAACCGCCACCACCAATTGGCACAATCACGTTATCGACATCATAGATATCTTcc                                                                                                                                            >  2:697808/1‑150 (MQ=255)
       agaagCCGCCATGCCGTGAACGTTTTCAGACTGTACGCCAATAACACGAATGGTCGGGTTAATAGATTTAATTGCCACCGCAATACCAGCAATTAAACCGCCACCACCAATTGGCACAATCACGTTATCGACATCATAGATATCTTCCAt                                                                                                                                          >  1:2362986/1‑150 (MQ=255)
        gaagCCGCCATGCCGTGAACGTTTTCAGACTGTACGCCAATAACACGAATGGTCGGGTTAATAGATTTAATTGCCACCGCAATACCAGCAATTAAACCGCCACCACCAATTGGCACAATCACGTTATCGACATCATAGATATCTTCCATa                                                                                                                                         >  2:502580/1‑150 (MQ=255)
          agCCGCCATGCCGTGAACGTTTTCAGACTGTACGCCAATAACACGAATGGTCGGGTTAATAGATTTAATTGCCACCGCAATACCAGCAATTAAACCGCCACCACCAATTGGCACAATCACGTTATCGACATCATAGATATCTTCCATAAt                                                                                                                                       >  1:5022535/1‑150 (MQ=255)
                cATGCCGTGAACGTTTTCAGACTGTACGCCAATAACACGAATGGTCGGGTTAATAGATTTAATTGCCACCGCAATACCAGCAATTAAACCGCCACCACCAATTGGCACAATCACGTTATCGACATCATAGATATCTTCCATAATTTCCAg                                                                                                                                 <  2:2003432/150‑1 (MQ=255)
                   gCCGTGAACGTTTTCAGACTGTACGCCAATAACACGAATGGTCGGGTTAATAGATTTAATTGCCACCGCAATACCAGCAATTAAACCGCCACCACCAATTGGCACAATCACGTTATCGACATCATAGATATCTTCCATAATTTCCAGAcc                                                                                                                              >  1:1577833/1‑150 (MQ=255)
                   gCCGTGAACGTTTTCAGACTGTACGCCAATAACACGAATGGTCGGGTTAATAGATTTAATTGCCACCGCAATACCAGCAATTAAACCGCCACCACCAATTGGCACAATCACGTTATCGACATCATAGATATCTTCCATAATTTCCAGAcc                                                                                                                              >  1:788735/1‑150 (MQ=255)
                     cGTGAACGTTTTCAGACTGTACGCCAATAACACGAATGGTCGGGTTAATAGATTTAATTGCCACCGCAATACCAGCAATTAAACCGCCACCACCAATTGGCACAATCACGTTATCGACATCATAGATATCTTCCATAATTTCCAGACCaa                                                                                                                            >  1:2448759/1‑150 (MQ=255)
                      gTGAACGTTTTCAGACTGTACGCCAATAACACGAATGGTCGGGTTAATAGATTTAATTGCCACCGCAATACCAGCAATTAAACCGCCACCACCAATTGGCACAATCACGTTATCGACATCATAGATATCTTCCATAATTTCCAGACCAAt                                                                                                                           >  2:2330402/1‑150 (MQ=255)
                      gTGAACGTTTTCAGACTGTACGCCAATAACACGAATGGTCGGGTTAATAGATTTAATTGCCACCGCAATACCAGCAATTAAACCGCCACCACCAATTGGCACAATCACGTTATCGACATCATAGATATCTTCCATAATTTCCAGACCAAt                                                                                                                           >  1:3769347/1‑150 (MQ=255)
                       tGAACGTTTTCAGACTGTACGCCAAGAACACGAATGGTCGGGTTAATAGATTTAATTGCCACCGCAATACCAGCAATTAAACCGCCACCACCAATTGGCACAATCACGTTATCGACATCATAGATATCTTCCATAATTTCCAGACCAATc                                                                                                                          <  2:2040477/150‑1 (MQ=255)
                         aaCGTTTTCAGACTGTACGCCAATAACACGAATGGTCGGGTTAATAGATTTAATTGCCACCGCAATACCAGCAATTAAACCGCCACCACCAATTGGCACAATCACGTTATCGACATCATAGATATCTTCCATAATTTCCAGACCAATCGt                                                                                                                        >  1:4861772/1‑150 (MQ=255)
                          aCGTTTTCAGACTGTACGCCAATAACACGAATGGTCGGGTTAATAGATTTAATTGCCACCGCAATACCAGCAATTAAACCGCCACCACCAATTGGCACAATCACGTTATCGACATCATAGATATCTTCCATAATTTCCAGACCAATCGtt                                                                                                                       <  2:158226/150‑1 (MQ=255)
                             ttttCAGACTGTACGCCAATAACACGAATGGTCGGGTTAATAGATTTAATTGCCACCGCAATACCAGCAATTAAACCGCCACCACCAATTGGCACAATCACGTTATCGACATCATAGATATCTTCCATACTTTCCAGACCAATCGTTccc                                                                                                                    >  2:3683639/1‑150 (MQ=255)
                              tttCAGACTGTACGCAAATAACACGAATGGTCGGGTTAATAGATTTAATTGCCACCGCAATACCAGCAATTAAACCGCCACCACCAATTGGCACAATCACGTTATCGACATCATAGATATCTTCCATAATTTCCAGACCAATCGTTCCCt                                                                                                                   <  1:2540747/150‑1 (MQ=255)
                               ttCAGACTGTACGCCAATAACACGAATGGTCGGGTTAATAGATTTAATTGCCACCGCAATACCAGCAATTAAACCGCCACCACCAATTGGCACAATCACGTTATCGACATCATAGATATCTTCCATAATTTCCAGACCAATCGTTCCCTg                                                                                                                  <  2:1612064/150‑1 (MQ=255)
                                   gACTGTACGCCAATAACACGAATGGTCGGGTTAATAGATTTAATTGCCACCGCAATACCAGCAATTAAACCGCCACCACCAATTGGCACAATCACGTTATCGACATCATAGATATCTTCCATAATTTCCAGACCAATCGTTCCCTGGCCa                                                                                                              >  2:4015243/1‑150 (MQ=255)
                                     cTGTACGCCAATAACACGAATGGTCGGGTTAATAGATTTAATTGCCACCGCAATACCAGCAATTAAACCGCCACCACCAATTGGCACAATCACGTTATCGACATCATAGATATCTTCCATAATTTCCAGACCAATCGTTCCCTGGCCAGc                                                                                                            >  1:2328773/1‑150 (MQ=255)
                                      tGTACGCCAATAACACGAATGGTCGGGTTAATAGATTTAATTGCCACCGCAATACCAGCAATTAAACCGCCACCACCAATTGGCACAATCACGTTATCGACATCATAGATATCTTCCATAATTTCCAGACCAATCGTTCCCTGGCCAGCa                                                                                                           >  2:2709692/1‑150 (MQ=255)
                                       gTACGCCAATAACACGAATGGTCGGGTTAATAGATTTAATTGCCACCGCAATACCAGCAATTAAACCGCCACCACCAATTGGCACAATCACGTTATCGACATCATAGATATCTTCCATAATTTCCAGACCAATCGTTCCCTGGCCAGCaa                                                                                                          <  1:185227/150‑1 (MQ=255)
                                        tACGCCAATAACACGAATGGTCGGGTTAATAGATTTAATTGCCACCGCAATACCAGCAATTAAACCGCCACCACCAATTGGCACAATCACGTTATCGACATCATAGATATCTTCCATAATTTCCAGACCAATCGTTCCCTGGCCAGCAAt                                                                                                         >  1:4815668/1‑150 (MQ=255)
                                          cGCCAATAACACGAATGGTCGGGTTAATAGATTTAATTGCCACCGCAATACCAGCAATTAAACCGCCACCACCAATTGGCACAATCACGTTATCGACATCATAGATATCTTCCATAATTTCCAGACCAATCGTTCCCTGGCCAGCaaatc                                                                                                       >  1:256893/1‑147 (MQ=255)
                                          cGCCAATAACACGAATGGTCGGGTTAATAGATTTAATTGCCACCGCAATACCAGCAATTAAACCGCCACCACCAATTGGCACAATCACGTTATCGACATCATAGATATCTTCCATAATTTCCAGACCAATCGTTCCCTGGCCAGCAATCa                                                                                                       >  2:3740976/1‑150 (MQ=255)
                                           gCCAATAACACGAATGGTCGGGTTAATAGATTTAATTGCCACCGCAATACCAGCAATTAAACCGCCACCACCAATTGGCACAATCACGTTATCGACATCATAGATATCTTCCATAATTTCCAGACCAATCGTTCCCTGGCCAGCAATCAc                                                                                                      >  2:3171047/1‑150 (MQ=255)
                                           gCCAATAACACGAATGGTCGGGTTAATAGATTTAATTGCCACCGCAATACCAGCAATTAAACCGCCACCACCAATTGGCACAATCACGTTATCGACATCATAGATATCTTCCATAATTTCCAGACCAATCGTTCCCTGGCCAGCAATCAc                                                                                                      >  1:4387450/1‑150 (MQ=255)
                                           gCCAATAACACGAATGGTCGGGTTAATAGATTTAATTGCCACCGCAATACCAGCAATTAAACCGCCACCACCAATTGGCACAATCACGTTATCGACATCATAGAGATCTTCCATAATTTCCAGACCAATCGTTCCCTGGCCAGCAATCAc                                                                                                      >  2:3027134/1‑150 (MQ=255)
                                            ccAATAACACGAATGGTCGGGTTAATAGATTTAATTGCCACCGCAATACCAGCAATTAAACCGCCACCACCAATTGGCACAATCACGTTATCGACATCATAGATATCTTCCATAATTTCCAGACCAATCGTTCCCTGGCCAGCAATCACt                                                                                                     >  2:4859325/1‑150 (MQ=255)
                                            ccAATAACACGAATGGTCGGGTTAATAGATTTAATTGCCACCGCAATACCAGCAATTAAACCGCCACCACCAATTGGCACAATCACGTTATCGACATCATAGATATCTTCCATAATTTCCAGACCAATCGTTCCCTGGCCAGCAATCACt                                                                                                     <  2:3927495/150‑1 (MQ=255)
                                             cAATAACACGAATGGTCGGGTTAATAGATTTAATTGCCACCGCAATACCAGCAATTAAACCGCCACCACCAATTGGCACAATCACGTTATCGACATCATAGATATCTTCCATAATTTCCAGACCAATCGTTCCCTGGCCAGCAATCACtt                                                                                                    >  1:1046920/1‑150 (MQ=255)
                                             cAATAACACGAATGGTCGGGTTAATAGATTTAATTGCCACCGCAATACCAGCAATTAAACCGCCACCACCAATTGGCACAATCACGTTATCGACATCATAGATATCTTCCATAATTTCCAGACCAATCGTTCCCTGGCCAGCAATCACtt                                                                                                    >  1:2181135/1‑150 (MQ=255)
                                              aaTAACACGAATGGTCGGGTTAATAGATTTAATTGCCACCGCAATACCAGCAATTAAACCGCCACCACCAATTGGCACAATCACGTTATCGACATCATAGAGATCTTCCATAATTTCCAGACCAATCGTTCCCTGGCCAGCAATCACttt                                                                                                   >  1:837763/1‑150 (MQ=255)
                                               aTAACACGAATGGTCGGGTTAATAGATTTAATTGCCACCGCAATACCAGCAATTAAACCGCCACCACCAATTGGCACAATCACGTTATCGACATCATAGATATCTTCCATAATTTCCAGACCAATCGTTCCCTGGCCAGCAATCACTTTc                                                                                                  <  1:4746447/150‑1 (MQ=255)
                                                tAACACGAATGGTCGGGTTAATAGATTTAATTGCCACCGCAATACCAGCAATTAAACCGCCACCACCAATTGGCACAATCACGTTATCGACATCATAGATATCTTCCATAATTTCCAGACCAATCGTTCCCTGGCCAGCAATCACTTTCg                                                                                                 >  1:1496632/1‑150 (MQ=255)
                                                 aaCACGAATGGTCGGGTTAATAGATTTAATTGCCACCGCAATACCAGCAATTAAACCGCCACCACCAATTGGCACAATCACGTTATCGACATCATAGATATCTTCCATACTTTCCAGACCAATCGTTCCCTGGCCAGCAATCACTTTCgg                                                                                                >  2:4185720/1‑150 (MQ=255)
                                                 aaCACGAATGGTCGGGTTAATAGATTTAATTGCCACCGCAATACCAGCAATTAAACCGCCACCACCAATTGGCACAATCACGTTATCGACATCATAGATATCTTCCATAATTTCCAGACCAATCGTTCCCTGGCCAGCAATCACTTTCgg                                                                                                <  1:3730483/150‑1 (MQ=255)
                                                       aaTGGTCGGGTTAATAGATTTAATTGCCACCGCAATACCAGCAATTAAACCGCCACCACCAATTGGCACAATCACGTTATCGACATCATAGATATCTTCCATAATTTCCAGACCAATCGTTCCCTGGCCAGCAATCACTTTCGGatcatc                                                                                          <  1:1273326/150‑1 (MQ=255)
                                                        aTGGTCGGGTTAATAGATTTAATTGCCACCGCAATACCAGCAATTAAACCGCCACCACCAATTGGCACAATCACGTTATCGACATCATAGATATCTTCCATAATTTCCAGACCAATCGTTCCCTGGCCAGCAATCACTTTCGGATCATCg                                                                                         <  1:200675/150‑1 (MQ=255)
                                                          ggTCGGGTTAATAGATTTAATTGCCACCGCAATACCAGCAATTAAACCGCCACCACCAATTGGCACAATCACGTTATCGACATCATAGATATCTTCCATAATTTCCAGACCAATCGTTCCCTGGCCAGCAATCACTTTCGGATCATCGTa                                                                                       >  2:404558/1‑150 (MQ=255)
                                                          ggTCGGGTTAATAGATTTAATTGCCACCGCAATACCAGCAATTAAACCGCCACCACCAATTGGCACAATCACGTTATCGACATCATAGATATCTTCCATAATTTCCAGACCAATCGTTCCCTGGCCAGCAATCACTTTCGGATCATCGTa                                                                                       <  1:3102937/150‑1 (MQ=255)
                                                            tCGGGTTAATAGATTTAATTGCCACCGCAATACCAGCAATTAAACCGCCACCACCAATTGGCACAATCACGTTATCGACATCATAGATATCTTCCATAATTTCCAGACCAATCGTTCCCTGGCCAGCAATCACTTTCGGATCATCGTAAg                                                                                     >  2:1359467/1‑150 (MQ=255)
                                                            tCGGGTTAATAGATTTAATTGCCACCGCAATACCAGCAATTAAACCGCCACCACCAATTGGCACAATCACGTTATCGACATCATAGATATCTTCCATAATTTCCAGACCAATCGTTCCCTGGCCAGCAATCACTTTCGGATCATCGTAAg                                                                                     >  2:2879225/1‑150 (MQ=255)
                                                            tCGGGTTAATAGATTTAATTGCCACCGCAATACCAGCAATTAAACCGCCACCACCAATTGGCACAATCACGTTATCGACATCATAGATATCTTCCATAATTTCCAGACCAATCGTTCCCTGGCCAGCAATCACTTTCGGATCATCGTAAg                                                                                     >  2:186676/1‑150 (MQ=255)
                                                            tCGGGTTAATAGATTTAATTGCCACCGCAATACCAGCAATTAAACCGCCACCACCAATTGGCACAATCACGTTATCGACATCATAGATATCTTCCATAATTTCCAGACCAATCGTTCCCTGGCCAGCAATCACTTTCGGATCATCGTAAg                                                                                     <  2:3859054/150‑1 (MQ=255)
                                                             cGGGTTAATAGATTTAATTGCCACCGCAATACCAGCAATTAAACCGCCACCACCAATTGGCACAATCACGTTATCGACATCATAGATATCTTCCATAATTTCCAGACCAATCGTTCCCTGGCCAGCAATCACTTTCGGATCATCGTAAgg                                                                                    >  1:3400118/1‑150 (MQ=255)
                                                                     tAGATTTAATTGCCACCGCAATACCAGCAATTAAACCGCCACCACCAATTGGCACAATCACGTTATCGACATCATAGATATCTTCCATAATTTCCAGACCAATCGTTCCCTGGCCAGCAATCACTTTCGGATCATCGTAAGGTGGGATaa                                                                            >  1:213209/1‑150 (MQ=255)
                                                                       gATTTAATTGCCACCGCAATACCAGCAATTAAACCGCCACCACCAATTGGCACAATCACGTTATCGACATCATAGATATCTTCCATAATTTCCAGACCAATCGTTCCCTGGCCAGCAATCACTTTCGGATCATCGTAAGGTGGGATaaaa                                                                          >  2:4615600/1‑150 (MQ=255)
                                                                          ttAATTGCCACCGCAATACCAGCAATTAAACCGCCACCACCAATTGGCACAATCACGTTATCGACATCATAGATATCTTCCATAATTTCCAGACCAATCGTTCCCTGGCCAGCAATCACTTTCGGATCCTCGTAAGGTGGGATAAAAATa                                                                       <  2:4698957/150‑1 (MQ=255)
                                                                           tAATTGCCACCGCAATACCAGCAATTAAACCGCCACCACCAATTGGCACAATCACGTTATCGACATCATAGATATCTTCCATAATTTCCAGACCAATCGTTCCCTGGCCAGCAATCACTTTCGGATCATCGTAAGGTGGGATAAAAATAc                                                                      >  2:4483791/1‑150 (MQ=255)
                                                                            aaTTGCCACCGCAATACCAGCAATTAAACCGCCACCACCAATTGGCACAATCACGTTATCGACATCATAGATATCTTCCATAATTTCCAGACCAATCGTTCCCTGGCCAGCAATCACTTTCGGATCATCGTAAGGTGGGATAAAAATACg                                                                     >  2:3413925/1‑150 (MQ=255)
                                                                              ttGCCACCGCAATACCAGCAATTAAACCGCCACCACCAATTGGCACAATCACGTTATCGACATCATAGATATCTTCCATAATTTCCAGACCAATCGTTCCCTGGCCAGCAATCACTTTCGGATCATCGTAAGGTGGGATAAAAATACGGc                                                                   <  2:1719168/150‑1 (MQ=255)
                                                                               tGCCACCGCAATACCAGCAATTAAACCGCCACCACCAATTGGCACAATCACGTTATCGACATCATAGATATCTTCCATAATTTCCAGACCAATCGTTCCCTGGCCAGCAATCACTTTCGGATCATCGTAAGGTGGGATAAAAATACGGcc                                                                  <  2:2469076/150‑1 (MQ=255)
                                                                                 ccACCGCAATACCAGCAATTAAACCGCCACCACCAATTGGCACAATCACGTTATCGACATCATAGATATCTTCCATAATTTCCAGACCAATCGTTCCCTGGCCAGCAATCACTNTCGGATCATCGTAAGGTGGGATAAAAATACGGCCtt                                                                <  1:2675958/150‑1 (MQ=255)
                                                                                  cACCGCAATACCAGCAATTAAACCGCCACCACCAATTGGCACAATCACGTTATCGACATCATAGATATCTTCCATAATTTCCAGACCAATCGTTCCCTGGCCAGCAATCACTTTCGGATCATCGTAAGGTGGGATAAAAATACGGCCTTc                                                               <  2:1266106/150‑1 (MQ=255)
                                                                                  cACCGCAATACCAGCAATTAAACCGCCACCACCAATTGGCACAATCACGTTATCGACATCATAGATATCTTCCATAATTTCCAGAACAATCGTTCCCTGGCCAGCACTCACTTTCGGATCATCGTAAGGTGGGATAAAAATACGGCCTTc                                                               >  1:5070562/1‑150 (MQ=255)
                                                                                   aCCGCAATACCAGCAATTAAACCGCCACCACCAATTGGCACAATCACGTTATCGACATCATAGATATCTTCCATAATTTCCAGACCAATCGTTCCCTGGCCAGCAATCACTTTCGGATCATCGTAAGGTGGGATAAAAATACGGCCTTcc                                                              <  1:3633874/150‑1 (MQ=255)
                                                                                   aCCGCAATACCAGCAATTAAACCGCCACCACCAATTGGCACAATCACGTTATCGACATCATAGATATCTTCCATAATTTCCAGACCAATCGTTCCCTGGCCAGCAATCACTTTCGGATCATCGTAAGGTGGGATAAAAATACGGCCTTcc                                                              <  1:4240381/150‑1 (MQ=255)
                                                                                     cGCAATACCAGCAATTAAACCGCCACCACCAATTGGCACAATCACGTTATCGACATCATAGATATCTTCCATAATTTCCAGACCAATCGTTCCCTGGCCAGCAATCACTTTCGGATCATCGTAAGGTGGGATAAAAATACGGCCTTCCAt                                                            >  2:1426100/1‑150 (MQ=255)
                                                                                     cGCAATACCAGCAATTAAACCGCCACCACCAATTGGCACAATCACGTTATCGACATCATAGATATCTTCCATAATTTCCAGACCAATCGTTCCCTGGCCAGCAATCACTTTCGGATCATCGTAAGGTGGGATAAAAATACGGCCTTCCAt                                                            >  2:606902/1‑150 (MQ=255)
                                                                                      gCAATACCAGCAATTAAACCGCCACCACCAATTGGCACAATCACGTTATCGACATCATAGATATCTTCCATAATTTCCAGACCAATCGTTCCCTGGCCAGCAATCACTTTCGGATCATCGTAAGGTGGGATAAAAATACGGCCTTCCAtt                                                           >  2:4352697/1‑150 (MQ=255)
                                                                                       cAATACCAGCAATTAAACCGCCACCACCAATTGGCACAATCACGTTATCGACATCATAGATATCTTCCATAATTTCCAGACCAATCGTTCCCTGGCCAGCAATCACTTTCGGATCATCGTAAGGTGGGATAAAAATACGGCCTTCCAttt                                                          <  1:4689823/150‑1 (MQ=255)
                                                                                         aTACCAGCAATTAAACCGCCACCACCAATTGGCACAATCACGTTATCGACATCATAGATATCTTCCATAATTTCCAGACCAATCGTTCCCTGGCCAGCAATCACTTTCGGATCATCGTAAGGTGGGATAAAAATACGGCCTTCCATTTCg                                                        >  2:4457676/1‑150 (MQ=255)
                                                                                          tACCAGCAATTAAACCGCCACCACCAATTGGCACAATCACGTTATCGACATCATAGATATCTTCCATAATTTCCAGACCAATCGTTCCCTGGCCAGCAATCACTTTCGGATCATCGTAAGGTGGTATAAAAATACGGCCTTCCATTTCGa                                                       >  2:2626141/1‑150 (MQ=255)
                                                                                             cAGCAATTAAACCGCCACCACCAATTGGCACAATCACGTTATCGACATCATAGATATCTTCCATAATTTCCAGACCAATCGTTCCCTGGCCAGCAATCACTTTCGGATCATCGTAAGGTGGGATAAAAATACGGCCTTCCATTTCGACaa                                                    >  1:4173478/1‑150 (MQ=255)
                                                                                              aGCAATTAAACCGCCACCACCAATTGGCACAATCACGTTATCGACATCATAGATATCTTCCATAATTTCCAGACCAATCGTTCCCTGGCCAGCAATCACTTTCGGATCATCGTAAGGTGGGATAAAAATACGGCCTTCCATTTCGACAAt                                                   >  2:784051/1‑150 (MQ=255)
                                                                                                   ttAAACCGCCACCACCAATTGGCACAATCACGTTATCGACATCATAGATATCTTCCATAATTTCCAGACCAATCGTTCCCTGGCCAGCAATCACTTTCGGATCATCGTAAGGTGGGATAAAAATACGGCCTTCCATTTCGACAATTTCGc                                              <  2:4499433/150‑1 (MQ=255)
                                                                                                     aaaCCGCCACCACCAATTGGCACAATCACGTTATCGACATCATAGATATCTTCCATAATTTCCAGACCAATCGTTCCCTGGCCAGCAATCACTTTCGGATCATCGTAAGGTGGGATAAAAATACGGCCTTCCATTTCGACAATTTCGCTc                                            <  2:1514084/150‑1 (MQ=255)
                                                                                                      aaCCGCCACCACCAATTGGCACAATCACGTTATCGACATCATAGATATCTTCCATAATTTCCAGACCAATCGTTCCCTGGCCAGCAATCACTTTCGGATCATCGTAAGGTGGGATAAAAATACGGCCTTCCATTTCGACAATTTCGCTCa                                           <  2:218835/150‑1 (MQ=255)
                                                                                                        ccGCCACCACCAATTGGCACAATCACGTTATCGACATCATAGATATCTTCCATAATTTCCAGACCAATCGTTCCCTGGCCAGCAATCACTTTCGGATCATCGTAAGGTGGGATAAAAATACGGCCTTCCATTTCGACAATTTCGCTCACt                                         <  2:4218372/150‑1 (MQ=255)
                                                                                                        ccGCCACCACCAATTGGCACAATCACGTTATCGACATCATAGATATCTTCCATAATTTCCAGACCAATCGTTCCCTGGCCAGCAATCACTTTCGGATCATCGTAAGGTGGGATAAAAATACGGCCTTCCATTTCGACAATTTCGCTCACt                                         >  2:2116993/1‑150 (MQ=255)
                                                                                                          gCCACCACCAATTGGCACAATCACGTTATCGACATCATAGATATCTTCCATAATTTCCAGACCACTCGTTCCCTGGCCAGCAATCACTTTCGGATCATCGTAAGGTGGGATAAAAATACGGCCTTCCATTTCGACAATTTCGCTCACttt                                       <  2:3251475/150‑1 (MQ=255)
                                                                                                          gCCACCACCAATTGGCACAATCACGTTATCGACATCATAGATATCTTCCATAATTTCCAGACCAATCGTTCCCTGGCCAGCAATCACTTTCGGATCATCGTAAGGTGGGATAAAAATACGGCCTTCCATTTCGACAATTTCGCTCACttt                                       >  2:738080/1‑150 (MQ=255)
                                                                                                          gCCACCACCAATTGGCACAATCACGTTATCGACATCATAGATATCTTCCATAATTTCCAGACCAATCGTTCCCTGGCCAGCAATCACTTTCGGATCATCGTAAGGTGGGATAAAAATACGGCCTTCCATTTCGACAATTTCGCTCACttt                                       >  2:4956021/1‑150 (MQ=255)
                                                                                                           ccaccaccaATTGGCACAATCACGTTATCGACATCATAGATATCTTCCATAATTTCCAGACCAATCGTTCCCTGGCCAGCAATCACTTTCGGATCATCGTAAGGTGGGATAAAAATACGGCCTTCCATTTCGACAATTTCGCTCACTTTa                                      <  1:3519542/150‑1 (MQ=255)
                                                                                                           ccaccaccaATTGGCACAATCACGTTATCGACATCATAGATATCTTCCATAATTTCCAGACCAATCGTTCCCTGGCCAGCAATCACTTTCGGATCATCGTAAGGTGGGATAAAAATACGGCCTTCCATTTCGACAATTTCGCTCACTTTa                                      >  1:304068/1‑150 (MQ=255)
                                                                                                           ccaccaccaATTGGCACAATCACGTTATCGACATCATAGATATCTTCCATAATTTCCAGACCAATCGTTCCCTGGCCAGCAATCACTTTCGGATCATCGTAAGGTGGGATAAAAATACGGCCTTCCATTTCGACAATTTCGCTCACTTTa                                      <  2:1042694/150‑1 (MQ=255)
                                                                                                           ccaccaccaATTGGCACAATCACGTTATCGACATCATAGATATCTTCCATAATTTCCAGACCAATCGGTCCCTGGCCAGCAATCACTTTCGGATCATCGTAAGGTGGGATAAAAATACGGCCTTCCATTTCGACAATTTCGCTCACTTTa                                      <  1:558225/150‑1 (MQ=255)
                                                                                                             accaccaATTGGCACAATCACGTTATCGACATCATAGATATCTTCCATAATTTCCAGACCAATCGTTCCCTGGCCAGCAATCACTTTCGGATCATCGTAAGGTGGGATAAAAATACGGCCTTCCATTTCGACAATTTCGCTCACTTTAGc                                    <  2:1039157/150‑1 (MQ=255)
                                                                                                               caccaATTGGCACAATCACGTTATCGACATCATAGATATCTTCCATAATTTCCAGACCAATCGTTCCCTGGCCAGCAATCACTTTCGGATCATCGTAAGGTGGGATAAAAATACGGCCTTCCATTTCGACAATTTCGCTCACTTTAGCGa                                  >  1:4216355/1‑150 (MQ=255)
                                                                                                                accaATTGGCACAATCACGTTATCGACATCATAGATATCTTCCATAATTTCCAGACCAATCGTTCCCTGGCCAGCAATCACTTTCGGATCATCGTAAGGTGGGATAAAAATACGGCCTTCCATTTCGACAATTTCGCTCACTTTAGCGAt                                 >  1:4651078/1‑150 (MQ=255)
                                                                                                                   aaTTGGCACAATCACGTTATCGACATCATAGATATCTTCCATAATTTCCAGACCAATCGTTCCCTGGCCAGCAATCACTTTCGGATCATCGTAAGGTGGGATAAAAATACGGCCTTCCATTTCGACAATTTCGCTCACTTTAGCGATAgt                              <  1:490656/150‑1 (MQ=255)
                                                                                                                         cacaATCACGTTATCGACATCATAGATATCTTCCATAATTTCCAGACCAATCGTTCCCTGGCCAGCAATCACTTTCGGATCATCGTAAGGTGGGATAAAAATACGGCCTTCCATTTCGACAATTTCGCTCACTTTAGCGATAGTGTCGtt                        <  1:5100933/150‑1 (MQ=255)
                                                                                                                         cacaATCACGTTATCGACATCATAGATATCTTCCATAATTTCCAGACCAATCGTTCCCTGGCCAGCAATCACTTTCGGATCATCGTAAGGTGGGATAAAAATACGGCCTTCCATTTCGACAATTTCGCTCACTTTAGCGATAGTGTCGtt                        <  2:150935/150‑1 (MQ=255)
                                                                                                                         cacaATCACGTTATCGACATCATAGATATCTTCCATAATTTCCAGACCAATCGTTCCCTGGCCAGCAATCACTTTCGGATCATCGTAAGGTGGGATAAAAATACGGCCTTCCATTTCGACAATTTCGCTCACTTTAGCGATAGTGTCGtt                        <  1:2709692/150‑1 (MQ=255)
                                                                                                                          acaATCACGTTATCGACATCATAGATATCTTCCATAATTTCCAGACCAATCGGTCCCTGGCCAGCAATCACTTTCGGATCATCGTAAGGTGGGATAAAAATACGGCCTTCCATTTCGACAATTTCGCTCACTTTAGCGATAGTGTCGTTg                       <  2:174674/150‑1 (MQ=255)
                                                                                                                            aaTCACGTTATCGACATCATAGATATCTTCCATAATTTCCAGACCAATCGTTCCCTGGCCAGCAATCACTTTCGGATCATCGTAAGGTGGGATAAAAATACGGCCTTCCATTTCGACAATTTCGCTCACTTTAGCGATAGTGTCGTTGaa                     >  1:2761236/1‑150 (MQ=255)
                                                                                                                            aaTCACGTTATCGACATCATAGATATCTTCCATAATTTCCAGACCAATCGTTCCCTGGCCAGCAATCACTTTCGGATCATCGTAAGGTGGGATAAAAATACGGCCTTCCATTTCGACAATTTCGCTCACTTTAGCGATAGTGTCGTTGaa                     <  2:4966131/150‑1 (MQ=255)
                                                                                                                            aaTCACGTTATCGACATCATAGATATCTTCCATAATTTCCAGACCAATCGTTCCCTGGCCAGCAATCACTTTCGGATCATCGTAAGGTGGGATAAAAATACGGCCTTCCATTTCGACAATTTCGCTCACTTTAGCGATAGTGTCGTTGaa                     >  1:1684552/1‑150 (MQ=255)
                                                                                                                            aaTCACGTTATCGACATCATAGATATCTTCCATAATTTCCAGACCAATCGTTCCCTGGCCAGCAATCACTTTCGGATCATCGTAAGGTGGGATAAAAATAAGGCCTTCCATTTCGACAATTTCGCTCACTTTAGAGATAGTGTCGTTGaa                     >  1:1988640/1‑150 (MQ=255)
                                                                                                                                aCGTTATCGACATCATAGATATCTTCCATAATTTCCAGACCAATCGTTCCCTGGCCAGCAATCACTTTCGGATCATCGTAAGGTGGGATAAAAATACGGCCTTCCATTTCGACAATTTCGCTCACTTTAGCGATAGTGTCGTTGAAGTTa                 <  1:4457676/150‑1 (MQ=255)
                                                                                                                                aCGTTATCGACATCATAGATATCTTCCATAATTTCCAGACCAATCGTTCCCTGGCCAGCAATCACTTTCGGATCATCGTAAGGTGGGATAAAAATACGGCCTTCCATTTCGACAATTTCGCTCACTTTAGCGATAGTGTCGTTGAAGTTa                 >  1:118228/1‑150 (MQ=255)
                                                                                                                                   ttATCGACATCATAGATATCTTCCATAATTTCCAGACCACTCGTTCCCTGGCCAGCAATCACTTTCGGATCATCGTAAGGTGGGATAAAAATACGGCCTTCCATTTCGACAATTTCGCTCACTTTAGCGATAGTGTCGTTGAAGTTATCa              <  1:3236246/150‑1 (MQ=255)
                                                                                                                                       cGACATCATAGATATCTTCCATAATTTCCAGACCAATCGTTCCCTGGCCAGCAATCACTTTCGGATCATCGTAAGGTGGGATAAAAATACGGCCTTCCATTTCGACAATTTCGCTCACTTTAGCGATAGTGTCGTTGAAGTTATCACCAt          <  2:24766/150‑1 (MQ=255)
                                                                                                                                       cGACATCATAGATATCTTCCATAATTTCCAGACCAATCGTTCCCTGGCCAGCAATCACTTTCGGATCATCGTAAGGTGGGATAAAAATACGGCCTTCCATTTCGACAATTTCGCTCACTTTAGCGATAGTGTCGTTGAAGTTATCACCAt          <  1:4961284/150‑1 (MQ=255)
                                                                                                                                       cGACATCATAGATATCTTCCATAATTTCCAGACCAATCGTTCCCTGGCCAGCAATCACTTTCGGATCATCGTAAGGTGGGATAAAAATACGGCCTTCCATTTCGACAATTTCGCTCACTTTAGCGATAGTGTCGTTGAAGTTATCACCAt          >  2:3432864/1‑150 (MQ=255)
                                                                                                                                        gACATCATAGATATCTTCCATAATTTCCAGACCAATCGTTCCCTGGCCAGCAATCACTTTCGGATCATCGTAAGGTGGGATAAAAATACGGCCTTCCATTTCGACAATTTCGCTCACTTTAGCGATAGTGTCGTTGAAGTTATCACCATg         <  2:4002214/150‑1 (MQ=255)
                                                                                                                                           atcatAGATATCTTCCATAATTTCCAGACCAATCGTTCCCTGGCCAGCAATCACTTTCGGATCATCGTAAGGTGGGATAAAAATACGGCCTTCCATTTCGACAATTTCGCTCACTTTAGCGATAGTGTCGTTGAAGTTATCACCATGCAg      <  2:4928338/150‑1 (MQ=255)
                                                                                                                                              atAGATATCTTCCATAATTTCCAGACCAATCGTTCCCTGGCCAGCAATCACTTTCGGATCATCGTAAGGTGGGATAAAAATACGGCCTTCCATTTCGACAATTTCGCTCACTTTAGCGATAGTGTCGTTGAAGTTATCACCATGCAGAAc   >  2:602370/1‑150 (MQ=255)
                                                                                                                                              atAGATATCTTCCATAATTTCCAGACCAATCGTTCCCTGGCCAGCAATCACTTTCGGATCATCGTAAGGTGGGATAAAAATACGGCCTTCCATTTCGACAATTTCGCTCACTTTAGCGATAGTGTCGTTGAAGTTATCACCATGCAGAAc   <  2:676138/150‑1 (MQ=255)
                                                                                                                                              atAGATATCTTCCATAATTTCCAGACCAATCGTTCCCTGGCCAGCAATCACTTTCGGATCATCGTAAGGTGGGATAAAAATACGGCCTTCCATTTCGACAATTTCGCTCACTTTAGCGATAGTGTCGTTGAAGTTATCACCATGCAGAAc   <  1:1359467/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
AGTGGAAAGAAGCCGCCATGCCGTGAACGTTTTCAGACTGTACGCCAATAACACGAATGGTCGGGTTAATAGATTTAATTGCCACCGCAATACCAGCAATTAAACCGCCACCACCAATTGGCACAATCACGTTATCGACATCATAGAGATCTTCCATAATTTCCAGACCAATCGTTCCCTGGCCAGCAATCACTTTCGGATCATCGTAAGGTGGGATAAAAATACGGCCTTCCATTTCGACAATTTCGCTCACTTTAGCGATAGTGTCGTTGAAGTTATCACCATGCAGAACG  >  CP009273/3258724‑3259016

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: