Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,337,533 G→T 87.4% A18E (GCA→GAA)  osmB ← lipoprotein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,337,5330GT87.4% 728.1 / 105.5 285A18E (GCA→GAA) osmBlipoprotein
Reads supporting (aligned to +/- strand):  ref base G (15/21);  new base T (74/174);  total (90/195)
Fisher's exact test for biased strand distribution p-value = 1.79e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

Only 100 of 286 total aligned reads displayed.

CCAATAACACCTCCGACAGCTGCACCACCTAATGTACCCAACGTACTGCCATCGGTCAGTACTGCACCGCCTAATGCCCCTGCACCCGCGCCGATTGCGGTGTTGCGGTCCCGTTTAGACCAGTTAGAACAGGCACTCAGAGACATTGCCAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTGCTCGTTACAAACATAATACTCTCTCCTGAATTTATGATTCACGGAAGTAAGCTCTCTATAACTATAATAGCTAAGAATAAGTCTGGTGAAATTATTCCGTGAAATCTGCTC  >  CP009273/1337386‑1337682
                                                                                                                                                   |                                                                                                                                                     
ccAATAACACCTCCGACAGCTGCACCACCTAATGTACCCAACGTACTGCCATCGGTCAGTACTGCACCGCCTAATGCCCCTGCACCCGCGCCGATTGCGGTGTTGCGGTCCCGTTTAGACCAGTTAGAACAGGCACTCAGAGACATTTcc                                                                                                                                                     <  2:981698/150‑1 (MQ=255)
 cAATAACACCTCCGACAGCTGCACCACCTAATGTACCCAACGTACTGCCATCGGTCAGTACTGCACCGCCTAATGCCCCTGCACCCGCGCCGATTGCGGTGTTGCGGTCCCGTTTAGACCAGTTAGAACAGGCACTCAGAGACATTTCCa                                                                                                                                                    <  1:4058295/150‑1 (MQ=255)
  aaTAACACCTCCGACAGCTGCACCACCTAATGTACCCAACGTACTGCCATCGGTCAGTACTGCACCGCCTAATGCCCCTGCACCCGCGCCGATTGCGGTGTTGCGGTCCCGTTTAGACCAGTTAGAACAGGCACTCAGAGACATTTCCaa                                                                                                                                                   >  2:3661375/1‑150 (MQ=255)
   aTAACACCTCCGACAGCTGCACCACCTAATGTACCCAACGTACTGCCATCGGTCAGTACTGCACCGCCTAATGCCCCTGCACCCGCGCCGCTTGCGGTGTTGCGGTCCCGTTTAGACCAGTTAGAACAGGCACTCAGAGACATTGCCaaa                                                                                                                                                  <  1:1381657/150‑1 (MQ=255)
   aTAACACCTCCGACAGCTGCACCACCTAATGTACCCAACGTACTGCCATCGGTCAGTACTGCACCGCCTAATGCCCCTGCACCCGCGCCGATTGCGGTGTTGCGGTCCCGTTTAGACCAGTTAGAACAGGCACTCAGAGACATTTCCaaa                                                                                                                                                  >  2:2121592/1‑150 (MQ=255)
      acacCTCCGACAGCTGCACCACCTAATGTACCCAACGTACTGCCATCGGTCAGTACTGCACCGCCTAATGCCCCTGCACCCGCGCCGATTGCGGTGTTGCGGTCCCGTTTAGACCAGTTAGAACAGGCACTCAGAGACATTTCCAAAGTa                                                                                                                                               >  1:3305585/1‑150 (MQ=255)
         ccTCCGACAGCTGCACCACCTAATGTACCCAACGTACTGCCATCGGTCAGTACTGCACCGCCTAATGCCCCTGCACCCGCGCCGATTGCGGTGTTGCGGTCCCGTTTAGACCAGTTAGAACAGGCACTCAGAGACATTTCCAAAGTAAtt                                                                                                                                            >  2:1000050/1‑150 (MQ=255)
          cTCCGACAGCTGCACCACCTAATGTACCCAACGTACTGCCATCGGTCAGTACTGCACCGCCTAATGCCCCTGCACCCGCGCCGATTGCGGTGTTGCGGTCCCGTTTAGACCAGTTAGAACAGGCACTCAGAGACATTTCCAAAGTAATTg                                                                                                                                           >  2:1178350/1‑150 (MQ=255)
           tCCGACAGCTGCACCACCTAATGTACCCAACGTACTGCCATCGGTCAGTACTGCACCGCCTAATGCCCCTGCACCCGCGCCGATTGCGGTGTTGCGGTCCCGTTTAGACCAGTTAGAACAGGCACTCAGAGACATTTCCAAAGTAATTGc                                                                                                                                          >  1:2833044/1‑150 (MQ=255)
             cGACAGCTGCACCACCTAATGTACCCAACGTACTGCCATCGGTCAGTACTGCACCGCCTAATGCCCCTGCACCCGCGCCGATTGCGGTGTTGCGGTCCCGTTTAGACCAGTTAGAACAGGCACTCAGAGACATTTCCAAAGTAATTGCCa                                                                                                                                        <  2:684869/150‑1 (MQ=255)
                cAGCTGCACCACCTAATGTACCCAACGTACTGCCATCGGTCAGTACTGCACCGCCTAATGCCCCTGCACCCGCGCCGATTGCGGTGTTGCGGTCCCGTTTAGACCAGTTAGAACAGGCACTCAGAGACATTTCCAAAGTAATTGCCAGaa                                                                                                                                     >  2:4454719/1‑150 (MQ=255)
                 aGCTGCACCACCTAATGTACCCAACGTACTGCCATCGGTCAGTACTGCACCGCCTAATGCCCCTGCACCCGCGCCGATTGCGGTGTTGCGGTCCCGTTTAGACCAGTTAGAACAGGCACTCAGAGACATTTCCAAAGTAATTGCCAGAAc                                                                                                                                    <  2:1904487/150‑1 (MQ=255)
                      caccacCTAATGTACCCAACGTACTGCCATCGGTCAGTACTGCACCGCCTAATGCCCCTGCACCCGCGCCGATTGCGGTGTTGCGGTCCCGTTTAGACCAGTTAGAACAGGCACTCAGAGACATTTCCAAAGTAATTGCCAGAACAGCcg                                                                                                                               >  2:2168204/1‑150 (MQ=255)
                         caccTAATGTACCCAACGTACTGCCATCGGTCAGTACTGCACCGCCTAATGCCCCTGCACCCGCGCCGATTGCGGTGTTGCGGTCCCGTTTAGACCAGTTAGAACAGGCACTCAGAGACATTTCCAAAGTAATTGCCAGAACAGCCGCgg                                                                                                                            >  1:3747087/1‑150 (MQ=255)
                            cTAATGTACCCAACGTACTGCCATCGGTCAGTACTGCACCGCCTAATGCCCCTGCACCCGCGCCGATTGCGGTGTTGCGGTCCCGTTTAGACCAGTTAGAACAGGCACTCAGAGACATTTCCAAAGTAATTGCCAGAACAGCCGCGGTCa                                                                                                                         <  1:1543603/150‑1 (MQ=255)
                            cTAATGTACCCAACGTACTGCCATCGGTCAGTACTGCACCGCCTAATGCCCCTGCACCCGCGCCGATTGCGGTGTTGCGGTCCCGTTTAGACCAGTTAGAACAGGCACTCAGAGACATTTCCAAAGTAATTGCCAGAACAGCCGCGGTCa                                                                                                                         >  2:1541097/1‑150 (MQ=255)
                            cTAATGTACCCAACGTACTGCCATCGGTCAGTACTGCACCGCCTAATGCCCCTGCACCCGCGCCGATTGCGGTGTTGCGGTCCCGTTTAGACCAGTTAGAACAGGCACTCAGAGACATTTCCAAAGTAATTGCCAGAACAGCCGCGGTCa                                                                                                                         >  2:543442/1‑150 (MQ=255)
                            cTAATGTACCCAACGTACTGCCATCGGTCAGTACTGCACCGCCTAATGCCCCTGCACCCGCGCCGATTGCGGTGTTGCGGTCCCGTTTAGACCAGTTAGAACAGGCACTCAGAGACATTTCCAAAGTAATTGCCAGAACAGCCGCGGTCa                                                                                                                         <  2:3590028/150‑1 (MQ=255)
                              aaTGTACCCAACGTACTGCCATCGGTCAGTACTGCACCGCCTAATGCCCCTGCACCCGCGCCGATTGCGGTGTTGCGGTCCCGTTTAGACCAGTTAGAACAGGCACTCAGAGACATTTCCAAAGTAATTGCCAGAACAGCCGCGGTCAtt                                                                                                                       >  1:1889852/1‑150 (MQ=255)
                              aaTGTACCCAACGTACTGCCATCGGTCAGTACTGCACCGCCTAATGCCCCTGCACCCGCGCCGATTGCGGTGTTGCGGTCCCGTTTAGACCAGTTAGAACAGGCACTCAGAGACATTTCCAAAGTAATTGCCAGAACAGCCGCGGTCAtt                                                                                                                       <  2:2833352/150‑1 (MQ=255)
                               aTGTACCCAACGTACTGCCATCGGTCAGTACTGCACCGCCTAATGCCCCTGCACCCGCGCCGATTGCGGTGTTGCGGTCCCGTTTAGACCAGTTAGAACAGGCACTCAGAGACATTTCCAAATTAATTGCCAGAACAGCCGCGGTCAttt                                                                                                                      <  2:1437949/150‑1 (MQ=255)
                               aTGTACCCAACGTACTGCCATCGGTCAGTACTGCACCGCCTAATGCCCCTGCACCCGCGCCGATTGCGGTGTTGCGGTCCCGTTTAGACCAGTTAGAACAGGCACTCAGAGACATTTCCAAAGTAATTGCCAGAACAGCCGCGGTCAttt                                                                                                                      <  2:4082485/150‑1 (MQ=255)
                                tggaCCCAACGTACTGCCATCGGTCAGTACTGCACCGCCTAATGCCCCTGCACCCGCGCCGATTGCGGTGTTGCGGTCCCGTTTAGACCAGTTAGAACAGGCACTCAGAGACATTTCCAAAGTAATTGCCAGAACAGCCGCGGTCAtttt                                                                                                                     <  1:2365651/147‑1 (MQ=255)
                                tGTACCCAACGTACTGCCATCGGTCAGTACTGCACCGCCTAATGCCCCTGCACCCGCGCCGATTGCGGTGTTGCGGTCCCGTTTAGACCAGTTAGAACAGGCACTCAGAGACATTGCCAAAGTAATTGCCAGAACAGCCGCGGTCAtttt                                                                                                                     <  1:98277/150‑1 (MQ=255)
                                     ccAACGTACTGCCATCGGTCAGTACTGCACCGCCTAATGCCCCTGCACCCGCGCCGATTGCGGTGTTGCGGTCCCGTTTAGACCAGTTAGAACAGGCACTCAGAGACATTTCCAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTGc                                                                                                                <  2:2642205/150‑1 (MQ=255)
                                      cAACGTACTGCCATCGGTCAGTACTGCACCGCCTAATGCCCCTGCACCCGCGCCGATTGCGGTGTTGCGGTCCCGTTTAGACCAGTTAGAACAGGCACTCAGAGACATTTCCAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTGCt                                                                                                               >  1:4140043/1‑150 (MQ=255)
                                      cAACGTACTGCCATCGGTCAGTACTGCACCGCCTAATGCCCCTGCACCCGCGCCGATTGCGGTGTTGCGGTCCAGTTTAGACCAGTTAGAACAGGCACTCAGAGACATTGCCAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTGCt                                                                                                               <  1:583055/150‑1 (MQ=255)
                                       aaCGTACTGCCATCGGTCAGTACTGCACCGCCTAATGCCCCTGCACCCGCGCCGATTGCGGTGTTGCGGTCCCGTTTAGACCAGTTAGAACAGGCACTCAGAGACATTGCCAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTGCTc                                                                                                              >  1:1004323/1‑150 (MQ=255)
                                           tACTGCCATCGGTCAGTACTGCACCGCCTAATGCCCCTGCACCCGCGCCGATTGCGGTGTTGCGGTCCCGTTTAGACCAGTTAGAACAGGCACTCAGAGACATTTCCAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTGCTCGTTa                                                                                                          <  1:4168553/150‑1 (MQ=255)
                                           tACTGCCATCGGTCAGTACTGCACCGCCTAATGCCCCTGCACCCGCGCCGATTGCGGTGTTGCGGTCCCGTTTAGACCAGTTAGAACAGGCACTCAGAGACATTTCCAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTGCTCGGTa                                                                                                          <  2:3980746/150‑1 (MQ=255)
                                            aCTGCCATCGGTCAGTACTGCACCGCCTAATGCCCCTGCACCCGCGCCGATTGCGGTGTTGCGGTCCCGTTTAGACCAGTTAGAACAGGCACTCAGAGACATTGCCAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTGCTCGTTAc                                                                                                         >  2:4161591/1‑150 (MQ=255)
                                            aCTGCCATCGGTCAGTACTGCACCGCCTAATGCCCCTGCACCCGCGCCGATTGCGGTGTTGCGGTCCCGTTTAGACCAGTTAGAACAGGCACTCAGAGACATTGCCAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTGCTCGTTAc                                                                                                         >  2:1713449/1‑150 (MQ=255)
                                             cTGCCATCGGTCAGTACTGCACCGCCTAATGCCCCTGCACCCGCGCCGATTGCGGTGTTGCGGTCCCGTTTAGACCAGTTAGAACAGGCACTCAGAGACATTTCCAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTGCTCGTTACa                                                                                                        <  1:1796390/150‑1 (MQ=255)
                                                ccATCGGTCAGTACTGCACCGCCTAATGCCCCTGCACCCGCGCCGATTGCGGTGTTGCGGTCCCGTTTAGACCAGTTAGAACAGGCACTCAGAGACATTTCCAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTGCTCGTTACAAAc                                                                                                     >  1:1287925/1‑150 (MQ=255)
                                                ccATCGGTCAGTACTGCACCGCCTAATGCCCCTGCACCCGCGCCGATTGCGGTGTTGCGGTCCCGTTTAGACCAGTTAGAACAGGCACTCAGAGACATTTCCAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTGCTCGTTACAAAc                                                                                                     >  2:2482855/1‑150 (MQ=255)
                                                 cATCGGTCAGTACTGCACCGCCTAATGCCCCTGCACCCGCGCCGATTGCGGTGTTGCGGTCCCGTTTAGACCAGTTAGAACAGGCACTCAGAGACATTTCCAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTGCTCGTTACAAACa                                                                                                    <  2:3245692/150‑1 (MQ=255)
                                                 cATCGGTCAGTACTGCACCGCCTAATGCCCCTGCACCCGCGCCGATTGCGGTGTTGCGGTCCCGTTTAGACCAGTTAGAACAGGCACTCAGAGACATTGCCAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTGCTCGTTACAAACa                                                                                                    <  2:3132422/150‑1 (MQ=255)
                                                   tCGGTCAGTACTGCACCGCCTAATGCCCCTGCACCCGCGCCGATTGCGGTGTTGCGGTCCCGTTTAGACCAGTTAGAACAGGCACTCAGAGACATTTCCAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTGCTCGTTACAAACata                                                                                                  <  1:3686012/150‑1 (MQ=255)
                                                    cGGTCAGTACTGCACCGCCTAATGCCCCTGCACCCGCGCCGATTGCGGTGTTGCGGTCCCGTTTAGACCAGTTAGAACAGGCACTCAGAGACATTTCCAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTGCTCGTTACAAACataa                                                                                                 >  2:2191481/1‑150 (MQ=255)
                                                    cGGTCAGTACTGCACCGCCTAATGCCCCTGCACCCGCGCCGATTGCGGTGTTGCGGTCCCGTTTAGACCAGTTAGAACAGGCACTCAGAGACATTTCCAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTGCTCGTTACAAACataa                                                                                                 <  2:459719/150‑1 (MQ=255)
                                                    cGGTCAGTACTGCACCGCCTAATGCCCCTGCACCCGCGCCGATTGCGGTGTTGCGGTCCCGTTTAGACCAGTTAGAACAGGCACTCAGAGACATTTCCAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTGCTCGTTACAAACataa                                                                                                 <  1:116609/150‑1 (MQ=255)
                                                         aGTACTGCACCGCCTAATGCCCCTGCACCCGCGCCGATTGCGGTGTTGCGGTCCCGTTTAGACCAGTTAGAACAGGCACTCAGAGACATTTCCAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTGCTCGTTACAAACATAATActc                                                                                            <  1:807169/150‑1 (MQ=255)
                                                           tACTGCACCGCCTAATGCCCCTGCACCCGCGCCGATTGCGGTGTTGCGGTCCCGTTTAGACCAGTTAGAACAGGCACTCAGAGACATTTCCAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTGCTCGTTACAAACATAATActctc                                                                                          >  1:2103794/1‑150 (MQ=255)
                                                              tGCACCGCCTAATGCCCCTGCACCCGCGCCGATTGCGGTGTTGCGGTCCCGTTTAGACCAGTTAGAACAGGCACTCAGAGACATTTCCAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTGCTCGTTACAAACATAATACTCTCTcc                                                                                       <  2:887042/150‑1 (MQ=255)
                                                              tGCACCGCCTAATGCCCCTGCACCCGCGCCGATTGCGGTGTTGCGGTCCCGTTTAGACCAGTTAGAACAGGCACTCAGAGACATTTCCAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTGCTCGTTACAAACATAATACTCTCTcc                                                                                       <  1:4347898/150‑1 (MQ=255)
                                                               gCACCGCCTAATGCCCCTGCACCCGCGCCGATTGCGGTGTTGCGGTCCCGTTTAGACCAGTTAGAACAGGCACTCAGAGACATTTCCAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTGCTCGTTACAAACATAATACTCTCTCCt                                                                                      >  1:317312/1‑150 (MQ=255)
                                                               gCACCGCCTAATGCCCCTGCACCCGCGCCGATTGCGGTGTTGCGGTCCCGTTTAGACCAGTTAGAACAGGCACTCAGAGACATTGCCAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTGCTCGTTACAAACATAATACTCTCTCCt                                                                                      >  1:730935/1‑150 (MQ=255)
                                                                cACCGCCTAATGCCCCTGCACCCGCGCCGATTGCGGTGTTGCGGTCCCGTTTAGACCAGTTAGAACAGGCACTCAGAGACATTTCCAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTGCTCGTTACAAACATAATACTCTCTCCTg                                                                                     <  2:569399/150‑1 (MQ=255)
                                                                   cGCCTAATGCCCCTGCACCCGCGCCGATTGCGGTGTTGCGGTCCCGTTTAGACCAGTTAGAACAGGCACTCAGAGACATTGCCAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTGCTCGTTACAAACATAATACTCTCTCCTGAAt                                                                                  >  2:2994531/1‑150 (MQ=255)
                                                                   cGCCTAATGCCCCTGCACCCGCGCCGATTGCGGTGTTGCGGTCCCGTTTAGACCAGTTAGAACAGGCACTCAGAGACATTGCCAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTGCTCGTTACAAACATAATACTCTCTCCTGAAt                                                                                  >  2:1307919/1‑150 (MQ=255)
                                                                     ccTAATGCCCCTGCACCCGCGCCGATTGCGGTGTTGCGGTCCCGTTTAGACCAGTTAGAACAGGCACTCAGAGACATTTCCAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTGCTCGTTACAAACATAATACTCTCTCCTGAAttt                                                                                >  2:2060763/1‑150 (MQ=255)
                                                                     ccTAATGCCCCTGCACCCGCGCCGATTGCGGTGTTGCGGTCCCGTTTAGACCAGTTAGAACAGGCACTCAGAGACATTTCCAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTGCTCGTTACAAACATAATACTCTCTCCTGAAttt                                                                                <  1:413596/150‑1 (MQ=255)
                                                                        aaTGCCCCTGCACCCGCGCCGATTGCGGTGTTGCGGTCCCGTTTAGACCAGTTAGAACAGGCACTCAGAGACATTTCCAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTGCTCGTTACAAACATAATACTCTCTCCTGAATTTATg                                                                             >  1:1591830/1‑150 (MQ=255)
                                                                         aTGCCCCTGCACCCGCGCCGATTGCGGTGTTGCGGTCCCGTTTAGACCAGTTAGAACAGGCACTCAGAGACATTTCCAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTGCTCGTTACAAACATAATACTCTCTCCTGAATTTATGa                                                                            >  1:939528/1‑150 (MQ=255)
                                                                            ccccTGCACCCGCGCCGATTGCGGTGTTGCGGTCCCGTTTAGACCAGTTAGAACAGGCACTCAGAGACATTTCCAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTGCTCGTTACAAACATAATACTCTCTCCTGAATTTATGATTc                                                                         >  2:1288619/1‑150 (MQ=255)
                                                                            ccccTGCACCCGCGCCGATTGCGGTGTTGCGGTCCCGTTTAGACCAGTTAGAACAGGCACTCAGAGACATTTCCAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTGCTCGTTACAAACATAATACTCTCTCCTGAATTTATGATTc                                                                         >  1:599839/1‑150 (MQ=255)
                                                                             cccTGCACCCGCGCCGATTGCGGTGTTGCGGTCCCGTTTAGACCAGTTAGAACAGGCACTCAGAGACATTTCCAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTGCTCGTTACAAACATAATACTCTCTCCTGAATTTATGATTCa                                                                        <  2:3317278/150‑1 (MQ=255)
                                                                                 gCACCCGCGCCGATTGCGGTGTTGCGGTCCCGTTTAGACCAGTTAGAACAGGCACTCAGAGACATTTCCAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTGCTCGTTACAAACATAATACTCTCTCCTGAATTTATGATTCACGGa                                                                    <  2:2592264/150‑1 (MQ=255)
                                                                                  cACCCGCGCCGATTGCGGTGTTGCGGTCCCTTATAGACCAGTTAGAACAGGCACTCAGAGACATTTCCAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTGCTCGTTACAAACATAATACTCTCTCCTGAATTTATGATTCACGGaa                                                                   <  1:1594022/150‑1 (MQ=255)
                                                                                   aCCCGCGCCGATTGCGGTGTTGCGGTCCCGTTTAGACCAGTTAGAACAGGCACTCAGAGACATTTCCAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTGCTCGTTACAAACATAATACTCTCTCCTGAATTTATGATTCACGGAAg                                                                  <  1:1968162/150‑1 (MQ=255)
                                                                                     ccGCGCCGATTGCGGTGTTGCGGTCCCGTTTAGACCAGTTAGAACAGGCACTCAGAGACATTTCCAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTGCTCGTTACAAACATAATACTCTCTCCTGAATTTATGATTCACGGAAGTa                                                                <  1:4300964/150‑1 (MQ=255)
                                                                                      cgcgCCGATTGCGGTGTTGCGGTCCCGTTTAGACCAGTTAGAACAGGCACTCAGAGACATTTCCAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTGCTCGTTACAAACATAATACTCTCTCCTGAATTTATGATTCACGGAAGTaa                                                               <  1:3148490/150‑1 (MQ=255)
                                                                                      cgcgCCGATTGCGGTGTTGCGGTCCCGTTTAGACCAGTTAGAACAGGCACTCAGAGACATTTCCAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTGCTCGTTACAAACATAATACTCTCTCCTGAATTTATGATTCACGGAAGTaa                                                               <  1:3148535/150‑1 (MQ=255)
                                                                                      cgcgCCGATTGCGGTGTTGCGGTCCCGTTTAGACCAGTTAGAACAGGCACTCAGAGACATTTCCAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTGCTCGTTACAAACATAATACTCTCTCCTGAATTTATGATTCACGGAAGTaa                                                               <  1:95244/150‑1 (MQ=255)
                                                                                           cGATTGCGGTGTTGCGGTCCCGTTTAGACCAGTTAGAACAGGCACTCAGAGACATTTCCAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTGCTCGTTACAAACATAATACTCTCTCCTGAATTTATGATTCACGGAAGTAAGctct                                                          <  2:2321555/150‑1 (MQ=255)
                                                                                           cGATTGCGGTGTTGCGGTCCCGTTTAGACCAGTTAGAACAGGCACTCAGAGACATTTCCAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTGCTCGTTACAAACATAATACTCTCTCCTGAATTTATGATTCACGGAAGTAAGctct                                                          <  1:1837264/150‑1 (MQ=255)
                                                                                             aTTGCGGTGTTGCGGTCCCGTTTAGACCAGTTAGAACAGGCACTCAGAGACATTTCCAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTGCTCGTTACAAACATAATACTCTCTCCTGAATTTATGATTCACGGAAGTAAGctctct                                                        <  2:553760/150‑1 (MQ=255)
                                                                                              ttGCGGTGTTGCGGTCCCGTTTAGACCAGTTAGAACAGGCACTCAGAGACATTTCCAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTGCTCGTTACAAACATAATACTCTCTCCTGAATTTATGATTCACGGAAGTAAGCTCTcta                                                       <  1:2241903/150‑1 (MQ=255)
                                                                                              ttGCGGTGTTGCGGTCCCGTTTAGACCAGTTAGAACAGGCACTCAGAGACATTTCCAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTGCTCGTTACAAACATAATACTCTCTCCTGAATTTATGATTCACGGAAGTAAGCTCTcta                                                       <  1:4478395/150‑1 (MQ=255)
                                                                                                gCGGTGTTGCGGTCCCGTTTAGACCAGTTAGAACAGGCACTCAGAGACATTTCCAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTGCTCGTTACAAACATAATACTCTCTCCTGAATTTATGATTCACGGAAGTAAGCTCTctata                                                     <  2:2218877/150‑1 (MQ=255)
                                                                                                 cGGTGTTGCGGTCCCGTTTAGACCAGTTAGAACAGGCACTCAGAGACATTTCCAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTGCTCGTTACAAACATAATACTCTCTCCTGAATTTATGATTCACGGAAGTAAGCTCTctataa                                                    <  1:1535646/150‑1 (MQ=255)
                                                                                                  ggTGTTGCGGTCCCGTTTAGACCAGTTAGAACAGGCACTCAGAGACATTTCCAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTGCTCGTTACAAACATAATACTCTCTCCTGAATTTATGATTCACGGAAGTAAGCTCTctataac                                                   <  1:1375733/150‑1 (MQ=255)
                                                                                                  ggTGTTGCGGTCCCGTTTAGACCAGTTAGAACAGGCACTCAGAGACATTTCCAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTGCTCGTTACAAACATAATACTCTCTCCTGAATTTATGATTCACGGAAGTAAGCTCTctataac                                                   <  1:4021189/150‑1 (MQ=255)
                                                                                                        gCGGTCCCGTTTAGACCAGTTAGAACAGTCACTCAGAGACATTTCCAAAGTAAGTGCCAGAACAGCCGCGGTCATTTTTTTGCTCGTTACAAACATAATACTCTCTCCTGAATTTATGATTCACGGAAGTAAGCTCTCTATAACTataat                                             <  1:3682677/150‑1 (MQ=255)
                                                                                                          ggTCCCGTTTAGACCAGTTAGAACAGGCACTCAGAGACATTTCCAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTGCTCGTTACAAACATAATACTCTCTCCTGAATTTATGATTCACGGAAGTAAGCTCTCTATAACTATAATAg                                           <  2:623550/150‑1 (MQ=255)
                                                                                                          ggTCCCGTTTAGACCAGTTAGAACAGGCACTCAGAGACATTTCCAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTGCTCGTTACAAACATAATACTCTCTCCTGAATTTATGATTCACGGAAGTAAGCTCTCTATAACTATAATAg                                           <  2:1867417/150‑1 (MQ=255)
                                                                                                          ggTCCCGTTTAGACCAGTTAGAACAGGCACTCAGAGACATTGCCAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTGCTCGTTACAAACATAATACTCTCTCCTGAATTTATGATTCACGGAAGTAAGCTCTCTATAACTATAATAg                                           <  2:4223630/150‑1 (MQ=255)
                                                                                                               cGTTTAGACCAGTTAGAACAGGCACTCAGAGACATTTCCAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTGCTCGTTACAAACATAATACTCTCTCCTGAATTTATGATTCACGGAAGTAAGCTCTCTATAACTATAATAGCTAAg                                      <  1:1642934/150‑1 (MQ=255)
                                                                                                               cGTTTAGACCAGTTAGAACAGGCACTCAGAGACATTTCCAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTGCTCGTTACAAACATAATACTCTCTCCTGAATTTATGATTCACGGAAGTAAGCTCTCTATAACTATAATAGCTAAg                                      >  1:4542465/1‑150 (MQ=255)
                                                                                                                gTTTAGACCAGTTAGAACAGGAACTCAGAGACATTTACATAGTAATTGCGAGAACAGCCGCGGTCATTTTTTTGCTCGTTACAAACATAATACTCTCTCCTGAATTTATGATTCTCGGAAGTAAGCTCTCTATAACTATAATAGCTAAGa                                     <  1:3949968/150‑1 (MQ=255)
                                                                                                                   tAGACCAGTTAGAACAGGCACTCAGAGACATTTCCAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTGCTCGTTACAAACATAATACTCTCTCCTGAATTTATGATTCACGGAAGTAAGCTCTCTATAACTATAATAGCTAAGAATa                                  <  2:260171/150‑1 (MQ=255)
                                                                                                                   tAGACCAGTTAGAACAGGCACTCAGAGACATTTCCAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTGCTCGTTACAAACATAATACTCTCTCCTGAATTTATGATTCACGGAAGTAAGCTCTCTATAACTATAATAGCTAAGAATa                                  >  1:2603777/1‑150 (MQ=255)
                                                                                                                    aGACCAGTTAGAACAGGCACTCAGAGACATTTCCAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTGCTCGTTACAAACATAATACTCTCTCCTGAATTTATGATTCACGGAAGTAAGCTCTCTATAACTATAATAGCTAAGAATaa                                 <  2:774759/150‑1 (MQ=255)
                                                                                                                      aCCAGTTAGAACAGGCACTCAGAGACATTGCCAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTGCTCGTTACAAACATAATACTCTCTCCTGAATTTATGATTCACGGAAGTAAGCTCTCTATAACTATAATAGCTAAGAATAAGt                               >  2:1855197/1‑150 (MQ=255)
                                                                                                                        cAGTTAGAACAGGCACTCAGAGACATTTCCAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTGCTCGTTACAAACATAATACTCTCTCCTGAATTTATGATTCACGGAAGTAAGCTCTCTATAACTATAATAGCTAAGAATAAGTCt                             <  1:1100690/150‑1 (MQ=255)
                                                                                                                        cAGTTAGAACAGGCACTCAGAGACATTTCCAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTGCTCGTTACAAACATAATACTCTCTCCTGAATTTATGATTCACGGAAGTAAGCTCTCTATAACTATAATAGCTAAGAATAAGTCt                             <  1:344317/150‑1 (MQ=255)
                                                                                                                          gTTAGAACAGGCACTCAGAGACATTTCCAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTGCTCGTTACAAACATAATACTCTCTCCTGAATTTATGATTCACGGAAGTAAGCTCTCTATAACTATAATAGCTAAGAATAAGTCTgg                           <  1:175101/150‑1 (MQ=255)
                                                                                                                          gTTAGAACAGGCACTCAGAGACATTTCCAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTGCTCGTTACAAACATAATACTCTCTCCTGAATTTATGATTCACGGAAGTAAGCTCTCTATAACTATAATAGCTAAGAATAAGTCTgg                           >  1:250013/1‑150 (MQ=255)
                                                                                                                           ttAGAACAGGCACTCAGAGACATTTCCAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTGCTCGTTACAAACATAATACTCTCTCCTGAATTTATGATTCACGGAAGTAAGCTCTCTATAACTATAATAGCTAAGAATAAGTCTGGt                          <  1:2697209/150‑1 (MQ=255)
                                                                                                                              gAACAGGCACTCAGAGACATTTCCAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTGCTCGTTACAAACATAATACTCTCTCCTGAATTTATGATTCACGGAAGTAAGCTCTCTATAACTATAATAGCTAAGAATAAGTCTGGTGaa                       >  1:2309677/1‑150 (MQ=255)
                                                                                                                                   ggCACTCAGAGACATTGCCAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTGCTCGTTACAAACATAATACTCTCTCCTGAATTTATGATTCACGGAAGTAAGCTCTCTATAACTATAATAGCTAAGAATAAGTCTGGTGAAattat                  <  2:426804/150‑1 (MQ=255)
                                                                                                                                    gCACTCAGAGACATTTCCAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTGCTCGTTACAAACATAATACTCTCTCCTGAATTTATGATTCACGGAAGTAAGCTCTCTATAACTATAATAGCTAAGAATAAGTCTGGTGAAattatt                 <  1:2121592/150‑1 (MQ=255)
                                                                                                                                      aCTCAGAGACATTGCCAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTGCTCGTTACAAACATAATACTCTCTCCTGAATTTATGATTCACGGAAGTAAGCTCTCTATAACTATAATAGCTAAGAATAAGTCTGGTGAAATTATTcc               <  2:805663/150‑1 (MQ=255)
                                                                                                                                          agagACATTTCCAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTGCTCGTTACAAACATAATACTCTCTCCTGAATTTATGATTCACGGAAGTAAGCTCTCTATAACTATAATAGCTAAGAATAAGTCTGGTGAAATTATTCCGTGa           >  1:3009781/1‑150 (MQ=255)
                                                                                                                                            agaCATTTCCAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTGCTCGTTACAAACATAATACTCTCTCCTGAATTTATGATTCACGGAAGTAAGCTCTCTATAACTATAATAGCTAAGAATAAGTCTGGTGAAATTATTCCGTGaaa         <  1:1178350/150‑1 (MQ=255)
                                                                                                                                             gaCATTTCCAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTGCTCGTTACAAACATAATACTCTCTCCTGAATTTATGATTCACGGAAGTAAGCTCTCTATAACTATAATAGCTAAGAATAAGTCTGGTGAAATTATTCCGTGAAAt        >  1:2175938/1‑150 (MQ=255)
                                                                                                                                              aCATTTCCAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTGCTCGTTACAAACATAATACTCTCTCCTGAATTTATGATTCACGGAAGTAAGCTCTCTATAACTATAATAGCTAAGAATAAGTCTGGTGAAATTATTCCGTGAAATc       <  2:4128654/150‑1 (MQ=255)
                                                                                                                                               cATTTCCAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTGCTCGTTACAAACATAATACTCTCTCCTGAATTTATGATTCACGGAAGTAAGCTCTCTATAACTATAATAGCTAAGAATAAGTCTGGTGAAATTATTCCGTGAAATCt      <  2:827315/150‑1 (MQ=255)
                                                                                                                                                aTTGCCAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTGCTCGTTACAAACATAATACTCTCTCCTGAATTTATGATTCACGGAAGTAAGCTCTCTATAACTATAATAGCTAAGAATAAGTCTGGTGAAATTATTCCGTGAAATCTg     >  1:552354/1‑150 (MQ=255)
                                                                                                                                                aTTGCCAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTGCTCGTTACAAACATAATACTCTCTCCTGAATTTATGATTCACGGAAGTAAGCTCTCTATAACTATAATAGCTAAGAATAAGTCTGGTGAAATTATTCCGTGAAATCTg     >  1:552661/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
CCAATAACACCTCCGACAGCTGCACCACCTAATGTACCCAACGTACTGCCATCGGTCAGTACTGCACCGCCTAATGCCCCTGCACCCGCGCCGATTGCGGTGTTGCGGTCCCGTTTAGACCAGTTAGAACAGGCACTCAGAGACATTGCCAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTGCTCGTTACAAACATAATACTCTCTCCTGAATTTATGATTCACGGAAGTAAGCTCTCTATAACTATAATAGCTAAGAATAAGTCTGGTGAAATTATTCCGTGAAATCTGCTC  >  CP009273/1337386‑1337682

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: ‑