Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,106,685 G→T 100% S51Y (TCT→TAT)  yghG ← secretin (GspDbeta) OM localization lipoprotein pilotin

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,106,6850GT98.2% 1181.7 / 6.2 377S51Y (TCT→TAT) yghGsecretin (GspDbeta) OM localization lipoprotein pilotin
Reads supporting (aligned to +/- strand):  ref base G (5/4);  new base T (185/183);  total (190/187)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 380 total aligned reads displayed.

TTCACCGTCGGGTCAGCGCACATTTGTCGTTGATAGCTGGTTAAAAAGGCGTCAGGCGTCTGCGTGGTTTGTGTACCCGCTTCGCTGATAATGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCAGATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTTAATAATCCTGTAACGCTCAACAATAGCGATATTAATACCCTCCCTGGCATTTGTTTTATC  >  CP009273/3106538‑3106831
                                                                                                                                                   |                                                                                                                                                  
ttCACCGTCGGGTCAGCGCACATTTGTCGTTGATAGCTGGTTAAAAAGGCGTCAGGCGTCTGCGTGGTTTGTGTACCCGCTTCGCTGATAATGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCatat                                                                                                                                                  <  1:4826312/150‑1 (MQ=255)
   aCCGTCGGGTCAGCGCACATTTGTCGTTGATAGCTGGTTAAAAAGGCGTCAGGCGTCTGCGTGGTTTGTGTACCCGCTTCGCTGATAATGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCATATTTa                                                                                                                                               <  1:117447/150‑1 (MQ=255)
    ccGTCGGGTCAGCGCACATTTGTCGTTGATAGCTGGTTAAAAAGGCGTCAGGCGTCTGCGTGGTTTGTGTACCCGCTTCGCTGATAATGGTCATTTTTACCGTTTTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCATATTTaa                                                                                                                                              >  1:3731819/1‑150 (MQ=255)
    ccGTCGGGTCAGCGCACATTTGTCGTTGATAGCTGGTTAAAAAGGCGTCAGGCGTCTGCGTGGTTTGTGTACCCGCTTCGCTGATAATGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCATATTTaa                                                                                                                                              <  2:280164/150‑1 (MQ=255)
     cGTCGGGTCAGCGCACATTTGTCGTTGATAGCTGGTTAAAAAGGCGTCAGGCGTCTGCGTGGTTTGTGTACCCGCTTCGCTGATAATGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCATATTTAAt                                                                                                                                             <  1:4238316/150‑1 (MQ=255)
      gTCGGGTCAGCGCACATTTGTCGTTGATAGCTGGTTAAAAAGGCGTCAGGCGTCTGCGTGGTTTGTGTACCCGCTTCGCTGATAATGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCATATTTAATc                                                                                                                                            >  1:5053552/1‑150 (MQ=255)
          ggTCAGCGCACATTTGTCGTTGATAGCTGGTTAAAAAGGCGTCAGGCGTCTGCGTGGTTTGTGTACCCGCTTCGCTGATAATGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCATATTTAATCGGCa                                                                                                                                        <  2:155577/150‑1 (MQ=255)
          ggTCAGCGCACATTTGTCGTTGATAGCTGGTTAAAAAGGCGTCAGGCGTCTGCGTGGTTTGTGTACCCGCTTCGCTGATAATGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCATATTTAATCGGCa                                                                                                                                        >  2:2026894/1‑150 (MQ=255)
            tCAGCGCACATTTGTCGTTGATAGCTGGTTAAAAAGGCGTCAGGCGTCTGCGTGGTTTGTGTACCCGCTTCGCTGATAATGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCATATTTAATCGGCAgg                                                                                                                                      >  1:3017895/1‑150 (MQ=255)
              aGCGCACATTTGTCGTTGATAGCTGGTTAAAAAGCCGACAGGCGTCTGCGTGTTTTGTGTACACGCTTCGCTGATAATGGTAATTTTTACAGTTGTGCCTCTACTTTGCGACAGCACTAAGGTATATCCAGCATATTTAATCGGCAGGtt                                                                                                                                    >  2:3486094/1‑150 (MQ=255)
               gcgcACATTTGTCGTTGATAGCTGGTTAAAAAGGCGTCAGGCGTCTGCGTGGTTTGTGTACCCGCTTCGCTGATAATGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCATATTTAATCGGCAGGttt                                                                                                                                   <  1:88217/150‑1 (MQ=255)
               gcgcACATTTGTCGTTGATAGCTGGTTAAAAAGGCGTCAGGCGTCTGCGTGGTTTGTGTACCCGCTTCGCTGATAATGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCATATTTAATCGGCAGGttt                                                                                                                                   <  2:6259/150‑1 (MQ=255)
                cgcACATTTGTCGTTGATAGCTGGTTAAAAAGGCGTCAGGCGTCTGCGTGGTTTGTGTACCCGCTTCGCTGATAATGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCATATTTAATCGGCAGGtttt                                                                                                                                  >  2:4163296/1‑150 (MQ=255)
                cgcACATTTGTCGTTGATAGCTGGTTAAAAAGGCGTCAGGCGTCTGCGTGGTTTGTGTACCCGCTTAGCTGATAATGGTCATTTTTACCGTTGTGACGCTACTTTTCGCCAGCACTAAGGTATATCCCGCATATTTAATAGGAATGtttt                                                                                                                                  >  1:4018742/1‑150 (MQ=255)
                  cacaTTTGTCGTTGATAGCTGGTTAAAAAGGCGTCAGGCGTCTGCGTGGTTTGTGTACCCGCTTCGCTGATAATGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTACGGTATATCCCGCATATTTAATCGGCAGGTTTTgg                                                                                                                                <  1:1855591/150‑1 (MQ=255)
                   acaTTTGTCGTTGATAGCTGGTTAAAAAGGCGTCAGGCGTCTGCGTGGTTTGTGTACCCGCTTCGCTGATAATGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCATATTTAATCGGCAGGTTTTGGc                                                                                                                               >  2:4099609/1‑150 (MQ=255)
                      tttGTCGTTGATAGCTGGTTAAAAAGGCGTCAGGCGTCTGCGTGGTTTGTGTACCCGCTTCGCTGATAATGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCATATTTAATCGGCAGGTTTTGGCTGa                                                                                                                            >  2:5284844/1‑150 (MQ=255)
                       ttGTCGTTGATCGCTGGTTAAAAAGGCGTCAGGCGTCTGCGTGGTTTGTGTACCCGCTTCGCTGATAATGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGTTATTTCCCGCATATTTAATCGGCAGGTTTTGGCTGat                                                                                                                           <  1:1402798/150‑1 (MQ=255)
                        tGTCGTTGATAGCTGGTTAAAAAGGCGTCAGGCGTCTGCGTGGTTTGTGTACCCGCTTCGCTGATAATGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCATATTTAATCGGCAGGTTTTGGCTGata                                                                                                                          >  2:3999952/1‑150 (MQ=255)
                            gTTGATAGCTGGTTAAAAAGGCGTCAGGCGTCTGCGTGGTTTGTGTACCCGCTTCGCTGATAATGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCATATTTAATCGGCAGGTTTTGGCTGATAtttt                                                                                                                      <  1:689087/150‑1 (MQ=255)
                            gTTGATAGCTGGTTAAAAAGGCGTCAGGCGTCTGCGTGGTTTGTGTACCCGCTTCGCTGATAATGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCATATTTAATCGGCAGGTTTTGGCTGATAtttt                                                                                                                      >  1:1022724/1‑150 (MQ=255)
                             ttGATAGCTGGTTAAAAAGGCGTCAGGCGTCTGCGTGGTTTGTGTACCCGCTTCGCTGATAATGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCATATTTAATCGGCAGGTTTTGGCTGATATTTTg                                                                                                                     >  2:2695224/1‑150 (MQ=255)
                                 tAGCTGGTTAAAAAGGCGTCAGGCGTCTGCGTGGTTTGTGTACCCGCTTCGCTGATAATGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCATATTTAATCGGCAGGTTTTGGCTGATATTTTGCGcc                                                                                                                 >  2:1244764/1‑150 (MQ=255)
                                  aGCTGGTTAAAAAGGCGTCAGGCGTCTGCGTGGTTTGTGTACCCGCTTCGCTGATAATGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCATATTTAATCGGAAGGTTTTGGCTGATATTTTGAGCCt                                                                                                                >  2:1878371/1‑150 (MQ=255)
                                    cTGGTTAAAAAGGCGTCAGGCGTCTGCGTGGTTTGTGTACCCGCTTCGCTGATAATGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCATATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGt                                                                                                              >  2:4966950/1‑150 (MQ=255)
                                     tGGTTAAAAAGGCGTCAGGCGTCTGCGTGGTTTGTGTACCCGCTTCGCTGATAATGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCATATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGtt                                                                                                             >  1:4408480/1‑150 (MQ=255)
                                      ggTTAAAAAGGCGTCAGGCGTCTGCGTGGTTTGTGTACCCGCTTCGCTGATAATGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCATATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGttt                                                                                                            <  2:3644356/150‑1 (MQ=255)
                                       gTTAAAAAGGCGTCAGGCGTCTGCGTGGTTTGTGTACCCGCTTCGCTGATAATGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCATATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGtttt                                                                                                           >  2:3900965/1‑150 (MQ=255)
                                        ttAAAAAGGCGTCAGGCGTCTGCGTGGTTTGTGTACCCGCTTCGCTGATAATGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCATATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGttttt                                                                                                          <  2:2247868/150‑1 (MQ=255)
                                               ggcgtcaggcgtcTGCGTGGTTTGTGTACCCGCTTCGCTGATAATGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCATATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAg                                                                                                   <  2:1697217/150‑1 (MQ=255)
                                                  gtcaggcgtcTGCGTGGTTTGTGTACCCGCTTCGCTGATAATGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCATATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTaa                                                                                                >  1:1273697/1‑150 (MQ=255)
                                                  gtcaggcgtcTGCGTGGTTTGTGTACCCGCTTCGCTGATAATGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCATATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTaa                                                                                                >  1:410041/1‑150 (MQ=255)
                                                      ggcgtcTGCGTGGTTTGTGTACCCGCTTCGCTGATAATGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCATATTTAATCGGCAGGGTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTg                                                                                            >  2:1295420/1‑150 (MQ=255)
                                                        cgtcTGCGTGGTTTGTGTACCCGCTTCGCTGATAATGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCATATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGc                                                                                          >  1:2900049/1‑150 (MQ=255)
                                                         gtcTGCGTGGTTTGTGTACCCGCTTCGCTGATAATGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCATATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCa                                                                                         >  1:2212171/1‑150 (MQ=255)
                                                            tGCGTGGTTTGTGTACCCGCTTCGCTGATAATGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCATATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCAttt                                                                                      >  1:3359787/1‑150 (MQ=255)
                                                              cGTGGTTTGTGTACCCGCTTCGCTGATAATGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCATATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTc                                                                                    >  2:2129257/1‑150 (MQ=255)
                                                              cGTGGTTTGTGAACCCGCTTCGCTGATAATGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCATATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTc                                                                                    >  1:2734575/1‑150 (MQ=255)
                                                               gTGGTTTGTGTACCCGCTTCGCTGATAATGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCATATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCa                                                                                   >  2:181794/1‑150 (MQ=255)
                                                               gTGGTTTGTGTACCCGCTTCGCTGATAATGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCATATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCa                                                                                   >  1:3519763/1‑150 (MQ=255)
                                                               gTGGTTTGTGTACCCGCTTCGCTGATAATGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCATATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCa                                                                                   >  1:3096766/1‑150 (MQ=255)
                                                                tGGTTTGTGTACCCGCTTCGCTGATAATGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCATATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCAt                                                                                  <  1:1482920/150‑1 (MQ=255)
                                                                   tttGTGTACCCGCTTCGCTGATAATGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCATATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTAt                                                                               <  1:3062958/150‑1 (MQ=255)
                                                                      gtgtACCCGCTTCGCTGATAATGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCATATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTAtggc                                                                            >  2:236827/1‑150 (MQ=255)
                                                                       tgtCCCCGCTTCGCTGATAATGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCATATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTAtggct                                                                           <  2:1104249/150‑1 (MQ=255)
                                                                         tACCCGCTTCGCTGATAATGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCATATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTAtggctgg                                                                         >  2:886146/1‑150 (MQ=255)
                                                                          aCCCGCTTCGCTGATAATGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCATATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTAtggctggc                                                                        >  2:4900718/1‑150 (MQ=255)
                                                                              gCTTCGCTGATAATGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCATATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAg                                                                    >  2:4350422/1‑150 (MQ=255)
                                                                              gCTTCGCTGATAATGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCATATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAg                                                                    >  1:5164413/1‑150 (MQ=255)
                                                                               cTTCGCTGATAATGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCATATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGc                                                                   >  2:2698211/1‑150 (MQ=255)
                                                                                  cGCTGATAATGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCAGATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCAc                                                                <  2:1021421/150‑1 (MQ=255)
                                                                                   gCTGATAATGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCATATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACt                                                               <  1:1344066/150‑1 (MQ=255)
                                                                                    cTGATAATGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTTTATCCCGCATTTTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACtt                                                              <  1:4348665/150‑1 (MQ=255)
                                                                                          aTGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCATATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTTaataat                                                        >  2:4668982/1‑150 (MQ=255)
                                                                                             gTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCATATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTTAATAATCCt                                                     >  2:736423/1‑150 (MQ=255)
                                                                                              tCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCATATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTTAATAATCCTg                                                    >  1:4692841/1‑150 (MQ=255)
                                                                                               cATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCATATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTTAATAATCCTGt                                                   <  2:3149312/150‑1 (MQ=255)
                                                                                                    ttACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCATATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTTAATAATCCTGTAACGc                                              <  1:2653141/150‑1 (MQ=255)
                                                                                                    ttACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCATATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTTAATAATCCTGTAACGc                                              >  2:5097339/1‑150 (MQ=255)
                                                                                                        cGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCATATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTTAATAATCCTGTAACGCTcaa                                          >  2:2982510/1‑150 (MQ=255)
                                                                                                          ttGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCATATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTTAATAATCCTGTAACGCTcaaca                                        <  1:4899644/150‑1 (MQ=255)
                                                                                                            gtgCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCATATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTTAATAATCCTGTAACGCTCAACAAt                                      <  2:288069/150‑1 (MQ=255)
                                                                                                               ccGCTACTTTGCGCCAGCACTAAGGTATATCCCGCATATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTTAATAATCCTGTAACGCTCAACAATAGc                                   <  2:1056919/150‑1 (MQ=255)
                                                                                                               ccGCTACTTTGCGCCAGCACTAAGGTATATCCCGCATATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTTAATAATCCTGTAACGCTCAACAATAGc                                   >  2:3605584/1‑150 (MQ=255)
                                                                                                               ccGCTACTTTGCGCCAGCACTAAGGTATATCCCGCATATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTTAATAATCCTGTAACGCTCAACAATAGc                                   <  1:3667753/150‑1 (MQ=255)
                                                                                                               ccGCTACTTTGCGCCAGCACTAAGGTATATCCCGCATATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTTAATAATCCTGTAACGCTCAACAATAGc                                   <  2:3731819/150‑1 (MQ=255)
                                                                                                               ccGCTACTTTGCGCCAGCACTAAGGTATATCCCGCATATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTTAATAATCCTGTAACGCTCAACAATAGc                                   <  1:2021530/150‑1 (MQ=255)
                                                                                                                cGCTACTTTGCGCCAGCACTAAGGTATATCCCGCATATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACCGCCACTTAATAATCCTGTAACGCTCAACAATAGCg                                  >  2:2334995/1‑150 (MQ=255)
                                                                                                                 gCTACTTTGCGCCAGCACTAAGGTATATCCCGCATATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTTAATAATCCTGTAACGCTCAACAATAGCGa                                 >  2:2535303/1‑150 (MQ=255)
                                                                                                                   tACTTTGCGCCAGCACTAAGGTATATCCCGCATATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTTAATAATCCTGTAACGCTCAACAATAGCGata                               <  1:4460959/150‑1 (MQ=255)
                                                                                                                   tACTTTGCGCCAGCACTAAGGTATATCCCGCATATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTTAATAATCCTGTAACGCTCAACAATAGCGata                               >  2:1417481/1‑150 (MQ=255)
                                                                                                                   tACTTTGCGCCAGCACTAAGGTATATCCCGCATATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTTAATAATCCTGTAACGCTCAACAATAGCGata                               <  2:3359787/150‑1 (MQ=255)
                                                                                                                    aCTTTGCGCCAGCACTAAGGTATATCCCGCATATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTTAATAATCCTGTAACGCTCAACAATAGCGatat                              <  1:4932917/150‑1 (MQ=255)
                                                                                                                     cTTTGCGCCAGCACTAAGGTATATCCCGCATATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCGACTTAATAATCCTGTAACGCTCAACAATAGCGATAtt                             >  2:1340591/1‑150 (MQ=255)
                                                                                                                        tGCGCCAGCACTAAGGTATATCCCGCATATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTTAATAATCCTGTAACGCTCAACAATAGCGATATTAAt                          <  2:3849669/150‑1 (MQ=255)
                                                                                                                        tGCGCCAGCACTAAGGTATATCCCGCATATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTTAATAATCCTGTAACGCTCAACAATAGCGATATTAAt                          <  2:3237155/150‑1 (MQ=255)
                                                                                                                              aGCACTAAGGTATATCCCGCATATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTTAATAATCCTGTAACGCTCAACAATAGCGATATTAATAccctc                    <  1:5284844/150‑1 (MQ=255)
                                                                                                                                 aCTAAGGTATATCCCGCATATTTTATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTTAATAATCCTGTAACGCTCAACAATAGCGATATTAATAccctccct                 <  1:1097388/150‑1 (MQ=255)
                                                                                                                                 aCTAAGGTATATCCCGCATATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTTAATAATCCTGTAACGCTCAACAATAGCGATATTAATAccctccct                 <  1:3218629/150‑1 (MQ=255)
                                                                                                                                 aCTAAGGTATATCCCGCATATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTTAATAATCCTGTAACGCTCAACAATAGCGATATTAATAccctccct                 <  2:3432575/150‑1 (MQ=255)
                                                                                                                                 aCTAAGGTAGACCCCGCATATTTCATCGGCAGGTGTGGGCTGATATTGTGCGCCTGTTGTTTCGCCAGTCACCTGGCATTTTCAGTATGGCTGGCCCAGCCACTTAATAATCCTGTAACGCTCAACAATAGCGATATTAATAccctccct                 <  2:4834894/150‑1 (MQ=255)
                                                                                                                                  cTAAGGTATATCCCGCATATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTTAATAATCCTGTAACGCTCAACAATAGCGATATTAATACCCTCCCTg                >  2:1506432/1‑150 (MQ=255)
                                                                                                                                  cTAAGGTATATCCCGCATATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTTAATAATCCTGTAACGCTCAACAATAGCGATATTAATACCCTCCCTg                >  2:1958820/1‑150 (MQ=255)
                                                                                                                                    aaGGTATATCCCGCATATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTTAATAATCCTGTAACGCTCAACAATAGCGATATTAATACCCTCCCTGGc              <  1:3380731/150‑1 (MQ=255)
                                                                                                                                       gTATATCCCGCATATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTTAATAATCCTGTAACGCTCAACAATAGCGATATTAATACCCTCCCTGGCAtt           <  1:161860/150‑1 (MQ=255)
                                                                                                                                        tataTCCCGCATATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTTAATAATCCTGTAACGCTCAACAATAGCGATATTAATACCCTCCCTGGCAttt          >  1:2480327/1‑150 (MQ=255)
                                                                                                                                        tataTCCCGCATATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTTAATAATCCTGTAACGCTCAACAATAGCGATATTAATACCCTCCCTGGCAttt          >  1:2376321/1‑150 (MQ=255)
                                                                                                                                        tataTCCCGCATATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTTAATAATCCTGTAACGCTCAACAATAGCGATATTAATACCCTCCCTGGCAttt          <  1:2308906/150‑1 (MQ=255)
                                                                                                                                         atatCCCGCATATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTTAATAATCCTGTAACGCTCAACAATAGCGATATTAATACCCTCCCTGGCATTTg         <  1:3860353/150‑1 (MQ=255)
                                                                                                                                          tatCCCGCATATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTTAATAATCCTGTAACGCTCAACAATAGCGATATTAATACCCTCCCTGGCATTTGt        <  1:4282013/150‑1 (MQ=255)
                                                                                                                                           atCCCGCATATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTTAATAATCCTGTAACGCTCAACAATAGCGATATTAATACCCTCCCTGGCATTTGtt       >  2:4496485/1‑150 (MQ=255)
                                                                                                                                           atCCCGCATATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTTAATAATCCTGTAACGCTCAACAATAGCGATATTAATACCCTCCCTGGCATTTGtt       <  1:4163296/150‑1 (MQ=255)
                                                                                                                                           atCCCGCATATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTTAATAATCCCGTAACGCTCAACAATAGCGATATTAATACCCTCCCTGGCATTTGtt       <  2:4039042/150‑1 (MQ=255)
                                                                                                                                            tCCCGCATATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTTAATAATCCTGTAACGCTCAACAATAGCGATATTAATACCCTCCCTGGCATTTGttt      >  1:714260/1‑150 (MQ=255)
                                                                                                                                            tCCCGCATATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTTAATAATCCTGTAACGCTCAACAATAGCGATATTAATACCCTCCCTGGCATTTGttt      >  1:4067333/1‑150 (MQ=255)
                                                                                                                                             cccGCATATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTTAATAATCCTGTAACGCTCAACAATAGCGATATTAATACCCTCCCTGGCATTTGtttt     <  1:1687359/150‑1 (MQ=255)
                                                                                                                                             cccGCATATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTTAATAATCCTGTAACGCTCAACAATAGCGATATTAATACCCTCCCTGGCATTTGtttt     <  1:1738096/150‑1 (MQ=255)
                                                                                                                                               cGCATATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTTAATAATCCTGTAACGCTCAACAATAGCGATATTAATACCCTCCCTGGCATTTGTTTTAt   >  1:2444537/1‑150 (MQ=255)
                                                                                                                                                gCATATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTTAATAATCCTGTAACGCTCAACAATAGCGATATTAATACCCTCCCTGGCATTTGTTTTATc  <  1:3976128/150‑1 (MQ=255)
                                                                                                                                                gCATATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTTAATAATCCTGTAACGCTCAACAATAGCGATATTAATACCCTCCCTGGCATTTGTTTTATc  <  2:2629312/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
TTCACCGTCGGGTCAGCGCACATTTGTCGTTGATAGCTGGTTAAAAAGGCGTCAGGCGTCTGCGTGGTTTGTGTACCCGCTTCGCTGATAATGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCAGATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTTAATAATCCTGTAACGCTCAACAATAGCGATATTAATACCCTCCCTGGCATTTGTTTTATC  >  CP009273/3106538‑3106831

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: