Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,000,841 G→T 6.3% A206A (GCG→GCT pyrD → dihydro‑orotate oxidase, FMN‑linked

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,000,8410GT 6.3% 1083.0 / 31.6 283A206A (GCG→GCTpyrDdihydro‑orotate oxidase, FMN‑linked
Reads supporting (aligned to +/- strand):  ref base G (144/121);  new base T (9/9);  total (153/130)
Fisher's exact test for biased strand distribution p-value = 8.09e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.33e-01

Only 100 of 294 total aligned reads displayed.

ATTTGTATGGAAAAAATCTATGCCTATGCGGGATATATCGCCATCAATATTTCATCGCCGAATACCCCAGGATTACGCACGCTGCAATATGGTGAAGCGCTGGATGATCTCTTAACCGCGATTAAAAATAAGCAAAATGATTTGCAAGCGATGCACCATAAATATGTGCCGATCGCAGTGAAGATCGCGCCGGATCTTTCTGAAGAAGAATTGATCCAGGTTGCCGATAGTTTAGTTCGCCATAATATTGATGGCGTTATTGCAACCAATACCACACTCGATCGTTCTCTTGTTCAGGG  >  CP009273/1000692‑1000990
                                                                                                                                                     |                                                                                                                                                     
        ggAAAAAATCTATGCCTATGCGGGATATATCGCCATCAATATTTCATCGCCGAATACCCCAGGATTACGCACGCTGCAATATGGTGAAGCGCTGGATGATCTCTTAACCGCGATTAAAAATAAGCAAAATGATTTGCAAGCGATGCACCa                                                                                                                                               >  1:1790846/1‑150 (MQ=255)
        ggAAAAAATCTATGCCTATGCGGGATATATCGCCATCAATATTTCATCGCCGAATACCCCAGGATTACGCACGCTGCAATATGGTGAAGCGCTGGATGATCTCTTAACCGCGATTAAAAATAAGCAAAATGATTTGCAAGCGATGCACCa                                                                                                                                               >  1:3662449/1‑150 (MQ=255)
            aaaaTCTATGCCTATGCGGGATATATCGCCATCAATATTTCATCGCCGAATACCCCAGGATTACGCACGCTGCAATATGGTGAAGCGCTGGATGATCTCTTAACCGCGATTAAAAATAAGCAAAATGATTTGCAAGCGATGCACCATaaa                                                                                                                                           >  2:2183025/1‑150 (MQ=255)
              aaTCTATGCCTATGCGGGATATATCGCCATCAATATTTCATCGCCGAATACCCCAGGATTACGCACGCTGCAATATGGTGAAGCGCTGGATGATCTCTTAACCGCGATTAAAAATAAGCAAAATGATTTGCAAGCGATGCACCATAAata                                                                                                                                         >  2:2453112/1‑150 (MQ=255)
               aTCTATGCCTATGCGGGATATATCGCCATCAATATTTCATCGCCGAATACCCCAGGATTACGCACGCTGCAATATGGTGAAGCGCTGGATGATCTCTTAACCGCGATTAAAAATAAGCAAAATGATTTGCAAGCGATGCACCATAAatat                                                                                                                                        >  1:3289430/1‑150 (MQ=255)
                tCTATGCCTATGCGGGATATATCGCCATCAATATTTCATCGCCGAATACCCCAGGATTACGCACGCTGCAATATGGTGAAGCGCTGGATGATCTCTTAACCGCGATTAAAAATAAGCAAAATGATTTGCAAGCGATGCACCATAAATAtg                                                                                                                                       >  1:3861199/1‑150 (MQ=255)
                 ctatgcctatgcGGGATATATCGACATCAATATTTCATCGCCGAATACCCCAGGATTACGCACGCTGCAATATGGTGAAGCGCTGGATGATCTCTTAACCGCGATTAAAAATAAGCAAAATGATTTGCAAGCGATGCACCATAAATAtgt                                                                                                                                      <  2:2902101/150‑1 (MQ=255)
                  tatgcctatgcGGGATATATCGCCATCAATATTTCATCGCCGAATACCCCAGGATTACGCACGCTGCAATATGGTGAAGCGCTGGATGATCTCTTAACCGCGATTAAAAATAAGCAAAATGATTTGCAAGCTATGCACCATAAATAtgtg                                                                                                                                     <  1:4476201/150‑1 (MQ=255)
                  tatgcctatgcGGGATATATCGCCATCAATATTTCATCGCCGAATACCCCAGGATTACGCACGCTGCAATATGGTGAAGCGCTGGATGATCTCTTAACCGCGATTAAAAATAAGCAAAATGATTTGCAAGCGATGCACCATAAATAtgtg                                                                                                                                     >  2:1797002/1‑150 (MQ=255)
                   atgcctatgcGGGATATATCGCCATCAATATTTCATCGCCGAATACCCCAGGATTACGCACGCTGCAATATGGTGAAGCGCTGGATGATCTCTTAACCGCGATTAAAAATAAGCAAAATGATTTGCAAGCGATGCACCATAAATATGTGc                                                                                                                                    >  2:429954/1‑150 (MQ=255)
                    tgcctatgcGGGATATATCGCCATCAATATTTCATCGCCGAATACCCCAGGATTACGCACGCTGCAATATGGTGAAGCGCTGGATGATCTCTTAACCGCGATTAAAAATAAGCAAAATGATTTGCAAGCGATGCACCATAAATATGTGcc                                                                                                                                   <  2:3018583/150‑1 (MQ=255)
                          tgcGGGATATATCGCCATCAATATTTCATCGCCGAATACCCCAGGATTACGCACGCTGCAATATGGTGAAGCGCTGGATGATCTCTTAACCGCGATTAAAAATAAGCAAAATGATTTGCAAGCGATGCACCATAAATATGTGCCGATCGc                                                                                                                             >  2:4366493/1‑150 (MQ=255)
                            cGGGATATATCGCCATCAATATTTCATCGCCGAATACCCCAGGATTACGCACGCTGCAATATGGTGAAGCGCTGGATGATCTCTTAACCGCGATTAAAAATAAGCAAAATGATTTGCAAGCGATGCACCATAAATATGTGCCGATCGCAg                                                                                                                           <  1:2051546/150‑1 (MQ=255)
                            cGGGATATATCGCCATCAATATTTCATCGCCGAATACCCCAGGATTACGCACGCTGCAATATGGTGAAGCGCTGGATGATCTCTTAACCGCGATTAAAAATAAGCAAAATGATTTGCAAGCGATGCACCATAAATATGTGCCGATCGCAg                                                                                                                           <  1:4134530/150‑1 (MQ=255)
                             gggATATATCGCCATCAATATTTCATCGCCGAATACCCCAGGATTACGCACGCTGCAATATGGTGAAGCGCTGGATGATCTCTTAACCGCGATTAAAAATAAGCAAAATGATTTGCAAGCGATGCACCATAAATATGTGCCGATCGCAGt                                                                                                                          >  2:1429893/1‑150 (MQ=255)
                                 tatatCGCCATCAATATTTCATCGCCGAATACCCCAGGATTACGCACGCTGCAATATGGTGAAGCGCTGGATGATCTCTTAACCGCGATTAAAAATAAGCAAAATGATTTGCAAGCGATGCACCATAAATATGTGCCGATCGCAGTGAAg                                                                                                                      >  2:2468109/1‑150 (MQ=255)
                                  atatCGCCATCAATATTTCATCGCCGAATACCCCAGGATTACGCACGCTGCAATATGGTGAAGCGCTGGATGATCTCTTAACCGCGATTAAAAATAAGCAAAATGATTTGCAAGCGATGCACCATAAATATGTGCCGATCGCAGTGAAGa                                                                                                                     >  1:2145666/1‑150 (MQ=255)
                                   tatCGCCATCAATATTTCATCGCCGAATACCCCAGGATTACGCACGCTGCAATATGGTGAAGCGCTGGATGATCTCTTAACCGCGATTAAAAATAAGCAAAATGATTTGCAAGCTATGCACCATAAATATGTGCCGATCGCAGTGAAGAt                                                                                                                    >  1:2889613/1‑150 (MQ=255)
                                    atCGCCATCAATATTTCATCGCCGAATACCCCAGGATTACGCACGCTGCAATATGGTGAAGCGCTGGATGATCTCTTAACCGCGATTAAAAATAAGCAAAATGATTTGCAAGCTATGCACCATAAATATGTGCCGATCGCAGTGAAGATc                                                                                                                   <  2:1754554/150‑1 (MQ=255)
                                    atCGCCATCAATATTTCATCGCCGAATACCCCAGGATTACGCACGCTGCAATATGGTGAAGCGCTGGATGATCTCTTAACCGCGATTAAAAATAAGCAAAATGATTTGCAAGCGATGCACCATAAATATGTGCCGATCGCAGTGAAGATc                                                                                                                   >  1:2828237/1‑150 (MQ=255)
                                    atCGCCATCAATATTTCATCGCCGAATACCCCAGGATTACGCACGCTGCAATATGGTGAAGCGCTGGATGATCTCTTAACCGCGATTAAAAATAAGCAAAATGATTTGCAAGCGATGCACCATAAATATGTGCCGATCGCAGTGAAGATc                                                                                                                   <  2:596429/150‑1 (MQ=255)
                                    atCGCCATCAATATTTCATCGCCGAATACCCCAGGATTACGCACGCTGCAATATGGTGAAGCGCTGGATGATCTCTTAACCGCGATTAAAAATAAGCAAAATGATTTGCAAGCGATGCACCATAAATATGTGCCGATCGCAGTGAAGATc                                                                                                                   >  2:4294881/1‑150 (MQ=255)
                                         cATCAATATTTCATCGCCGAATACCCCAGGATTACGCACGCTGCAATATGGTGAAGCGCTGGATGATCTCTTAACCGCGATTAAAAATAAGCAAAATGATTTGCAAGCGATGCACCATAAATATGTGCCGATCGCAGTGAAGATCGCGcc                                                                                                              <  1:2204019/150‑1 (MQ=255)
                                         cATCAATATTTCATCGCCGAATACCCCAGGATTACGCACGCTGCAATATGGTGAAGCGCTGGATGATCTCTTAACCGCGATTAAAAATAAGCAAAATGATTTGCAAGCGATGCACCATAAATATGTGCCGATCGCAGTGAAGATCGCGcc                                                                                                              >  1:2270396/1‑150 (MQ=255)
                                            cAATATTTCATCGCCGAATACCCCAGGATTACGCACGCTGCAATATGGTGAAGCGCTGGATGATCTCTTAACCGCGATTAAAAATAAGCAAAATGATTTGCAAGCGATGCACCATAAATATGTGCCGATCGCAGTGAAGATCGCGCCGGa                                                                                                           >  1:559900/1‑150 (MQ=255)
                                            cAATATTTCATCGCCGAATACCCCAGGATTACGCACGCTGCAATATGGTGAAGCGCTGGATGATCTCTTAACCGCGATTAAAAATAAGCAAAATGATTTGCAAGCGATGCACCATAAATATGTGCCGATCGCAGTGAAGATCGCGCCGGa                                                                                                           >  2:3825433/1‑150 (MQ=255)
                                                atTTCATCGCCGAATACCCCAGGATTACGCACGCTGCAATATGGTGAAGCGCTGGATGATCTCTTAACCGCGATTAAAAATAAGCAAAATGATTTGCAAGCGATGCACCATAAATATGTGCCGATCGCAGTGAAGATCGCGCCGGATCtt                                                                                                       <  1:94947/150‑1 (MQ=255)
                                                  ttCATCGCCGAATACCCCAGGATTACGCACGCTGCAATATGGTGAAGCGCTGGATGATCTCTTAACCGCGATTAAAAATAAGCAAAATGATTTGCAAGCGATGCACCATAAATATGTGCCGATCGCAGTGAAGATCGCGCCGGATCTTTc                                                                                                     >  1:3253398/1‑150 (MQ=255)
                                                   tCATCGCCGAATACCCCAGGATTACGCACGCTGCAATATGGTGAAGCGCTGGATGATCTCTTAACCGCGATTAAAAATAAGCAAAATGATTTGCAAGCGATGCACCATAAATATGTGCCGATCGCAGTGAAGATCGCGCCGGATCTTTCt                                                                                                    <  2:3727702/150‑1 (MQ=255)
                                                    cATCGCCGAATACCCCAGGATTACGCACGCTGCAATATGGTGAAGCGCTGGATGATCTCTTAACCGCGATTAAAAATAAGCAAAATGATTTGCAAGCGATGCACCATAAATATGTGCCGATCGCAGTGAAGATCGCGCCGGATCTTTCTg                                                                                                   >  2:3920009/1‑150 (MQ=255)
                                                       cGCCGAATACCCCAGGATTACGCACGCTGCAATATGGTGAAGCGCTGGATGATCTCTTAACCGCGATTAAAAATAAGCAAAATGATTTGCAAGCGATGCACCATAAATATGTGCCGATCGCAGTGAAGATCGCGCCGGATCTTTCTgaag                                                                                                >  1:1264901/1‑150 (MQ=255)
                                                         ccGAATACCCCAGGATTACGCACGCTGCAATATGGTGAAGCGCTGGATGATCTCTTAACCGCGATTAAAAATAAGCAAAATGATTTGCAAGCGATGCACCATAAATATGTGCCGATCGCAGTGAAGATCGCGCCGGATCTTTCTgaagaa                                                                                              >  2:1688654/1‑150 (MQ=255)
                                                         ccGAATACCCCAGGATTACGCACGCTGCAATATGGTGAAGCGCTGGATGATCTCTTAACCGCGATTAAAAATAAGCAAAATGATTTGCAAGCGATGCACCATAAATATGTGCCGATCGCAGTGAAGATCGCGCCGGATCTTTCTgaagaa                                                                                              >  2:758950/1‑150 (MQ=255)
                                                           gAATACCCCAGGATTACGCACGCTGCAATATGGTGAAGCGCTGGATGATCTCTTAACCGCGATTAAAAATAAGCAAAATGATTTGCAAGCGATGCACCATAAATATGTGCCGATCGCAGTGAAGATCGCGCCGGATCTTTCTgaagaaga                                                                                            >  2:2359684/1‑150 (MQ=255)
                                                            aaTACCCCAGGATTACGCACGCTGCAATATGGTGAAGCGCTGGATGATCTCTTAACCGCGATTAAAAATAAGCAAAATGATTTGCAAGCGATGCACCATAAATATGTGCCGATCGCAGTGAAGATCGCGCCGGATCTTTCTgaagaagaa                                                                                           >  2:3606595/1‑150 (MQ=255)
                                                             aTACCCCAGGATTACGCACGCTGCAATATGGTGAAGCGCTGGATGATCTCTTAACCGCGATTAAAAATAAGCAAAATGATTTGCAAGCGATGCACCATAAATATGTGCCGATCGCAGTGAAGATCGCGCCGGATCTTTCTGAAGAAGAAt                                                                                          >  1:349627/1‑150 (MQ=255)
                                                                 cccAGGATTACGCACGCTGCAATATGGTGAAGCGCTGGATGATCTCTTAACCGCGATTAAAAATAAGCAAAATGATTTGCAAGCGATGCACCATAAATATGTGCCGATCGCAGTGAAGATCGCGCCGGATCTTTCTGAAGAAGAATTGAt                                                                                      <  1:3781914/150‑1 (MQ=255)
                                                                 cccAGGATTACGCACGCTGCAATATGGTGAAGCGCTGGATGATCTCTTAACCGCGATTAAAAATAAGCAAAATGATTTGCAAGCGATGCACCATAAATATGTGCCGATCGCAGTGAAGATCGCGCCGGATCTTTCTGAAGAAGAATTGAt                                                                                      <  1:3048962/150‑1 (MQ=255)
                                                                  ccAGGATTACGCACGCTGCAATATGGTGAAGCGCTGGATGATCTCTTAACCGCGATTAAAAATAAGCAAAATGATTTGCAAGCGATGCACCATAAATATGTGCCGATCGCAGTGAAGATCGCGCCGGATCTTTCTGAAGAAGAATTGATc                                                                                     >  1:735727/1‑150 (MQ=255)
                                                                  ccAGGATTACGCACGCTGCAATATGGTGAAGCGCTGGATGATCTCTTAACCGCGATTAAAAATAAGCAAAATGATTTGCAAGCGATGCACCATAAATATGTGCCGATCGCAGTGAAGATCGCGCCGGATCTTTCTGAAGAAGAATTGATc                                                                                     >  1:4056215/1‑150 (MQ=255)
                                                                    aGGATTACGCACGCTGCAATATGGTGAAGCGCTGGATGATCTCTTAACCGCGATTAAAAATAAGCAAAATGATTTGCAAGCGATGCACCATAAATATGTGCCGATCGCAGTGAAGATCGCGCCGGATCTTTCTGAAGAAGAATTGATCCa                                                                                   >  2:4258581/1‑150 (MQ=255)
                                                                        ttACGCACGCTGCAATATGGTGAAGCGCTGGATGATCTCTTAACCGCGATTAAAAATAAGCAAAATGATTTGCAAGCGATGCACCATAAATATGTGCCGATCGCAGTGAAGATCGCGCCGGATCTTTCTGAAGAAGAATTGATCCAGGtt                                                                               <  2:2025441/150‑1 (MQ=255)
                                                                          acgcacgcTGCAATATGGTGAAGCGCTGGATGATCTCTTAACCGCGATTAAAAATAAGCAAAATGATTTGCAAGCGATGCACCATAAATATGTGCCGATCGCAGTGAAGATCGCGCCGGATCTTTCTGAAGAAGAATTGATCCAGGTTGc                                                                             <  1:4544812/150‑1 (MQ=255)
                                                                          acgcacgcTGCAATATGGTGAAGCGCTGGATGATCTCTTAACCGCGATTAAAAATAAGCAAAATGATTTGCAAGCGATGCACCATAAATATGTGCCGATCGCAGTGAAGATCGCGCCGGATCTTTCTGAAGAAGAATTGATCCAGGTTGc                                                                             <  1:2650469/150‑1 (MQ=255)
                                                                            gcacgcTGCAATATGGTGAAGCGCTGGATGATCTCTTAACCGCGATTAAAAATAAGCAAAATGATTTGCAAGCTATGCACCATAAATATGTGCCGATCGCAGTGAAGATCGCGCCGGATCTTTCTGAAGAAGAATTGATCCAGGTTGCCg                                                                           <  2:3671134/150‑1 (MQ=255)
                                                                            gcacgcTGCAATATGGTGAAGCGCTGGATGATCTCTTAACCGCGATTAAAAATAAGCAAAATGATTTGCAAGCGATGCACCATAAATATGTGCCGATCGCAGTGAAGATCGCGCCGGATCTTTCTGAAGAAGAATTGATCCAGGTTGCCg                                                                           <  2:3937855/150‑1 (MQ=255)
                                                                            gcacgcTGCAATATGGTGAAGCGCTGGATGATCTCTTAACCGCGATTAAAAATAAGCAAAATGATTTGCAAGCGATGCACCATAAATATGTGCCGATCGCAGTGAAGATCGCGCCGGATCTTTCTGAAGAAGAATTGATCCAGGTTGCCg                                                                           >  2:2687306/1‑150 (MQ=255)
                                                                                 cTGCAATATGGTGAAGCGCTGGATGATCTCTTAACCGCGATTAAAAATAAGCAAAATGATTTGCAAGCGATGCACCATAAATATGTGCCGATCGCAGTGAAGATCGCGCCGGATCTTTCTGAAGAAGAATTGATCCAGGTTGCCGAtagt                                                                      >  2:402203/1‑150 (MQ=255)
                                                                                  tGCAATATGGTGAAGCGCTGGATGATCTCTTAACCGCGATTAAAAATAAGCAAAATGATTTGCAAGCGATGCACCATAAATATGTGCCGATCGCAGTGAAGATCGCGCCGGATCTTTCTGAAGAAGAATTGATCCAGGTTGCCGAtagtt                                                                     >  1:3983122/1‑150 (MQ=255)
                                                                                    cAATATGGTGAAGCGCTGGATGATCTCTTAACCGCGATTAAAAATAAGCAAAATGATTTGCAAGCGATGCACCATAAATATGTGCCGATCGCAGTGAAGATCGCGCCGGATCTTTCTGAAGAAGAATTGATCCAGGTTGCCGAtagttta                                                                   >  2:1899904/1‑150 (MQ=255)
                                                                                    cAATATGGTGAAGCGCTGGATGATCTCTTAACCGCGATTAAAAATAAGCAAAATGATTTGCAAGCGATGCACCATAAATATGTGCCGATCGCAGTGAAGATCGCGCCGGATCTTTCTGAAGAAGAATTGATCCAGGTTGCCGAtagttta                                                                   <  1:973410/150‑1 (MQ=255)
                                                                                     aaTATGGTGAAGCGCTGGATGATCTCTTAACCGCGATTAAAAATAAGCAAAATGATTTGCAAGCGATGCACCATAAATATGTGCCGATCGCAGTGAAGATCGCGCCGGATCTTTCTGAAGAAGAATTGATCCAGGTTGCCGAtagtttag                                                                  <  1:1449204/150‑1 (MQ=255)
                                                                                      atatGGTGAAGCGCTGGATGATCTCTTAACCGCGATTAAAAATAAGCAAAATGATTTGCAAGCGATGCACCATAAATATGTGCCGATCGCAGTGAAGATCGCGCCGGATCTTTCTGAAGAAGAATTGATCCAGGTTGCCGAtagtttagt                                                                 <  1:2760934/150‑1 (MQ=255)
                                                                                      atatGGTGAAGCGCTGGATGATCTCTTAACCGCGATTAAAAATAAGCAAAATGATTTGCAAGCGATGCACCATAAATATGTGCCGATCGCAGTGAAGATCGCGCCGGATCTTTCTGAAGAAGAATTGATCCAGGTTGCCGAtaggttagt                                                                 >  1:3897033/1‑150 (MQ=255)
                                                                                       tatGGTGAAGCGCTGGATGATCTCTTAACCGCGATTAAAAATAAGCAAAATGATTTGCAAGCGATGCACCATAAATATGTGCCGATCGCAGTGAAGATCGCGCCGGATCTTTCTGAAGAAGAATTGATCCAGGTTGCCGAtagtttagtt                                                                <  1:4063068/150‑1 (MQ=255)
                                                                                        atGGTGAAGCGCTGGATGATCTCTTAACCGCGATTAAAAATAAGCAAAATGCTTTGCAAGCGATGCACCATAAATATGTGCCGATCGCAGTGAAGATCGCGCCGGATCTTTCTGAAGAAGAATTGATCCAGGTTGCCGATAGTTTAGTTc                                                               <  2:1722248/150‑1 (MQ=255)
                                                                                         tGGTGAAGCGCTGGATGATCTCTTAACCGCGATTAAAAATAAGCAAAATGATTTGCAAGCGATGCACCATAAATATGTGCCGATCGCAGTGAAGATCGCGCCGGATCTTTCTGAAGAAGAATTGATCCAGGTTGCCGATAGTTTAGTTCg                                                              <  1:4312217/150‑1 (MQ=255)
                                                                                          ggTGAAGCGCTGGATGATCTCTTAACCGCGATTAAAAATAAGCAAAATGATTTGCAAGCGATGCACCATAAATATGTGCCGATCGCAGTGAAGATCGCGCCGGATCTTTCTGAAGAAGAATTGATCCAGGTTGCCGATAGTTTAGTTCGc                                                             >  2:2010405/1‑150 (MQ=255)
                                                                                              aaGCGCTGGATGATCTCTTAACCGCGATTAAAAATAAGCAAAATGATTTGCAAGCGATGCACCATAAATATGTGCCGATCGCAGTGAAGATCGCGCCGGATCTTTCTGAAGAAGAATTGATCCAGGTTGCCGATAGTTTAGTTCGCCata                                                         <  1:3683978/150‑1 (MQ=255)
                                                                                               aGCGCTGGATGATCTCTTAACCGCGATTAAAAATAAGCAAAATGATTTGCAAGCGATGCACCATAAATATGTGCCGATCGCAGTGAAGATCGCGCCGGATCTTTCTGAAGAAGAATTGATCCAGGTTGCCGATAGTTTAGTTCGCCataa                                                        <  2:2140451/150‑1 (MQ=255)
                                                                                                gcgcTGGATGATCTCTTAACCGCGATTAAAAATAAGCAAAATGATTTGCAAGCGATGCACCATAAATATGTGCCGATCGCAGTGAAGATCGCGCCGGATCTTTCTGAAGAAGAATTGATCCAGGTTGCCGATAGTTTAGTTCGCCataat                                                       <  1:1962925/150‑1 (MQ=255)
                                                                                                  gcTGGATGATCTCTTAACCGCGATTAAAAATAAGCAAAATGATTTGCAAGCGATGCACCATAAATATGTGCCGATCGCAGTGAAGATCGCGCCGGATCTTTCTGAAGAAGAATTGATCCAGGTTGCCGATAGTTTAGTTCGCCATAatat                                                     <  1:510709/150‑1 (MQ=255)
                                                                                                   cTGGATGATCTCTTAACCGCGATTAAAAATAAGCAAAATGATTTGCAAGCGATGCACCATAAATATGTGCCGATCGCAGTGAAGATCGCGCCGGATCTTTCTGAAGAAGAATTGATCCAGGTTGCCGATAGTTTAGTTCGCCATAATAtt                                                    <  2:1926844/150‑1 (MQ=255)
                                                                                                   cTGGATGATCTCTTAACCGCGATTAAAAATAAGCAAAATGATTTGCAAGCGATGCACCATAAATATGTGCCGATCGCAGTGAAGATCGCGCCGGATCTTTCTGAAGAAGAATTGATCCAGGTTGCCGATAGTTTAGTTCGCCATAATAtt                                                    <  2:2141907/150‑1 (MQ=255)
                                                                                                     ggATGATCTCTTAACCGCGATTAAAAATAAGCAAAATGATTTGCAAGCGATGCACCATAAATATGTGCCGATCGCAGTGAAGATCGCGCCGGATCTTTCTGAAGAAGAATTGATCCAGGTTGCCGATAGTTTAGTTCGCCATAATATTGa                                                  >  1:2916470/1‑150 (MQ=255)
                                                                                                      gatgatCTCTTAACCGCGATTAAAAATAAGCAAAATGATTTGCAAGCTATGCACCATAAATATGTGCCGATCGCAGTGAAGATCGCGCCGGATCTTTCTGAAGAAGAATTGATCCAGGTTGCCGATAGTTTAGTTCGCCATAATATTGAt                                                 >  2:2282198/1‑150 (MQ=255)
                                                                                                      gatgatCTCTTAACCGCGATTAAAAATAAGCAAAATGATTTGCAAGCGATGCACCATAAATATGTGCCGATCGCAGTGAAGATCGCGCCGGATCTTTCTGAAGAAGAATTGATCCAGGTTGCCGATAGTTTAGTTCGCCATAATATTGAt                                                 <  2:3316983/150‑1 (MQ=255)
                                                                                                       atgatCTCTTAACCGCGATTAAAAATAAGCAAAATGATTTGCAAGCGATGCACCATAAATATGTGCCGATCGCAGTGAAGATCGCGCCGGATCTTTCTGAAGAAGAATTGATCCAGGTTGCCGATAGTTTAGTTCGCCATAATATTGATg                                                >  1:4267681/1‑150 (MQ=255)
                                                                                                       atgatCTCTTAACCGCGATTAAAAATAAGCAAAATGATTTGCAAGCGATGCACCATAAATATGTGCCGATCGCAGTGAAGATCGCGCCGGATCTTTCTGAAGAAGAATTGATCCAGGTTGCCGATAGTTTAGTTCGCCATAATATTGATg                                                >  1:60719/1‑150 (MQ=255)
                                                                                                        tgatCTCTTAACCGCGATTAAAAATAAGCAAAATGATTTGCAAGCGATGCACCATAAATATGTGCCGATCGCAGTGAAGATCGCGCCGGATCTTTCTGAAGAAGAATTGATCCAGGTTGCCGATAGTTTAGTTCGCCATAATATTGATgg                                               <  2:1268060/150‑1 (MQ=255)
                                                                                                        tgatCTCTTAACCGCGATTAAAAATAAGCAAAATGATTTGCAAGCGATGCACCATAAATATGTGCCGATCGCAGTGAAGATCGCGCCGGATCTTTCTGAAGAAGAATTGATCCAGGTTGCCGATAGTTTAGTTCGCCATAATATTGATgg                                               <  2:4214689/150‑1 (MQ=255)
                                                                                                          atCTCTTAACCGCGATTAAAAATAAGCAAAATGATTTGCAAGCGATGCACCATAAATATGTGCCGATCGCAGTGAAGATCGCGCCGGATCTTTCTGAAGAAGAATTGATCCAGGTTGCCGATAGTTTAGTTCGCCATAATATTGATGGCg                                             <  1:4162312/150‑1 (MQ=255)
                                                                                                            ctctTAACCGCGATTAAAAATAAGCAAAATGATTTGCAAGCGATGCACCATAAATATGTGCCGATCGCAGTGAAGATCGCGCCGGATCTTTCTGAAGAAGAATTGATCCAGGTTGCCGATAGTTTAGTTCGCCATAATATTGATGGCGtt                                           >  2:4561720/1‑150 (MQ=255)
                                                                                                             tctTAACCGCGATTAAAAATAAGCAAAATGATTTGCAAGCGATGCACCATAAATATGTGCCGATCGCAGTGAAGATCGCGCCGGATCTTTCTGAAGAAGAATTGATCCAGGTTGCCGATAGTTTAGTTCGCCATAATATTGATGGCGTTa                                          >  1:1009987/1‑150 (MQ=255)
                                                                                                              ctTAACCGCGATTAAAAATAAGCAAAATGATTTGCAAGCGATGCACCATAAATATGTGCCGATCGCAGTGAAGATCGCGCCGGATCTTTCTGAAGAAGAATTGATCCAGGTTGCCGATAGTTTAGTTCGCCATAATATTGATGGCGTTAt                                         <  1:158857/150‑1 (MQ=255)
                                                                                                                 aaCCGCGATTAAAAATAAGCAAAATGATTTGCAAGCGATGCACCATAAATATGTGCCGATCGCAGTGAAGATCGCGCCGGATCTTTCTGAAGAAGAATTGATCCAGGTTGCCGATAGTTTAGTTCGCCATAATATTGATGGCGTTATTGc                                      >  1:799454/1‑150 (MQ=255)
                                                                                                                 aaCCGCGATTAAAAATAAGCAAAATGATTTGCAAGCGATGCACCATAAATATGTGCCGATCGCAGTGAAGATCGCGCCGGATCTTTCTGAAGAAGAATTGATCCAGGTTGCCGATAGTTTAGTTCGCCATAATATTGATGGCGTTATTGc                                      >  1:4181518/1‑150 (MQ=255)
                                                                                                                      cgATTAAAAATAAGCAAAATGATTTGCAAGCGATGCACCATAAATATGTGCCGATCGCAGTGAAGATCGCGCCGGATCTTTCTGAAGAAGAATTGATCCAGGTTGCCGATAGTTTAGTTCGCCATAATATTGATGGCGTTATTGCAACCa                                 >  1:3423038/1‑150 (MQ=255)
                                                                                                                       gATTAAAAATAAGCAAAATGATTTGCAAGCGATGCACCATAAATATGTGCCGATCGCAGTGAAGATCGCGCCGGATCTTTCTGAAGAAGAATTGATCCAGGTTGCCGATAGTTTAGTTCGCCATAATATTGATGGCGTTATTGCAACCaa                                <  2:101451/150‑1 (MQ=255)
                                                                                                                        aTTAAAAATAAGCAAAATGATTTGCAAGCGATGCACCATAAATATGTGCCGATCGCAGTGAAGATCGCGCCGGATCTTTCTGAAGAAGAATTGATCCAGGTTGCCGATAGTTTAGTTCGCCATAATATTGATGGCGTTATTGCAACCAAt                               <  1:139338/150‑1 (MQ=255)
                                                                                                                          tAAAAATAAGCAAAATGATTTGCAAGCGATGCACCATAAATATGTGCCGATCGCAGTGAAGATCGCGCCGGATCTTTCTGAAGAAGAATTGATCCAGGTTGCCGATAGTTTAGTTCGCCATAATATTGATGGCGTTATTGCAACCAATAc                             <  1:2442900/150‑1 (MQ=255)
                                                                                                                          tAAAAATAAGCAAAATGATTTGCAAGCGATGCACCATAAATATGTGCCGATCGCAGTGAAGATCGCGCCGGATCTTTCTGAAGAAGAATTGATCCAGGTTGCCGATAGTTTAGTTCGCCATAATATTGATGGCGTTATTGCAACCAATAc                             >  2:3139721/1‑150 (MQ=255)
                                                                                                                          tAAAAATAAGCAAAATGATTTGCAAGCGATGCACCATAAATATGTGCCGATCGCAGTGAAGATCGCGCCGGATCTTTCTGAAGAAGAATTGATCCAGGTTGCCGATAGTTTAGTTCGCCATAATATTGATGGCGTTATTGCAACCAATAc                             >  1:47055/1‑150 (MQ=255)
                                                                                                                             aaaTAAGCAAAATGATTTGCAAGCGATGCACCATAAATATGTGCCGATCGCAGTGAAGATCGCGCCGGATCTTTCTGAAGAAGAATTGATCCAGGTTGCCGATAGTTTAGTTCGCCATAATATTGATGGCGTTATTGCAACCAATACcac                          >  1:3204787/1‑150 (MQ=255)
                                                                                                                               aTAAGCAAAATGATTTGCAAGCGATGCACCATAAATATGTGCCGATCGCAGTGAAGATCGCGCCGGATCTTTCTGAAGAAGAATTGATCCAGGTTGCCGATAGTTTAGTTCGCCATAATATTGATGGCGTTATTGCAACCAATACCacac                        <  1:1979908/150‑1 (MQ=255)
                                                                                                                                tAAGCAAAATGATTTGCAAGCTATGCACCATAAATATGTGCCGATCGCAGTGAAGATCGCGCCGGATCTTTCTGAAGAAGAATTGATCCAGGTTGCCGATAGTTTAGTTCGCCATAATATTGATGGCGTTATTGCAACCAATACCACACt                       <  1:2579503/150‑1 (MQ=255)
                                                                                                                                tAAGCAAAATGATTTGCAAGCGATGCACCATAAATATGTGCCGATCGCAGTGAAGATCGCGCCGGATCTTTCTGAAGAAGAATTGATCCAGGTTGCCGATAGTTTAGTTCGCCATAATATTGATGGCGTTATTGCAACCAATACCACACt                       <  1:1001741/150‑1 (MQ=255)
                                                                                                                                  aGCAAAATGATTTGCAAGCGATGCACCATAAATATGTGCCGATCGCAGTGAAGATCGCGCCGGATCTTTCTGAAGAAGAATTGATCCAGGTTGCCGATAGTTTAGTTCGCCATAATATTGATGGCGTTATTGCAACCAATACCACACTCg                     <  2:3497930/150‑1 (MQ=255)
                                                                                                                                     aaaaTGATTTGCAAGCGATGCACCATAAATATGTGCCGATCGCAGTGAAGATCGCGCCGGATCTTTCTGAAGAAGAATTGATCCAGGTTGCCGATAGTTTAGTTCGCCATAATATTGATGGCGTTATTGCAACCAATACCACACTCGATc                  <  2:4436683/150‑1 (MQ=255)
                                                                                                                                      aaaTGATTTGCAAGCGATGCACCATAAATATGTGCCGATCGCAGTGAAGATCGCGCCGGATCTTTCTGAAGAAGAATTGATCCAGGTTGCCGATAGTTTAGTTCGCCATAATATTGATGGCGTTATTGCAACCAATACCACACTCGATCg                 <  2:2631814/150‑1 (MQ=255)
                                                                                                                                      aaaTGATTTGCAAGCGATGCACCATAAATATGTGCCGATCGCAGTGAAGATCGCGCCGGATCTTTCTGAAGAAGAATTGATCCAGGTTGCCGAGAGTTTAGTTCGCCATAATATTGATGGCGTTATTGCAACCAATACCACACTCGATCg                 >  2:2808848/1‑150 (MQ=255)
                                                                                                                                        aTGATTTGCAAGCGATGCACCATAAATATGTGCCGATCGCAGTGAAGATCGCGCCGGATCTTTCTGAAGAAGAATTGATCCAGGTTGCCGATAGTTTAGTTCGCCATAATATTGATGGCGTTATTGCAACCAATACCACACTCGATCGtt               >  2:2975704/1‑150 (MQ=255)
                                                                                                                                               gCAAGCGATGCACCATAAATATGTGCCGATCGCAGTGAAGATCGCGCCGGATCTTTCTGAAGAAGAATTGATCCAGGTTGCCGATAGTTTAGTTCGCCATAATATTGATGGCGTTATTGCAACCAATACCACACTCGATCGTTCTCTTGt        <  1:1081032/150‑1 (MQ=255)
                                                                                                                                                cAAGCGATGCACCATAAATATGTGCCGATCGCAGTGAAGATCGCGCCGGATCTTTCTGAAGAAGAATTGATCCAGGTTGCCGATAGTTTAGTTCGCCATAATATTGATGGCGTTATTGCAACCAATACCACACTCGATCGTTCTCTTGtt       <  1:3158594/150‑1 (MQ=255)
                                                                                                                                                 aaGCGATGCACCATAAATATGTGCCGATCGCAGTGAAGATCGCGCCGGATCTTTCTGAAGAAGAATTGATCCAGGTTGCCGATAGTTTAGTTCGCCATAATATTGATGGCGTTATTGCAACCAATACCACACTCGATCGTTCTCTTGTTc      <  2:4110843/150‑1 (MQ=255)
                                                                                                                                                 aaGCGATGCACCATAAATATGTGCCGATCGCAGTGAAGATCGCGCCGGATCTTTCTGAAGAAGAATTGATCCAGGTTGCCGATAGTTTAGTTCGCCATAATATTGATGGCGTTATTGCAACCAATACCACACTCGATCGTTCTCTTGTTc      <  2:1207880/150‑1 (MQ=255)
                                                                                                                                                  aGCGATGCACCATAAATATGTGCCGATCGCAGTGAAGATCGCGCCGGATCTTTCTGAAGAAGAATTGATCCAGGTTGCCGATAGTTTAGTTCGCCATAATATTGATGGCGTTATTGCAACCAATACCACACTCGATCGTTCTCTTGTTCa     <  2:1109109/150‑1 (MQ=255)
                                                                                                                                                   gCGATGCACCATAAATATGTGCCGATCGCAGTGAAGATCGCGCCGGATCTTTCTGAAGAAGAATTGATCCAGGTTGCCGATAGTTTAGTTCGCCATAATATTGATGGCGTTATTGCAACCAATACCACACTCGATCGTTCTCTTGTTCAg    <  2:1061595/150‑1 (MQ=255)
                                                                                                                                                     gATGCACCATAAATATGTGCCGATCGCAGTGAAGATCGCGCCGGATCTTTCTGAAGAAGAATTGATCCAGGTTGCCGATAGTTTAGTTCGCCATAATATTGATGGCGTTATTGCAACCAATACCACACTCGATCGTTCTCTTGTTCAggg  >  1:43855/1‑150 (MQ=255)
                                                                                                                                                     gATGCACCATAAATATGTGCCGATCGCAGTGAAGATCGCGCCGGATCTTTCTGAAGAAGAATTGATCCAGGTTGCCGATAGTTTAGTTCGCCATAATATTGATGGCGTTATTGCAACCAATACCACACTCGATCGTTCTCTTGTTCAggg  >  1:3182070/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
ATTTGTATGGAAAAAATCTATGCCTATGCGGGATATATCGCCATCAATATTTCATCGCCGAATACCCCAGGATTACGCACGCTGCAATATGGTGAAGCGCTGGATGATCTCTTAACCGCGATTAAAAATAAGCAAAATGATTTGCAAGCGATGCACCATAAATATGTGCCGATCGCAGTGAAGATCGCGCCGGATCTTTCTGAAGAAGAATTGATCCAGGTTGCCGATAGTTTAGTTCGCCATAATATTGATGGCGTTATTGCAACCAATACCACACTCGATCGTTCTCTTGTTCAGGG  >  CP009273/1000692‑1000990

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: