Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,359,673 G→T 58.7% M89I (ATG→ATT arnB → uridine 5'‑(beta‑1‑threo‑pentapyranosyl‑4‑ulose diphosphate) aminotransferase, PLP‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,359,6730GT58.7% 113.5 / inf 303M89I (ATG→ATTarnBuridine 5'‑(beta‑1‑threo‑pentapyranosyl‑4‑ulose diphosphate) aminotransferase, PLP‑dependent
Reads supporting (aligned to +/- strand):  ref base G (73/52);  new base T (92/86);  total (165/138)
Fisher's exact test for biased strand distribution p-value = 2.92e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.87e-01

Only 100 of 303 total aligned reads displayed.

AAGCTTTTTGCCAGTTGACGGGAAATCAGCATGCCATCGCGGTCAGTTCAGCCACCGCCGGAATGCATATCACGCTAATGGCGTTGAAAATTGGCAAGGGCGATGAAGTGATTACGCCTTCCCTGACCTGGGTTTCAACCCTCAATATGATTTCCTTGTTGGGTGCAACGCCGGTAATGGTGGATGTCGACCGCGATACGCTGATGGTCACGCCTGAAGCTATCGAGTCAGCCATTACGCCACGCACTAAAGCCATCATTCCGGTGCATTATGCCGGTGCGCCAGCAGATATTGACG  >  CP009273/2359525‑2359821
                                                                                                                                                    |                                                                                                                                                    
aaGCTTTTTGCCAGTTGACGGGAAATCAGCATGACATCGCGGTCAGTTCAGCCACCGCCGGAATGCATATCACGCTAATGGCGTTGAAAATTGGCAAGGGCGATGAAGTGATTACGCCTTCCCTGACCTGGGTTTCAACCCTCAATATGa                                                                                                                                                     >  1:1291353/1‑150 (MQ=255)
    tttttGCCAGTTGACGGGAAATCAGCATGCCATCGCGGTCAGTTCAGCCACCGCCGGAATGCATATCACGCTAATGGCGTTGAAAATTGGCAAGGGCGATGAAGTGATTACGCCTTCCCTGACCTGGGTTTCAACCCTCAATATGATTTc                                                                                                                                                 >  1:4154467/1‑150 (MQ=255)
       ttGCCAGTTGACGGGAAATCAGCATGCCATCGCGGTCAGTTCAGCCACCGCCGGAATGCATATCACGCTAATGGCGTTGAAAATTGGCAAGGGCGATGAAGTGATTACGCCTTCCCTGACCTGGGTTTCAACCCTCAATATTATTTCCtt                                                                                                                                              >  1:4377713/1‑150 (MQ=255)
       ttGCCAGTTGACGGGAAATCAGCATGCCATCGCGGTCAGTTCAGCCACCGCCGGAATGCATATCACGCTAATGGCGTTGAAAATTGGCAAGGGCGATGAAGTGATTACGCCTTCCCTGACCTGGGTTTCAACCCTCAATATTATTTCCtt                                                                                                                                              <  2:3038554/150‑1 (MQ=255)
       ttGCCAGTTGACGGGAAATCAGCATGCCATCGCGGTCAGTTCAGCCACCGCCGGAATGCATATCACGCTAATGGCGTTGAAAATTGGCAAGGGCGATGAAGTGATTACGCCTTCCCTGACCTGGGTTTCAACCCTCAATATTATTTCCtt                                                                                                                                              <  1:2022392/150‑1 (MQ=255)
        tGCCAGTTGACGGGAAATCAGCATGCCATCGCGGTCAGTTCAGCCACCGCCGGAATGCATATCACGCTAATGGCGTTGAAAATTGGCAAGGGCGATGAAGTGATTACGCCTTCCCTGACCTGGGTTTCAACCCTCAATATTATTTCCttg                                                                                                                                             <  2:729958/150‑1 (MQ=255)
        tGCCAGTTGACGGGAAATCAGCATGCCATCGCGGTCAGTTCAGCCACCGCCGGAATGCATATCACGCTAATGGCGTTGAAAATTGGCAAGGGCGATGAAGTGATTACGCCTTCCCTGACCTGGGTTTCAACCCTCAATATTATTTCCttg                                                                                                                                             >  1:2712560/1‑150 (MQ=255)
          ccAGTTGACGGGAAATCAGCATGCCATCGCGGTCAGTTCAGCCACCGCCGGAATGCATATCACGCTAATGGCGTTGAAAATTGGCAAGGGCGATGAAGTGATTACGCCTTCCCTGACCTGGGTTTCAACCCTCAATATGATTTCCttgtt                                                                                                                                           <  1:2613253/150‑1 (MQ=255)
          ccAGTTGACGGGAAATCAGCATGCCATCGCGGTCAGTTCAGCCACCGCCGGAATGCATATCACGCTAATGGCGTTGAAAATTGGCAAGGGCGATGAAGTGATTACGCCTTCCCTGACCTGGGTTTCAACCCTCAATATGATTTCCttgtt                                                                                                                                           >  2:3466579/1‑150 (MQ=255)
            aGTTGACGGGAAATCAGCATGCCATCGCGGTCAGTTCAGCCACCGCCGGAATGCATATCACGCTAATGGCGTTGAAAATTGGCAAGGGCGATGAAGTGATTACGCCTTCCCTGACCTGGGTTTCAACCCTCAATATTATTTCCTTGTTgg                                                                                                                                         >  1:3766779/1‑150 (MQ=255)
            aGTTGACGGGAAATCAGCATGCCATCGCGGTCAGTTCAGCCACCGCCGGAATGCATATCACGCTAATGGCGTTGAAAATTGGCAAGGGCGATGAAGTGATTACGCCTTCCCTGACCTGGGTTTCAACCCTCAATATGATTTCCTTGTTgg                                                                                                                                         >  1:2992828/1‑150 (MQ=255)
               tGACGGGAAATCAGCATGCCATCGCGGTCAGTTCAGCCACCGCCGGAATGCATATCACGCTAATGGCGTTGAAAATTGGCAAGGGCGATGAAGTGATTACGCCTTCCCTGACCTGGGTTTCAACCCTCAATATGATTTCCTTGTTGGGTg                                                                                                                                      <  2:3992168/150‑1 (MQ=255)
                gACGGGAAATCAGCATGCCATCGCGGTCAGTTCAGCCACCGCCGGAATGCATATCACGCTAATGGCGTTGAAAATTGGCAAGGGCGATGAAGTGATTACGCCTTCCCTGACCTGGGTTTCAACCCTCAATATTATTTCCTTGTTGGGTGc                                                                                                                                     >  2:2278349/1‑150 (MQ=255)
                 aCGGGAAATCAGCATGCCATCGCGGTCAGTTCAGCCACCGCCGGAATGCATATCACGCTAATGGCGTTGAAAATTGGCAAGGGCGATGAAGTGATTACGCCTTCCCTGACCTGGGTTTCAACCCTCAATATTATTTCCTTGTTGGGTGCa                                                                                                                                    <  1:4195192/150‑1 (MQ=255)
                 aCGGGAAATCAGCATGCCATCGCGGTCAGTTCAGCCACCGCCGGAATGCATATCACGCTAATGGCGTTGAAAATTGGCAAGGGCGATGAAGTGATTACGCCTTCCCTGACCTGGGTTTCAACCCTCAATATGATTTCCTTGTTGGGTGCa                                                                                                                                    <  2:654678/150‑1 (MQ=255)
                   gggAAATCAGCATGCCATCGCGGTCAGTTCAGCCACCGCCGGAATGCATATCACGCTAATGGCGTTGAAAATTGGCAAGGGCGATGAAGTGATTACGCCTTCCCTGACCTGGGTTTCAACCCTCAATATTATTTCCTTGTTGGGTGCAAc                                                                                                                                  >  1:3708145/1‑150 (MQ=255)
                    ggAAATCAGCATGCCATCGCGGTCAGTTCAGCCACCGCCGGAATGCATATCACGCTAATGGCGTTGAAAATTGGCAAGGGCGATGAAGTGATTACGCCTTCCCTGACCTGGGTTTCAACCCTCAATATTATTTCCTTGTTGGGTGCAACg                                                                                                                                 <  1:2151828/150‑1 (MQ=255)
                    ggAAATCAGCATGCCATCGCGGTCAGTTCAGCCACCGCCGGAATGCATATCACGCTAATGGCGTTGAAAATTGGCAAGGGCGATGAAGTGATTACGCCTTCCCTGACCTGGGTTTCAACCCTCAATATGATTTCCTTGTTGGGTGCAACg                                                                                                                                 <  1:1196439/150‑1 (MQ=255)
                       aaTCAGCATGCCATCGCGGTCAGTTCAGCCACCGCCGGAATGCATATCACGCTAATGGCGTTGAAAATTGGCAAGGGCGATGAAGTGATTACGCCTTCCCTGACCTGGGTTTCAACCCTCAATATTATTTCCTTGTTGGGTGCAACGCCg                                                                                                                              <  1:4726474/150‑1 (MQ=255)
                          cAGGATGCCATCGCGGTCTGTTCAGCCACCGCCGGTATGCTTATCACGCTAATGGCGTTGAAATTTGGCATGGGCGATGATGTGATTACGCCTGCCCTGACCTGGGTTTCAACCCTCAATCTGATTTCCTTGTTGGGTGCAACGCCGGtt                                                                                                                           <  2:2684818/150‑2 (MQ=255)
                           aGCATGCCATCGCGGTCAGTTCAGCCACCGCCGGAATGCATATCACGCTAATGGCGTTGAAAATTGGCAAGGGCGATGAAGTGATTACGCCTTCCCTGACCTGGGTTTCAACCCTCAATATGATTTCCTTGTTGGGTGCAACGCCGGTaa                                                                                                                          >  1:2368358/1‑150 (MQ=255)
                            gCATGCCATCGCGGTCAGTTCAGCCACCGCCGGAATGCATATCACGCTAATGGCGTTGAAAATTGGCAAGGGCGATGAAGTGATTACGCCTTCCCTGACCTGGGTTTCAACCCTCAATATTATTTCCTTGTTGGGTGCAACGCCGGTAAt                                                                                                                         >  2:1470210/1‑150 (MQ=255)
                              aTGCCATCGCGGTCAGTTCAGCCACCGCCGGAATGCATATCACGCTAATGGCGTTGAAAATTGGCAAGGGCGATGAAGTGATTACGCCTTCCCTGACCTGGGTTTCAACCCTCAATATTATTTCCTTGTTGGGTGCAACGCCGGTAAtgg                                                                                                                       >  2:4319178/1‑150 (MQ=255)
                               tGCCATCGCGGTCAGTTCAGCCACCGCCGGAATGCATATCACGCTAATGGCGTTGAAAATTGGCAAGGGCGATGAAGTGATTACGCCTTCCCTGACCTGGGTTTCAACCCTCAATATTATTTCCTTGTTGGGTGCAACGCCGGTAAtggt                                                                                                                      <  2:3241709/150‑1 (MQ=255)
                                 ccATCGCGGTCAGTTCAGCCACCGCCGGAATGCATATCACGCTAATGGCGTTGAAAATTGGCAAGGGCGATGAAGTGATTACGCCTTCCCTGACCTGGGTTTCAACCCTCAATATTATTTCCTTGTTGGGTGCAACGCCGGTAAtggtgg                                                                                                                    <  2:829625/150‑1 (MQ=255)
                                 ccATCGCGGTCAGTTCAGCCACCGCCGGAATGCATATCACGCTAATGGCGTTGAAAATTGGCAAGGGCGATGAAGTGATTACGCCTTCACTGACCTGGGTTTCAACCCTCAATATGATTTCCTTGTTGGGTGCAACGCCGGTCAtggtgg                                                                                                                    >  2:2982410/1‑150 (MQ=255)
                                  cATCGCGGTCAGTTCAGCCACCGCCGGAATGCATATCACGCTAATGGCGTTGAAAATTGGCAAGGGCGATGAAGTGATTACGCCTTCCCTGACCTGGGTTTCAACCCTCAATATTATTTCCTTGTTGGGTGCAACGCCGGTAATGGTGGa                                                                                                                   >  1:4875034/1‑150 (MQ=255)
                                  cATCGCGGTCAGTTCAGCCACCGCCGGAATGCATATCACGCTAATGGCGTTGAAAATTGGCAAGGGCGATGAAGTGATTACGCCTTCCCTGACCTGGGTTTCAACCCTCAATATGATTTCCTTGTTGGGTGCAACGCCGGTAATGGTGGa                                                                                                                   >  2:1735916/1‑150 (MQ=255)
                                     cgcgGTCAGTTCAGCCACCGCCGGAATGCATATCACGCTAATGGCGTTGAAAATTGGCAAGGGCGATGAAGTGATTACGCCTTCCCTGACCTGGGTTTCAACCCTCAATATGATTTCCTTGTTGGGTGCAACGCCGGTAATGGTGGATGt                                                                                                                <  1:3616153/150‑1 (MQ=255)
                                      gcgGTCAGTTCAGCCACCGCCGGAATGCATATCACGCTAATGGCGTTGAAAATTGGCAAGGGCGATGAAGTGATTACGCCTTCCCTGACCTGGGTTTCAACCCTCAATATTATTTCCTTGTTGGGTGCAACGCCGGTAATGGTGGATGTc                                                                                                               <  2:1805095/150‑1 (MQ=255)
                                       cgGTCAGTTCAGCCACCGCCGGAATGCATATCACGCTAATGGCGTTGAAAATTGGCAAGGGCGATGAAGTGATTACGCCTTCCCTGACCTGGGTTTCAACCCTCAATATTATTTCCTTGTTGGGTGCAACGCCGGTAATGGTGGATGTCg                                                                                                              >  2:1939512/1‑150 (MQ=255)
                                       cgGTCAGTTCAGCCACCGCCGGAATGCATATCACGCTAATGGCGTTGAAAATTGGCAAGGGCGATGAAGTGATTACGCCTTCCCTGACCTGGGTTTCAACCCTCAATATGATTTCCTTGTTGGGTGCAACGCCGGTAATGGTGGATGTCg                                                                                                              <  2:1559772/150‑1 (MQ=255)
                                        ggTCAGTTCAGCCACCGCCGGAATGCATATCACGCTAATGGCGTTGAAAATTGGCAAGGGCGATGAAGTGATTACGCCTTCCCTGACCTGGGTTTCAACCCTCAATATTATTTCCTTGTTGGGTGCAACGCCGGTAATGGTGGATGTCGa                                                                                                             <  2:4337028/150‑1 (MQ=255)
                                        ggTCAGTTCAGCCACCGCCGGAATGCATATCACGCTAATGGCGTTGAAAATTGGCAAGGGCGATGAAGTGATTACGCCTTCCCTGACCTGGGTTTCAACCCTCAATATGATTTCCTTGTTGGGTGCAACGCCGGTAATGGTGGATGTCGa                                                                                                             >  1:2211598/1‑150 (MQ=255)
                                        ggTCAGTTCAGCCACCGCCGGAATGCATATCACGCTAATGGCGTTGAAAATTGGCAAGGGCGATGAAGTGATTACGCCTTCCCTGACCTGGGTTTCAACCCTCAATATGATTTCCTTGTTGGGTGCAACGCCGGTAATGGTGGATGTCGa                                                                                                             >  1:4663676/1‑150 (MQ=255)
                                            aGTTCAGCCACCGCCGGAATGCATATCACGCTAATGGCGTTGAAAATTGGCAAGGGCGATGAAGTGATTACGCCTTCCCTGACCTGGGTTTCAACCCTCAATATTATTTCCTTGTTGGGTGCAACGCCGGTAATGGTGGATGTCGACcgc                                                                                                         >  2:1665684/1‑150 (MQ=255)
                                            aGTTCAGCCACCGCCGGAATGCATATCACGCTAATGGCGTTGAAAATTGGCAAGGGCGATGAAGTGATTACGCCTTCCCTGACCTGGGTTTCAACCCTCAATATGATTTCCTTGTTGGGTGCAACGCCGGTAATGGTGGATGTCGACcgc                                                                                                         >  1:4024854/1‑150 (MQ=255)
                                              ttCAGCCACCGCCGGAATGCATATCACGCTAATGGCGTTGAAAATTGGCAAGGGCGATGAAGTGATTACGCCTTCCCTGACCTGGGTTTCAACCCTCAATATTATTTCCTTGTTGGGTGCAACGCCGGTAATGGTGGATGTCGACCGCGa                                                                                                       <  2:2183796/150‑1 (MQ=255)
                                               tCAGCCACCGCCGGAATGCATATCACGCTAATGGCGTTGAAAATTGGCAAGGGCGATGAAGTGATTACGCCTTCCCTGACCTGGGTTTCAACCCTCAATATTATTTCCTTGTTGGGTGCAACGCCGGTAATGGTGGATGTCGACCGCGAt                                                                                                      >  1:3316371/1‑150 (MQ=255)
                                                cAGCCACCGCCGGAATGCATATCACGCTAATGGCGTTGAAAATTGGCAAGGGCGATGAAGTGATTACGCCTTCCCTGACCTGGGTTTCAACCCTCAATATGATTTCCTTGTTGGGTGCAACGCCGGTAATGGTGGATGTCGACCGCGATa                                                                                                     >  2:4671152/1‑150 (MQ=255)
                                                cAGCCACCGCCGGAATGCATATCACGCTAATGGCGTTGAAAATTGGCAAGGGCGATGAAGTGATTACGCCTTCCCTGACCTGGGTTTCAACCCTCAATATGATTTCCTTGTTGGGTGCAACGCCGGTAATGGTGGATGTCGACCGCGATa                                                                                                     >  2:641807/1‑150 (MQ=255)
                                                 aGCCACCGCCGGAATGCATATCACGCTAATGGCGTTGAAAATTGGCAAGGGCGATGAAGTGATTACGCCTTCCCTGACCTGGGTTTCAACCCTCAATATTATTTCCTTGTTGGGTGCAACGCCGGTAATGGTGGATGTCGACCGCGATAc                                                                                                    >  2:1386435/1‑150 (MQ=255)
                                                 aGCCACCGCCGGAATGCATATCACGCTAATGGCGTTGAAAATTGGCAAGGGCGATGAAGTGATTACGCCTTCCCTGACCTGGGTTTCAACCCTCAATATGATTTCCTTGTTGGGTGCAACGCCGGTAATGGTGGATGTCGACCGCGATAc                                                                                                    <  1:1069776/150‑1 (MQ=255)
                                                  gCCACCGCCGGAATGCATATCACGCTAATGGCGTTGAAAATTGGCAAGGGCGATGAAGTGATTACGCCTTCCCTGACCTGGGTTTCAACCCTCAATATTATTTCCTTGTTGGGTGCAACGCCGGTAATGGTGGATGTCGACCGCGATACg                                                                                                   <  1:300390/150‑1 (MQ=255)
                                                  gCCACCGCCGGAATGCATATCACGCTAATGGCGTTGAAAATTGGCAAGGGCGATGAAGTGATTACGCCTTCCCTGACCTGGGTTTCAACCCTCAATATGATTTCCTTGTTGGGTGCAACGCCGGTAATGGTGGATGTCGACCGCGATACg                                                                                                   >  1:1946885/1‑150 (MQ=255)
                                                   ccACCGCCGGAATGCATATCACGCTAAAGGCGTTGAAAATTGGCAAGGGCGATGAAGTGATTACGCCTTCCCTGACCTGGGTTTCAACCCTCAATATTATTTCCTTGTTGGGTGCAAAGCCGGTAATGGTGGATGTCGACCGCGATACGc                                                                                                  >  2:2851259/1‑150 (MQ=255)
                                                        gccgGAATGCATATCACGCTAATGGCGTTGAAAATTGGCAAGGGCGATGAAGTGATTACGCCTTCCCTGACCTGGGTTTCAACCCTCAATATTATTTCCTTGTTGGGTGCAACGCCGGTAATGGTGGATGTCGACCGCGATACGCTGATg                                                                                             >  1:3232291/1‑150 (MQ=255)
                                                         ccgGAATGCATATCACGCTAATGGCGTTGAAAATTGGCAAGGGCGATGAAGTGATTACGCCTTCCCTGACCTGGGTTTCAACCCTCAATATGATTTCCTTGTTGGGTGCAACGCCGGTAATGGTGGATGTCGACCGCGATACGCTGATgg                                                                                            <  2:3917665/150‑1 (MQ=255)
                                                          cgGAATGCATATCACGCTAATGGCGTTGAAAATTGGCAAGGGCGATGAAGTGATTACGCCTTCCCTGACCTGGGTTTCAACCCTCAATATGATTTCCTTGTTGGGTGCAACGCCGGTAATGGTGGATGTCGACCGCGATACGCTGATGGt                                                                                           <  2:609494/150‑1 (MQ=255)
                                                           ggAATGCATATCACGCTAATGGCGTTGAAAATTGGCAAGGGCGATGAAGTGATTACGCCTTCCCTGACCTGGGTTTCAACCCTCAATATTATTTCCTTGTTGGGTGCAACGCCGGTAATGGTGGATGTCGACCGCGATACGCTGATGGTc                                                                                          >  1:4814614/1‑150 (MQ=255)
                                                           ggAATGCATATCACGCTAATGGCGTTGAAAATTGGCAAGGGCGATGAAGTGATTACGCCTTCCCTGACCTGGGTTTCAACCCTCAATATGATTTCCTTGTTGGGTGCAACGCCGGTAATGGTGGATGTCGACCGCGATACGCTGATGGTc                                                                                          >  2:1089193/1‑150 (MQ=255)
                                                                           tAATGGCGTTGAAAATTGGCAAGGGCGATGAAGTGATTACGCCTTCCCTGACCTGGGTTTCAACCCTCAATATTATTTCCTTGTTGGGTGCAACGCCGGTAATGGTGGATGTCGACCGCGATACGCTGATGGTCACGCCTGAAGCTATCg                                                                          >  2:3357728/1‑150 (MQ=255)
                                                                             aTGGCGTTGAAAATTGGCAAGGGCGATGAAGTGATTACGCCTTCCCTGACCTGGGTTTCAACCCTCAATATTATTTCCTTGTTGGGTGCAACGCCGGTAATGGTGGATGTCGACCGCGATACGCTGATGGTCACGCCTGAAGCTATCGAg                                                                        <  1:113959/150‑1 (MQ=255)
                                                                              tGGCGTTGAAAATTGGCAAGGGCGATGAAGTGATTACGCCTTCCCTGACCTGGGTTTCAACCCTCAATATTATTTCCTTGTTGGGTGCAACGCCGGTAATGGTGGATGTCGACCGCGATACGCTGATGGTCACGCCTGAAGCTATCGAGt                                                                       <  1:2550242/150‑1 (MQ=255)
                                                                               ggCGTTGAAAATTGGCAAGGGCGATGAAGTGATTACGCCTTCCCTGACCTGGGTTTCAACCCTCAATATTATTTCCTTGTTGGGTGCAACGCCGGTAATGGTGGATGTCGACCGCGATACGCTGATGGTCACGCCTGAAGCTATCGAGTc                                                                      >  2:4750248/1‑150 (MQ=255)
                                                                                  gTTGAAAATTGGCAAGGGCGATGAAGTGATTACGCCTTCCCTGACCTGGGTTTCAACCCTCAATATGATTTCCTTGTTGGGTGCAACGCCGGTAATGGTGGATGTCGACCGCGATACGCTGATGGTCACGCCTGAAGCTATCGAGTCAGc                                                                   >  1:147857/1‑150 (MQ=255)
                                                                                  gTTGAAAATTGGCAAGGGCGATGAAGTGATTACGCCTTCCCTGACCTGGGTTTCAACCCTCAATATGATTTCCTTGTTGGGTGCAACGCCGGTAATGGTGGATGTCGACCGCGATACGCTGATGGTCACGCCTGAAGCTATCGAGTCAGc                                                                   <  1:4667812/150‑1 (MQ=255)
                                                                                     gAAAATTGGCAAGGGCGATGAAGTGATTACGCCTTCCCTGACCTGGGTTTCAACCCTCAATATTATTTCCTTGTTGGGTGCAACGCCGGTAATGGTGGATGTCGACCGCGATACGCTGATGGTCACGCCTGAAGCTATCGAGTCAGCCAt                                                                <  1:2906906/150‑1 (MQ=255)
                                                                                      aaaaTTGGCAAGGGCGATGAAGTGATTACGCCTTCCCTGACCTGGGTTTCAACCCTCAATATTATTTCCTTGTTGGGTGCAACGCCGGTAATGGTGGATGTCGACCGCGATACGCTGATGGTCACGCCTGAAGCTATCGAGTCAGCCAtt                                                               <  1:1033438/150‑1 (MQ=255)
                                                                                           tGGCAAGGGCGATGAAGTGATTACGCCTTCCCTGACCTGGGTTTCAACCCTCAATATTATTTCCTTGTTGGGTGCAACGCCGGTAATGGTGGATGTCGACCGCGATACGCTGATGGTCACGCCTGAAGCTATCGAGTCAGCCATTACGcc                                                          <  1:3118913/150‑1 (MQ=255)
                                                                                             gCAAGGGCGATGAAGTGATTACGCCTTCCCTGACCTGGGTTTCAACCCTCAATATGATTTCCTTGTTGGGTGCAACGCCGGTAATGGTGGATGTCGACCGCGATACGCTGATGGTCACGCCTGAAGCTATCGAGTCAGCCATTACGCCAc                                                        <  1:1832688/150‑1 (MQ=255)
                                                                                              cAAGGGCGATGAAGTGATTACGCCTTCCCTGACCTGGGTTTCAACCCTCAATATTATTTCCTTGTTGGGTGCAACGCCGGTAATGGTGGATGTCGACCGCGATACGCTGATGGTCACGCCTGAAGCTATCGAGTCAGCCATTACGCCACg                                                       >  2:2565611/1‑150 (MQ=255)
                                                                                              cAAGGGCGATGAAGTGATTACGCCTTCCCTGACCTGGGTTTCAACCCTCAATATGATTTCCTTGTTGGGTGCAACGCCGGTAATGGTGGATGTCGACCGCGATACGCTGATGGTCACGCCTGAAGCTATCGAGTCAGCCATTACGCCACg                                                       <  2:1448615/150‑1 (MQ=255)
                                                                                              cAAGGGCGATGAAGTGATTACGCCTTCCCTGACCTGGGTTTCAACCCTCAATATGATTTCCTTGTTGGGTGCAACGCCGGTAATGGTGGATGTCGACCGCGATACGCTGATGGTCACGCCTGAAGCTATCGAGTCAGCCATTACGCCACg                                                       >  1:2556902/1‑150 (MQ=255)
                                                                                                    cGATGAAGTGATTACGCCTTCCCTGACCTGGGTTTCAACCCTCAATATGATTTCCTTGTTGGGTGCAACGCCGGTAATGGTGGATGTCGACCGCGATACGCTGATGGTCACGCCTGAAGCTATCGAGTCAGCCATTACGCCACGCACTaa                                                 <  2:3872192/150‑1 (MQ=255)
                                                                                                      aTGAAGTGATTACGCCTTCCCTGACCTGGGTTTCAACCCTCAATATTATTTCCTTGTTGGGTGCAACGCCGGTAATGGTGGATGTCGACCGCGATACGCTGATGGTCACGCCTGAAGCTATCGAGTCAGCCATTACGCCACGCACTAAAg                                               <  1:571219/150‑1 (MQ=255)
                                                                                                        gAAGTGATTACGCCTTCCCTGACCTGGGTTTCAACCCTCAATATTATTTCCTTGTTGGGTGCAACGCCGGTAATGGTGGATGTCGACCGCGATACGCTGATGGTCACGCCTGAAGCTATCGAGTCAGCCATTACGCCACGCACTAAAGcc                                             >  1:746841/1‑150 (MQ=255)
                                                                                                          aGTGATTACGCCTTCCCTGACCTGGGTTTCAACCCTCAATATTATTTCCTTGTTGGGTGCAACGCCGGTAATGGTGGATGTCGACCGCGATACGCTGATGGTCACGCCTGAAGCTATCGAGTCAGCCATTACGCCACGCACTAAAGCcat                                           >  1:1146971/1‑150 (MQ=255)
                                                                                                          aGTGATTACGCCTTCCCTGACCTGGGTTTCAACCCTCAATATGATTTCCTTGTTGGGTGCAACGCCGGTAATGGTGGATGTCGACCGCGATACGCTGATGGTCACGCCTGAAGCTATCGAGTCAGCCATTACGCCACGCACTAAAGCcat                                           <  2:4816308/150‑1 (MQ=255)
                                                                                                           gTGATTACGCCTTCCCTGACCTGGGTTTCAACCCTCAATATGATTTCCTTGTTGGGTGCAACGCCGGTAATGGTGGATGTCGACCGCGATACGCTGATGGTCACGCCTGAAGCTATCGAGTCAGCCATTACGCCACGCACTAAAGCcatc                                          <  2:1354210/150‑1 (MQ=255)
                                                                                                             gATTACGCCTTCCCTGACCTGGGTTTCAACCCTCAATATGATTTCCTTGTTGGGTGCAACGCCGGTAATGGTGGATGTCGACCGCGATACGCTGATGGTCACGCCTGAAGCTATCGAGTCAGCCATTACGCCACGCACTAAAGCcatcat                                        <  1:2073330/150‑1 (MQ=255)
                                                                                                              aTTACGCCTTCCCTGACCTGGGTTTCAACCCTCAATATTATTTCCTTGTTGGGTGCAACGCCGGTAATGGTGGATGTCGACCGCGATACGCTGATGGTCACGCCTGAAGCTATCGAGTCAGCCATTACGCCACGCACTAAAGCCATCAtt                                       <  2:314936/150‑1 (MQ=255)
                                                                                                                   gCCTTCCCTGACCTGGGTTTCAACCCTCAATATTATTTCCTTGTTGGGTGCAACGCCGGTAATGGTGGATGTCGACCGCGATACGCTGATGGTCACGCCTGAAGCTATCGAGTCAGCCATTACGCCACGCACTAAAGCCATCATTCCGGt                                  <  1:209037/150‑1 (MQ=255)
                                                                                                                   gCCTTCCCTGACCTGGGTTTCAACCCTCAATATTATTTCCTTGTTGGGTGCAACGCCGGTAATGGTGGATGTCGACCGCGATACGCTGATGGTCACGCCTGAAGCTATCGAGTCAGCCATTACGCCACGCACTAAAGCCATCATTCCGGt                                  <  2:3180499/150‑1 (MQ=255)
                                                                                                                   gCCTTCCCTGACCTGGGTTTCAACCCTCAATATGATTTCCTTGTTGGGTGCAACGCCGGTAATGGTGGATGTCGACCGCGATACGCTGATGGTCACGCCTGAAGCTATCGAGTCAGCCATTACGCCACGCACTAAAGCCATCATTCCGGt                                  >  1:3859742/1‑150 (MQ=255)
                                                                                                                     cTTCCCTGACCTGGGTTTCAACCCTCAATATGATTTCCTTGTTGGGTGCAACGCCGGTAATGGTGGATGTCGACCGCGATACGCTGATGGTCACGCCTGAAGCTATCGAGTCAGCCATTACGCCACGCACTAAAGCCATCATTCCGGTGc                                <  2:1986739/150‑1 (MQ=255)
                                                                                                                      ttCCCTGACCTGGGTTTCAACCCTCAATATTATTTCCTTGTTGGGTGCAACGCCGGTAATGGTGGATGTCGACCGCGATACGCTGATGGTCACGCCTGAAGCTATCGAGTCAGCCATTACGCCACGCACTAAAGCCATCATTCCGGTGCa                               >  2:1998005/1‑150 (MQ=255)
                                                                                                                       tCCCTGACCTGGGTTTCAACCCTCAATATTATTTCCTTGTTGGGTGCAACGCCGGTAATGGTGGATGTCGACCGCGATACGCTGATGGTCACGCCTGAAGCTATCGAGTCAGCCATTACGCCACGCACTAAAGCCATCATTCCGGTGCAt                              >  2:4362341/1‑150 (MQ=255)
                                                                                                                        cccTGACCTGGGTTTCAACCCTCAATATGATTTCCTTGTTGGGTGCAACGCCGGTAATGGTGGATGTCGACCGCGATACGCTGATGGTCACGCCTGAAGCTATCGAGTCAGCCATTACGCCACGCACTAAAGCCATCATTCCGGTGCAtt                             <  2:1291353/150‑1 (MQ=255)
                                                                                                                        cccTGACCTGGGTTTCAACCCTCAATATGATTTCCTTGTTGGGTGCAACGCCGGTAATGGTGGATGTCGACCGCGATACGCTGATGGTCACGCCTGAAGCTATCGAGTCAGCCATTACGCCACGCACTAAAGCCATCATTCCGGTGCAtt                             <  2:3753101/150‑1 (MQ=255)
                                                                                                                          cTGACCTGGGTTTCAACCCTCAATATTATTTCCTTGTTGGGTGCAACGCCGGTAATGGTGGATGTCGACCGCGATACGCTGATGGTCACGCCTGAAGCTATCGAGTCAGCCATTACGCCACGCACTAAAGCCATCATTCCGGTGCATTAt                           >  1:4466215/1‑150 (MQ=255)
                                                                                                                           tGACCTGGGTTTCAACCCTCAATATGATTTCCTTGTTGGGTGCAACGCCGGTAATGGTGGATGTCGACCGCGATACGCTGATGGTCACGCCTGAAGCTATCGAGTCAGCCATTACGCCACGCACTAAAGCCATCATTCCGGTGCATTATg                          >  1:797020/1‑150 (MQ=255)
                                                                                                                             aCCTGGGTTTCAACCCTCAATATTATTTCCTTGTTGGGTGCAACGCCGGTAATGGTGGATGTCGACCGCGATACGCTGATGGTCACGCCTGAAGCTATCGAGTCAGCCATTACGCCACGCACTAAAGCCATCATTCCGGTGCATTATGcc                        >  1:2005587/1‑150 (MQ=255)
                                                                                                                             aCCTGGGTTTCAACCCTCAATATTATTTCCTTGTTGGGTGCAACGCCGGTAATGGTGGATGTCGACCGCGATACGCTGATGGTCACGCCTGAAGCTATCGAGTCAGCCATTACGCCACGCACTAAAGCCATCATTCCGGTGCATTATGcc                        >  2:1144674/1‑150 (MQ=255)
                                                                                                                             aCCTGGGTTTCAACCCTCAATATTATTTCCTTGTTGGGTGCAACGCCGGTAATGGTGGATGTCGACCGCGATACGCTGATGGTCACGCCTGAAGCTATCGAGTCAGCCATTACGCCACGCACTAAAGCCATCATTCCGGTGCATTATGcc                        <  2:491802/150‑1 (MQ=255)
                                                                                                                              ccTGGGTTTCAACCCTCAATATGATTTCCTTGTTGGGTGCAACGCCGGTAATGGTGGATGTCGACCGCGATACGCTGATGGTCACGCCTGAAGCTATCGAGTCAGCCATTACGCCACGCACTAAAGCCATCATTCCGGTGCATTATGCCg                       <  2:2407443/150‑1 (MQ=255)
                                                                                                                               cTGGGTTTCAACCCTCAATATGATTTCCTTGTTGGGTGCAACGCCGGTAATGGTGGATGTCGACCGCGATACGCTGATGGTCACGCCTGAAGCTATCGAGTCAGCCATTACGCCACGCACTAAAGCCATCATTCCGGTGCATTATGCCgg                      >  1:3430394/1‑150 (MQ=255)
                                                                                                                                tGGGTTTCAACCCTCAATATGATTTCCTTGTTGGGTGCAACGCCGGTAATGGTGGATGTCGACCGCGATACGCTGATGGTCACGCCTGAAGCTATCGAGTCAGCCATTACGCCACGCACTAAAGCCATCATTCCGGTGCATTATGCCGGt                     >  2:323560/1‑150 (MQ=255)
                                                                                                                                   gTTTCAACCCTCAATATGATTTCCTTGTTGGGTGCAACGCCGGTAATGGTGGATGTCGACCGCGATACGCTGATGGTCACGCCTGAAGCTATCGAGTCAGCCATTACGCCACGCACTAAAGCCATCATTCCGGTGCATTATGCCGGTgcg                  <  2:1904255/150‑1 (MQ=255)
                                                                                                                                     ttCAACCCTCAATATTATTTCCTTGTTGGGTGCAACGCCGGTAATGGTGGATGTCGACCGCGATACGCTGATGGTCACGCCTGAAGCTATCGAGTCAGCCATTACGCCACGCACTAAAGCCATCATTCCGGTGCATTATGCCGGTGCGcc                <  1:1429040/150‑1 (MQ=255)
                                                                                                                                     ttCAACCCTCAATATGATTTCCTTGTTGGGTGCAACGCCGGTAATGGTGGATGTCGACCGCGATACGCTGATGGTCACGCCTGAAGCTATCGAGTCAGCCATTACGCCACGCACTAAAGCCATCATTCCGGTGCATTATGCCGGTGCGcc                >  1:3341363/1‑150 (MQ=255)
                                                                                                                                       cAACCCTCAATATTATTTCCTTGTTGGGTGCAACGCCGGTAATGGTGGATGTCGACCGCGATACGCTGATGGTCACGCCTGAAGCTATCGAGTCAGCCATTACGCCACGCACTAAAGCCATCATTCCGGTGCATTATGCCGGTGCGCcag              >  1:3170618/1‑150 (MQ=255)
                                                                                                                                        aaCCCTCAATATTATTTCCTTGTTGGGTGCAACGCCGGTAATGGTGGATGTCGACCGCGATACGCTGATGGTCACGCCTGAAGCTATCGAGTCAGCCATTACGCCACGCACTAAAGCCATCATTCCGGTGCATTATGCCGGTGCGCcagc             >  2:402810/1‑150 (MQ=255)
                                                                                                                                         aCCCTCAATATGATTTCCTTGTTGGGTGCAACGCCGGTAATGGTGGATGTCGACCGCGATACGCTGATGGTCACGCCTGAAGCTATCGAGTCAGCCATTACGCCACGCACTAAAGCCATCATTCCGGTGCATTATGCCGGTGCGCcagca            <  2:565514/150‑1 (MQ=255)
                                                                                                                                           ccTCAATATTATTTCCTTGTTGGGTGCAACGCCGGTAATGGTGGATGTCGACCGCGATACGCTGATGGTCACGCCTGAAGCTATCGAGTCAGCCATTACGCCACGCACTAAAGCCATCATTCCGGTGCATTATGCCGGTGCGCCAGCAGa          <  2:219896/150‑1 (MQ=255)
                                                                                                                                               aaTATTATTTCCTTGTTGGGTGCAACGCCGGTAATGGTGGATGTCGACCGCGATACGCTGATGGTCACGCCTGAAGCTATCGAGTCAGCCATTACGCCACGCACTAAAGCCATCATTCCGGTGCATTATGCCGGTGCGCCAGCAGATAtt      >  2:4196937/1‑150 (MQ=255)
                                                                                                                                               aaTATGATTTCCTTGTTGGGTGCAACGCCGGTAATGGTGGATGTCGACCGCGATACGCTGATGGTCACGCCTGAAGCTATCGAGTCAGCCATTACGCCACGCACTAAAGCCATCATTCCGGTGCATTATGCCGGTGCGCCAGCAGATAtt      <  2:2994275/150‑1 (MQ=255)
                                                                                                                                                atatTATTTCCTTGTTGGGTGCAACGCCGGTAATGGTGGATGTCGACCGCGATACGCTGATGGTCACGCCTGAAGCTATCGAGTCAGCCATTACGCCACGCACTAAAGCCATCATTCCGGTGCATTATGCCGGTGCGCCAGCAGATATTg     >  1:1706712/1‑150 (MQ=255)
                                                                                                                                                   tGATTTCCTTGTTGGGTGCAACGCCGGTAATGGTGGATGTCGACCGCGATACGCTGATGGTCACGCCTGAAGCTATCGAGTCAGCCATTACGCCACGCACTAAAGCCATCATTCCGGTGCATTATGCCGGTGCGCCAGCAGATATTGACg  >  2:2065670/1‑150 (MQ=255)
                                                                                                                                                   tGATTTCCTTGTTGGGTGCAACGCCGGTAATGGTGGATGTCGACCGCGATACGCTGATGGTCACGCCTGAAGCTATCGAGTCAGCCATTACGCCACGCACTAAAGCCATCATTCCGGTGCATTATGCCGGTGCGCCAGCAGATATTGACg  >  1:296410/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
AAGCTTTTTGCCAGTTGACGGGAAATCAGCATGCCATCGCGGTCAGTTCAGCCACCGCCGGAATGCATATCACGCTAATGGCGTTGAAAATTGGCAAGGGCGATGAAGTGATTACGCCTTCCCTGACCTGGGTTTCAACCCTCAATATGATTTCCTTGTTGGGTGCAACGCCGGTAATGGTGGATGTCGACCGCGATACGCTGATGGTCACGCCTGAAGCTATCGAGTCAGCCATTACGCCACGCACTAAAGCCATCATTCCGGTGCATTATGCCGGTGCGCCAGCAGATATTGACG  >  CP009273/2359525‑2359821

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: