Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,106,685 G→T 88.7% 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,6850GT88.7% 762.3 / 92.9 299S51Y (TCT→TAT) yghGsecretin (GspDbeta) OM localization lipoprotein pilotin
Reads supporting (aligned to +/- strand):  ref base G (16/18);  new base T (120/145);  total (136/163)
Fisher's exact test for biased strand distribution p-value = 8.57e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 301 total aligned reads displayed.

TCACCGTCGGGTCAGCGCACATTTGTCGTTGATAGCTGGTTAAAAAGGCGTCAGGCGTCTGCGTGGTTTGTGTACCCGCTTCGCTGATAATGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCAGATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTTAATAATCCTGTAACGCTCAACAATAGCGATATTAATACCCTCCCTGGCATTTGTTTTATCG  >  CP009273/3106539‑3106832
                                                                                                                                                  |                                                                                                                                                   
tCACCGTCGGGTCAGCGCACATTTGTCGTTGATAGCTGGTTAAAAAGGCGTCAGGCGTCTGCGTGGTTTGTGTACCCGCTTCGCTGATAATGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCatat                                                                                                                                                   <  2:2448890/149‑1 (MQ=255)
tCACCGTCGGGTCAGCGCACATTTGTCGTTGATAGCTGGTTAAAAAGGCGTCAGGCGTCTGCGTGGTTTGTGTACCCGCTTCGCTGATAATGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCatat                                                                                                                                                   >  1:3590871/1‑149 (MQ=255)
 cACCGTCGGGTCAGCGCACATTTGTCGTTGATAGCTGGTTAAAAAGGCGTCAGGCGTCTGCGTGGTTTGTGTACCCGCTTCGCTGATAATGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCATAtt                                                                                                                                                  <  1:4437174/149‑1 (MQ=255)
 cACCGTCGGGTCAGCGCACATTTGTCGTTGATAGCTGGTTAAAAAGGCGTCAGGCGTCTGCGTGGTTTGTGTACCCGCTTCGCTGATAATGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCATAtt                                                                                                                                                  <  1:1633538/149‑1 (MQ=255)
 cACCGTCGGGTCAGCGCACATTTGTCGTTGATAGCTGGTTAAAAAGGCGTCAGGCGTCTGCGTGGTTTGTGTACCCGCTTCGCTGATAATGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCATAtt                                                                                                                                                  >  2:177829/1‑149 (MQ=255)
    cGTCGGGTCAGCGCACATTTGTCGTTGATAGCTGGTTAAAAAGGCGTCAGGCGTCTGCGTGGTTTGTGTACCCGCTTCGCTGATAATGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCATATTTaa                                                                                                                                               <  2:2313431/149‑1 (MQ=255)
      tCGGGTCAGCGCACATTTGTCGTTGATAGCTGGTTAAAAAGGCGTCAGGCGTCTGCGTGGTTTGTGTACCCGCTTCGCTGATAATGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCATATTTAATc                                                                                                                                             <  2:1537075/149‑1 (MQ=255)
      tCGGGTCAGCGCACATTTGTCGTTGATAGCTGGTTAAAAAGGCGTCAGGCGTCTGCGTGGTTTGTGTACCCGCTTCGCTGATAATGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCATATTTAATc                                                                                                                                             <  2:1613010/149‑1 (MQ=255)
        gggTCAGCGCACATTTGTCGTTGATAGCTGGTTAAAAAGGCGTCAGGCGTCTGCGTGGTTTGTGTACCCGCTTCGCTGATAATGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCATATTTAATCgg                                                                                                                                           <  2:792952/149‑1 (MQ=255)
        gggTCAGCGCACATTTGTCGTTGATAGCTGGTTAAAAAGGCGTCAGGCGTCTGCGTGGTTTGTGTACCCGCTTCGCTGATAATGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCATATTTAATCgg                                                                                                                                           >  1:3260729/1‑149 (MQ=255)
         ggTCAGCGCACATTTGTCGTTGATAGCTGGTTAAAAAGGCGTCAGGCGTCTGCGTGGTTTGTGTACCCGCTTCGCTGATAATGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCATATTTAATCGGc                                                                                                                                          <  2:3144537/149‑1 (MQ=255)
         ggTCAGCGCACATTTGTCGTTGATAGCTGGTTAAAAAGGCGTCAGGCGTCTGCGTGGTTTGTGTACCCGCTTCGCTGATAATGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCATATTTAATCGGc                                                                                                                                          >  2:3593039/1‑149 (MQ=255)
             aGCGCACATTTGTCGTTGATAGCTGGTTAAAAAGGCGTCAGGCGTCTGCGTGGTTTGTGTACCCGCTTCGCTGATAATGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCAGATTTAATCGGCAGGt                                                                                                                                      <  2:3866738/149‑1 (MQ=255)
              gcgcACATTTGTCGTTGATAGCTGGTTAAAAAGGCGTCAGGCGTCTGCGTGGTTTGTGTACCCGCTTCGCTGATAATGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCATATTTAATCGGCAGGtt                                                                                                                                     <  2:533498/149‑1 (MQ=255)
                          cGTTGATAGCTGGTTAAAAAGGCGTCAGGCGTCTGCGTGGTTTGTGTACCCGCTTCGCTGATAATGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCATATTTAATCGGCAGGTTTTGGCTGATAtt                                                                                                                         <  1:2617493/149‑1 (MQ=255)
                          cGTTGATAGCTGGTTAAAAAGGCGTCAGGCGTCTGCGTGGTTTGTGTACCCGCTTCGCTGATAATGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCATATTTAATCGGCAGGTTTTGGCTGATAtt                                                                                                                         <  2:2748231/149‑1 (MQ=255)
                            ttGATAGCTGGTTAAAAAGGCGTCAGGCGTCTGCGTGGTTTGTGGACCCGCTTCGCTGATACTGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCAGATTTAATCGGCAGGTTTTGGCTGATAtttt                                                                                                                       <  1:748052/149‑1 (MQ=255)
                             tGATAGCTGGTTAAAAAGGCGTCAGGCGTCTGCGTGGTTTGTGTACCCGCTTCGCTGATAATGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCATATTTAATCGGCAGGTTTTGGCTGATATTTTg                                                                                                                      <  1:4459265/149‑1 (MQ=255)
                                 aGCTGGTTAAAAAGGCGTCAGGCGTCTGCGTGGTTTGTGTACCCGCTTCGCTGATAATGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCATATTTAATCGGCAGGTTTTGGCTGATATTTTGCGcc                                                                                                                  >  1:1022338/1‑149 (MQ=255)
                                 aGCTGGTTAAAAAGGCGTCAGGCGTCTGCGTGGTTTGTGTACCCGCTTCGCTGATAATGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCATATTTAATCGGCAGGTTTTGGCTGATATTTTGCGcc                                                                                                                  >  2:334189/1‑149 (MQ=255)
                                 aGCTGGTTAAAAAGGCGTCAGGCGTCTGCGTGGTTTGTGTACCCGCTTCGCTGATAATGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCATATTTAATCGGCAGGTTTTGGCTGATATTTTGCGcc                                                                                                                  >  2:732794/1‑149 (MQ=255)
                                 aGCTGGTTAAAAAGGCGTCAGGCGTCTGCGTGGTTTGTGTACCCGCTTCGCTGATAATGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCATATTTAATCGGCAGGTTTTGGCTGATATTTTGCGcc                                                                                                                  >  1:2005026/1‑149 (MQ=255)
                                  gCTGGTTAAAAAGGCGTCAGGCGTCTGCGTGGTTTGTGTACCCGCTTCGCTGATAATGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCATATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCt                                                                                                                 <  1:1823537/149‑1 (MQ=255)
                                    tGGTTAAAAAGGCGTCAGGCGTCTGCGTGGTTTGTGTACCCGCTTCGCTGATAATGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCATATTTAATCGGCAGGTTTTGGCTGATAGTTTGCGCCTGt                                                                                                               >  1:2257624/1‑149 (MQ=255)
                                     ggTTAAAAAGGCGTCAGGCGTCTGCGTGGTTTGTGTACCCGCTTCGCTGATAATGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCATATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGtt                                                                                                              <  1:3841724/149‑1 (MQ=255)
                                      gTTAAAAAGGCGTCAGGCGTCTGCGTGGTTTGTGTACCCGCTTCGCTGATAATGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCATATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGttt                                                                                                             <  1:536535/149‑1 (MQ=255)
                                      gTTAAAAAGGCGTCAGGCGTCTGCGTGGTTTGTGTACCCGCTTCGCTGATAATGGTCATTTTTACCGTTGTGCCGCTACCTTGCGCCAGCACTAAGGTATATCCCGCATATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGttt                                                                                                             <  2:3937361/149‑1 (MQ=255)
                                       ttAAAAAGGCGTCAGGCGTCTGCGTGGTTTGTGTACCCGCTTCGCTGATAATGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCATATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGtttt                                                                                                            >  2:256873/1‑149 (MQ=255)
                                       ttAAAAAGGCGTCAGGCGTCTGCGTGGTTTGTGTACCCGCTTCGCTGATAATGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCATATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGtttt                                                                                                            >  2:1048547/1‑149 (MQ=255)
                                           aaaGGCGTCAGGCGTCTGCGTGGTTTGTGTACCCGCTTCGCAGATAATGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCATATGTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCg                                                                                                        >  2:4373632/1‑149 (MQ=255)
                                            aaGGCGTCAGGCGTCTGCGTGGTTTGTGTACCCGCTTCGCTGATAATGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCAGATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGc                                                                                                       <  1:420506/149‑1 (MQ=255)
                                             aggcgtcaggcgtcTGCGTGGTTTGTGTACCCGCTTCGCTGATAATGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCATATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGcc                                                                                                      >  2:4493808/1‑149 (MQ=255)
                                              ggcgtcaggcgtcTGCGTGGTTTGTGTACCCGCTTCGCTGATAATGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCATATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCa                                                                                                     <  1:4535395/149‑1 (MQ=255)
                                              ggcgtcaggcgtcTGCGTGGTTTGTGTACCCGCTTCGCTGATAATGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCATATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCa                                                                                                     <  2:4276790/149‑1 (MQ=255)
                                                     ggcgtcTGCGTGGTTTGTGTACCCGCTTCGCTGATAATGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCAGATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACt                                                                                              <  1:645238/149‑1 (MQ=255)
                                                      gcgtcTGCGTGGTTTGTGTACCCGCTTCGCTGATAATGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCATATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTg                                                                                             >  2:4716260/1‑149 (MQ=255)
                                                       cgtcTGCGTGGTTTGTGTACCCGCTTCGCTGATAATGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCATATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTgg                                                                                            <  2:3730316/149‑1 (MQ=255)
                                                        gtcTGCGTGGTTTGTGTACCCGCTTCGCTGATAATGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCATATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGc                                                                                           >  2:1216690/1‑149 (MQ=255)
                                                        gtcTGCGTGGTTTGTGTACCCGCTTCGCTGATAATGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCATATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGc                                                                                           <  2:3506996/149‑1 (MQ=255)
                                                        gtcTGCGTGGTTTGTGTACCCGCTTCGCTGATAATGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCATATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGc                                                                                           <  1:2316551/149‑1 (MQ=255)
                                                          cTGCGTGGTTTGTGTACCCGCTTCGCTGATAATGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCATATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCAt                                                                                         <  1:4634156/149‑1 (MQ=255)
                                                             cGTGGTTTGTGTACCCGCTTCGCTGATAATGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCATATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCAtttt                                                                                      <  2:4764729/149‑1 (MQ=255)
                                                                ggTTTGTGTACCCGCTTCGCTGATAATGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCATATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCAt                                                                                   >  2:399089/1‑149 (MQ=255)
                                                                ggTTTGTGGACCCTCTTCGTTGATTCTGGTCATTTTTAACGTTGTGCAGCTACTTTGCGCCAGCACTAATGTATATCACTCATTTTTATTCGGAAGGTTTTGGCTGATATTTTGCGCCTGTATTTTCTCCAGTAATCTGGCATTTTCAt                                                                                   <  2:2123178/149‑1 (MQ=255)
                                                                 gTTTGTGTACCCGCTTCGCTGATAATGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCAGATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCAtt                                                                                  >  2:2727913/1‑149 (MQ=255)
                                                                       gtACCCGCTTCGCTGATAATGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCATATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTAtggct                                                                            <  1:148883/149‑1 (MQ=255)
                                                                          cccGCTTCGCTGATAATGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCAGATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTAtggctggc                                                                         <  1:3336279/149‑1 (MQ=255)
                                                                            cGCTTCGCTGATAATGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCATATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGcac                                                                       <  2:2510527/149‑1 (MQ=255)
                                                                             gCTTCGCTGATAATGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCATATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGcaca                                                                      <  1:1940810/149‑1 (MQ=255)
                                                                              cTTCGCTGATAATGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCATATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAg                                                                     >  1:296274/1‑149 (MQ=255)
                                                                              cTTCGCTGATAATGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCAGATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAg                                                                     >  2:1401613/1‑149 (MQ=255)
                                                                                  gCTGATAATGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCATATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCAc                                                                 <  1:2857723/149‑1 (MQ=255)
                                                                                  gCTGATAATGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCATATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCAc                                                                 >  1:1401482/1‑149 (MQ=255)
                                                                                   cTGATAATGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCATATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACt                                                                >  1:2572840/1‑149 (MQ=255)
                                                                                   cTGATAATGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCATATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACt                                                                >  1:578189/1‑149 (MQ=255)
                                                                                     gATAATGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCATATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTta                                                              >  2:3651735/1‑149 (MQ=255)
                                                                                     gATAATGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCATATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTta                                                              >  2:2432585/1‑149 (MQ=255)
                                                                                      aTAATGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCATATTTAATCGTCAGGTTTTGGCTGATATTTTGCGCATGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTtaa                                                             >  2:19568/1‑149 (MQ=255)
                                                                                      aTAATGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCATATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTtaa                                                             >  2:4211942/1‑149 (MQ=255)
                                                                                       tAATGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCATATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATTGCTGGCACAGCCACTtaat                                                            >  2:1190093/1‑149 (MQ=255)
                                                                                       tAATGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCATATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTtaat                                                            <  1:3661516/149‑1 (MQ=255)
                                                                                        aaTGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCATATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTtaata                                                           <  1:1782862/149‑1 (MQ=255)
                                                                                         aTGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCATATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTtaataa                                                          >  2:1841607/1‑149 (MQ=255)
                                                                                          tGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCATATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTTaataat                                                         >  2:4546646/1‑149 (MQ=255)
                                                                                              cATTTGTACCGTTGTGCCGCTACTGTGCGCCAGCACTAAGGTATATCCCGCATATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTTAATAAGCCTg                                                     <  2:932773/149‑1 (MQ=255)
                                                                                               aTTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCAGATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTTAATAATCCTGt                                                    >  1:599920/1‑149 (MQ=255)
                                                                                                   ttACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCATATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTTAATAATCCTGTAACg                                                <  1:3887102/149‑1 (MQ=255)
                                                                                                    tACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCATATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTTAATAATCCTGTAACGc                                               >  1:1769538/1‑149 (MQ=255)
                                                                                                    tACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCATATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTTAATAATCCTGTAACGc                                               >  1:1186332/1‑149 (MQ=255)
                                                                                                    tACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCATATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTTAATAATCCTGTAACGc                                               >  1:4787403/1‑149 (MQ=255)
                                                                                                    tACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCAGATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTTAATAATCCTGTAACGc                                               >  2:2177178/1‑149 (MQ=255)
                                                                                                     aCCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCATATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTTAATAATCCTGTAACGCt                                              <  2:517996/149‑1 (MQ=255)
                                                                                                     aCCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCATATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTTAATAATCCTGTAACGCt                                              >  2:2829180/1‑149 (MQ=255)
                                                                                                       cGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCATATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTTAATAATCCTGTAACGCTca                                            <  1:875619/149‑1 (MQ=255)
                                                                                                       cGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCATATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTTAATAATCCTGTAACGCTca                                            <  1:2154493/149‑1 (MQ=255)
                                                                                                        gTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCATATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTTAATAATCCTGTAACGCTcaa                                           <  1:3979789/149‑1 (MQ=255)
                                                                                                         ttGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCATATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTTAATAATCCTGTAACGCTcaac                                          <  2:4596637/149‑1 (MQ=255)
                                                                                                            tgCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCATATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCGGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTTAATAATCCTGTAACGCTCAACAAt                                       <  2:3306896/149‑1 (MQ=255)
                                                                                                             gCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCATATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTTAATAATCCTGTAACGCTCAACAATa                                      <  2:3260729/149‑1 (MQ=255)
                                                                                                              ccGCTACTTTGCGCCAGCACTAAGGTATATCCCGCAGATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTTAATAATCCTGTAACGCTCAACAATAg                                     <  2:2603688/149‑1 (MQ=255)
                                                                                                                 cTACTTTGCGCCAGCACTAAGGTATATCCCGCAGATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTTAATAATCCTGTAACGCTCAACAATAGCGa                                  >  1:2429597/1‑149 (MQ=255)
                                                                                                                   aCTTTGCGCCAGCACTAAGGTATATCCCGCATATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTTAATAATCCTGTAACGCTCAACAATAGCGata                                >  2:3188465/1‑149 (MQ=255)
                                                                                                                       tGCGCCAGCACTAAGGTATATCCCGCATATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTTAATAATCCTGTAACGCTCAACAATAGCGATATTaa                            <  2:221400/149‑1 (MQ=255)
                                                                                                                       tGCGCCAGCACTAAGGTATATCCCGCAGATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTTAATAATCCTGTAACGCTCAACAATAGCGATATTaa                            <  1:3165034/149‑1 (MQ=255)
                                                                                                                        gcgcCAGCACTAAGGTATATCCCGCATATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTTAATAATCCTGTAACGCTCAACAATAGCGATATTAAt                           <  1:1105028/149‑1 (MQ=255)
                                                                                                                         cgcCAGCACTAAGGTATATCCCGCATATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTTAATAATCCTGTAACGCTCAACAATAGCGATATTAATa                          >  1:2831314/1‑149 (MQ=255)
                                                                                                                           ccAGCACTAAGGTATATCCCGCATATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTTAATAATCCTGTAACGCTCAACAATAGCGATATTAATAcc                        >  2:3073036/1‑149 (MQ=255)
                                                                                                                           ccAGCACTAAGGTATATCCCGCATATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTTAATAATCCTGTAACGCTCAACAATAGCGATATTAATAcc                        >  1:3889650/1‑149 (MQ=255)
                                                                                                                            cAGCACTAAGGTATATCCCGCATATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTTAATAATCCTGTAACGCTCAACAATAGCGATATTAATAccc                       <  2:4521926/149‑1 (MQ=255)
                                                                                                                              gCACTAAGGTATATCCCGCATATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTTAATAATCCTGTAACGCTCAACAATAGCGATATTAATAccctc                     <  1:2282597/149‑1 (MQ=255)
                                                                                                                              gCACTAAGGTATATCCCGCATATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTTAATAATCCTGTAACGCTCAACAATAGCGATATTAATAccctc                     <  1:988283/149‑1 (MQ=255)
                                                                                                                              gCACTAAGGTATATCCCGCATATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTTAATAATCCTGTAACGCTCAACAATAGCGATATTAATAccctc                     >  2:4648098/1‑149 (MQ=255)
                                                                                                                                   aaGGTATATCCCGCATATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTTAATAATCCTGTAACGCTCAACAATAGCGATATTAATACCCTCCCTgg                <  1:1665051/149‑1 (MQ=255)
                                                                                                                                       tataTCCCGCATATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTTAATAATCCTGTAACGCTCAACAATAGCGATATTAATACCCTCCCTGGCAtt            <  1:101717/149‑1 (MQ=255)
                                                                                                                                       tataTCCCGCATATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTTAATAATCCTGTAACGCTCAACAATAGCGATATTAATACCCTCCCTGGCAtt            >  2:607114/1‑149 (MQ=255)
                                                                                                                                       tataTCCCGCATATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTTAATAATCCTGTAACGCTCAACAATAGCGATATTAATACCCTCCCTGGCAtt            <  2:2913987/149‑1 (MQ=255)
                                                                                                                                       tataTCCCGCATATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTTAATAATCCTGTAACGCTCAACAATAGCGATATTAATACCCTCCCTGGCAtt            <  1:2494998/149‑1 (MQ=255)
                                                                                                                                           tCCCGCATATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTTAATAATCCTGTAACGCTCAACAATAGCGATATTAATACCCTCCCTGGCATTTGtt        >  2:563471/1‑149 (MQ=255)
                                                                                                                                              cGCATATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTTAATAATCCTGTAACGCTCAACAATAGCGATATTAATACCCTCCCTGGCATTTGTTTTa     >  1:328799/1‑149 (MQ=255)
                                                                                                                                                 aGATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTTAATAATCCTGTAACGCTCAACAATAGCGATATTAATACCCTCCCTGGCATTTGTTTTATCg  >  1:1367570/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
TCACCGTCGGGTCAGCGCACATTTGTCGTTGATAGCTGGTTAAAAAGGCGTCAGGCGTCTGCGTGGTTTGTGTACCCGCTTCGCTGATAATGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCAGATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTTAATAATCCTGTAACGCTCAACAATAGCGATATTAATACCCTCCCTGGCATTTGTTTTATCG  >  CP009273/3106539‑3106832

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: