Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,109,544 C→A 5.8% T976T (ACG→ACT yghJ ← DUF4092 family putative lipoprotein peptidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,109,5440CA 5.8% 1067.4 / 25.0 283T976T (ACG→ACTyghJDUF4092 family putative lipoprotein peptidase
Reads supporting (aligned to +/- strand):  ref base C (139/127);  new base A (7/10);  total (146/137)
Fisher's exact test for biased strand distribution p-value = 4.56e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 290 total aligned reads displayed.

AGCGGGTAGCTCGGGTTCATCATGCCCGCTTTGCTGCTACCGTCACCGTAGATCATGTTGTTATCGACCAGCGATTTTTTCAGATCTGCAGAACACTTGGTGCCGCCTGCATAGGCATCATTGGCGTAGCAGTTCAGGAACTCGGTGAACGTTTTAAAGCCCAGCTCGTCATTTTTGCCTTCTTCATAACGATAGCTCGTATCGTTCCACAGCCAGACCGACATGTTCTGGTACAGACGTTCCAGATCGACGCTGCTCAGACGCTCACCTTTGCGACCATTGGTCCGGAAGTAGAGCTC  >  CP009273/3109395‑3109693
                                                                                                                                                     |                                                                                                                                                     
aGCGGGTAGCTCGGGTTCATCATGCCCGCTTTGCTGCTACCGTCACCGTAGATCATGTTGTTATCGACCAGCGATTTTTTCAGATCTGCAGAACACTTGGTGCCGCCTGCATAGGCATCATTGGCGTAGCAGTTCAGGAACTCGGTGAAc                                                                                                                                                       <  2:1133039/150‑1 (MQ=255)
aGCGGGTAGCACGGGATCATCATGCCCGCTTAGCTGCTACCGTCACCGTAGATCATGTTGTTATCGAACAGCGATTTTTTCAGATCTGCAGAACACTTGGTGCCGCCTGCATAGGCATCATTGGCGTAGCAGTTCAGGAACACGGTGAAc                                                                                                                                                       >  2:2842998/1‑150 (MQ=255)
   gggTAGCTCGGGTTCATCATGCCCGCTTTGCTGCTACCGTCACCGTAGATCATGTTGTTATCGACCAGCGATTTTTTCAGATCTGCAGAACACTTGGTGCCGCCTGCATAGGCATCATTGGCGTAGCAGTTCAGGAACTCGGTGAACGtt                                                                                                                                                    >  2:2476292/1‑150 (MQ=255)
        gCTCGTGTTCATCATGCCCGCTTTGCTGCTACCGTCACCGTAGATCATGTTGTTATCGACCAGCGATTTTTTCAGATCTGCAGAACACTTGGTGCCGCCTGCATAGGCATCATTGGCGTAGCAGTTCAGGAACTCGGTGAACGTTTTaaa                                                                                                                                               >  1:2877648/1‑150 (MQ=255)
          tCGGGTTCATCATGCCCGCTTTGCTGCTACCGTCACCGTAGATCATGTTGTTATCGACCAGCGATTTTTTCAGATCTGCAGAACACTTGGTGCCGCCTGCATAGGCATCATTGGCGTAGCAGTTCAGGAACTCGGTGAACGTTTTAAAGc                                                                                                                                             <  2:863601/150‑1 (MQ=255)
          tCGGGTTCATCATGCCCGCTTTGCTGCTACCGTCACCGTAGATCATGTTGTTATCGACCAGCGATTTTTTCAGATCTGCAGAACACTTGGTGCCGCCTGCATAGGCATCATTGGCGTAGCAGTTCAGGAACTCGGTGAACGTTTTAAAGc                                                                                                                                             >  1:4325572/1‑150 (MQ=255)
             ggTTCATCATGCCCGCTTTGCTGCTACCGTCACCGTAGATCATGTTGTTATCGACCAGCGATTTTTTCAGATCTGCAGAACACTTGGTGCCGCCTGCATAGGCATCATTGGCGTAGCAGTTCAGGAACTCGGTGAACGTTTTAAAGCCCa                                                                                                                                          >  1:715836/1‑150 (MQ=255)
                  atcatGCCCGCTTTGCTGCTACCGTCACCGTAGATCATGTTGTTATCGACCAGCGATTTTTTCAGATCTGCAGAACACTTGGTGCCGCCTGCATAGGCATCATTGGCGTAGCAGTTCAGGAACTCGGTGAACGTTTTAAAGCCCAGCTCg                                                                                                                                     >  1:631133/1‑150 (MQ=255)
                  atcatGCCCGCTTTGCTGCTACCGTCACCGTAGATCATGTTGTTATCGACCAGCGATTTTTTCAGATCTGCAGAACACTTGGTGCCGCCTGCATAGGCATCATTGGCGTAGCAGTTCAGGAACTCGGTGAACGTTTTAAAGCCCAGCTCg                                                                                                                                     >  1:2307671/1‑150 (MQ=255)
                   tcatGCCCGCTTTGCTGCTACCGTCACCGTAGATCATGTTGTTATCGACCAGCGATTTTTTCAGATCTGCAGAACACTTGGTGCCGCCTGCATAGGCATCATTGGCGTAGCAGTTCAGGAACTCGGTGAAAGTTTTAAAGCCCAGCTCGt                                                                                                                                    <  1:1572774/150‑1 (MQ=255)
                      tGCCCGCTTTGCTGCTACCGTCACCGTAGATCATGTTGTTATCGACCAGCGATTTTTTCAGATCTGCAGAACACTTGGTGCCGCCTGCATAGGCATCATTGGCGTAGCAGTTCAGGAACTCGGTGAACGTTTTAAAGCCCAGCTCGTCAt                                                                                                                                 >  1:758164/1‑150 (MQ=255)
                         ccGCTTTGCTGCTACCGTCACCGTAGATCATGTTGTTATCGACCAGCGATTTTTTCAGATCTGCAGAACACTTGGTGCCGCCTGCATAGGCATCATTGGCGTAGCAGTTCAGGAACTCGGTGAACGTTTTAAAGCCCAGCTCGTCAtttt                                                                                                                              <  2:3781573/150‑1 (MQ=255)
                          cGCTTTGCTGCTACCGTCACCGTAGATCATGTTGTTATCGACCAGCGATTTTTTCAGATCTGCAGAACACTTGGTGCCGCCTGCATAGGCATCATTGGCGTAGCAGTTCAGGAACTCGGTGAACGTTTTAAAGCCCAGCTCGTCAttttt                                                                                                                             >  1:2020078/1‑150 (MQ=255)
                          cGCTTTGCTGCTACCGTCACCGTAGATCATGTTGTTATCGACCAGCGATTTTTTCAGATCTGCAGAACACTTGGTGCCGCCTGCATAGGCATCATTGGCGTAGCAGTTCAGGAACTCGGTGAACGTTTTAAAGCCCAGCTCGTCAttttt                                                                                                                             <  1:529991/150‑1 (MQ=255)
                          cGCTTTGCTGCTACCGTCACCGTAGATCATGTTGTTATCGACCAGCGATTTTTTCAGATCTGCAGAACACTTGGTGCCGCCTGCATAGGCATCATTGGCGTAGCAGTTCAGGAACTCGGTGAACGTTTTAAAGCCCAGCTCGTCAttttt                                                                                                                             <  1:2757326/150‑1 (MQ=255)
                          cGCTTTGCTGCTACCGTCACCGTAGATCATGTTGTTATCGACCAGCGATTTTTTCAGATCTGCAGAACACTTGGTGCCGCCTGCATAGGCATCATTGGCGTAGCAGTTCAGGAACTCGGTGAAAGTTTTAAAGCCCAGCTCGTCAttttt                                                                                                                             <  2:470270/150‑1 (MQ=255)
                          cGCTTTGCTGCTACCGTCACCGTAGATCATGTTGTTATCGACCAGCGATTTTTTCAGATCTGCAGAACACTTGGTGCCGCCTGCATAGGCATCATTGGCGTAGCAGTTCAGGAACTCGGTGAAAGTTTTAAAGCCCAGCTCGTCAttttt                                                                                                                             <  1:3493903/150‑1 (MQ=255)
                           gCTTTGCTGCTACCGTCACCGTAGATCATGTTGTTATCGACCAGCGATTTTTTCAGATCTGCAGAACACTTGGTGCCGCCTGCATAGGCATCATTGGCGTAGCAGTTCAGGAACTCGGTGAACGTTTTAAAGCCCAGCTCGTCATTTTTg                                                                                                                            <  2:3191157/150‑1 (MQ=255)
                              ttGCTGCTACCGTCACCGTAGATCATGTTGTTATCGACCAGCGATTTTTTCAGATCTGCAGAACACTTGGTGCCGCCTGCATAGGCATCATTGGCGTAGCAGTTCAGGAACTCGGTGAACGTTTTAAAGCCCAGCTCGTCATTTTTGCct                                                                                                                         >  1:3355246/1‑150 (MQ=255)
                              gtgctgcTACCGTCACAGTAGATCATGTTGTTCTCGACCAGCGATTTTTTCAGATCTGCAGTACTGTTGGTGCCGCCTGCATAGGCATCATTGGCGTAGCAGTTCAGGATCTCGGTGAACGTTTTAAAGCCCAGCTCGTCATTTTTGCct                                                                                                                         <  2:593849/149‑1 (MQ=255)
                                 ctgctACCGTCACCGTAGATCATGTTGTTATCGACCAGCGATTTTTTCAGATCTGCAGAACACTTGGTGCCGCCTGCATAGGCATCATTGGCGTAGCAGTTCAGGAACTCGGTGAACGTTTTAAAGCCCAGCTCGTCATTTTTGCcttct                                                                                                                      >  2:628224/1‑150 (MQ=255)
                                    ctACCGTCACCGTAGATCATGTTGTTATCGACCAGCGATTTTTTCAGATCTGCAGAACACTTGGTGCCGCCTGCATAGGCATCATTGGCGTAGCAGTTCAGGAACTCGGTGAACGTTTTAAAGCCCAGCTCGTCATTTTTGCCTTCTTCa                                                                                                                   >  1:2839780/1‑150 (MQ=255)
                                       ccGTCACCGTAGATCATGTTGTTATCGACCAGCGATTTTTTCAGATATGCAGAACACTTGGTGCCGCCTGCATAGGCATCATTGGCGTAGCAGTTCAGGAACTCGGTGACCGTTTTAAAGCCCAGCTCTTCATTTTTGCCTTCTTCATaa                                                                                                                >  2:4407443/1‑150 (MQ=255)
                                        cGTCACCGTAGATCATGTTGTTATCGACCAGCGATTTTTTCAGATCTGCAGAACACTTGGTGCCGCCTGCATAGGCATCATTGGCGTAGCAGTTCAGGAACTCGGTGAACGTTTTAAAGCCCAGCTCGTCATTTTTGCCTTCTTCATAAc                                                                                                               <  2:3239797/150‑1 (MQ=255)
                                            aCCGTAGATCATGTTGTTATCGACCAGCGATTTTTTCAGATCTGCAGAACACTTGGTGCCGCCTGCATAGGCATCATTGGCGTAGCAGTTCAGGAACTCGGTGAACGTTTTAAAGCCCAGCTCGTCATTTTTGCCTTCTTCATAACGATa                                                                                                           >  2:148688/1‑150 (MQ=255)
                                             ccGTAGATCATGTTGTTATCGACCAGCGATTTTTTCAGATCTGCAGAACACTTGGTGCCGCCTGCATAGGCATCATTGGCGTAGCAGTTCAGGAACTCGGTGAACGTTTTAAAGCCCAGCTCGTCATTTTTGCCTTCTTCATAACGATAg                                                                                                          <  1:1808721/150‑1 (MQ=255)
                                             ccGTAGATCATGTTGTTATCGACCAGCGATTTTTTCAGATCTGCAGAACACTTGGTGCCGCCTGCATAGGCATCATTGGCGTAGCAGTTCAGGAACTCGGTGAACGTTTTAAAGCCCAGCTCGTCATTTTTGCCTTCTTCATAACGATAg                                                                                                          >  1:914614/1‑150 (MQ=255)
                                             ccGTAGATCATGTTGTTATCGACCAGCGATTTTTTCAGATCTGCAGAACACTTGGTGCCGCCTGCATAGGCATCATTGGCGTAGCAGTTCAGGAACTCGGTGAACGTTTTAAAGCCCAGCTCGTCATTTTTGCCTTCTTCATAACGATAg                                                                                                          <  1:4061883/150‑1 (MQ=255)
                                             ccGTAGATCATGTTGTTATCGACCAGCGATTTTTTCAGATCTGCAGAACACTTGGTGCCGCCTGCATAGGCATCATTGGCGTAGCAGTTCAGGAACTCGGTGAACGTTTTAAAGCCCAGCTCGTCATTTTTGCCTTCTTCATAACGATAg                                                                                                          <  1:526622/150‑1 (MQ=255)
                                             ccGTAGATCATGTTGTTATCGACCAGCGATTTTTTCAGATCTGCAGAACACTTGGTGCCGCCTGCATAGGCATCATTGGCGTAGCAGTTCAGGAACTCGGTGAACGTTTTAAAGCCCAGCTCGTCATTTTTGCCTTCTTCATAACGATAg                                                                                                          >  1:3728785/1‑150 (MQ=255)
                                              cGTAGATCATGTTGTTATCGACCAGCGATTTTTTCAGATCTGCAGAACACTTGGTGCCGCCTGCATAGGCATCATTGGCGTAGCAGTTCAGGAACTCGGTGAACGTTTTAAAGCCCAGCTCGTCATTTTTGCCTTCTTCATAACGATAGc                                                                                                         >  2:489350/1‑150 (MQ=255)
                                               gTAGATCATGTTGTTATCGACCAGCGATTTTTTCAGATCTGCAGAACACTTGGTGCCGCCTGCATAGGCATCATTGGCGTAGCAGTTCAGGAACTCGGTGAACGTTTTAAAGCCCAGCTCGTCATTTTTGCCTTCTTCATAACGATAGCt                                                                                                        <  1:4276739/150‑1 (MQ=255)
                                                tAGATCATGTTGTTATCGACCAGCGATTTTTTCAGATCTGCAGAACACTTGGTGCCGCCTGCATAGGCATCATTGGCGTAGCAGTTCAGGAACTCGGTGAACGTTTTAAAGCCCAGCTCGTCATTTTTGCCTTCTTCATAACGATAGCTc                                                                                                       >  1:1632310/1‑150 (MQ=255)
                                                    tCATGTTGTTATCGACCAGCGATTTTTTCAGATCTGCAGAACACTTGGTGCCGCCTGCATAGGCATCATTGGCGTAGCAGTTCAGGAACTCGGTGAACGTTTTAAAGCCCAGCTCGTCATTTTTGCCTTCTTCATAACGATAGCTCGTAt                                                                                                   <  2:2527157/150‑1 (MQ=255)
                                                    tCATGTTGTTATCGACCAGCGATTTTTTCAGATCTGCAGAACACTTGGTGCCGCCTGCATAGGCATCATTGGCGTAGCAGTTCAGGAACTCGGTGAACGTTTTAAAGCCCAGCTCGTCATTTTTGCCTTCTTCATAACGATAGCTCGTAt                                                                                                   <  2:434630/150‑1 (MQ=255)
                                                    tCATGTTGTTATCGACCAGCGATTTTTTCAGATCTGCAGAACACTTGGTGCCGCCTGCATAGGCATCATTGGCGTAGCAGTTCAGGAACTCGGTGAACGTTTTAAAGCCCAGCTCGTCATTTTTGCCTTCTTCATAACGATAGCTCGTAt                                                                                                   >  1:3215073/1‑150 (MQ=255)
                                                     cATGTTGTTATCGACCAGCGATTTTTTCAGATCTGCAGAACACTTGGTGCCGCCTGCATAGGCATCATTGGCGTAGCAGTTCAGGAACTCGGTGAACGTTTTAAAGCCCAGCTCGTCATTTTTGCCTTCTTCATAACGATAGCTCGTATc                                                                                                  <  2:4492623/150‑1 (MQ=255)
                                                        gttggtATCGACCAGCGATTTTTTCAGATCTGCAGAACACTTGGTGCCGCCTGCATAGGCATCATTGGCGTAGCAGTTCAGGAACTCGGTGAACGTTTTAAAGCCCAGCTCGTCATTTTTGCCTTCTTCATAACGATAGCTCGTATCGtt                                                                                               <  2:4208495/150‑1 (MQ=255)
                                                         ttgttATCGACCAGCGTTTTTTTCAGATCTGCAGAACACTTGGTGCCGCCTGCATAGGCATCATTGGCGTAGCAGTTCAGGAACTCGGTGAACGTTTTAAAGCCCAGCTCGTCATTTTTGCCTTCTTCATAACGATAGCTCGTATCGTTc                                                                                              <  2:2190177/150‑1 (MQ=255)
                                                                cGACCAGCGATTTTTTCAGATCTGCAGAACACTTGGTGCCGCCTGCATAGGCATCATTGGCGTAGCAGTTCAGGAACTCGGTGAACGTTTTAAAGCCCAGCTCGTCATTTTTGCCTTCTTCATAACGATAGCTCGTATCGTTCCACAGcc                                                                                       <  2:401197/150‑1 (MQ=255)
                                                                 gACCAGCGATTTTTTCAGATCTGCAGAACACTTGGTGCCGCCTGCATAGGCATCATTGGCGTAGCAGTTCAGGAACTCGGTGAACGTTTTAAAGCCCAGCTCGTCATTTTTGCCTTCTTCATAACGATAGCTCGTATCGTTCCACAGCCa                                                                                      <  2:1606857/150‑1 (MQ=255)
                                                                 gACCAGCGATTTTTTCAGATCTGCAGAACACTTGGTGCCGCCTGCATAGGCATCATTGGCGTAGCAGTTCAGGAACTCGGTGAACGTTTTAAAGCCCAGCTCGTCATTTTTGCCTTCTTCATAACGATAGCTCGTATCGTTCCACAGCCa                                                                                      <  2:3866473/150‑1 (MQ=255)
                                                                   ccAGCGATTTTTTCAGATCTGCAGAACACTTGGTGCCGCCTGCATAGGCATCATTGGCGTAGCAGTTCAGGAACTCGGTGAACGTTTTAAAGCCCAGCTCGTCATTTTTGCCTTCTTCATAACGATAGCTCGTATCGTTCCACAGCCAGa                                                                                    >  1:2530683/1‑150 (MQ=255)
                                                                   ccAGCGATTTTTTCAGATCTGCAGAACACTTGGTGCCGCCTGCATAGGCATCATTGGCGTAGCAGTTCAGGAACTCGGTGAACGTTTTAAAGCCCAGCTCGTCATTTTTGCCTTCTTCATAACGATAGCTCGTATCGTTCCACAGCCAGa                                                                                    >  1:454496/1‑150 (MQ=255)
                                                                    cAGCGATTTTTTCAGATCTGCAGAACACTTGGTGCCGCCTGCATAGGCATCATTGGCGTAGCAGTTCAGGAACTCGGTGAACGTTTTAAAGCCCAGCTCGTCATTTTTGCCTTCTTCATAACGATAGCTCGTATCGTTCCACAGCCAGAc                                                                                   >  1:2138517/1‑150 (MQ=255)
                                                                       cGATTTTTTCAGATCTGCAGAACACTTGGTGCCGCCTGCATAGGCATCATTGGCGTAGCAGTTCAGGAACTCGGTGAACGTTTTAAAGCCCAGCTCGTCATTTTTGCCTTCTTCATAACGATAGCTCGTATCGTTCCACAGCCAGACCGa                                                                                <  1:1239635/150‑1 (MQ=255)
                                                                        gATTTTTTCAGATCTGCAGAACACTTGGTGCCGCCTGCATAGGCATCATTGGCGTAGCAGTTCAGGAACTCGGTGAACGTTTTAAAGCCCAGCTCGTCATTTTTGCCTTCTTCATAACGATAGCTCGTATCGTTCCACAGCCAGACCGAc                                                                               >  1:3149224/1‑150 (MQ=255)
                                                                        gATTTTTTCAGATCTGCAGAACACTTGGTGCCGCCTGCATAGGCATCATTGGCGTAGCAGTTCAGGAACTCGGTGAACGTTTTAAAGCCCAGCTCGTCATTTTTGCCTTCTTCATAACGATAGCTCGTATCGTTCCACAGCCAGACCGAc                                                                               >  2:1384306/1‑150 (MQ=255)
                                                                              ttCAGATCTGCAGAACACTTGGAGCCGCCTGCATAGGCATCATTGGCGTAGCAGATCAGGAACTCGGTGAACTTTTTAAAGCCCAGCTCGTCATTTTTGCCTTCTTCATAACGAAAGCTCGTATCGATCCACAGCCAGACCGACATGTTc                                                                         >  2:2636041/1‑150 (MQ=255)
                                                                                cAGATCTGCAGAACACTTGGTGCCGCCTGCATAGGCATCATTGGCGTAGCAGTTCAGGAACTCGGTGAACGTTTTAAAGCCCAGCTCGTCATTTTTGCCTTCTTCATAACGATAGCTCGTATCGTTCCACAGCCAGACCGACATGTTCTg                                                                       >  1:3100833/1‑150 (MQ=255)
                                                                                 aGATCTGCAGAACACTTGGTGCCGCCTGCATAGGCATCATTGGCGCAGCAGTTCAGGAACTCGGTGAACGTTTTAAAGCCCAGCTCGTCATTTTTGCCTTCTTCATAACGATAGCTCGTATCGTTCCACAGCCAGACCGACATGTTCTgg                                                                      >  1:2795588/1‑150 (MQ=255)
                                                                                    tCTGCAGAACACTTGGTGCCGCCTGCATAGGCATCATTGGCGTAGCAGTTCAGGAACTCGGTGAACGTTTTAAAGCCCAGCTCGTCATTTTTGCCTTCTTCATAACGATAGCTCGTATCGTTCCACAGCCAGACCGACATGTTCTGGTAc                                                                   <  2:1753100/150‑1 (MQ=255)
                                                                                        cAGAACACTTGGTGCCGCCTGCATAGGCATCATTGGCGTAGCAGTTCAGGAACTCGGTGAACGTTTTAAAGCCCAGCTCGTCATTTTTGCCTTCTTCATAACGATAGCTCGTATCGTTCCACAGCCAGACCGACATGTTCTGGTACAGAc                                                               <  2:2236473/150‑1 (MQ=255)
                                                                                        cAGAACACTTGGTGCCGCCTGCATAGGCATCATTGGCGTAGCAGTTCAGGAACTCGGTGAACGTTTTAAAGCCCAGCTCGTCATTTTTGCCTTCTTCATAACGATAGCTCGTATCGTTCCACAGCCAGACCGACATGTTCTGGTACAGAc                                                               >  1:1052362/1‑150 (MQ=255)
                                                                                             cacTTGGTGCCGCCTGCATAGGCATCATTGGCGTAGCAGTTCAGGAACTCGGTGAACGTTTTAAAGCCCAGCTCGTCATTTTTGCCTTCTTCATAACGATAGCTCGTATCGTTCCACAGCCAGACCGACATGTTCTGGTACAGACGTTcc                                                          >  2:1709658/1‑150 (MQ=255)
                                                                                              acTTGGTGCCGCCTGCATAGGCATCATTGGCGTAGCAGTTCAGGAACTCGGTGAACGTTTTAAAGCCCAGCTCGTCATTTTTGCCTTCTTCATAACTATAGCTCGTATCGTTCCACAGCCAGACCGACATGTTCTGGTACAGACGTTCCa                                                         >  2:4291059/1‑150 (MQ=255)
                                                                                              acTTGGTGCCGCCTGCATAGGCATCATTGGCGTAGCAGTTCAGGAACTCGGTGAACGTTTTAAAGCCCAGCTCGTCATTTTTGCCTTCTTCATAACGATAGCTCGTATCGTTCCACAGCCAGACCGACATGTTCTGGTACAGACGTTCCa                                                         >  2:2324454/1‑150 (MQ=255)
                                                                                               cTTGGTGCCGCCTGCATAGGCATCATTGGCGTAGCAGTTCAGGAACTCGGTGAACGTTTTAAAGCCCAGCTCGTCATTTTTGCCTTCTTCATAACGATAGCTCGTATCGTTCCACAGCCAGACCGACATGTTCTGGTACAGACGTTCCAg                                                        >  2:3379747/1‑150 (MQ=255)
                                                                                                 tGGTGCCGCCTGCATAGGCATCATTGGCGTAGCAGTTCAGGAACTCGGTGAACGTTTTAAAGCCCAGCTCGTCATTTTTGCCTTCTTCATAACGATAGCTCGTATCGTTCCACAGCCAGACCGACATGTTCTGGTACAGACGTTCCAGAt                                                      <  1:2233999/150‑1 (MQ=255)
                                                                                                 tGGTGCCGCCTGCATAGGCATCATTGGCGTAGCAGTTCAGGAACTCGGTGAACGTTTTAAAGCCCAGCTCGTCATTTTTGCCTTCTTCATAACGATAGCTCGTATCGTTCCACAGCCAGACCGACATGTTCTGGTACAGACGTTCCAGAt                                                      <  1:4101431/150‑1 (MQ=255)
                                                                                                   gTGCCGCCTGCATAGGCATCATTGGCGTAGCAGTTCAGGAACTCGGTGAACGTTTTAAAGCCCAGCTCGTCATTTTTGCCTTCTTCATAACGATAGCTCGTATCGTTCCACAGCCAGACCGACATGTTCTGGTACAGACGTTCCAGATCg                                                    >  2:4073658/1‑150 (MQ=255)
                                                                                                    tGCCGCCTGCATAGGCATCATTGGCGTAGCAGTTCAGGAACTCGGTGAACGTTTTAAAGCCCAGCTCGTCATTTTTGCCTTCTTCATAACGATAGCTCGTATCGTTCCACAGCCAGACCAACATGTTCTGGTACAGACGTTCCAGATCGa                                                   >  2:1224893/1‑150 (MQ=255)
                                                                                                        gccTGCATAGGCATCATTGGCGTAGCAGTTCAGGAACTCGGTGAACGTTTTAACGCCCAGCTCGTCATTTTTGCCTTCTTCATAACGATAGCTCGTATCGTTCCACAGCCAGACCGACATGTTCTGGTACAGACGTTCCAGATCGACgct                                               <  2:4312862/150‑1 (MQ=255)
                                                                                                         ccTGCATAGGCATCATTGGCGTAGCAGTTCAGGAACTCGGTGAACGTTTTAAAGCCCAGCTCGTCATTTTTGCCTTCTTCATCACGATAGCTCGTATCGTTCCACAGCCAGACCGACATGTTCTGGTACAGACGTTCCAGATCGACgctg                                              >  2:1800305/1‑150 (MQ=255)
                                                                                                          cTGCATAGGCATCATTGGCGTAGCAGTTCAGGAACTCGGTGAACGTTTTAAAGCCCAGCTCGTCATTTTTGCCTTCTTCATAACGATAGCTCGTATCGTTCCACAGCCAGACCGACATGTTCTGGTACAGACGTTCCAGATCGACgctgc                                             <  1:3870037/150‑1 (MQ=255)
                                                                                                             cATAGGCATCATTGGCGTAGCAGTTCAGGAACTCGGTGAACGTTTTAAAGCCCAGCTCGTCATTTTTGCCTTCTTCATAACGATAGCTCGTATCGTTCCACAGCCAGACCGACATGTTCTGGTACAGACGTTCCAGATCGACGCTGCTCa                                          <  2:3954151/150‑1 (MQ=255)
                                                                                                             cATAGGCATCATTGGCGTAGCAGTTCAGGAACTCGGTGAACGTTTTAAAGCCCAGCTCGTCATTTTTGCCTTCTTCATAACGATAGCTCGTATCGTTCCACAGCCAGACCGACATGTTCTGGTACAGACGTTCCAGATCGACGCTGCTCa                                          <  1:4495089/150‑1 (MQ=255)
                                                                                                              aTAGGCATCATTGGCGTAGCAGTTCAGGAACTCGGTGAACGTTTTAAAGCCCAGCTCGTCATTTTTGCCTTCTTCATAACGATAGCTCGTATCGTTCCACAGCCAGACCGACATGTTCTGGTACAGACGTTCCAGATCGACGCTGCTCAg                                         >  1:307/1‑150 (MQ=255)
                                                                                                                cggCATCAGTGGCGTAGCAGTACAGGAACTCGGTGAACGTTTTAAAGCCCAGCTCGTCATTTTTGCCTTCTTCATAACGATAGCTCGGATCGTTCCACAGCCAGACCGACCTGTTCTGGTACAGACGTTCCAGATCGACGCTGCTCAGAc                                       <  2:2710611/149‑1 (MQ=255)
                                                                                                                 ggCATCATTGGCGTAGCAGTTCAGGAACTCGGTGAACGTTTTAAAGCCCAGCTCGTCATTTTTGCCTTCTTCATAACGATAGCTCGTATCGTTCCACAGCCAGACCGACATGTTCTGGTACAGACGTTCCAGATCGACGCTGCTCAGACg                                      <  1:3494651/150‑1 (MQ=255)
                                                                                                                   catcatTGGCGTAGCAGTTCAGGAACTCGGTGAACGTTTTAAAGCCCAGCTCGTCATTTTTGCCTTCTTCATAACGATAGCTCGTATCGTTCCACAGCCAGACCGACATGTTCTGGTACAGACGTTCCAGATCGACGCTGCTCAGACGCt                                    >  2:37077/1‑150 (MQ=255)
                                                                                                                    atcatTGGCGTAGCAGTTCAGGAACTCGGTGAACGTTTTAAAGCCCAGCTCGTCATTTTTGCCTTCTTCATAACGATAGCTCGTATCGTTCCACAGCCAGACCGACATGTTCTGGTACAGACGTTCCAGATCGACGCTGCTCAGACGCTc                                   >  1:379633/1‑150 (MQ=255)
                                                                                                                    atcatTGGCGTAGCAGTTCAGGAACTCGGTGAACGTTTTAAAGCCCAGCTCGTCATTTTTGCCTTCTTCATAACGATAGCTCGTATCGTTCCACAGCCAGACCGACATGTTCTGGTACAGACGTTCCAGATCGACGCTGCTCAGACGCTc                                   >  2:1737777/1‑150 (MQ=255)
                                                                                                                     tcatTGGCGTAGCAGTTCGGAACTCGGTGAACGTTTTAAAGCCCAGCTCGTCATTTTTGCCTTCTTCATAACGATAGCTCGTATCGTTCCACAGCCAGACCGACATGTTCTGGTACAGACGTTCCAGATCGACGCTGCTCAGACGCTCAc                                 <  2:3559019/150‑1 (MQ=255)
                                                                                                                      catTGGCGTAGCAGTTCAGGAACTCGGTGAACGTTTTAAAGCCCAGCTCGTCATTTTTGCCTTCTTCATAACGATAGCTCGTATCGTTCCACAGCCAGACCGACATGTTCTGGTACAGACGTTCCAGATCGACGCTGCTCAGACGCTCAc                                 <  2:3558955/150‑1 (MQ=255)
                                                                                                                      catTGGCGTAGCAGTTCAGGAACTCGGTGAACGTTTTAAAGCCCAGCTCGTCATTTTTGCCTTCTTCATAACGATAGCTCGTATCGTTCCACAGCCAGACCGACATGTTCTGGTACAGACGTTCCAGATCGACGCTGCTCAGACGCTCAc                                 >  2:1532887/1‑150 (MQ=255)
                                                                                                                       atTGGCGTAGCAGTTCAGGAACTCGGTGAACGTTTTAAAGCCCAGCTCGTCATTTTTGCCTTCTTCATAACGATAGCTCGTATCGTTCCACAGCCAGACCGACATGTTCTGGTACAGACGTTCCAGATCGACGCTGCTCAGACGCTCAcc                                >  2:1867736/1‑150 (MQ=255)
                                                                                                                       atTGGCGTAGCAGTTCAGGAACTCGGTGAACGTTTTAAAGCCCAGCTCGTCATTTTTGCCTTCTTCATAACGATAGCTCGTATCGTTCCACAGCCAGACCGACATGTTCTGGTACAGACGTTCCAGATCGACGCTGCTCAGACGCTCAcc                                >  2:304172/1‑150 (MQ=255)
                                                                                                                       atTGGCGTAGCAGTTCAGGAACTCGGTGAACGTTTTAAAGCCCAGCTCGTCATTTTTGCCTTCTTCATAACGATAGCTCGTATCGTTCCACAGCCAGACCGACATGTTCTGGTACAGACGTTCCAGATCGACGCTGCTCAGACGCTCAca                                >  2:1053218/1‑149 (MQ=255)
                                                                                                                          ggCGTAGCAGTTCAGGAACTCGGTGAACGTTTTAAAGCCCAGCTCGTCATTTTTGCCTTCTTCATAACGATAGCTCGTATCGTTCCACAGCCAGACCGACATGTTCTGGTACAGACGTTCCAGATCGACGCTGCTCAGACGCTCACCttt                             >  1:4126899/1‑150 (MQ=255)
                                                                                                                          ggCGTAGCAGTTCAGGAACTCGGTGAACGTTTTAAAGCCCAGCTCGTCATTTTTGCCTTCTTCATAACGATAGCTCGTATCGTTCCACAGCCAGACCGACATGTTCTGGTACAGACGTTCCAGATCGACGCTGCTCAGACGCTCACCttt                             <  1:2560381/150‑1 (MQ=255)
                                                                                                                                 cAGTTCAGGAACTCGGTGAACGTTTTAAAGCCCAGCTCGTCATTTTTGCCTTCTTCATAACGATAGCTCGTATCGTTCCACAGCCAGACCGACATGTTCTGGTACAGACGTTCCAGATCGACGCTGCTCAGACGCTCACCTTTGCGACCa                      >  1:3533039/1‑150 (MQ=255)
                                                                                                                                 cAGTTCAGGAACTCGGTGAACGTTTTAAAGCCCAGCTCGTCATTTTTGCCTTCTTCATAACGATAGCTCGTATCGTTCCACAGCCAGACCGACATGTTCTGGTACAGACGTTCCAGATCGACGCTGCTCAGACGCTCACCTTTGCGACCa                      >  1:1031280/1‑150 (MQ=255)
                                                                                                                                  aGTTCAGGAACTCGGTGAACGTTTTAAAGCCCAGCTCGTCATTTTTGCCTTCTTCATAACGATAGCTCGTATCGTTCCACAGCCAGACCGACATGTTCTGGTACAGACGTTCCAGATCGACGCTGCTCAGACGCTCACCTTTGCGACCAt                     >  1:3576294/1‑150 (MQ=255)
                                                                                                                                   gTTCAGGAACTCGGTGAACGTTTTAAAGCCCAGCTCGTCATTTTTGCCTTCTTCATAACGATAGCTCGTATCGTTCCACAGCCAGACCGACATGTTCTGGTACAGACGTTCCAGATCGACGCTGCTCAGACGCTCACCTTTGCGACCCtt                    >  2:3682540/1‑150 (MQ=255)
                                                                                                                                   gTTCAGGAACTCGGTGAACGTTTTAAAGCCCAGCTCGTCATTTTTGCCTTCTTCATAACGATAGCTCGTATCGTTCCACAGCCAGACCGACATGTTCTGGTACAGACGTTCCAGATCGACGCTGCTCAGACGCTCACCTTTGCGACCAtt                    >  1:1891114/1‑150 (MQ=255)
                                                                                                                                    ttCAGGAACTCGGTGAACGTTTTAAAGCCCAGCTCGTCATTTTTGCCTTCTTCATAACGATAGCTCGTATCGTTCCACAGCCAGACCGACATGTTCTGGTACAGACGTTCCAGATCGACGCTGCTCAGACGCTCACCTTTGCGACCATTg                   <  1:4516772/150‑1 (MQ=255)
                                                                                                                                       aGGAACTCGGTGAACGTTTTAAAGCCCAGCTCGTCATTTTTGCCTTCTTCATAACGATAGCTCGTATCGTTCCACAGCCAGACCGACATGTTCTGGTACAGACGTTCCAGATCGACGCTGCTCAGACGCTCACCTTTGCGACCATTGGTc                >  1:2100234/1‑150 (MQ=255)
                                                                                                                                       aGGAACTCGGTGAACGTTTTAAAGCCCAGCTCGTCATTTTTGCCTTCTTCATAACGATAGCTCGTATCGTTCCACAGCCAGACCGACATGTTCTGGTACAGACGTTCCAGATCGACGCTGCTCAGACGCTCACCATTGCGACCATTGGTc                >  2:2922174/1‑150 (MQ=255)
                                                                                                                                        ggAACTCGGTGAACGTTTTAAAGCCCAGCTCGTCATTTTTGCCTTCTTCATAACGATAGCTCGTATCGTTCCACAGCCAGACCGACATGTTCTGGTACAGACGTTCCAGATCGACGCTGCTCAGACGCTCACCTTTGCGACCATTGGTcc               <  2:2020078/150‑1 (MQ=255)
                                                                                                                                              cGGTGAACGTTTTAAAGCCCAGCTCGTCATTTTTGCCTTCTTCATAACGATAGCTCGTATCGTTCCACAGCCAGACCGACATGTTCTGGTACAGACGTTCCAGATCGACGCTGCTCAGACGCTCACCTTTGCGACCATTGGTCCGGAAGt         >  1:64542/1‑150 (MQ=255)
                                                                                                                                               ggTGAACGTTTTAAAGCCCAGCTCGTCATTTTTGCCTTCTTCATAACGATAGCTCGTATCGTTCCACAGCCAGACCGACATGTTCTGGTACAGACGTTCCAGATCGACGCTGCTCAGACGCTCACCTTTGCGACCATTGGTCCGGAAGTa        <  1:1438466/150‑1 (MQ=255)
                                                                                                                                               ggTGAACGTTTTAAAGCCCAGCTCGTCATTTTTGCCTTCTTCATAACGATAGCTCGTATCGTTCCACAGCCAGACCGACATGTTCTGGTACAGACGTTCCAGATCGACGCTGCTCAGACGCTCACCTTTGCGACCATTGGTCCGGAAGTa        >  1:2500376/1‑150 (MQ=255)
                                                                                                                                                gTGAACGTTTTAAAGCCCAGCTCGTCATTTTTGCCTTCTTCATAACGATAGCTCGTATCGTTCCACAGCCAGACCGACATGTTCTGGTACAGACGTTCCAGATCGACGCTGCTCAGACGCTCACCTTTGCGACCATTGGTCCGGAAGTag       <  1:498758/150‑1 (MQ=255)
                                                                                                                                                 tGAACGTTTTAAAGCCCAGCTCGTCATTTTTGCCTTCTTCATAACGATAGCTCGTATCGTTCCACAGCCAGACCGACATGTTCTGGTACAGACGTTCCAGATCGACGCTGCTCAGACGCTCACCTTTGCGACCATTGGTCCGGAAGTaga      >  1:3381727/1‑150 (MQ=255)
                                                                                                                                                  gAAAGTTTTAAAGCCCAGCTCGTCATTTTTGCCTTCTTCATAACGATAGCTCGTATCGTTCCACAGCCAGACCGACATGTTCTGGTACAGACGTTCCAGATCGACGCTGCTCAGACGCTCACCTTTGCGACCATTGGTCCGGAAGTagag     >  2:1758719/1‑150 (MQ=255)
                                                                                                                                                   aaCGTTTTAAAGCCCAGCTCGTCATTTTTGCCTTCTTCATAACGATAGCTCGTATCGTTCCACAGCCAGACCGACATGTTCTGGTACAGACGTTCCAGATCGACGCTGCTCAGACGCTCACCTTTGCGACCATTGGTCCGGAAGTAGAGc    <  1:628224/150‑1 (MQ=255)
                                                                                                                                                    aCGTTTTAAAGCCCAGCTCGTCATTTTTGCCTTCTTCATAACGATAGCTCGTATCGTTCCACAGCCAGACCGACATGTTCTGGTACAGACGTTCCAGATCGACGCTGCTCAGACGCTCACCTTTGCGACCATTGGTCCGGAAGTAGAGCt   <  1:3712260/150‑1 (MQ=255)
                                                                                                                                                    aCGTTTTAAAGCCCAGCTCGTCATTTTTGCCTTCTTCATAACGATAGCTCGTATCGTTCCACAGCCAGACCGACATGTTCTGGTACAGACGTTCCAGATCGACGCTGCTCAGACGCTCACCTTTGCGACCATTGGTCCGGAAGTAGAGCt   <  1:4188505/150‑1 (MQ=255)
                                                                                                                                                     cGTTTTAAAGCCCAGCTCGTCATTTTTGCCTTCTTCATAACGATAGCTCGTATCGTTCCACAGCCAGACCGACATGTTCTGGTACAGACGTTCCAGATCGACGCTGCTCAGACGCTCACCTTTGCGACCATTGGTCCGGAAGTAGAGCTc  <  1:4538658/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
AGCGGGTAGCTCGGGTTCATCATGCCCGCTTTGCTGCTACCGTCACCGTAGATCATGTTGTTATCGACCAGCGATTTTTTCAGATCTGCAGAACACTTGGTGCCGCCTGCATAGGCATCATTGGCGTAGCAGTTCAGGAACTCGGTGAACGTTTTAAAGCCCAGCTCGTCATTTTTGCCTTCTTCATAACGATAGCTCGTATCGTTCCACAGCCAGACCGACATGTTCTGGTACAGACGTTCCAGATCGACGCTGCTCAGACGCTCACCTTTGCGACCATTGGTCCGGAAGTAGAGCTC  >  CP009273/3109395‑3109693

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: