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,6850GT96.2% 859.7 / 19.8 286S51Y (TCT→TAT) yghGsecretin (GspDbeta) OM localization lipoprotein pilotin
Reads supporting (aligned to +/- strand):  ref base G (7/4);  new base T (138/137);  total (145/141)
Fisher's exact test for biased strand distribution p-value = 5.41e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 289 total aligned reads displayed.

CACCGTCGGGTCAGCGCACATTTGTCGTTGATAGCTGGTTAAAAAGGCGTCAGGCGTCTGCGTGGTTTGTGTACCCGCTTCGCTGATAATGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCAGATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTTAATAATCCTGTAACGCTCAACAATAGCGATATTAATACCCTCCCTGGCATTTGTTTTA  >  CP009273/3106540‑3106829
                                                                                                                                                 |                                                                                                                                                
cACCGTCGGGTCAGCGCACATTTGTCGTTGATAGCTGGTTAAAAAGGCGTCAGGCGTCTGCGTGGTTTGTGTACCCGCTTCGCTGATAATGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCatat                                                                                                                                                >  1:983036/1‑148 (MQ=255)
    gTCGGGTCAGCGCACATTTGTCGTTGATAGCTGGTTAAAAAGGCGTCAGGCGTCTGCGTGGTTTGTGTACCCGCTTCGCTGATAATGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCATATTTaa                                                                                                                                            <  2:1564717/148‑1 (MQ=255)
            aGCGCACATTTGTCGTTGATAGCTGGTTAAAAAGGCGTCAGGCGTCTGCGTGGTTTGTGTACCCGCTTCGCTGATAATGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCATATTTAATCGGCAgg                                                                                                                                    <  2:2581238/148‑1 (MQ=255)
            aGCGCACATTTGTCGTTGATAGCTGGTTAAAAAGGCGTCAGGCGTCTGCGTGGTTTGTGTACCCGCTTCGCTGATAATGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCATATTTAATCGGCAgg                                                                                                                                    >  2:3347569/1‑148 (MQ=255)
             gcgcACATTTGTCGTTGATAGCTGGTTAAAAAGGCGTCAGGCGTCTGCGTGGTTTGTGTACCCGCTTCGCTGATAATGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCATATTTAATCGGCAGGt                                                                                                                                   >  1:349632/1‑148 (MQ=255)
                cacaTTTGTCGTTGATAGCTGGTTAAAAAGGCGTCAGGCGTCTGCGTGGTTTGTGTACCCGCTTCGCTGATAATGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCATATTTAATCGGCAGGtttt                                                                                                                                <  2:2349097/148‑1 (MQ=255)
                cacaTTTGTCGTTGATAGCTGGTTAAAAAGGCGTCAGGCGTCTGCGTGGTTTGTGTACCCGCTTCGCTGATAATGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCATATTTAATCGGCAGGtttt                                                                                                                                >  2:3017965/1‑148 (MQ=255)
                 acaTTTGTCGTTGATAGCTGGTTAAAAAGGCGTCAGGCGTCTGCGTGGTTTGTGTACCCGCTTCGCTGATAATGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCATATTTAATCGGCAGGTTTTg                                                                                                                               >  1:44522/1‑148 (MQ=255)
                 acaTTTGTCGTTGATAGCTGGTTAAAAAGGCGTCAGGCGTCTGCGTGGTTTGTGTACCCGCTTCGCTGATAATGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCATATTTAATCGGCAGGTTTTg                                                                                                                               >  1:1569621/1‑148 (MQ=255)
                  caTTTGTCGTTGATAGCTGGTTAAAAAGGCGTCAGGCGTCTGCGTGGTTTGTGTACCCGCTTCGCTGATAATGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCATATTTAATCGGCAGGTTTTgg                                                                                                                              >  2:3197541/1‑148 (MQ=255)
                  caTTTGTCGTTGATAGCTGGTTAAAAAGGCGTCAGGCGTCTGCGTGGTTTGTGTACCCGCTTCGCTGATAATGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCATATTTAATCGGCAGGTTTTgg                                                                                                                              <  1:561720/148‑1 (MQ=255)
                  caTTTGTCGTTGATAGCTGGTTAAAAAGGCGTCAGGCGTCTGCGTGGTTTGTGTACCCGCTTCGCTGATAATGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCATATTTAATCGGCAGGTTTTgg                                                                                                                              >  2:1864411/1‑148 (MQ=255)
                         cGTTGATAGCTGGTTAAAAAGGCGTCAGGCGTCTGCGTGGTTTGTGTACCCGCTTCGCTGATAATGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCATATTTAATCGGCAGGTTTTGGCTGatat                                                                                                                       <  1:3224828/148‑1 (MQ=255)
                          gTTGATCGCTGGTTAAAAAGGCGTCAGGCGTCTGCGTGGTTTGTGTACCCGCTTCGCTGATAATGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCATATTTAATCGGCAGGTTTTGGCTGATAtt                                                                                                                      <  1:3564466/148‑1 (MQ=255)
                           ttGATAGCTGGTTAAAAAGGCGTCAGGCGTCTGCGTGGTTTGTGTACCCGCTTCGCTGATAATGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCATATTTAATCGGCAGGTTTTGGCTGATAttt                                                                                                                     <  1:2358980/148‑1 (MQ=255)
                               tAGCTGGTTAAAAAGGCGTCAGGCGTCTGCGTGGTTTGTGTACCCGCTTCGCTGATAATGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCATATTTAATCGGCAGGTTTTGGCTGATATTTTgcg                                                                                                                 >  1:515018/1‑148 (MQ=255)
                                aGCTGGTTAAAAAGGCGTCAGGCGTCTGCGTGGTTTGTGTACCCGCTTCGCTGATAATGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCATATTTAATCGGCAGGTTTTGGCTGATATTTTgcgc                                                                                                                >  1:139575/1‑148 (MQ=255)
                                 gCTGGTTAAAAAGGCGTCAGGCGTCTGCGTGGTTTGTGTACCCGCTTCGCTGATAATGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCATATTTAATCGGCAGGTTTTGGCTGATATTTTGCGcc                                                                                                               >  2:2109377/1‑148 (MQ=255)
                                     gTTAAAAAGGCGTCAGGCGTCTGCGTGGTTTGTGTACCCGCTTCGCTGATAATGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCATATTTAATCGGCAGGTTTTGGCTGTTATTTTGCGCCTGtt                                                                                                           >  1:936921/1‑148 (MQ=255)
                                       tAAAAAGGCGTCAGGCGTCTGCGTGGTTTGTGTACCCGCTTCGCTGATAATGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCATATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGtttt                                                                                                         <  2:73433/148‑1 (MQ=255)
                                        aaaaaGGCGTCAGGCGTCTGCGTGGTTTGTGTACCCGCTTCGCTGATAATGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCATATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGttttt                                                                                                        <  1:2955228/148‑1 (MQ=255)
                                          aaaGGCGTCAGGCGTCTGCGTGGTTTGTGTACCCGCTTCGCTGATAATGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCATATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTc                                                                                                      >  2:1051711/1‑148 (MQ=255)
                                              gcgtcaggcgtcTGCGTGGTTTGTGTACCCGCTTCGCTGATAATGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCATATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCa                                                                                                  >  1:810790/1‑148 (MQ=255)
                                              gcgtcaggcgtcTGCGTGGTTTGTGTACCCGCTTCGCTGATAATGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCATATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCa                                                                                                  >  1:2597556/1‑148 (MQ=255)
                                              gcgtcaggcgtcTGCGTGGTTTGTGTACCCGCTTCGCTGATAATGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCATATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCa                                                                                                  >  1:3062848/1‑148 (MQ=255)
                                                 tcaggcgtcTGCGTGGTTTGTGTACCCGCTTCGCTGATAATGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCATATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTa                                                                                               >  1:1137212/1‑148 (MQ=255)
                                                  caggcgtcTGCGTGGTTTGTGTACCCGCTTCGCTGATAATGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCATATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTaa                                                                                              >  2:442430/1‑148 (MQ=255)
                                                  caggcgtcTGCGTGGTTTGTGTACCCGCTTCGCTGATAATGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCATATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTaa                                                                                              >  1:661832/1‑148 (MQ=255)
                                                    ggcgtcTGCGTGGGTTGTGTACCCGCTTCGCTGATAATGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCATATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAAc                                                                                            >  2:3437969/1‑148 (MQ=255)
                                                     gcgtcTGCGTGGTTTGTGTACCCGCTTCGCTGATAATGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCATATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACt                                                                                           <  1:1608205/148‑1 (MQ=255)
                                                     gcgtcTGCGTGGTTTGTGTACCCGCTTCGCTGATAATGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCATATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACt                                                                                           >  1:2826500/1‑148 (MQ=255)
                                                      cgtcTGCGTGGTTTGTGTACCCGCTTCGCTGATAATGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCATATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTg                                                                                          <  2:145136/148‑1 (MQ=255)
                                                      cgtcTGCGTGGTTTGTGTACCCGCTTCGCTGATAATGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCATATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTg                                                                                          >  1:3280021/1‑148 (MQ=255)
                                                         cTGCGTGGTTTGTGTACCCGCTTCGCTGATAATGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCATATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCa                                                                                       <  1:3098093/148‑1 (MQ=255)
                                                          tGCGTGGTGTGTGTACCCGCTTCGCTGATAATGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCATATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCAt                                                                                      <  1:1954994/148‑1 (MQ=255)
                                                              tGGTTTGTGTACCCGCTTCGCTGATAATGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCATATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTc                                                                                  >  2:2826199/1‑148 (MQ=255)
                                                              tGGTTTGTGTACCCGCTTCGCTGATAATGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCATATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTc                                                                                  <  2:2447285/148‑1 (MQ=255)
                                                               ggTTTGTGTACCCGCTTCGCTGATAATGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCATATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCa                                                                                 <  2:3039692/148‑1 (MQ=255)
                                                               ggTTTGTGTACCCGCTTCGCTGATAATGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCATATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCa                                                                                 <  1:2188947/148‑1 (MQ=255)
                                                                 tttGTGTACCCGCTTCGCTGATAATGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCATATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCAtt                                                                               <  1:683062/148‑1 (MQ=255)
                                                                 tttGTGTACCCGCTTCGCTGATAATGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCATATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCAtt                                                                               >  2:2150462/1‑148 (MQ=255)
                                                                  ttGTGTACCCGCTTCGCTGATAATGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCATATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTa                                                                              <  1:1017509/148‑1 (MQ=255)
                                                                    gtgtACCCGCTTCGCTGATAATGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCATATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTAtg                                                                            <  2:1492676/148‑1 (MQ=255)
                                                                    gtgtACCCGCTTCGCTGATAATGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCATATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTAtg                                                                            <  2:2384288/148‑1 (MQ=255)
                                                                     tgtACCCGCTTCGCTGATAATGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCATATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTAtgg                                                                           >  1:172703/1‑148 (MQ=255)
                                                                       tACCCGCTTCGCTGATAATGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCATATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTAtggct                                                                         >  1:2513421/1‑148 (MQ=255)
                                                                         cccGCTTCGCTGATAATGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCATATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTAtggctgg                                                                       >  1:2567680/1‑148 (MQ=255)
                                                                          ccGCTTCGCTGATAATGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCATATTTAATCGTCAGGTTTTGGCTGAAATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTAtggctggc                                                                      >  2:192989/1‑148 (MQ=255)
                                                                          ccGCTTCGCTGATAATGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCATATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTAtggctggc                                                                      >  2:578999/1‑148 (MQ=255)
                                                                           cGCTTCGCTGATAATGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCATATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGca                                                                     >  1:2052334/1‑148 (MQ=255)
                                                                            gCTTCGCTGATAATGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCATATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGcac                                                                    >  2:2242824/1‑148 (MQ=255)
                                                                             cTTCGCTGATAATGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCATATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGcaca                                                                   >  1:337718/1‑148 (MQ=255)
                                                                              ttCGCTGATAATGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCATATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAg                                                                  <  2:525024/148‑1 (MQ=255)
                                                                               tCGCTGATAATGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCATATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGc                                                                 <  2:765902/148‑1 (MQ=255)
                                                                               tCGCTGATAATGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCATATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGc                                                                 <  1:744526/148‑1 (MQ=255)
                                                                                cGCTGATAATGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCATATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGcc                                                                >  1:2097095/1‑148 (MQ=255)
                                                                                  cTGATAATGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCATATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCAc                                                              >  1:120824/1‑148 (MQ=255)
                                                                                  cTGATAATGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCATATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCAc                                                              >  1:1923016/1‑148 (MQ=255)
                                                                                     aTAATGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCATATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTta                                                           >  1:1063500/1‑148 (MQ=255)
                                                                                     aTAATGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCAGATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTta                                                           >  2:1427872/1‑148 (MQ=255)
                                                                                      tAATGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCATATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTtaa                                                          <  1:1292193/148‑1 (MQ=255)
                                                                                       aaTGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCATATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTtaat                                                         <  1:1640975/148‑1 (MQ=255)
                                                                                        aTGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCATATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTtaata                                                        >  1:329805/1‑148 (MQ=255)
                                                                                            tCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCATATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCATGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTTAATAATcc                                                    <  1:3391299/148‑1 (MQ=255)
                                                                                              aTTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCATATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTTAATAATCCTg                                                  <  2:3460046/148‑1 (MQ=255)
                                                                                                ttttACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCATATTTAATCGGCAGGTTTTGGCTGATCTTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTTAATAATCCTGTa                                                >  1:223224/1‑148 (MQ=255)
                                                                                                   tACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCATATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTTAATAATCCTGTAACg                                             <  2:676278/148‑1 (MQ=255)
                                                                                                   tACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCATATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTTAATAATCCTGTAACg                                             <  1:2015607/148‑1 (MQ=255)
                                                                                                     ccGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCATATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTTAATAATCCTGTAACGCt                                           <  1:125648/148‑1 (MQ=255)
                                                                                                     ccGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCATATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTTAATAATCCTGTAACGCt                                           >  1:1440788/1‑148 (MQ=255)
                                                                                                      cGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCATATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTTAATAATCCTGTAACGCTc                                          <  1:1528110/148‑1 (MQ=255)
                                                                                                      cGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCATATTTAATCGGCAGGTTTTGGCGGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTTAATAATCCTGTAACGCTc                                          >  1:414042/1‑148 (MQ=255)
                                                                                                       gTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCATATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTTAATAATCCTGTAACGCTca                                         <  2:2138839/148‑1 (MQ=255)
                                                                                                       gTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCATATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTTAATAATCCTGTAACGCTca                                         <  2:1801444/148‑1 (MQ=255)
                                                                                                        ttGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCATATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTTAATAATCCTGTAACGCTcaa                                        <  2:2455202/148‑1 (MQ=255)
                                                                                                          gtgCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCATATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTTAATAATCCTGTAACGCTcaaca                                      <  1:2678155/148‑1 (MQ=255)
                                                                                                           tgCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCATATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTTAATAATCCTGTAACGCTcaacaa                                     <  2:2517732/148‑1 (MQ=255)
                                                                                                             ccGCTACTTTGCGCCAGCACTAAGGTATATCCCGCATATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTTAATAATCCTGTAACGCTCAACAATa                                   <  2:2723338/148‑1 (MQ=255)
                                                                                                              cGCTACTTTGCGCCAGCACTAAGGTATATCCCGCATATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTTAATAATCCTGTAACGCTCAACAATAg                                  <  1:1790467/148‑1 (MQ=255)
                                                                                                                cTACTTTGCGCCAGCACTAAGGTATATCCCGCATATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTTAATAATCCTGTAACGCTCAACAATAGCg                                <  2:1706974/148‑1 (MQ=255)
                                                                                                                  aCTTTGCGCCAGCACTAAGGTATATCCCGCATATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTTAATAATCCTGTAACGCTCAACAATAGCGat                              >  1:277376/1‑148 (MQ=255)
                                                                                                                  aCTTTGCGCCAGCACTAAGGTATATCCCGCATATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTTAATAATCCTGTAACGCTCAACAATAGCGat                              >  1:697020/1‑148 (MQ=255)
                                                                                                                   cTTTGCGCCAGCACTAAGGTATATCCCGCATATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTTAATAATCCTGTAACGCTCAACAATAGCGata                             >  1:880240/1‑148 (MQ=255)
                                                                                                                     ttGCGCCAGCACTAAGGTATATCCCGCATATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTTAATAATCCTGTAACGCTCAACAATAGCGATAtt                           >  1:2298450/1‑148 (MQ=255)
                                                                                                                          ccAGCACTAAGGTATATCCCGCATATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTTAATAATCCTGTAACGCTCAACAATAGCGATATTAATAc                      <  2:2488496/148‑1 (MQ=255)
                                                                                                                              cACTAAGGTATATCCCGCATATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTTAATAATCCTGTAACGCTCAACAATAGCGATATTAATAccctc                  >  2:2076383/1‑148 (MQ=255)
                                                                                                                                 tAAGGTATATCCCGCATATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTTAATAATCCTGTAACGCTCAACAATAGCGATATTAATAccctccct               <  1:1662538/148‑1 (MQ=255)
                                                                                                                                 tAAGGTATATCCCGCATATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTTAATAATCCTGTAACGCTCAACAATAGCGATATTAATAccctccct               <  2:3290401/148‑1 (MQ=255)
                                                                                                                                  aaGGTATATCCCGCATATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTTAATAATCCTGTAACGCTCAACAATAGCGATATTAATACCCTCCCTg              >  1:2887285/1‑148 (MQ=255)
                                                                                                                                    ggTATATCCCGCATATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTTAATAATCCTGTAACGCTCAACAATAGCGATATTAATACCCTCCCTGGc            >  2:3225997/1‑148 (MQ=255)
                                                                                                                                      tataTCCCGCATATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTTAATAATCCTGTAACGCTCAACAATAGCGATATTAATACCCTCCCTGGCAt          <  2:1094550/148‑1 (MQ=255)
                                                                                                                                      tataTCCCGCATATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTTAATAATCCTGTAACGCTCAACAATAGCGATATTAATACCCTCCCTGGCAt          <  2:1150975/148‑1 (MQ=255)
                                                                                                                                       atatCCCGCATATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTTAATAATCCTGTAACGCTCAACAATAGCGATATTAATACCCTCCCTGGCAtt         <  1:2109377/148‑1 (MQ=255)
                                                                                                                                           cccGCATATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTTAATAATCCTGTAACGCTCAACAATAGCGATATTAATACCCTCCCTGGCATTTGtt     <  1:182436/148‑1 (MQ=255)
                                                                                                                                           cccGCATATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTTAATAATCCTGTAACGCTCAACAATAGCGATATTAATACCCTCCCTGGCATTTGtt     >  1:2124611/1‑148 (MQ=255)
                                                                                                                                           cccGCATATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTTAATAATCCTGTAACGCTCAACAATAGCGATATTAATACCCTCCCTGGCATTTGtt     <  2:120824/148‑1 (MQ=255)
                                                                                                                                             cGCATATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTTAATAATCCTGTAACGCTCAACAATAGCGATATTAATACCCTCCCTGGCATTTGtttt   <  1:3103611/148‑1 (MQ=255)
                                                                                                                                              gCATATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTTAATAATCCTGTAACGCTCAACAATAGCGATATTAATACCCTCCCTGGCATTTGTTTTa  >  2:2936734/1‑148 (MQ=255)
                                                                                                                                              gCATATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTTAATAATCCTGTAACGCTCAACAATAGCGATATTAATACCCTCCCTGGCATTTGTTTTa  >  1:589202/1‑148 (MQ=255)
                                                                                                                                              gCATATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTTAATAATCCTGTAACGCTCAACAATAGCGATATTAATACCCTCCCTGGCATTTGTTTTa  <  2:2614748/148‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                
CACCGTCGGGTCAGCGCACATTTGTCGTTGATAGCTGGTTAAAAAGGCGTCAGGCGTCTGCGTGGTTTGTGTACCCGCTTCGCTGATAATGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCAGATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTTAATAATCCTGTAACGCTCAACAATAGCGATATTAATACCCTCCCTGGCATTTGTTTTA  >  CP009273/3106540‑3106829

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: