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

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,337,5330GT33.8% 313.3 / 212.2 207A18E (GCA→GAA) osmBlipoprotein
Reads supporting (aligned to +/- strand):  ref base G (73/64);  new base T (26/44);  total (99/108)
Fisher's exact test for biased strand distribution p-value = 3.91e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.85e-01

Only 100 of 208 total aligned reads displayed.

CCAATAACACCTCCGACAGCTGCACCACCTAATGTACCCAACGTACTGCCATCGGTCAGTACTGCACCGCCTAATGCCCCTGCACCCGCGCCGATTGCGGTGTTGCGGTCCCGTTTAGACCAGTTAGAACAGGCACTCAGAGACATTGCCAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTGCTCGTTACAAACATAATACTCTCTCCTGAATTTATGATTCACGGAAGTAAGCTCTCTATAACTATAATAGCTAAGAATAAGTCTGGTGAAATTATTCCGTGAAATCTGCTC  >  CP009273/1337386‑1337682
                                                                                                                                                   |                                                                                                                                                     
 cAATAACACCTCCGACAGCTGCACCACCTAATGTACCCAACGTACTGCCATCGGTCAGTACTGCACCGCCTAATGCCCCTGCACCCGCGCCGATTGCGGTGTTGCGGTCCCGTTTAGACCAGTTAGAACAGGCACTCAGAGACATTGCCa                                                                                                                                                    >  1:502285/1‑150 (MQ=255)
   aTAACACCTCCGACAGCTGCACCACCTAATGTACCCAACGTACTGCCATCGGTCAGTACTGCACCGCCTAATGCCCCTGCACCCGCGCCGATTGCGGTGTTGCGGTCCCGTTTAGACCAGTTAGAACAGGCACTCAGAGACATTTCCaaa                                                                                                                                                  >  2:2837944/1‑150 (MQ=255)
    tAACACCTCCGACAGCTGCACCACCTAATGTACCCAACGTACTGCCATCGGTCAGTACTGCACCGCCTAATGCCCCTGCACCCGCGCCGATTGCGGTGTTGCGGTCCCGTTTAGACCAGTTAGAACAGGCACTCAGAGACATTGCCAAAg                                                                                                                                                 >  1:1346828/1‑150 (MQ=255)
     aaCACCTCCGACAGCTGCACCACCTAATGTACCCAACGTACTGCCATCGGTCAGTACTGCACCGCCTAATGCCCCTGCACCCGCGCCGATTGCGGTGTTGCGGTCCCGTTTAGACCAGTTAGAACAGGCACTCAGAGACATTGCCAAAGt                                                                                                                                                >  2:1156014/1‑150 (MQ=255)
      acacCTCCGACAGCTGCACCACCTAATGTACCCAACGTACTGCCATCGGTCAGTACTGCACCGCCTAATGCCCCTGCACCCGCGCCGATTGCGGTGTTGCGGTCCCGTTTAGACCAGTTAGAACAGGCACTCAGAGACATTTCCAAAGTa                                                                                                                                               <  1:3975879/150‑1 (MQ=255)
      acacCTCCGACAGCTGCACCACCTAATGTACCCAACGTACTGCCATCGGTCAGTACTGCACCGCCTAATGCCCCTGCACCCGCGCCGATTGCGGTGTTGCGGTCCCGTTTAGACCAGTTAGAACAGGCACTCAGAGACATTGCCAAAGTa                                                                                                                                               >  1:3442471/1‑150 (MQ=255)
       cacCTCCGACAGCTGCACCACCTAATGTACCCAACGTACTGCCATCGGTCAGTACTGCACCGCCTAATGCCCCTGCACCCGCGCCGATTGCGGTGTTGCGGTCCCGTTTAGACCAGTTAGAACAGGCACTCAGAGACATTGCCAAAGTaa                                                                                                                                              <  2:3139695/150‑1 (MQ=255)
        acCTCCGACAGCTGCACCACCTAATGTACCCAACGTACTGCCATCGGTCAGTACTGCACCGCCTAATGCCCCTGCACCCGCGCCGATTGCGGTGTTGCGGTCCCGTTTAGACCAGTTAGAACAGGCACTCAGAGACATTGCCAAAGTAAt                                                                                                                                             <  2:1251393/150‑1 (MQ=255)
              gACAGCTGCACCACCTAATGTACCCAACGTACTGCCATCGGTCAGTACTGCACCGCCTAATGCCCCTGCACCCGCGCCGATTGCGGTGTTGCGGTCCCGTTTAGACCAGTTAGAACAGGCACTCAGAGACATTGCCAAAGTAATTGCCAg                                                                                                                                       >  2:3744861/1‑150 (MQ=255)
                      caccacCTAATGTACCCAACGTACTGCCATCGGTCAGTACTGCACCGCCTAATGCCCCTGCACCCGCGCCGATTGCGGTGTTGCGGTCCCGTTTAGACCAGTTAGAACAGGCACTCAGAGACATTGCCAAAGTAATTGCCAGAACAGCcg                                                                                                                               <  1:782896/150‑1 (MQ=255)
                          accTAATGTACCCAACGTACTGCCATCGGTCAGTACTGCACCGCCTAATGCCCCTGCACCCGCGCCGATTGCGGTGTTGCGGTCCCGTTTAGACCAGTTAGAACAGGCACTCAGAGACATTGCCAAAGTAATTGCCAGAACAGCCGCGGt                                                                                                                           <  2:4587326/150‑1 (MQ=255)
                          accTAAGGTACCCAACGTACTGCCATCGGTCAGTACTGCACCGCCTAATGCCCCTGCACCCGCGCCGATTGCGGTGTTGCGGTCCCGTTTAGACCAGTTAGAACAGGCACTCAGAGACATTGCCAAAGTAATTGCCAGAACAGCCGCGGt                                                                                                                           <  1:3231661/150‑1 (MQ=255)
                            cTAATGTACCCAACGTACTGCCATCGGTCAGTACTGCACCGCCTAATGCCCCTGCACCCGCGCCGATTGCGGTGTTGCGGTCCCGTTTAGACCAGTTAGAACAGGCACTCAGAGACATTTCCAAAGTAATTGCCAGAACAGCCGCGGTCa                                                                                                                         <  1:2471082/150‑1 (MQ=255)
                            cTAATGTACCCAACGTACTGCCATCGGTCAGTACTGCACCGCCTAATGCCCCTGCACCCGCGCCGATTGCGGTGTTGCGGTCCCGTTTAGACCAGTTAGAACAGGCACTCAGAGACATTGCCAAAGTAATTGCCAGAACAGCCGCGGTCa                                                                                                                         >  2:3383925/1‑150 (MQ=255)
                              aaTGTACCCAACGTACTGCCATCGGTCAGTACTGCACCGCCTAATGCCCCTGCACCCGCGCCGATTGCGGTGTTGCGGTCCCGTTTAGACCAGTTAGAACAGGCACTCAGAGACATTGCCAAAGTAATTGCCAGAACAGCCGCGGTCAtt                                                                                                                       <  2:2221586/150‑1 (MQ=255)
                               aTGTACCCAACGTACTGCCATCGGTCAGTACTGCACCGCCTAATGCCCCTGCACCCGCGCCGATTGCGGTGTTGCGGTCCCGTTTAGACCAGTTAGAACAGGCACTCAGAGACATTGCCAAAGTAATTGCCAGAACAGCCGCGGTCAttt                                                                                                                      >  2:710378/1‑150 (MQ=255)
                               aTGTACCCAACGTACTGCCATCGGTCAGTACTGCACCGCCTAATGCCCCTGCACCCGCGCCGATTGCGGTGTTGCGGTCCCGTTTAGACCAGTTAGAACAGGCACTCAGAGACATTGCCAAAGTAATTGCCAGAACAGCCGCGGTCAttt                                                                                                                      >  2:2275630/1‑150 (MQ=255)
                                tGTACCCAACGTACTGCCATCGGTCAGTACTGCACCGCCTAATGCCCCTGCACCCGCGCCGATTGCGGTGTTGCGGTCCCGTTTAGACCAGTTAGAACAGGCACTCAGAGACATTGCCAAAGTAATTGCCAGAACAGCCGCGGTCAtttt                                                                                                                     <  1:3528072/150‑1 (MQ=255)
                                tGTACCCAACGTACTGCCATCGGTCAGTACTGCACCGCCTAATGCCCCTGCACCCGCGCCGATTGCGGTGTTGCGGTCCCGTTTAGACCAGTTAGAACAGGCACTCAGAGACATTGCCAAAGTAATTGCCAGAACAGCCGCGGTCAtttt                                                                                                                     >  1:3714974/1‑150 (MQ=255)
                                 ggaCCCAACGTACTGCCATCGGTCAGTACTGCACCGCCTAATGCCCCTGCACCCGCGCCGATTGCGGTGTTGCGGTCCCGTTTAGACCAGTTAGAACAGGCACTCAGAGACATTGCCAAAGTAATTGCCAGAACAGCCGCGGTCAttttt                                                                                                                    <  2:5184921/148‑1 (MQ=255)
                                 gTACCCAACGTACTGCCATCGGTCAGTACTGCACCGCCTAATGCCCCTGCACCCGCGCCGATTGCGGTGTTGCGGTCCCGTTTAGACCAGTTAGAACAGGCACTCAGAGACATTTCCAAAGTAATTGCCAGAACAGCCGCGGTCAttttt                                                                                                                    <  2:5164816/150‑1 (MQ=255)
                                        tcGTACTGCCATCGGTCAGTTCTGCACAGCCTTATGCCCCTGCACCCGCGCCGAGTGGGGTGTTGCGGTCCCGTTTAGACCAGTTAGAACAGGCTCTCAGATACATTTCCAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTGCTCg                                                                                                             <  1:3062894/149‑1 (MQ=255)
                                          gTACTGCCATCGGTCAGTACTGCACCGCCTAATGCCCCTGCACCCGCGCCGATTGCGGTGTTGCGGTCCCGTTTAGACCAGTTAGAACAGGCACTCAGAGACATTTCCAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTGCTCGtt                                                                                                           <  2:1434703/150‑1 (MQ=255)
                                           tACTGCCATCGGTCAGTACTGCACCGCCTAATGCCCCTGCACCCGCGCCGATTGCGGTGTTGCGGTCCCGTTTAGACCAGTTAGAACAGGCACTCAGAGACATTTCCAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTGCTCGTTa                                                                                                          >  1:4939078/1‑150 (MQ=255)
                                           tACTGCCATCGGTCAGTACTGCACCGCCTAATGCCCCTGCACCCGCGCCGATTGCGGTGTTGCGGTCCCGTTTAGACCAGTTAGAACAGGCACTCAGAGACATTGCCAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTGCTCGTTa                                                                                                          <  1:2751135/150‑1 (MQ=255)
                                            aCTGCCATCGGTCAGTACTGCACCGCCTAATGCCCCTGCACCCGCGCCGATTGCGGTGTTGCGGTCCCGTTTAGACCAGTTAGAACAGGCACTCAGAGACATTTCCAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTGCTCGTTAc                                                                                                         <  1:4700597/150‑1 (MQ=255)
                                            aCTGCCATCGGTCAGTACTGCACCGCCTAAGGCCCCTGCACCCGCGCCGATTGCGGTGTTGCGGTCCCGTTTAGACCAGTTAGAACAGGCACTCAGAGACATTTCCAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTGCTCGTTAc                                                                                                         <  1:4700667/150‑1 (MQ=255)
                                             cTGCCATCGGTCAGTACTGCACCGCCTAATGCCCCTGCACCCGCGCCGATTGCGGTGTTGCGGTCCCGTTTAGACCAGTTAGAACAGGCACTCAGAGACATTTCCAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTGCTCGTTACa                                                                                                        <  2:2733143/150‑1 (MQ=255)
                                             cTGCCATCGGTCAGTACTGCACCGCCTAATGCCCCTGCACCCGCGCCGATTGCGGTGTTGCGGTCCCGTTTAGACCAGTTAGAACAGGCACTCAGAGACATTGCCAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTGCTCGTTACa                                                                                                        <  2:74418/150‑1 (MQ=255)
                                               gCCATCGGTCAGTACTGCACCGCCTAATGCCCCTGCACCCGCGCCGATTGCGGTGTTGCGGTCCCGTTTAGACCAGTTAGAACAGGCACTCAGAGACATTTCCAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTGCTCGTTACaaa                                                                                                      <  1:1932243/150‑1 (MQ=255)
                                                   tCGGTCAGTACTGCACCGCCTAATGCCCCTGCACCCGCGCCGATTGCGGTGTTGCGGTCCCGTTTAGACCAGTTAGAACAGGCACTCAGAGACATTTCCAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTGCTCGTTACAAACata                                                                                                  <  2:1644782/150‑1 (MQ=255)
                                                        cAGTACTGCACCGCCTAATGCCCCTGCACCCGCGCCGATTGCGGTGTTGCGGTCCCGTTTAGACCAGTTAGAACAGGCACTCAGAGACATTGCCAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTGCTCGTTACAAACATAATAct                                                                                             <  1:2644263/150‑1 (MQ=255)
                                                          gTACTGCACCGCCTAATGCCCCTGCACCCGCGCCGATTGCGGTGTTGCGGTCCCGTTTAGACCAGTTAGAACAGGCACTCAGAGACATTGCCAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTGCTCGTTACAAACATAATActct                                                                                           >  1:4996011/1‑150 (MQ=255)
                                                          gTACTGCACCGCCTAATGCCCCTGCACCCGCGCCGATTGCGGTGTTGCGGTCCCGTTTAGACCAGTTAGAACAGGCACTCAGAGACATTGCCAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTGCTCGTTACAAACATAATActct                                                                                           >  1:4523428/1‑150 (MQ=255)
                                                           tACTGCACCGCCTAATGCCCCTGCACCCGCGCCGATTGCGGTGTTGCGGTCCCGTTTAGACCAGTTAGAACAGGCACTCAGAGACATTTCCAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTGCTCGTTACAAACATAATActctc                                                                                          >  1:3070373/1‑150 (MQ=255)
                                                            aCTGCACCGCCTAAGGCCCCTGCACCCGCGCCGATTGCGGTGTTGCGGTCCCGTTTAGACCAGTTAGAACAGGCACTCAGAGTCATTGCCTAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTGCTCGTTACAAACATAATActctct                                                                                         <  2:4173830/150‑1 (MQ=255)
                                                             cTGCACCGCCTAATGCCCCTGCACCCGCGCCGATTGCGGTGTTGCGGTCCCGTTTAGACCAGTTAGAACAGGCACTCAGAGACATTTCCAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTGCTCGTTACAAACATAATACtctctc                                                                                        >  1:2851279/1‑150 (MQ=255)
                                                             cTGCACCGCCTAATGCCCCTGCACCCGCGCCGATTGCGGTGTTGCGGTCCCGTTTAGACCAGTTAGAACAGGCACTCAGAGACATTGCCAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTGCTCGTTACAAACATAATACtctctc                                                                                        >  1:1839624/1‑150 (MQ=255)
                                                                     ccTAATGCCCCTGCACCCGCGCCGATTGCGGTGTTGCGGTCCCGTTTAGACCAGTTAGAACAGGCCCTCAGAGACATTGCCAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTGCTCGTTACAAACATAATACTCTCTCCTGAAttt                                                                                >  2:4015586/1‑150 (MQ=255)
                                                                      cTAATGCCCCTGCACCCGCGCCGAGTGCGGTGTTGCGGTCCCGTTTAGACCAGTTAGAACAGGCACTCAGAGACATTGCCAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTGCTCGTTACAAACATAATACTCTCTCCTGAATTTa                                                                               >  1:3751421/1‑150 (MQ=255)
                                                                        aaTGCCCCTGCACCCGCGCCGATTGCGGTGTTGCGGTCCCGTTTAGACCAGTTAGAACAGGCACTCAGAGACATTGCCAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTGCTCGTTACAAACATAATACTCTCTCCTGAATTTATg                                                                             <  1:1566175/150‑1 (MQ=255)
                                                                        aaTGCCCCTGCACCCGCGCCGATTGCGGTGTTGCGGTCCCGTTTAGACCAGTTAGAACAGGCACTCAGAGACATTGCCAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTGCTCGTTACAAACATAATACTCTCTCCTGAATTTATg                                                                             >  1:575695/1‑150 (MQ=255)
                                                                        aaTGCCCCTGCACCCGCGCCGATTGCGGTGTTGCGGTCCCGTTTAGACCAGTTAGAACAGGCACTCAGAGACATTGACAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTGCTCGTTACAAACATAATACTCTCTCCTGAATTTATg                                                                             <  1:4781887/150‑1 (MQ=255)
                                                                         aTGCCCCTGCACCCGCGCCGATTGCGGTGTTGCGGTCCCGTTTAGACCAGTTAGAACAGGCACTCAGAGACATTGCCAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTGCTCGTTACAAACATAATACTCTCTCCTGAATTTATGa                                                                            >  2:3143310/1‑150 (MQ=255)
                                                                          tGCCCCTGCACCCGCGCCGATTGCGGTGTTGCGGTCCCGTTTAGACCAGTTAGAACAGGCACTCAGAGACATTTCCAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTGCTCGTTACAAACATAATACTCTCTCCTGAATTTATGAt                                                                           <  2:3595643/150‑1 (MQ=255)
                                                                          tGCCCCTGCACCCGCGCCGATTGCGGTGTTGCGGTCCCGTTTAGACCAGTTAGAACAGGCACTCAGAGACATTGCCAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTGCTCGTTACAAACATAATACTCTCTCCTGAATTTATGAt                                                                           >  1:3875667/1‑150 (MQ=255)
                                                                           gCCCCTGCACCCGCGCCGATTGCGGTGTTGCGGTCCCGTTTAGACCAGTTAGAACAGGCACTCAGAGACATTTCGAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTGCTCGTTACAAACATAATACTCTCTCCTGAATTTATGAtt                                                                          <  2:4915316/150‑1 (MQ=255)
                                                                           gCCCCTGCACCCGCGCCGATTGCGGTGTTGCGGTCCCGTTTAGACCAGTTAGAACAGGCACTCAGAGACATTTCCAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTGCTCGTTACAAACATAATACTCTCTCCTGAATTTATGAtt                                                                          >  2:3179814/1‑150 (MQ=255)
                                                                            ccccTGCACCCGCGCCGATTGCGGTGTTGCGGTCCCGTTTAGACCAGTTAGAACAGGCACTCAGAGACATTGCCAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTGCTCGTTACAAACATAATACTCTCTCCTGAATTTATGATTc                                                                         >  1:4234730/1‑150 (MQ=255)
                                                                             cccTGCACCCGCGCCGATTGCGGTGTTGCGGTCCCGTTTAGACCAGTTAGAACAGGCACCCAGAGACATTTCCAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTGCTCGTTACAAACATAATACTCTCTCCTGAATTTATGATTCa                                                                        <  2:1814741/150‑1 (MQ=255)
                                                                             cccTGCACCCGCGCCGATTGCGGTGTTGCGGTCCCGTTTAGACCAGTTAGAACAGGCACCCAGAGACATTTCCAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTGCTCGTTACAAACATAATACTCTCTCCTGAATTTATGATTCa                                                                        <  2:1814293/150‑1 (MQ=255)
                                                                              ccTGCACCCGCGCCGATTGCGGTGTTGCGGTCCCGTTTAGACCAGTTAGAACAGGCACTCAGAGACATTTCCAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTGCTCGTTACAAACATAATACTCTCTCCTGAATTTATGATTCAc                                                                       >  1:2845391/1‑150 (MQ=255)
                                                                              ccTGCACCCGCGCCGATTGCGGTGTTGCGGTCCCGTTTAGACCAGTTAGAACAGGCACTCAGAGACATTGCCAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTGCTCGTTACAAACATAATACTCTCTCCTGAATTTATGATTCAc                                                                       >  1:2150649/1‑150 (MQ=255)
                                                                              ccTGCACCCGCGCCGAGTGCGGTGTTGCGGTCCCGTTTAGACCAGTTAGAACAGGCACTCAGAGACATTGCCAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTGCTCGTTACAAACATAATACTCTCTCCTGAATTTATGATTCAc                                                                       <  1:4418921/150‑1 (MQ=255)
                                                                               cTGCACCCGCGCCGATTGCGGTGTTGCGGTCCCGTTTAGACCAGTTAGAACAGGCACTCAGAGACATTGCCAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTGCTCGTTACAAACATAATACTCTCTCCTGAATTTATGATTCACg                                                                      <  1:4718748/150‑1 (MQ=255)
                                                                               cTGCACCCGCGCCGATTGCGGTGTTGCGGTCCCGTTTAGACCAGTTAGAACAGGCACTCAGAGACATTGCCAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTGCTCGTTACAAACATAATACTCTCTCCTGAATTTATGATTCACg                                                                      <  1:4137287/150‑1 (MQ=255)
                                                                                  cACCCGCGCCGATTGCGGTGTTGCGGTCCCGTTTAGACCAGTTAGAACAGGCACTCAGAGACATTTCCAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTGCTCGTTACAAACATAATACTCTCTCCTGAATTTATGATTCACGGaa                                                                   <  2:895177/150‑1 (MQ=255)
                                                                                   aCCCGCGCCGATTGCGGTGTTGCGGTCCCGTTTAGACCAGTTAGAACAGGCACTCAGAGACATTTCCAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTGCTCGTTACAAACATAATACTCTCTCCTGAATTTATGATTCACGGAAg                                                                  >  2:927919/1‑150 (MQ=255)
                                                                                   aCCCGCGCCGATTGCGGTGTTGCGGTCCCGTTTAGACCAGTTAGAACAGGCACTCAGAGACATTGCCAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTGCTCGTTACAAACATAATACTCTCTCCTGAATTTATGATTCACGGAAg                                                                  >  2:2482561/1‑150 (MQ=255)
                                                                                   aCCCGCGCCGATTGCGGTGTTGCGGTCCCGTTTAGACCAGTTAGAACAGGCACTCAGAGACATTGCCAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTGCTCGTTACAAACATAATACTCTCTCCTGAATTTATGATTCACGGAAg                                                                  <  2:845962/150‑1 (MQ=255)
                                                                                      cgcgCCGATTGCGGTGTTGCGGTCCCGTTTAGACCAGTTAGAACAGGCACTCAGAGACATTTCCAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTGCTCGTTACAAACATAATACTCTCTCCTGAATTTATGATTCACGGAAGTaa                                                               <  2:2284958/150‑1 (MQ=255)
                                                                                      cgcgCCGATTGCGGTGTTGCGGTCCCGTTTAGACCAGTTAGAACAGGCACTCAGAGACATTGCCAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTGCTCGTTACAAACATAATACTCTCTCCTGAATTTATGATTCACGGAAGTaa                                                               <  2:4035555/150‑1 (MQ=255)
                                                                                      cgcgCCGATTGCGGTGTTGCGGTCCCGTTTAGACCAGTTAGAACAGGCACTCAGAGACATTGCCAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTGCTCGTTACAAACATAATACTCTCTCCTGAATTTATGATTCACGGAAGTaa                                                               <  2:2510066/150‑1 (MQ=255)
                                                                                       gcgcCGATTGCGGTGTTGCGGTCCCGTTTAGACCAGTTAGAACAGGCACTCAGAGACATTGCCAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTGCTCGTTACAAACATAATACTCTCTCCTGAATTTATGATTCACGGAAGTAAg                                                              <  1:3169775/150‑1 (MQ=255)
                                                                                             aTTGCGGTGTTGCGGTCCCGTTTAGACCAGTTAGAACAGGCACTCAGAGACATTTCCAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTGCTCGTTACAAACATAATACTCTCTCCTGAATTTATGATTCACGGAAGTAAGctctct                                                        >  1:4065674/1‑150 (MQ=255)
                                                                                             aTTGCGGTGTTGCGGTCCCGTTTAGACCAGTTAGAACAGGCACTCAGAGACATTGCCAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTGCTCGTTACAAACATAATACTCTCTCCTGAATTTATGATTCACGGAAGTAAGctctct                                                        <  1:1256480/150‑1 (MQ=255)
                                                                                               tGCGGTGTTGCGGTCCCGTTTAGACCAGTTAGAACAGGCACTCAGAGACATTGCCAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTGCTCGTTACAAACATAATACTCTCTCCTGAATTTATGATTCACGGAAGTAAGCTCTctat                                                      <  1:1681315/150‑1 (MQ=255)
                                                                                               tGCGGTGTTGCGGTCCCGTTTAGACCAGTTAGAACAGGCACTCAGAGACATTGCCAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTGCTCGTTACAAACATAATACTCTCTCCTGAATTTATGATTCACGGAAGTAAGCTCTctat                                                      <  2:3672098/150‑1 (MQ=255)
                                                                                                gCGGTGTTGCGGTCCCGTTTAGACCAGTTAGAACAGGCACTCAGAGACATTTCCAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTGCTCGTTACAAACATAATACTCTCTCCTGAATTTATGATTCACGGAAGTAAGCTCTctata                                                     <  1:1570772/150‑1 (MQ=255)
                                                                                                  ggTGTTGCGGTCCCGTTTAGACCAGTTAGAACAGGCACTCAGAGACATTGCCAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTGCTCGTTACAAACATAATACTCTCTCCTGAATTTATGATTCACGGAAGTAAGCTCTctataac                                                   <  2:4362804/150‑1 (MQ=255)
                                                                                                     gtTGCGGTCCCGTTTAGACCAGTTAGAACAGGCACTCAGAGACATTTCCAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTGCTCGTTACAAACATAATACTCTCTCCTGAATTTATGATTCACGGAAGTAAGCTCTctataactat                                                >  2:1279484/1‑150 (MQ=255)
                                                                                                          ggTCCCGTTTAGACCAGTTAGAACAGGCACTCAGAGACATTTCCAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTGCTCGTTACAAACATAATACTCTCTCCTGAATTTATGATTCACGGAAGTAAGCTCTCTATAACTATAATAg                                           >  2:571937/1‑150 (MQ=255)
                                                                                                            tCCCGTTTAGACCAGTTAGAACAGGCACTCAGAGACATTGCCAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTGCTCGTTACAAACATAATACTCTCTCCTGACTTTATGATTCACGGAAGTAAGCTCTCTATAACTATAATAGCt                                         >  2:2690276/1‑150 (MQ=255)
                                                                                                             cccGTTTAGACCAGTTAGAACAGGCACTCAGAGACATTGCCAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTGCTCGTTACAAACATAATACTCTCTCCTGAATTTATGATTCACGGAAGTAAGCTCTCTATAACTATAATAGCTa                                        >  1:2029723/1‑150 (MQ=255)
                                                                                                                gTTTAGACCAGTTAGAACAGGCACTCAGAGACATTTCCAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTGCTCGTTACAAACATAATACTCTCTCCTGAATTTATGATTCACGGAAGTAAGCTCTCTATAACTATAAGAGCTAAGa                                     <  2:783584/150‑1 (MQ=255)
                                                                                                                    aGACCAGTTAGAACAGGCACTCAGAGACATTGCCAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTGCTCGTTACAAACATAATACTCTCTCCTGAATTTATGATTCACGGAAGTAAGCTCTCTATAACTATAATAGCTAAGAATaa                                 <  1:863480/150‑1 (MQ=255)
                                                                                                                        cAGTTAGAACAGGCACTCAGAGACATTGCCAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTGCTCGTTACAAACATAATACTCTCTCCTGAATTTATGATTCACGGAAGTAAGCTCTCTATAACTATAATAGCTAAGAATAAGTCt                             <  2:2571638/150‑1 (MQ=255)
                                                                                                                            tAGAACAGGCACTCAGAGACATTTCCAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTGCTCGTTACAAACATAATACTCTCTCCTGAATTTATGATTCACGGAAGTAAGCTCTCTATAACTATAATAGCTAAGAATAAGTCTGGTg                         >  2:3933485/1‑150 (MQ=255)
                                                                                                                            tAGAACAGGCACTCAGAGACATTGCCAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTGGCTCGTTACAAACATAATACTCTCTCCTGAATTTATGATTCACGGAAGTAAGCTCTCTATAACTATAATAGCTAAGAATAAGTCTGGTg                         <  2:1593073/150‑1 (MQ=255)
                                                                                                                             aGAACAGGCACTCAGAGACATTGCCAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTGCTCGTTACAAACATAATACTCTCTCCTGAATTTATGATTCACGGAAGTAAGCTCTCTATAACTATAATAGCTAAGAATAAGTCTGGTGa                        >  2:1947642/1‑150 (MQ=255)
                                                                                                                               aaCAGGCACTCAGAGACATTGCCAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTGCTCGTTACAAACATAATACTCTCTCCTGAATTTATGATTCACGGAAGTAAGCTCTCTATAACTATAATAGCTAAGAATAAGTCTGGTGaaa                      >  1:3145052/1‑150 (MQ=255)
                                                                                                                                aCAGGCACTCAGAGACATTGCCAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTGCTCGTTACAAACATAATACTCTCTCCTGAATTTATGATTCACGGAAGTAAGCTCTCTATAACTATAATAGCTAAGAATAAGTCTGGTGAAat                     <  1:3837884/150‑1 (MQ=255)
                                                                                                                                 cAGGCACTCAGAGACATTGCCAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTGCTCGTTACAAACATAATACTCTCTCCTGAATTTATGATTCACGGAAGTAAGCTCTCTATAACTATAATAGCTAAGAATAAGTCTGGTGAAatt                    >  1:3052848/1‑150 (MQ=255)
                                                                                                                                 cAGGCACTCAGAGACATTGCCAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTGCTCGTTACAAACATAATACTCTATCCTGAATTTATGATTCACGGAAGTAAGCTCTCTATAACTATAATAGCTAAGAATAAGTCTGGTGAAatt                    >  2:3145608/1‑150 (MQ=255)
                                                                                                                                  aGGCACTCAGAGACATTTCCAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTGCTCGTTACAAACATAATACTCTCTCCTGAATTTATGATTCACGGAAGTAAGCTCTATATAACTATAATAGCTAAGAATAAGTCTGGTGAAatta                   <  2:3060839/150‑1 (MQ=255)
                                                                                                                                   ggCACTCAGAGACATTGCCAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTGCTCGTTACAAACATAATACTCTCTCCTGAATTTATGATTCACGGAAGTAAGCTCTCTATAACTATAATAGCTAAGAATAAGTCTGGTGAAaatat                  >  2:4508831/1‑150 (MQ=255)
                                                                                                                                    gCACTCAGAGACATTGCCAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTGCTCGTTACAAACATAATACTCTCTCCTGAATTTATGATTCACGGAAGTAAGCTCTCTATAACTATAATAGCTAAGAATAAGTCTGGTGAAattatt                 >  2:2579347/1‑150 (MQ=255)
                                                                                                                                     cACTCAGAGACATTTCCAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTGCTCGTTACAAACATAATACTCTCTCCTGAATTTATGATTCACGGAAGTAAGCTCTCTATAACTATAATAGCTAAGAATAAGTCTGGTGAAATTATTc                <  2:3755349/150‑1 (MQ=255)
                                                                                                                                     cACTCAGAGACATTGCCAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTGCTCGTTACAAACATAATACTCTCTCCTGAATTTATGATTCACGGAAGTAAGCTCTCTATAACTATAATAGCTAAGAATAAGTCTGGTGAAATTATTc                >  1:4964142/1‑150 (MQ=255)
                                                                                                                                      aCTCAGAGACATTTCCAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTGCTCGTTACAAACATAATACTCTCTCCTGAATTTATTATTCACGGAAGTAAGCTCTCTATAACTATAATAGCAAAGAATAAGTCTGGTGAAattattac               >  1:2234508/1‑148 (MQ=255)
                                                                                                                                      aCTCAGAGACATTGCCAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTGCTCGTTACAAACATAATACTCTCTCCTGAATTTATGATTCACGGAAGTAAGCTCTCTATAACTATAATAGCTAAGAATAAGTCTGGTGAAATTATTcc               >  2:4704926/1‑150 (MQ=255)
                                                                                                                                      aCTCAGAGACATTGCCAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTGCTCGTTACAAACATAATACTCTCTCCTGAATTTATGATTCACGGAAGTAAGCTCTCTATAACTATAATAGCTAAGAATAAGTCTGGTGAAATTATTcc               >  2:2120045/1‑150 (MQ=255)
                                                                                                                                         cAGAGACATTGCCAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTGCTCGTTACAAACATAATACTCTCTCCTGAATTTATGATTCACGGAAGTAAGCTCTCTATAACTATAATAGCTAAGAATAAGTCTGGTGAAATTATTCCGTg            >  2:2169440/1‑150 (MQ=255)
                                                                                                                                          agagACATTTCCAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTGCTCGTTACAAACATAATACTCTCTCCTGAATTTATGATTCACGGAAGTAAGCTCTCTATAACTATAATAGCTAAGAATAAGTCTGGTGAAATTATTCCGTGa           <  2:2641615/150‑1 (MQ=255)
                                                                                                                                            agaCATTTCCAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTGCTCGTTACAAACATAATACTCTCTCCTGAATTTATGATTCACGGAAGTAAGCTCTCTATAACTATAATAGCTAAGAATAAGTCTGGTGAAATTATTCCGTGaaa         <  1:1083384/150‑1 (MQ=255)
                                                                                                                                              aCATTGCCAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTGCTCGTTACAAACATAATACTCTCTCCTGAATTTATGATTCACGGAAGTAAGCTCTCTATAACTATAATAGCTAAGAATAAGTCTGGTGAAATTATTCCGTGAAATc       >  1:3174948/1‑150 (MQ=255)
                                                                                                                                               cATTTCCAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTGCTCGTTACAAACATAATACTCTCTCCTGAATTTATGATTCACGGAAGTAAGCTCTCTATAACTATAATAGCTAAGAATAAGTCTGGTGAAATTATTCCGTGAAATCt      <  1:143127/150‑1 (MQ=255)
                                                                                                                                                aTTTCCAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTGCTCGTTACAAACATAATACTCTCTCCTGAATTTATGATTCACGGAAGTAAGCTCTCTATAACTATAATAGCTAAGAATAAGTCTGGTGAAATTATTCCGTGAAATCTg     >  2:4480693/1‑150 (MQ=255)
                                                                                                                                                  tGCCAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTGCTCGTTACAAACATAATACTCTCTCCTGAATTTATGATTCACGGAAGTAAGCTCTCTATAACTATAATAGCTAAGAATAAGTCTGGTGAAATTATTCCGTGAAATCTGCt   >  1:4421544/1‑150 (MQ=255)
                                                                                                                                                   gCCAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTGCTCGTTACAAACATAATACTCTCTCCTGAATTTATGATTCACGGAAGTAAGCTCTCTATAACTATAATAGCTAAGAATAAGTCTGGTGAAATTATTCCGTGAAATCTGCTc  <  1:1288153/150‑1 (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: ‑