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

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,337,5330GT92.1% 807.9 / 57.3 266A18E (GCA→GAA) osmBlipoprotein
Reads supporting (aligned to +/- strand):  ref base G (12/9);  new base T (107/138);  total (119/147)
Fisher's exact test for biased strand distribution p-value = 2.59e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.74e-01

Only 100 of 266 total aligned reads displayed.

CCAATAACACCTCCGACAGCTGCACCACCTAATGTACCCAACGTACTGCCATCGGTCAGTACTGCACCGCCTAATGCCCCTGCACCCGCGCCGATTGCGGTGTTGCGGTCCCGTTTAGACCAGTTAGAACAGGCACTCAGAGACATTGCCAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTGCTCGTTACAAACATAATACTCTCTCCTGAATTTATGATTCACGGAAGTAAGCTCTCTATAACTATAATAGCTAAGAATAAGTCTGGTGAAATTATTCCGTGAAATCTG  >  CP009273/1337386‑1337679
                                                                                                                                                   |                                                                                                                                                  
ccAATAACACCTCCGACAGCTGCACCACCTAATGTACCCAACGTACTGCCATCGGTCAGTACTGCACCGCCTAATGCCCCTGCACCCGCGCCGATTGCGGTGTTGCGGTCCCGTTTAGACCAGTTAGAACAGGCACTCAGAGACATTTcc                                                                                                                                                  <  1:4169735/150‑1 (MQ=255)
ccAATAACACCTCCGACAGCTGCACCACCTAATGTACCCAACGTACTGCCATCGGTCAGTACTGCACCGCCTAATGCCCCTGCACCCGCGCCGATTGCGGTGTTGCGGTCCCGTTTAGACCAGTTAGAACAGGCACTCAGAGACATTTcc                                                                                                                                                  <  2:2885731/150‑1 (MQ=255)
 cAATAACACCTCCGACAGCTGCACCACCTAATGTACCCAACGTACTGCCATCGGTCAGTACTGCACCGCCTAATGCCCCTGCACCCGCGCCGATTGCGGTGTTGCGGTCCCGTTTAGACCAGTTAGAACAGGCACTCAGAGACATTTCCa                                                                                                                                                 <  2:711786/150‑1 (MQ=255)
 cAATAACACCTCCGACAGCTGCACCACCTAATGTACCCAACGTACTGCCATCGGTCAGTACTGCACCGCCTAATGCCCCTGCACCCGCGCCGATTGCGGTGTTGCGGTCCCGTTTAGACCAGTTAGAACAGGCACTCAGAGACATTTCCa                                                                                                                                                 <  2:2067732/150‑1 (MQ=255)
 cAATAACACCTCCGACAGCTGCACCACCTAATGTACCCAACGTACTGCCATCGGTCAGTACTGCACCGCCTAATGCCCCTGCACCCGCGCCGATTGCGGTGTTGCGGTCCCGTTTAGACCAGTTAGAACAGGCACTCAGAGACATTTCCa                                                                                                                                                 >  2:2157091/1‑150 (MQ=255)
        acCTCCGACAGCTGCACCACCTACTGTACCCAACGTACTGCCATCGGTCAGTACTGCACCGCCTAATGCCCCTGCACCCGCGCCGATTGCGGTGTTGCGGTCCCGTTTAGACCAGTTAGAACAGGCACTCAGAGACATTTCCAAAGTAAt                                                                                                                                          <  2:127637/150‑1 (MQ=255)
        acCTCCGACAGCTGCACCACCTAATGTACCCAACGTACTGCCATCGGTCAGTACTGCACCGCCTAATGCCCCTGCACCCGCGCCGATTGCGGTGTTGCGGTCCCGTTTAGACCAGTTAGAACAGGCACTCAGAGACATTTCCAAAGTAAt                                                                                                                                          <  2:2956115/150‑1 (MQ=255)
         ccTCCGACAGCTGCACCACCTAATGTACCCAACGTACTGCCATCGGTCAGTACTGCACCGCCTAATGCCCCTGCACCCGCGCCGATTGCGGTGTTGCGGTCCCGTTTAGACCAGTTAGAACAGGCACTCAGAGACATTTCCAAAGTAAtt                                                                                                                                         <  1:1931685/150‑1 (MQ=255)
          cTCCGACAGCTGCACCACCTAATGTACCCAACGTACTGCCATCGGTCAGTACTGCACCGCCTAATGCCCCTGCACCCGCGCCGATTGCGGTGTTGCGGTCCCGTTTAGACCAGTTAGAACAGGCACTCAGAGACATTTCCAAAGTAATTg                                                                                                                                        >  1:5126702/1‑150 (MQ=255)
          cTCCGACAGCTGCACCACCTAATGTACCCAACGTACTGCCATCGGTCAGTACTGCACCGCCTAATGCCCCTGCACCCGCGCCGATTGCGGTGTTGCGGTCCCGTTTAGACCAGTTAGAACAGGCACTCAGAGACATTGCCAAAGTAATTg                                                                                                                                        >  1:216017/1‑150 (MQ=255)
             cGACAGCTGCACCACCTAATGTACCCAACGTACTGCCATCGGTCAGTACTGCACCGCCTAATGCCCCTGCACCCGCGCCGATTGCGGTGTTGCGGTCCCGTTTAGACCAGTTAGAACAGGCACTCAGAGACATTTCCAAAGTAATTGCCa                                                                                                                                     <  2:773873/150‑1 (MQ=255)
             cGACAGCTGCACCACCTAATGTACCCAACGTACTGCCATCGGTCAGTACTGCACCGCCTAATGCCCCTGCACCCGCGCCGATTGCGGTGTTGCGGTCCCGTTTAGACCAGTTAGAACAGGCACTCAGAGACATTTCCAAAGTAATTGCCa                                                                                                                                     >  2:1573286/1‑150 (MQ=255)
              gACAGCTGCACCACCTAATGTACCCAACGTACTGCCATCGGTCAGTACTGCACCGCCTAATGCCCCTGCACCCGCGCCGATTGCGGTGTTGCGGTCCCGTTTAGACCAGTTAGAACAGGCACTCAGAGACATTTCCAAAGTAATTGCCAg                                                                                                                                    >  2:2734657/1‑150 (MQ=255)
              gACAGCTGCACCACCTAATGTACCCAACGTACTGCCATCGGTCAGTACTGCACCGCCTAATGCCCCTGCACCCGCGCCGATTGCGGTGTTGCGGTCCCGTTTAGACCAGTTAGAACAGGCACTCAGAGACATTTCCAAAGTAATTGCCAg                                                                                                                                    >  1:2428045/1‑150 (MQ=255)
                   cTGCACCACCTAATGTACCCAACGTACTGCCATCGGTCAGTACTGCACCGCCTAATGCCCCTGCACCCGCGCCGATTGCGGTGTTGCGGTCCCGTTTAGACCAGTTAGAACAGGCACTCAGAGACATTTCCAAAGTAATTGCCAGAACAg                                                                                                                               >  2:1169177/1‑150 (MQ=255)
                     gCACCACCTAATGTACCCAACGTACTGCCATCGGTCAGTACTGCACCGCCTAATGCCCCTGCACCCGCGCCGATTGCGGTGTTGCGGTCCCGTTTAGACCAGTTAGAACAGGCACTCAGAGACATTTCCAAAGTAATTGCCAGAACAGcc                                                                                                                             <  2:4984266/150‑1 (MQ=255)
                     gCACCACCTAATGTACCCAACGTACTGCCATCGGTCAGTACTGCACCGCCTAATGCCCCTGCACCCGCGCCGATTGCGGTGTTGCGGTCCCGTTTAGACCAGTTAGAACAGGCACTCAGAGACATTGCCAAAGTAATTGCCAGAACAGcc                                                                                                                             <  1:2280460/150‑1 (MQ=255)
                        ccaccTAATGTACCCAACGTACTGCCATCGGTCAGTACTGCACCGCCTAATGCCCCTGCACCCGCGCCGATTGCGGTGTTGCGGTCCCGTTTAGACCAGTTAGAACAGGCACTCAGAGACATTGCCAAAGTAATTGCCAGAACAGCcgcg                                                                                                                          >  1:1554690/1‑150 (MQ=255)
                         caccTAATGTACCCAACGTACTGCCATCGGTCAGTACTGCACCGCCTAATGCCCCTGCACCCGCGCCGATTGCGGTGTTGCGGTCCCGTTTAGACCAGTTAGAACAGGCACTCAGAGACATTTCCAAAGTACTTGCCAGAACAGCCGCgg                                                                                                                         >  2:4191516/1‑150 (MQ=255)
                           ccTAATGTACCCAACGTACTGCCATCGGTCAGTACTGCACCGCCTAATGCCCCTGCACCCGCGCCGATTGCGGTGTTGCGGTCCCGTTTAGACCAGTTAGAACAGGCACTCAGAGACATTTCCAAAGTAATTGCCAGAACAGCCGCGGTc                                                                                                                       <  2:3970061/150‑1 (MQ=255)
                            cTAATGTACCCAACGTACTGCCATCGGTCAGTACTGCACCGCCTAATGCCCCTGCACCCGCGCCGATTGCGGTGTTGCGGTCCCGTTTAGACCAGTTAGAACAGGCACTCAGAGACATTTCCAAAGTAATTGCCAGAACAGCCGCGGTCa                                                                                                                      <  1:4773552/150‑1 (MQ=255)
                            cTAATGTACCCAACGTACTGCCATCGGTCAGTACTGCACCGCCTAATGCCCCTGCACCCGCGCCGATTGCGGTGTTGCGGTCCCGTTTAGACCAGTTAGAACAGGCACTCAGAGACATTTCCAAAGTAATTGCCAGAACAGCCGCGGTCa                                                                                                                      <  1:4006856/150‑1 (MQ=255)
                               aTGTACCCAACGTACTGCCATCGGTCAGTACTGCACCGCCTAATGCCCCTGCACCCGCGCCGATTGCGGTGTTGCGGTCCCGTTTAGACCAGTTAGAACAGGCACTCAGAGACATTTCCAAAGTAATTGCCAGAACAGCCGCGGTCAttt                                                                                                                   >  2:2435894/1‑150 (MQ=255)
                               aTGTACCCAACGTACTGCCATCGGTCAGTACTGCACCGCCTAATGCCCCTGCACCCGCGCCGATTGCGGTGTTGCGGTCCCGTTTAGACCAGTTAGAACAGGCACTCAGAGACATTTCCAAAGTAATTGCCAGAACAGCCGCGGTCAttt                                                                                                                   >  1:1177139/1‑150 (MQ=255)
                                 ggaCCCAACGTACTGCCATCGGTCAGTACTGCACCGCCTAATGCCCCTGCACCCGCGCCGATTGCGGTGTTGCGGTCCCGTTTAGACCAGTTAGAACAGGCACTCAGAGACATTTCCAAAGTAATTGCCAGAACAGCCGCGGTCAttttt                                                                                                                 <  2:3778210/148‑1 (MQ=255)
                                 gTACCCAACGTACTGCCATCGGTCAGTACTGCACCGCCTAATGCCCCTGCACCCGCGCCGATTGCGGTGTTGCGGTCCCGTTTAGACCAGTTAGAACAGGCACTCAGAGACATTTCCAAAGTAATTGCCAGAACAGCCGCGGTCAttttt                                                                                                                 <  1:769524/150‑1 (MQ=255)
                                 gTACCCAACGTACTGCCATCGGTCAGTACTGCACCGCCTAATGCCCCTGCACCCGCGCCGATTGCGGTGTTGCGGTCCCGTTTAGACCAGTTAGAACAGGCACTCAGAGACATTTCCAAAGTAATTGCCAGAACAGCCGCGGTCAttttt                                                                                                                 >  1:471435/1‑150 (MQ=255)
                                    cccAACGTACTGCCATCGGTCAGTACTGCACCGCCTAATGCCCCTGCACCCGCGCCGATTGCGGTGTTGCGGTCCCGTTTAGACCAGTTAGAACAGGCACTCAGAGACATTTCCAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTg                                                                                                              <  1:2953458/150‑1 (MQ=255)
                                    cccAACGTACTGCCATCGGTCAGTACTGCACCGCCTAATGCCCCTGCACCCGCGCCGATTGCGGTGTTGCGGTCCCGTTTAGACCAGTTAGAACAGGCACTCAGAGACATTTCCAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTg                                                                                                              <  1:226495/150‑1 (MQ=255)
                                    cccAACGTACTGCCATCGGTCAGTACTGCACCGCCTAATGCCCCTGCACCCGCGCCGATTGCGGTGTTGCGGTCCCGTTTAGACCAGTTAGAACAGGCACTCAGAGACATTTCCAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTg                                                                                                              >  1:4958146/1‑150 (MQ=255)
                                       aaCGTACTGCCATCGGTCAGTACTGCACCGCCTAATGCCCCTGCACCCGCGCCGATTGCGGTGTTGCGGTCCCGTTTAGACCAGTTAGAACAGGCACTCAGAGACATTTCCAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTGCTc                                                                                                           <  2:3406186/150‑1 (MQ=255)
                                       aaCGTACTGCCATCGGTCAGTACTGCACCGCCTAATGCCCCTGCACCCGCGCCGATTGCGGTGTTGCGGTCCCGTTTAGACCAGTTAGAACAGGCACTCAGAGACATTTCCAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTGCTc                                                                                                           >  1:2227797/1‑150 (MQ=255)
                                         cGTACTGCCATCGGTCAGTACTGCACCGCCTAATGCCCCTGCACCCGCGCCGATTGCGGTGTTGCGGTCCCGTTTAGACCAGTTAGAACAGGCACTCAGAGACATTTCCAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTGCTCGt                                                                                                         >  2:2621649/1‑150 (MQ=255)
                                            aCTGCCATCGGTCAGTACTGCACCGCCTAATGCCCCTGCACCCGCGCCGATTGCGGTGTTGCGGTCCCGTTTAGACCAGTTAGAACAGGCACTCAGAGACATTTCCAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTGCTCGTTAc                                                                                                      >  2:976380/1‑150 (MQ=255)
                                            aCTGCCATCGGTCAGTACTGCACCGCCTAATGCCCCTGCACCCGCGCCGATTGCGGTGTTGCGGTCCCGTTTAGACCAGTTAGAACAGGCACTCAGAGACATTTCCAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTGCTCGTTAc                                                                                                      >  1:3761709/1‑150 (MQ=255)
                                              tGCCATCGGTCAGTACTGCACCGCCTAATGCCCCTGCACCCGCGCCGATTGCGGTGTTGCGGTCCCGTTTAGACCAGTTAGAACAGGCACTCAGAGACATTTCCAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTGCTCGTTACaa                                                                                                    <  1:2858826/150‑1 (MQ=255)
                                                 cATCGGTCAGTACTGCACCGCCTAATGCCCCTGCACCCGCGCCGATTGCGGTGTTGCGGTCCCGTTTAGACCAGTTAGAACAGGCACTCAGAGACATTTCCAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTGCTCGTTACAAACa                                                                                                 <  2:409410/150‑1 (MQ=255)
                                                  aTCGGTCAGTACTGCACCGCCTAATGCCCCTGCACCCGCGCCGATTGCGGTGTTGCGGTCCCGTTTAGACCAGTTAGAACAGGCACTCAGAGACATTTCCAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTGCTCGTTACAAACat                                                                                                >  1:2327675/1‑150 (MQ=255)
                                                      gTCAGTACTGCACCGCCTAATGCCCCTGCACCCGCGCCGATTGCGGTGTTGCGGTCCCGTTTAGACCAGTTAGAACAGGCACTCAGAGACATTTCCAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTGCTCGTTACAAACataata                                                                                            <  2:4623367/150‑1 (MQ=255)
                                                       tCAGTACTGCACCGCCTAATGCCCCTGCACCCGCGCCGATTGCGGTGTTGCGGTCCCGTTTAGACCAGTTAGAACAGGCACTCAGAGACATTTCCAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTGCTCGTTACAAACATAATAc                                                                                           <  2:1660456/150‑1 (MQ=255)
                                                          gTACTGCACCGCCTAATGCCCCTGCACCCGCGCCGATTGCGGTGTTGCGGTCCCGTTTAGACCAGTTAGAACAGGCACTCAGAGACATTTCCAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTGCTCGTTACAAACATAATActct                                                                                        >  1:2412524/1‑150 (MQ=255)
                                                             cTGCACCGCCTAATGCCCCTGCACCCGCGCCGATTGCGGTGTTGCGGTCCCGTTTAGACCAGTTAGAACAGGCACTCAGAGACATTTCCAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTGCTCGTTACAAACATAATACtctctc                                                                                     >  1:1048446/1‑150 (MQ=255)
                                                             cTGCACCGCCTAATGCCCCTGCACCCGCGCCGATTGCGGTGTTGCGGTCCCGTTTAGACCAGTTAGAACAGGCACTCAGAGACATTTCCAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTGCTCGTTACAAACATAATACtctctc                                                                                     <  2:4934551/150‑1 (MQ=255)
                                                                cACCGCCTAATGCCCCTGCACCCGCGCCGATTGCGGTGTTGCGGTCCCGTTTAGACCAGTTAGAACAGGCACTCAGAGACATTTCCAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTGCTCGTTACAAACATAATACTCTCTCCTg                                                                                  <  1:4416941/150‑1 (MQ=255)
                                                                cACCGCCTAATGCCCCTGCACCCGCGCCGATTGCGGTGTTGCGGTCCCGTTTAGACCAGTTAGAACAGGCACTCAGAGACATTTCCAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTGCTCGTTACAAACATAATACTCTCTCCTg                                                                                  >  1:3597373/1‑150 (MQ=255)
                                                                 aCCGCCTAATGCCCCTGCACCCGCGCCGATTGCGGTGTTGCGGTCCCGTTTAGACCAGTTAGAACAGGCACTCAGAGACATTTCCAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTGCTCGTTACAAACATAATACTCTCTCCTGa                                                                                 >  2:1469445/1‑150 (MQ=255)
                                                                  ccGCCTAATGCCCCTGCACCCGCGCCGATTGCGGTGTTGCGGTCCCGTTTAGACCAGTTAGAACAGGCACTCAGAGACATTTCCAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTGCTCGTTACAAACATAATACTCTCTCCTGaa                                                                                >  2:4315268/1‑150 (MQ=255)
                                                                  ccGCCTAATGCCCCTGCACCCGCGCCGATTGCGGTGTTGCGGTCCCGTTTAGACCAGTTAGAACAGGCACTCAGAGACATTTCCAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTGCTCGTTACAAACATAATACTCTCTCCTGaa                                                                                >  1:584978/1‑150 (MQ=255)
                                                                   cGCCTAATGCCCCTGCACCCGCGCCGATTGCGGTGTTGCGGTCCCGTTTAGACCAGTTAGAACAGGCACTCAGAGACATTTCCAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTGCTCGTTACAAACATAATACTCTCTCCTGAAt                                                                               <  1:2726370/150‑1 (MQ=255)
                                                                   cGCCTAATGCCCCTGCACCCGCGCCGATTGCGGTGTTGCGGTCCCGTTTAGACCAGTTAGAACAGGCACTCAGAGACATTTCCAAAGTAATTGCCAGAACAGCCGAGGTCATTTTTTTGCTCGTTACAAACATAATACTCTCTCCTGAAt                                                                               <  2:4502427/150‑1 (MQ=255)
                                                                    gCCTAATGCCCCTGCACCCGCGCCGATTGCGGTGTTGCGGTCCCGTTTAGACCAGTTAGAACAGGCACTCAGAGACATTTCCAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTGCTCGTTACAAACATAATACTCTCTCCTGAAtt                                                                              >  1:1169532/1‑150 (MQ=255)
                                                                    gCCTAATGCCCCTGCACCCGCGCCGATTGCGGTGTTGCGGTCCCGTTTAGACCAGTTAGAACAGGCACTCAGAGACATTTCCAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTGCTCGTTACAAACATAATACTCTCTCCTGAAtt                                                                              <  1:1110471/150‑1 (MQ=255)
                                                                     ccTAATGCCCCTGCACCCGCGCCGATTGCGGTGTTGCGGTCCCGTTTAGACCAGTTAGAACAGGCACTCAGAGACATTTCCAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTGCTCGTTACAAACATAATACTCTCTCCTGAAttt                                                                             <  1:4331778/150‑1 (MQ=255)
                                                                       tAATGCCCCTGCACCCGCGCCGATTGCGGTGTTGCGGTCCCGTTTAGACCAGTTAGAACAGGCACTCAGAGACATTGCCAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTGCTCGTTACAAACATAATACTCTCTCCTGAATTTAt                                                                           >  1:4763884/1‑150 (MQ=255)
                                                                        aaTGCCCCTGCACCCGCGCCGATTGCGGTGTTGCGGTCCCGTTTAGACCAGTTAGAACAGGCACTCAGAGACATTGCCAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTGCTCGTTACAAACATAATACTCTCTCCTGAATTTATg                                                                          <  1:884330/150‑1 (MQ=255)
                                                                           gCCCCTGCACCCGCGCCGATTGCGGTGTTGCGGTCCCGTTTAGACCAGTTAGAACAGGCACTCAGAGACATTTCCAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTGCTCGTTACAAACATAATACTCTCTCCTGAATTTATGAtt                                                                       <  1:1588011/150‑1 (MQ=255)
                                                                             cccTGCACCCGCGCCGATTGCGGTGTTGCTGTCCCGTTTAGACCATTTAGAACAGTCACTAAGAGACATTTCCAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTGCTCGTTACAAACATAATACTCTCTCCTGAATTTATGATTCa                                                                     <  1:4163585/150‑1 (MQ=255)
                                                                              ccTGCACCCGCGCCGATTGCGGTGTTGCGGTCCCGTTTAGACCAGTTAGAACAGGCACTCAGAGACATTGCCAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTGCTCGTTACAAACATAATACTCTCTCCTGAATTTATGATTCAc                                                                    >  2:2679421/1‑150 (MQ=255)
                                                                               cTGCACCCGCGCCGATTGCGGTGTTGCGGTCCCGTTTAGACCAGTTAGAACAGGCACTCAGAGACATTTCCAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTGCTCGTTACAAACATAATACTCTCTCCTGAATTTATGATTCACg                                                                   <  2:3807762/150‑1 (MQ=255)
                                                                               cTGCACCCGCGCCGATTGCGGTGTTGCGGTCCCGTTTAGACCAGTTAGAACAGGCACTCAGAGACATTTCCAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTGCTCGTTACAAACATAATACTCTCTCCTGAATTTATGATTCACg                                                                   >  1:4253879/1‑150 (MQ=255)
                                                                               cTGCACCCGCGCCGATTGCGGTGTTGCGGTCCCGTTTAGACCAGTTAGAACAGGCACTCAGAGACATTTCCAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTGCTCGTTACAAACATAATACTCTCTCCTGAATTTATGATTCACg                                                                   <  2:2795740/150‑1 (MQ=255)
                                                                               cTGCACCCGCGCCGATTGCGGTGTTGCGGTCCCGTTTAGACCAGTTAGAACAGGCACTCAGAGACATTTCCAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTGCTCGTTACAAACATAATACTCTCTCCTGAATTTATGATTCACg                                                                   <  1:718104/150‑1 (MQ=255)
                                                                                tGCACCCGCGCCGATTGCGGTGTTGCGGTCCCGTTTAGACCAGTTAGAACAGGCACTCAGAGACATTTCCAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTGCTCGTTACAAACATAATACTCTCTCCTGAATTTATGATTCACgg                                                                  <  1:4830899/150‑1 (MQ=255)
                                                                                 gCACCCGCGCCGATTGCGGTGTTGCGGTCCCGTTTAGACCAGTTAGAACAGGCACTCAGAGACATTTCCAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTGCTCGTTACAAACATAATACTCTCTCCTGAATTTATGATTCACGGa                                                                 <  1:257614/150‑1 (MQ=255)
                                                                                  cACCCGCGCCGATTGCGGTGTTGCGGTCCCGTTTAGACCAGTTAGAACAGGCACTCAGAGACATTTCCAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTGCTCGTTACAAACATAATACTCTCTCCTGAATTTATGATTCACGGaa                                                                <  2:3720927/150‑1 (MQ=255)
                                                                                   aCCCGCGCCGATTGCGGTGTTGCGGTCCCGTTTAGACCAGTTAGAACAGGCACTCAGAGACATTTCCAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTGCTCGTTACAAACATAATACTCTCTCCTGAATTTATGATTCACGGAAg                                                               >  1:2492579/1‑150 (MQ=255)
                                                                                    cccGCGCCGATTGCGGTGTTGCGGTCCCGTTTAGACCAGTTAGAACAGGCACTCAGAGACATTTCCAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTGCTCGTTACAAACATAATACTCTCTCCTGAATTTATGATTCACGGAAGt                                                              >  2:3623149/1‑150 (MQ=255)
                                                                                     ccGCGCCGATTGCGGTGTTGCGGTCCCGTTTAGACCAGTTAGAACAGGCACTCAGAGACATTTCCAACGTAATTGCCAGAACAGCCGCGGTCATTTTTTTGCTCGTTACAAACATAATACTCTCTCCTGAATTTATGATTCACGGAAGTa                                                             <  2:3512468/150‑1 (MQ=255)
                                                                                     ccGCGCCGATTGCGGTGTTGCGGTCCCGTTTAGACCAGTTAGAACAGGCACTCAGAGACATTTCCAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTGCTCGTTACAAACATAATACTCTCTCCTGAATTTATGATTCACGGAAGTa                                                             >  2:858450/1‑150 (MQ=255)
                                                                                      cgcgCCGATTGCGGTGTTGCGGTCCCGTTTAGACCAGTTAGAACAGGCACTCAGAGACATTTCCAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTGCTCGTTACAAACATAATACTCTCTCCTGAATTTATGATTCACGGAAGTaa                                                            >  2:1083411/1‑150 (MQ=255)
                                                                                         gcCGATTGCGGTGTTGCGGTCCCGTTTAGACCAGTTAGAACAGGCACTCAGAGACATTTCCAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTGCTCGTTACAAACATAATACTCTCTCCTGAATTTATGATTCACGGAAGTAAGct                                                         <  2:4817013/150‑1 (MQ=255)
                                                                                          ccGATTGCGGTGTTGCGGTCCCGTTTAGACCAGTTAGAACAGGCACTCAGAGACATTTCCAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTGCTCGTTACAAACATAATACTCTCTCCTGAATTTATGATTCACGGAAGTAAGctc                                                        <  1:3025782/150‑1 (MQ=255)
                                                                                          ccGATTGCGGTGTTGCGGTCCCGTTTAGACCAGTTAGAACAGGCACTCAGAGACATTGCCAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTGCTCGTTACAAACATAATACTCTCTCCTGAATTTATGATTCACGGAAGTAAGctc                                                        <  2:3133567/150‑1 (MQ=255)
                                                                                            gATTGCGGTGTTGCGGTCCCGTTTAGACCAGTTAGAACAGGCACTCAGAGACATTTCCAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTGCTCGTTACAAACATAATACTCTCTCCTGAATTTATGATTCACGGAAGTAAGctctc                                                      <  1:3624478/150‑1 (MQ=255)
                                                                                               tGCGGTGTTGCGGTCCCGTTTAGACCAGTTAGAACAGGCACTCAGAGACATTTCCAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTGCTCGTTACAAACATAATACTCTCTCCTGAATTTATGATTCACGGAAGTAAGCTCTctat                                                   <  1:1981499/150‑1 (MQ=255)
                                                                                                gCGGTGTTGCGGTCCCGTTTAGACCAGTTAGAACAGGCACTCAGAGACATTTCCAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTGCTCGTTACAAACATAATACTCTCTCCTGAATTTATGATTCACGGAAGTAAGCTCTctata                                                  <  2:3657868/150‑1 (MQ=255)
                                                                                                   gtgtTGCGGTCCCGTTTAGACCAGTTAGAACAGGCACTCAGAGACATTTCCAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTGCTCGTTACAAACATAATACTCTCTCCTGAATTTATGATTCACGGAAGTAAGCTCTctataact                                               >  1:3921794/1‑150 (MQ=255)
                                                                                                      ttGCGGTCCCGTTTAGACCAGTTAGAACAGGCACTCAGAGACATTTCCAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTGCTCGTTACAAACATAATACTCTCTCCTGAATTTATGATTCACGGAAGTAAGCTCTctataactata                                            >  1:3826961/1‑150 (MQ=255)
                                                                                                      ttGCGGTCCCGTTTAGACCAGTTAGAACAGGCACTCAGAGACATTTCCAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTGCTCGTTACAAACATAATACTCTCTCCTGAATTTATGATTCACGGAAGTAAGCTCTctataactata                                            >  1:3517278/1‑150 (MQ=255)
                                                                                                         cGGTCCCGTTTAGACCAGTTAGAACAGGCACTCAGAGACATTTCCAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTGCTCGTTACAAACATAATACTCTCTCCTGAATTTATGATTCACGGAAGTAAGCTCTCTATAACTataata                                         <  2:2400204/150‑1 (MQ=255)
                                                                                                             cccGTTTAGACCAGTTAGAACAGGCACTCAGAGACATTTCCAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTGCTCGTTACAAACATAATACTCTCTCCTGAATTTATGATTCACGGAAGTAAGCTCTCTATAACTATAATAGCTa                                     >  2:4906203/1‑150 (MQ=255)
                                                                                                              ccGTTTAGACCAGTTAGAACAGGCACTCAGAGACATTTCCAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTGCTCGTTACAAACATAATACTCTCTCCTGAATTTATGATTCACGGAAGTAAGCTCTCTATAACTATAATAGCTaa                                    <  2:5026663/150‑1 (MQ=255)
                                                                                                              ccGTTTAGACCAGTGAGAACAGGCACTCAGAGACATTTCCAAACTAATTGCCAGAACAGCCGCGGTCATTTTTTTGCTCGTTACAAACATAATACTCTCTCCTGAATTTATGATTCACGGAAGTAAGCTCTCTATAACTATAATAGCTaa                                    <  1:3746906/150‑1 (MQ=255)
                                                                                                               cGTTTAGACCAGTTAGAACAGGCACTCAGAGACATTTCCAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTGCTCGTTACAAACATAATACTCTCTCCTGAATTTATGATTCACGGAAGTAAGCTCTCTATAACTATAATAGCTAAg                                   >  2:4698680/1‑150 (MQ=255)
                                                                                                                gTTTAGACCAGTTAGAACAGGCACTCAGAGACATTTCCAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTGCTCGTTACAAACATAATACTCTCTCCTGAATTTATGATTCACGGAAGTAAGCTCTCTATAACTATAATAGCTAAGa                                  >  2:4031958/1‑150 (MQ=255)
                                                                                                                 tttAGACCAGTTAGAACAGGCACTCAGAGACATTTCCAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTGCTCGTTACAAACATAATACTCTCTCCTGAATTTATGATTCACGGAAGTAAGCTCTCTATAACTATAATAGCTAAGaa                                 <  1:1335734/150‑1 (MQ=255)
                                                                                                                 tttAGACCAGTTAGAACAGGCACTCAGAGACATGTCCAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTGCTCGTTACAAACATAATACTCTCTCCTGAATTTATGATTCACGGAAGTAAGCTCTCTATAACTATAATAGCTAAGaa                                 <  1:205037/150‑1 (MQ=255)
                                                                                                                  ttAGACCAGTTAGAACAGGCACTCAGAGACATTTCCAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTGCTCGTTACAAACATAATACTCTCTCCTGAATTTATGATTCACGGAAGTAAGCTCTCTATAACTATAATAGCTAAGAAt                                >  1:4317155/1‑150 (MQ=255)
                                                                                                                    aGACCAGTTAGAACAGGCACTCAGAGACATTTCCAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTGCTCGTTACAAACATAATACTCTCTCCTGAATTTATGATTCACGGAAGTAAGCTCTCTATAACTATAATAGCTAAGAATaa                              >  2:2353158/1‑150 (MQ=255)
                                                                                                                           ttAGAACAGGCACTCAGAGACATTTCCAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTGCTCGTTACAAACATAATACTCTCTCCTGAATTTATGATTCACGGAAGTAAGCTCTCTATAACTATAATAGCTAAGAATAAGTCTGGt                       >  1:3689551/1‑150 (MQ=255)
                                                                                                                              gAACAGGCACTCAGAGACATTTCCAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTGCTCGTTACAAACATAATACTCTCTCCTGAATTTATGATTCACGGAAGTAAGCTCTCTATAACTATAATAGCTAAGAATAAGTCTGGTGaa                    <  2:1840243/150‑1 (MQ=255)
                                                                                                                                      aCTCAGAGACATTTCCAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTGCTCGTTCCAAACATAATACTCTCTCCTGAATTTATGATTCACGGAAGTAAGCTCTCTATAACTATAATAGCTAAGAATAAGTCTGGTGAAATTATTcc            >  2:1346883/1‑150 (MQ=255)
                                                                                                                                         cAGAGACATTTCCAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTGCTCGTTACAAACATAATACTCTCTCCTGAATTTATGATTCACGGAAGTAAGCTCTCTATAACTATAATAGCTAAGAATAAGTCTGGTGAAATTATTCCGTg         <  1:552386/150‑1 (MQ=255)
                                                                                                                                          agagACATTTCCAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTGCTCGTTACAAACATAATACTCTCTCCTGAATTTATGATTCACGGAAGTAAGCTCTCTATAACTATAATAGCTAAGAATAAGTCTGGTGAAATTATTCCGTGa        >  1:3360499/1‑150 (MQ=255)
                                                                                                                                          agagACATTTCCAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTGCTCGTTACAAACATAATACTCTCTCCTGAATTTATGATTCACGGAAGTAAGCTCTCTATAACTATAATAGCTAAGAATAAGTCTGGTGAAATTATTCCGTGa        >  2:1685603/1‑150 (MQ=255)
                                                                                                                                              aCATTTCCAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTGCTCGTTACAAACATAATACTCTCTCCTGAATTTATGATTCACGGAAGTAAGCTCTCTATAACTATAATAGCTAAGAATAAGTCTGGTGAAATTATTCCGTGAAATc    <  2:903937/150‑1 (MQ=255)
                                                                                                                                              aCATTTCCAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTGCTCGTTACAAACATAATACTCTCTCCTGAATTTATGATTCACGGAAGTAAGCTCTCTATAACTATAATAGCTAAGAATAAGTCTGGTGAAATTATTCCGTGAAATCt   <  1:220973/150‑1 (MQ=255)
                                                                                                                                               cATTTCCAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTGCTCGTTACAAACATAATACTCTCTCCTGAATTTATGATTCACGGAAGTAAGCTCTCTATAACTATAATAGCTAAGAATAAGTCTGGTGAAATTATTCCGTGAAATCt   <  2:4419624/150‑1 (MQ=255)
                                                                                                                                                aTTTCCAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTGCTCGTTACAAACATAATACTCTCTCCTGAATTTATGATTCACGGAAGTAAGCTCTCTATAACTATAATAGCTAAGAATAAGTCTGGTGAAATTATTCCGTGAAATCTg  <  1:179437/150‑1 (MQ=255)
                                                                                                                                                aTTTCCAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTGCTCGTTACAAACATAATACTCTCTCCTGAATTTATGATTCACGGAAGTAAGCTCTCTATAACTATAATAGCTAAGAATAAGTCTGGTGAAATTATTCCGTGAAATCTg  <  2:3761709/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
CCAATAACACCTCCGACAGCTGCACCACCTAATGTACCCAACGTACTGCCATCGGTCAGTACTGCACCGCCTAATGCCCCTGCACCCGCGCCGATTGCGGTGTTGCGGTCCCGTTTAGACCAGTTAGAACAGGCACTCAGAGACATTGCCAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTGCTCGTTACAAACATAATACTCTCTCCTGAATTTATGATTCACGGAAGTAAGCTCTCTATAACTATAATAGCTAAGAATAAGTCTGGTGAAATTATTCCGTGAAATCTG  >  CP009273/1337386‑1337679

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: