Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,375,260 G→T 91.0% D387Y (GAT→TAT)  degQ → serine endoprotease, periplasmic

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,375,2600GT91.0% 989.7 / 84.0 338D387Y (GAT→TAT) degQserine endoprotease, periplasmic
Reads supporting (aligned to +/- strand):  ref base G (12/19);  new base T (161/146);  total (173/165)
Fisher's exact test for biased strand distribution p-value = 1.87e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 338 total aligned reads displayed.

GCACGAAAGTGAAGCTTGGCCTGCTGCGTAACGGCAAACCACTGGAAGTAGAAGTGACGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAAGATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAGGCTGGCTTGCAAAAAGACGATGTGATCATTGGCGTCAACCGCGATCGGGTGAACTCGATTGCTGAAATGCGTAAAGTGCTGGCGGCAAAA  >  CP009273/3375114‑3375406
                                                                                                                                                  |                                                                                                                                                  
gCACGAAAGTGAATCTTGGCCTGCTGCGTAACGGCAAACCACTGGAAGTAGAAGTGACGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAATATg                                                                                                                                                 <  2:277268/150‑1 (MQ=255)
  aCGAAAGTGAAGCTTGGCCTGCTGCGTAACGGCAAACCACTGGAAGTAGAAGTGACGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAATATggc                                                                                                                                               <  1:727539/150‑1 (MQ=255)
   cGAACGTGAAGCTTGGCCTGCTGCGTAACGGCAAACCACTGGAAGTAGAAGTGACGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAATATggcg                                                                                                                                              >  2:1544107/1‑150 (MQ=255)
     aaaGTGAAGCTTGGCCTGCTGCGTAACGGCAAACCACTGGAAGTAGAAGTGACGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAATATggcggc                                                                                                                                            >  1:2493754/1‑150 (MQ=255)
     aaaGTGAAGCTTGGCCTGCTGCGTAACGGCAAACCACTGGAAGTAGAAGTGACGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAATATggcggc                                                                                                                                            <  2:2109884/150‑1 (MQ=255)
          gAAGCTTGGCCTGCTGCGTAACGGCAAACCACTGGAAGTAGAAGTGACGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAATATGGCGGCAAAgg                                                                                                                                       <  1:3935348/150‑1 (MQ=255)
               ttGGCCTGCTGCGTAACGGCAAACCACTGGAAGTAGAAGTGACGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAATATGGCGGCAAAGGTATTa                                                                                                                                  >  2:4503495/1‑150 (MQ=255)
                tGGCCTGCTGCGTAACGGCAAACCACTGGAAGTAGAAGTGACGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAATATGGCGGCAAAGGTATTaa                                                                                                                                 >  1:2799074/1‑150 (MQ=255)
                tGGCCTGCTGCGTAACGGCAAACCACTGGAAGTAGAAGTGACGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAATATGGCGGCAAAGGTATTaa                                                                                                                                 >  1:4407676/1‑150 (MQ=255)
                 ggCCTGCTGCGTAACGGCAAACCACTGGAAGTAGAAGTGACGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAAGATGGCGGCAAAGGTATTaaa                                                                                                                                <  2:4044256/150‑1 (MQ=255)
                   ccTGCTGCGTAACGGCAAACCACTGGAAGTAGAAGTGACGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAATATGGCGGCACAGGTCTTAAAAt                                                                                                                              >  1:3592641/1‑150 (MQ=255)
                   ccTGCTGCGTAACGGCAAACCACTGGAAGTAGAAGTGACGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAATATGGCGGCAAAGGTATTAAAAt                                                                                                                              <  1:2001388/150‑1 (MQ=255)
                         gcGTAACGGCAAACCACTGGAAGTAGAAGTGACGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAATATGGCGGCAAAGGTATTAAAATCGATGa                                                                                                                        <  2:3180536/150‑1 (MQ=255)
                         gcGTAACGGCAAACCACTGGAAGTAGAAGTGACGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAATATGGCGGCAAAGGTATTAAAATCGATGa                                                                                                                        <  1:388524/150‑1 (MQ=255)
                          cGTAACGGCAAACCACTGGAAGTAGAAGTGACGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAAGATGGCGGCAAAGGTCTTAAAATCGATGaa                                                                                                                       >  1:549405/1‑150 (MQ=255)
                          cGTAACGGCAAACCACTGGAAGTAGAAGTGACGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAAGATGGCGGCAAAGGTATTAAAATCGATGaa                                                                                                                       >  1:2296388/1‑150 (MQ=255)
                            tAACGGCAAACCACTGGAAGTAGAAGTGACGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAATATGGCGGCAAAGGTATTAAAATCGATGAAGt                                                                                                                     >  1:2961471/1‑150 (MQ=255)
                             aaCGGCAAACCACTGGAAGTAGAAGTGACGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAATATGGCGGCAAAGGTATTAAAATCGATGAAGtt                                                                                                                    <  2:1832496/150‑1 (MQ=255)
                              aCGGCAAACCACTGGAAGTAGAAGTGACGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAATATGGCGGCAAAGGTATTAAAATCGATGAAGTTg                                                                                                                   <  2:3005913/150‑1 (MQ=255)
                                   aaaCCACTGGAAGTAGAAGTGACGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAATATGGCGGCAAAGGTATTGAAATCGATGAAGTTGTCaaa                                                                                                              >  2:518544/1‑150 (MQ=255)
                                    aaCCACTGGAAGTAGAAGTGACGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAATATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAg                                                                                                             <  2:338519/150‑1 (MQ=255)
                                     aCCACTGGAAGTAGAAGTGACGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAATATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAgg                                                                                                            >  2:38341/1‑150 (MQ=255)
                                       cACTGGAAGTAGAAGTGACGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAATATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGaa                                                                                                          >  1:1732753/1‑150 (MQ=255)
                                        aCTGGAAGTAGAAGTGACGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAATATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAg                                                                                                         >  2:1879252/1‑150 (MQ=255)
                                           ggAAGTAGAAGTGACGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAATATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGccc                                                                                                      <  2:80747/150‑1 (MQ=255)
                                           ggAAGTAGAAGTGACGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAAGATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGccc                                                                                                      <  1:4284957/150‑1 (MQ=255)
                                             aaGTAGAAGTGACGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAATATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAg                                                                                                    <  2:3480023/150‑1 (MQ=255)
                                               gTAGAAGTGACGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAATATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAgct                                                                                                  <  2:4433393/150‑1 (MQ=255)
                                               gTAGAAGTGACGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAAGATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAgct                                                                                                  <  1:790463/150‑1 (MQ=255)
                                                 aGAAGTGACGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAATATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAgctgc                                                                                                >  2:686856/1‑150 (MQ=255)
                                                   aaGTGACGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAATATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTc                                                                                              >  2:1310420/1‑150 (MQ=255)
                                                    aGTGACGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAATATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCa                                                                                             >  1:2183719/1‑150 (MQ=255)
                                                       gACGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAATATGGCGGCAAAGGTCTTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAggc                                                                                          >  2:2301786/1‑150 (MQ=255)
                                                       gACGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAATATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAggc                                                                                          <  1:4609995/150‑1 (MQ=255)
                                                         cGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAATATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAggctg                                                                                        <  2:1562905/150‑1 (MQ=255)
                                                         cGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAATATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAggctg                                                                                        >  2:5015382/1‑150 (MQ=255)
                                                          gCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAAGATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAggctgg                                                                                       <  2:3146146/150‑1 (MQ=255)
                                                             cGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGACGGTGCAACGTTGAGCGATGGTCAGCTAAAATATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAGGCTGGCtt                                                                                    >  1:1665534/1‑150 (MQ=255)
                                                             cGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAATATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAGGCTGGCtt                                                                                    >  1:2333012/1‑150 (MQ=255)
                                                             cGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAATATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAGGCTGGCtt                                                                                    >  1:4848000/1‑150 (MQ=255)
                                                             cGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAATATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAGGCTGGCtt                                                                                    >  1:3193847/1‑150 (MQ=255)
                                                               aTACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAATATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAGGCTGGCTTGc                                                                                  >  2:1029904/1‑150 (MQ=255)
                                                               aTACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAATATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAGGCTGGCTTGc                                                                                  >  2:4132610/1‑150 (MQ=255)
                                                                tACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAATATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAGGCTGGCTTGCa                                                                                 <  1:3802112/150‑1 (MQ=255)
                                                                 aCCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAATATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAGGCTGGCTTGCaa                                                                                <  2:1391187/150‑1 (MQ=255)
                                                                 aCCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAATATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAGGCTGGCTTGCaa                                                                                >  1:1918244/1‑150 (MQ=255)
                                                                    aGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAATATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAGGCTGGCTTGCaaaaa                                                                             >  2:4562023/1‑150 (MQ=255)
                                                                        ccTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAATATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAGGCTGGCTTGCAAAAAGACg                                                                         <  2:2004246/150‑1 (MQ=255)
                                                                         ctctTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAATATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAGGCTGGCTTGCAAAAAGACGa                                                                        >  2:1584702/1‑150 (MQ=255)
                                                                         ctctTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAATATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAGGCTGGCTTGCAAAAAGACGa                                                                        >  1:1644077/1‑150 (MQ=255)
                                                                             tcgtcgGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAATATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAGGCTGGCTTGCAAAAAGACGAtgtg                                                                    <  2:3326019/150‑1 (MQ=255)
                                                                                tcgGCCCGCGCTGAAATGATCACGCCCGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAACATCTGGCGGCAAAGGTAGTCAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAGGCTGGCTTGCAAAAAGACGATGTGATc                                                                 >  2:4905244/1‑150 (MQ=255)
                                                                                tcgGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAATATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAGGCTGGCTTGCACAAAGACGATGTGATc                                                                 >  2:2632736/1‑150 (MQ=255)
                                                                                tcgGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAATATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAGGCTGGCTTGCAAAAAGACGATGTGATc                                                                 >  1:3486761/1‑150 (MQ=255)
                                                                                 cgGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAATATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAGGCTGGCTTGCAAAAAGACGATGTGATCa                                                                <  1:2403363/150‑1 (MQ=255)
                                                                                   gCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAATATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAGGCTGGCTTGCAAAAAGACGATGTGATCAtt                                                              >  1:3290533/1‑150 (MQ=255)
                                                                                    ccAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAATATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAGGCTGGCTTGCAAAAAGACGATGTGATCATTg                                                             >  2:4217140/1‑150 (MQ=255)
                                                                                     cAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAATATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAGGCTGGCTTGCAAAAAGACGATGTGATCATTgg                                                            >  2:4776308/1‑150 (MQ=255)
                                                                                      aGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAATATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAGGCTGGCTTGCAAAAAGACGATGTGATCATTGGc                                                           >  1:1061/1‑150 (MQ=255)
                                                                                        cgcTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAATATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAGGCTGGCTTGCAAAAAGACGATGTGATCATTGGCGt                                                         <  1:4543402/150‑1 (MQ=255)
                                                                                         gcTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAAGATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAGGCTGGCTTGCAAAAAGACGATGTGATCATTGGCGTc                                                        >  1:2706198/1‑150 (MQ=255)
                                                                                             aaaTGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAATATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAGGCTGGCTTGCAAAAAGACGATGTGATCATTGGCGTCAAcc                                                    <  1:1560436/150‑1 (MQ=255)
                                                                                               aTGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAATATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAGGCTGGCTTGCAAAAAGACGATGTGATCATTGGCGTCAACcgc                                                  >  1:3694822/1‑150 (MQ=255)
                                                                                                 gATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAATATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAGGCTGGCTTGCAAAAAGACGATGTGATCATTGGCGTCAACCGCGa                                                <  1:5099804/150‑1 (MQ=255)
                                                                                                 gATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAATATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAGGCTGGCTTGCAAAAAGACGATGTGATCATTGGCGTCAACCGCGa                                                >  1:1406766/1‑150 (MQ=255)
                                                                                                  aTCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAATATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAGGCTGGCTTGCAAAAAGACGATGTGATCATTGGCGTCAACCGCGAt                                               >  1:4665117/1‑150 (MQ=255)
                                                                                                    cACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAATATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAGGCTGGCTTGCAAAAAGACGATGTGATCATTGGCGTCAACCGCGATCg                                             <  2:3910382/150‑1 (MQ=255)
                                                                                                    cACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAATATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAGGCTGGCTTGCAAAAAGACGATGTGATCATTGGCGTCAACCGCGATCg                                             <  1:1423113/150‑1 (MQ=255)
                                                                                                      cGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAATATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAGGCTGGCTTGCAAAAAGACGATGTGATCATTGGCGTCAACCGCGATCggg                                           >  2:2079129/1‑150 (MQ=255)
                                                                                                        ccAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAATATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAGGCTGGCTTGCAAAAAGACGATGTGATCATTGGCGTCAACCGCGATCGGGTg                                         >  1:898655/1‑150 (MQ=255)
                                                                                                        ccAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAATATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAGGCTGGCTTGCAAAAAGACGATGTGATCATTGGCGTCAACCGCGATCGGGTg                                         <  1:2434965/150‑1 (MQ=255)
                                                                                                         cAGCGCGGGAAGGTGCAACGTTGAGCGATGGTCAGCTCAAATATGGCGGCACAGGTATTCAAATCGAGGAAGTTGTCAAAGGAAGCCCAGCTGCTCAGGCTGGCTTGCAAAAAGACGATGTGATCATTGGCGTCAACCGCGATCGGGTGa                                        >  2:2362229/1‑150 (MQ=255)
                                                                                                            cgcTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAATATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAGGCTGGCTTGCAAAAAGACGATGTGATCATTGGCGTCAACCGCGATCGGGTGAACt                                     >  2:4366022/1‑150 (MQ=255)
                                                                                                            cgcTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAATATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAGGCTGGCTTGCAAAAAGACGATGTGATCATTGGCGTCAACCGCGATCGGGTGAACt                                     <  1:3114896/150‑1 (MQ=255)
                                                                                                            cgcTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAAGATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAGGCTGGCTTGCAAAAAGACGATGTGATCATTGGCGTCAACCGCGATCGGGTGAACt                                     <  2:1226699/150‑1 (MQ=255)
                                                                                                             gcTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAATATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAGGCTGGCTTGCAAAAAGACGATGTGATCATTGGCGTCAACCGCGATCGGGTGAACTc                                    <  2:2021705/150‑1 (MQ=255)
                                                                                                               tGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAATATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAGGCTGGCTTGCAAAAAGACGATGTGATCATTGGCGTCAACCGCGATCGGGTGAACTCGa                                  <  1:3711420/150‑1 (MQ=255)
                                                                                                                 gAAGGTGCAACGTTGAGCGATGGTCAGCTAAAATATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAGGCTGGCTTGCAAAAAGACGATGTGATCATTGGCGTCAACCGCGATCGGGTGAACTCGAtt                                <  1:3861995/150‑1 (MQ=255)
                                                                                                                 gAAGGTGCAACGTTGAGCGATGGTCAGCTAAAATATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAGGCTGGCTTGCAAAAAGACGATGTGATCATTGGCGTCAACCGCGATCGGGTGAACTCGAtt                                >  2:123587/1‑150 (MQ=255)
                                                                                                                 gAAGGTGCAACGTTGAGCGATGGTCAGCTAAAAGATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAGGCTGGCTTGCAAAAAGACGATGTGATCATTGGCGTCAACCGCGATCGGGTGAACTCGAtt                                <  1:4997655/150‑1 (MQ=255)
                                                                                                                  aaGGTGCAACGTTGAGCGATGGTCAGCTAAAATATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAGGCTGGCTTGCAAAAAGACGATGTGATCATTGGCGTCAACCGCGATCGGGTGAACTCGATTg                               >  1:67261/1‑150 (MQ=255)
                                                                                                                   aGGTGCAACGTTGAGCGATGGTCAGCTAAAATATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAGGCTGGCTTGCAAAAAGACGATGTGATCATTGGCGTCAACCGCGATCGGGTGAACTCGATTGc                              <  1:4374252/150‑1 (MQ=255)
                                                                                                                    ggTGCAACGTTGAGCGATGGTCAGCTAAAATATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAGGCTGGCTTGCAAAAAGACGATGTGATCATTGGCGTCAACCGCGATCGGGTGAACTCGATTGCt                             <  1:1840187/150‑1 (MQ=255)
                                                                                                                       gCAACGTTGAGCGATGGTCAGCTAAAATATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAGGCTGGCTTGCAAAAAGACGATGTGATCATTGGCGTCAACCGCGATCGGGTGAACTCGATTGCTGaa                          <  2:3816652/150‑1 (MQ=255)
                                                                                                                       gCAACGTTGAGCGATGGTCAGCTAAAATATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAGGCTGGCTTGCAAAAAGACGATGTGATCATTGGCGTCAACCGCGATCGGGTGAACTCGATTGCTGaa                          >  1:4185618/1‑150 (MQ=255)
                                                                                                                        cAACGTTGAGCGATGGTCAGCTAAAATATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAGGCTGGCTTGCAAAAAGACGATGTGATCATTGGCGTCAACCGCGATCGGGTGAACTCGATTGCTGaaa                         >  1:963288/1‑150 (MQ=255)
                                                                                                                        cAACGTTGAGCGATGGTCAGCTAAAATATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAGGCTGGCTTGCAAAAAGACGATGTGATCATTGGCGTCAACCGCGATCGGGTGAACTCGATTGCTGaaa                         <  1:2131928/150‑1 (MQ=255)
                                                                                                                         aaCGTTGAGCGATGGTCAGCTAAAATATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAGGCTGGCTTGCAAAAAGACGATGTGATCATTGGCGTCAACCGCGATCGGGTGAACTCGATTGCTGAAAt                        <  1:5044023/150‑1 (MQ=255)
                                                                                                                         aaCGTTGAGCGATGGTCAGCTAAAAGATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAGGCTGGCTTGCAAAAAGACGATGTGATCATTGGCGTCAACCGCGATCGGGTGAACTCGATTGCTGAAAt                        >  2:4668273/1‑150 (MQ=255)
                                                                                                                               gAGCGATGGTCAGCTAAAATATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAGGCTGGCTTGCAAAAAGACGATGTGATCATTGGCGTCAACCGCGATCGGGTGAACTCGATTGCTGAAATGCGTaa                  >  1:1871172/1‑150 (MQ=255)
                                                                                                                                    aTGGTCAGCTAAAATATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAGGCTGGCTTGCAAAAAGACGATGTGATCATTGGCGTCAACCGCGATCGGGTGAACTCGATTGCTGAAATGCGTAAAGTGc             >  1:1016220/1‑150 (MQ=255)
                                                                                                                                    aTGGTCAGCTAAAATATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAGGCTGGCTTGCAAAAAGACGATGTGATCATTGGCGTCAACCGCGATCGGGTGAACTCGATTGCTGAAATGCGTAAAGTGc             >  1:1095271/1‑150 (MQ=255)
                                                                                                                                    aTGGTCAGCTAAAAGATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAGGCTGGCTTGCAAAAAGACGATGTGATCATTGGCGTCAACCGCGATCGGGTGAACTCGATTGCTGAAATGCGTAAAGTGc             <  2:2879573/150‑1 (MQ=255)
                                                                                                                                      ggTCAGCTAAAATATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAGGCTGGCTTGCAAAAAGACGATGTGATCATTGGCGTCAACCGCGATCGGGTGAACTCGATTGCTGAAATGCGTAAAGTGCTg           >  2:2673930/1‑150 (MQ=255)
                                                                                                                                       gTCAGCTAAAATATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAGGCTGGCTTGCAAAAAGACGATGTGATCATTGGCGTCAACCGCGATCGGGTGAACTCGATTGCTGAAATGCGTAAAGTGCTgg          <  2:1558542/150‑1 (MQ=255)
                                                                                                                                       gTCAGCTAAAATATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAGGCTGGCTTGCAAAAAGACGATGTGATCATTGGCGTCAACCGCGATCGGGTGAACTCGATTGCTGAAATGCGTAAAGTGCTgg          >  1:1340552/1‑150 (MQ=255)
                                                                                                                                             tAAAATATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAGGCTGGCTTGCAAAAAGACGATGTGATCATTGGCGTCAACCGCGATCGGGTGAACTCGATTGCTGAAATGCGTAAAGTGCTGGCGGCaa    <  2:848333/150‑1 (MQ=255)
                                                                                                                                             tAAAAGATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAGGCTGGCTTCCAAAAAGACGATGTGATCATTGGCGTCAACCGCGATCGGGTGAACTCGATTGCTGAAATGCGTAAAGTGCTGGCGGCaa    >  1:5126873/1‑150 (MQ=255)
                                                                                                                                              aaaaTATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAGGCTGGCTTGCAAAAAGACGATGTGATCATTGGCGTCAACCGCGATCGGGTGAACTCGATTGCTGAAATGCGTAAAGTGCTGGCGGCaaa   <  1:634655/150‑1 (MQ=255)
                                                                                                                                              aaaaTATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAGGCTGGCTTGCAAAAAGACGATGTGATCATAGGCGTCAACCGCGATCGGGTGAACTCGATTGCTGAAATGCGTAAAGTGCTGGCGGCaaa   >  1:3941551/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
GCACGAAAGTGAAGCTTGGCCTGCTGCGTAACGGCAAACCACTGGAAGTAGAAGTGACGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAAGATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAGGCTGGCTTGCAAAAAGACGATGTGATCATTGGCGTCAACCGCGATCGGGTGAACTCGATTGCTGAAATGCGTAAAGTGCTGGCGGCAAAA  >  CP009273/3375114‑3375406

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: