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

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,337,5330GT93.9% 663.0 / 32.2 213A18E (GCA→GAA) osmBlipoprotein
Reads supporting (aligned to +/- strand):  ref base G (4/9);  new base T (88/112);  total (92/121)
Fisher's exact test for biased strand distribution p-value = 4.01e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

Only 100 of 213 total aligned reads displayed.

AATAACACCTCCGACAGCTGCACCACCTAATGTACCCAACGTACTGCCATCGGTCAGTACTGCACCGCCTAATGCCCCTGCACCCGCGCCGATTGCGGTGTTGCGGTCCCGTTTAGACCAGTTAGAACAGGCACTCAGAGACATTGCCAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTGCTCGTTACAAACATAATACTCTCTCCTGAATTTATGATTCACGGAAGTAAGCTCTCTATAACTATAATAGCTAAGAATAAGTCTGGTGAAATTATTCCGTGAAATCTG  >  CP009273/1337388‑1337679
                                                                                                                                                 |                                                                                                                                                  
aaTAACACCTCCGACAGCTGCACCACCTAATGTACCCAACGTACTGCCATCGGTCAGTACTGCACCGCCTAATGCCCCTGCACCCGCGCCGATTGCGGTGTTGCGGTCCCGTTTAGACCAGTTAGAACAGGCACTCAGAGACATTTCCaa                                                                                                                                                >  1:1709975/1‑150 (MQ=255)
  tAACACCTCCGACAGCTGCACCACCTAATGTACCCAACGTACTGCCATCGGTCAGTACTGCACCGCCTAATGCCCCTGCACCCGCGCCGATTGCGGTGTTGCGGTCCCGTTTAGACCAGTTAGAACAGGCACTCAGAGACATTTCCAAAg                                                                                                                                              >  2:2094736/1‑150 (MQ=255)
   aaCACCTCCGACAGCTGCACCACCTACTGTACCCAACGTACTGCCATCGGTCAGTACTGCACCGCCTAATGCCCCTGCACCCGCGCCGATTGCGGTGTTGCGGTCCCGTTTAGACCAGTTAGAACAGGCACTCAGAGACATTTCCAAAGt                                                                                                                                             <  1:1586669/150‑1 (MQ=255)
       ccTCCGACAGCTGCACCACCTAATGTACCCAACGTACTGCCATCGGTCAGTACTGCACCGCCTAATGCCCCTGCACCCGCGCCGATTGCGGTGTTGCGGTCCCGTTTAGACCAGTTAGAACAGGCACTCAGAGACATTTCCAAAGTAAtt                                                                                                                                         <  2:3581228/150‑1 (MQ=255)
          ccGACAGCTGCACCACCTAATGTACCCAACGTACTGCCATCGGTCAGTACTGCACCGCCTAATGCCCCTGCACCCGCGCCGATTGCGGTGTTGCGGTCCCGTTTAGACCAGTTAGAACAGGCACTCAGAGACATTTCCAAAGTAATTGcc                                                                                                                                      >  2:3308506/1‑150 (MQ=255)
           cGACAGCTGCACCACCTAATGTACCCAACGTACTGCCATCGGTCAGTACTGCACCGCCTAATGCCCCTGCACCCGCGCCGATTGCGGTGTTGCGGTCCCGTTTAGACCAGTTAGAACAGGCACTCAGAGACATTTCCAAAGTAATTGCCa                                                                                                                                     <  2:3253143/150‑1 (MQ=255)
              cAGCTGCACCACCTAATGTACCCAACGTACTGCCATCGGTCAGTACTGCACCGCCTAATGCCCCTGCACCCGCGCCGATTGCGGTGTTGCGGTCCCGTTTAGACCAGTTAGAACAGGCACTCAGAGACATTTCCAAAGTAATTGCCAGaa                                                                                                                                  >  1:500279/1‑150 (MQ=255)
                 cTGCACCACCTAATGTACCCAACGTACTGCCATCGGTCAGTACTGCACCGCCTAATGCCCCTGCACCCGCGCCGATTGCGGTGTTGCGGTCCCGTTTAGACCAGTTAGAACAGGCACTCAGAGACATTTCCAAAGTAATTGCCAGAACAg                                                                                                                               <  2:115322/150‑1 (MQ=255)
                  tGCACCACCTAATGTACCCAACGTACTGCCATCGGTCAGTACTGCACCGCCTAATGCCCCTGCACCCGCGCCGATTGCGGTGTTGCGGTCCCGTTTAGACCAGTTAGAACAGGCACTCAGAGACATTTCCAAAGTAATTGCCAGAACAGc                                                                                                                              <  1:760653/150‑1 (MQ=255)
                   gCACCACCTAATGTACCCAACGTACTGCCATCGGTCAGTACTGCACCGCCTAATGCCCCTGCACCCGCGCCGATTGCGGTGTTGCGGTCCCGTTTAGACCAGTTAGAACAGGCACTCAGAGACATTTCCAAAGTAATTGCCAGAACAGcc                                                                                                                             <  2:782823/150‑1 (MQ=255)
                   gCACCACCTAATGTACCCAACGTACTGCCATCGGTCAGTACTGCACCGCCTAATGCCCCTGCACCCGCGCCGATTGCGGTGTTGCGGTCCCGTTTAGACCAGTTAGAACAGGCACTCAGAGACATTTCCAAAGTAATTGCCAGAACAGcc                                                                                                                             <  2:2871121/150‑1 (MQ=255)
                   gCACCACCTAATGTACCCAACGTACTGCCATCGGTCAGTACTGCACCGCCTAATGCCCCTGCACCCGCGCCGATTGCGGTGTTGCGGTCCCGTTTAGACCAGTTAGAACAGGCACTCAGAGACATTTCCAAAGTAATTGCCAGAACAGcc                                                                                                                             >  1:2412289/1‑150 (MQ=255)
                      ccaccTAATGTACCCAACGTACTGCCATCGGTCAGTACTGCACCGCCTAATGCCCCTGCACCCGCGCCGATTGCGGTGTTGCGGTCCCGTTTAGACCAGTTAGAACAGGCACTCAGAGACATTTCCAAAGTAATTGCCAGAACAGCcgcg                                                                                                                          <  1:2020580/150‑1 (MQ=255)
                               gTACCCAACGTACTGCCATCGGTCAGTACTGCACCGCCTAATGCCCCTGCACCCGCGCCGATTGCGGTGTTGCGGTCCCGTTTAGACCAGTTAGAACAGGCACTCAGAGACATTTCCAAAGTAATTGCCAGAACAGCCGCGGTCAttttt                                                                                                                 >  2:3850636/1‑150 (MQ=255)
                               gTACCCAACGTACTGCCATCGGTCAGTACTGCACCGCCTAATGCCCCTGCACCCGCGCCGATTGCGGTGTTGCGGTCCCGTTTAGACCAGTTAGAACAGGCACTCAGAGACATTTCCAAAGTAATTGCCAGAACAGCCGCGGTCAttttt                                                                                                                 >  1:18042/1‑150 (MQ=255)
                                 aCCCAACGTACTGCCATCGGTCAGTACTGCACCGCCTAATGCCCCTGCACCCGCGCCGATTGCGGTGTTGCGGTCCCGTTTAGACCAGTTAGAACAGGCACTCAGAGACATTTCCAAAGTAATTGCCAGAACAGCCGCGGTCAttttttt                                                                                                               <  2:2935190/150‑1 (MQ=255)
                                 aCCCAACGTACTGCCATCGGTCAGTACTGCACCGCCTAATGCCCCTGCACCCGCGCCGATTGCGGTGTTGCGGTCCCGTTTAGACCAGTTAGAACAGGCACTCAGAGACATTTCCAAAGTAATTGCCAGAACAGCCGCGGTCAttttttt                                                                                                               <  2:3460121/150‑1 (MQ=255)
                                  cccAACGTACTGCCATCGGTCAGTACTGCACCGCCTAATGCCCCTGCACCCGCGCCTATTGCGGTGTTGCGGTCCCGTTTAGACCAGTTAGAACAGGCACTCAGAGACATTGACAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTg                                                                                                              <  1:1623605/150‑1 (MQ=255)
                                  cccAACGTACTGCCATCGGTCAGTACTGCACCGCCTAATGCCCCTGCACCCGCGCCGATTGCGGTGTTGCGGTCCCGTTTAGACCAGTTAGAACAGGCACTCAGAGACATTTCCAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTg                                                                                                              >  1:3429431/1‑150 (MQ=255)
                                   ccAACGTACTGCCATCGGTCAGTACTGCACCGCCTAATGCCCCTGCACCCGCGCCGATTGCGGTGTTGCGGTCCCGTTTAGACCAGTTAGAACAGGCACTCAGAGACATTTCCAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTGc                                                                                                             <  2:3817423/150‑1 (MQ=255)
                                   ccAACGTACTGCCATCGGTCAGTACTGCACCGCCTAATGCCCCTGCACCCGCGCCGATTGCGGTGTTGCGGTCCCGTTTAGACCAGTTAGAACAGGCACTCAGAGACATTGCCAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTGc                                                                                                             >  1:4178158/1‑150 (MQ=255)
                                     aaCGTACTGCCATCGGTCAGTACTGCACCGCCTAATGCCCCTGCACCCGCGCCGATTGCGGTGTTGCGGTCCCTTTTAGACCAGTTAGAACAGGCACTCAGAGACATTTCCAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTGCTc                                                                                                           >  1:3080040/1‑150 (MQ=255)
                                      aCGTACTGCCATCGGTCAGTACTGCACCGCCTAATGCCCCTGCACCCGCGCCGATTGCGGTGTTGCGGTCCCGTTTAGACCAGTTAGAACAGGCACTCAGAGACATTTCCAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTGCTCg                                                                                                          <  1:707408/150‑1 (MQ=255)
                                       cGTACTGCCATCGGTCAGTACTGCACCGCCTAATGCCCCTGCACCCGCGCCGATTGCGGTGTTGCGGTCCCGTTTAGACCAGTTAGAACAGGCACTCAGAGACATTTCCAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTGCTCGt                                                                                                         <  1:676586/150‑1 (MQ=255)
                                        gTACTGCCATCGGTCAGTACTGCACCGCCTAATGCCCCTGCACCCGCGCCGATTGCGGTGTTGCGGTCCCGTTTAGACCAGTTAGAACAGGCACTCAGAGACATTTCCAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTGCTCGtt                                                                                                        <  2:2530633/150‑1 (MQ=255)
                                        gTACTGCCATCGGTCAGTACTGCACCGCCTAATGCCCCTGCACCCGCGCCGATTGCGGTGTTGCGGTCCCGTTTAGACCAGTTAGAACAGGCACTCAGAGACATTTCCAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTGCTCGtt                                                                                                        <  2:1243330/150‑1 (MQ=255)
                                          aCTGCCATCGGTCAGTACTGCACCGCCTAATGCCCCTGCACCCGCGCCGATTGCGGTGTTGCGGTCCCTTTTAGACCAGTTAGAACAGGCACTCAGAGACATTTCCAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTGCTCGTTAc                                                                                                      >  2:3077608/1‑150 (MQ=255)
                                            tGCCATCGGTCAGTACTGCACCGCCTAATGCCCCTGCACCCGCGCCGATTGCGGTGTTGCGGTCCCGTTTAGACCAGTTAGAACAGGCACTCAGAGACATTTCCAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTGCTCGTTACaa                                                                                                    >  2:3474784/1‑150 (MQ=255)
                                            tGCCATCGGTCAGTACTGCACCGCCTAATGCCCCTGCACCCGCGCCGATTGCGGTGTTGCGGTCCCGTTTAGACCAGTTAGAACAGGCACTCAGAGACATTTCCAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTGCTCGTTACaa                                                                                                    >  1:1941393/1‑150 (MQ=255)
                                              ccATCGGTCAGTACTGCACCGCCTAATGCCCCTGCACCCGCGCCGATTGCGGTGTTGCGGTCCCGTTTAGACCAGTTAGAACAGGCACTCAGAGACATTTCCAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTGCTCGTTACAAAc                                                                                                  <  2:1870199/150‑1 (MQ=255)
                                              ccATCGGTCAGTACTGCACCGCCTAATGCCCCTGCACCCGCGCCGATTGCGGTGTTGCGGTCCCGTTTAGACCAGTTAGAACAGGCACTCAGAGACATTGCCAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTGCTCGTTACAAAc                                                                                                  <  1:1715606/150‑1 (MQ=255)
                                               cATCGGTCAGTACTGCACCGCCTAATGCCCCTGCACCCGCGCCGATTGCGGTGTTGCGGTCCCTTTTAGACCAGTTAGAACAGGCACTCAGAGACATTTCCAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTGCTCGTTACAAACa                                                                                                 <  2:1290625/150‑1 (MQ=255)
                                               cATCGGTCAGTACTGCACCGCCTAATGCCCCTGCACCCGCGCCGATTGCGGTGTTGCGGTCCCGTTTAGACCAGTTAGAACAGGCACTCAGAGCCATTTCCAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTGCTCGTTACAAACa                                                                                                 <  1:1055459/150‑1 (MQ=255)
                                               cATCGGTCAGTACTGCACCGCCTAATGCCCCTGCACCCGCGCCGATTGCGGTGTTGCGGTCCCGTTTAGACCAGTTAGAACAGGCACTCAGAGACATTTCCCAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTGCTCGTTACAAACa                                                                                                 >  2:1471235/1‑150 (MQ=255)
                                               cATCGGTCAGTACTGCACCGCCTAATGCCCCTGCACCCGCGCCGATTGCGGTGTTGCGGTCCCGTTTAGACCAGTTAGAACAGGCACTCAGAGACATTTCCAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTGCTCGTTACAAACa                                                                                                 >  2:1471199/1‑150 (MQ=255)
                                               cATCGGTCAGTACTGCACCGCCTAATGCCCCTGCACCCGCGCCGATTGCGGTGTTGCGGTCCCGTTTAGACCAGTTAGAACAGGCACTCAGAGACATTTCCAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTGCTCGTTACAAACa                                                                                                 >  2:140072/1‑150 (MQ=255)
                                                aTCGGTCAGTACTGCACCGCCTAATGCCCCTGCACCCGCGCCGATTGCGGTGTTGCGGTCCCGTTTAGACCAGTTAGAACAGGCACTCAGAGACATTTCCAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTGCTCGTTACAAACat                                                                                                >  2:269427/1‑150 (MQ=255)
                                                 tCGGTCAGTACTGCACCGCCTAATGCCCCTGCACCCGCGCCGATTGCGGTGTTGCGGTCCCGTTTAGACCAGTTAGAACAGGCACTCAGAGACATTTCCAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTGCTCGTTACAAACata                                                                                               >  2:1558892/1‑150 (MQ=255)
                                                 tCGGTCAGTACTGCACCGCCTAATGCCCCTGCACCCGCGCCGATTGCGGTGTTGCGGTCCCGTTTAGACCAGTTAGAACAGGCACTCAGAGACATTTCCAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTGCTCGTTACAAACata                                                                                               <  1:2710786/150‑1 (MQ=255)
                                                  cGGTCAGTACTGCACCGCCTAATGCCCCTGCACCCGCGCCGATTGCGGTGTTGCGGTCCCTTTTAGACCAGTTAGAACAGGCACTCAGAGACATTTCCAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTGCTCGTTACAAACataa                                                                                              >  2:170538/1‑150 (MQ=255)
                                                  cGGTCAGTACTGCACCGCCTAATGCCCCTGCACCCGCGCCGATTGCGGTGTTGCGGTCCCGTTTAGACCAGTTAGAACAGGCACTCAGAGACATTTCCAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTGCTCGTTACAAACataa                                                                                              >  1:4059668/1‑150 (MQ=255)
                                                    gTCAGTACTGCACCGCCTAATGCCCCTGCACCCGCGCCGATTGCGGTGTTGCGGTCCCGTTTAGACCAGTTAGAACAGGCACTCAGAGACATTTCCAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTGCTCGTTACAAACataata                                                                                            >  2:3244052/1‑150 (MQ=255)
                                                      cAGTACTGCACCGCCTAATGCCCCTGCACCCGCGCCGATTGCGGTGTTGCGGTCCCGTTTAGACCAGTTAGAACAGGCACTCAGAGACATTTCCAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTGCTCGTTACAAACATAATAct                                                                                          >  2:3086069/1‑150 (MQ=255)
                                                      cAGTACTGCACCGCCTAATGCCCCTGCACCCGCGCCGATTGCGGTGTTGCGGTCCCGTTTAGACCAGTTAGAACAGGCACTCAGAGACATTTCCAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTGCTCGTTACAAACATAATAct                                                                                          >  2:3393711/1‑150 (MQ=255)
                                                      cAGTACTGCACCGCCTAATGCCCCTGCACCCGCGCCGATTGCGGTGTTGCGGTCCCGTTTAGACCAGTTAGAACAGGCACTCAGAGACATTTCCAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTGCTCGTTACAAACATAATAct                                                                                          <  1:824003/150‑1 (MQ=255)
                                                       aGTACTGCACCGCCTAATGCCCCTGCACCCGCGCCGATTGCGGTGTTGCGGTCCCGTTTAGACCAGTTAGAACAGGCACTCAGAGACATTTCCAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTGCTCGTTACAAACATAATActc                                                                                         >  2:437843/1‑150 (MQ=255)
                                                        gTACTGCACCGCCTAATGCCCCTGCACCCGCGCCGATTGCGGTGTTGCGGTCCCGTTTAGACCAGTTAGAACAGGCACTCAGAGACATTTCCAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTGCTCGTTACAAACATAATActct                                                                                        <  2:3714436/150‑1 (MQ=255)
                                                         tACTGCACCGCCTAATGCCCCTGCACCCGCGCCGAGTGCGGTGTTGCGGTCCCGTTTAGACCAGTTAGAACAGGCACTCAGAGACATTTCCAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTGCTCGTTACAAACATAATActctc                                                                                       <  1:3534786/150‑1 (MQ=255)
                                                          aCTGCACCGCCTAATGCCCCTGCACCCGCGCCGATTGCGGTGTTGCGGTCCCGTTTAGACCAGTTAGAACAGGCACTCAGAGACATTTCCAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTGCTCGTTACAAACATAATActctct                                                                                      <  1:2168049/150‑1 (MQ=255)
                                                           cTGCACCGCCTAATGCCCCTGCACCCGCGCCGATTGCGGTGTTGCGGTCCCGTTTAGACCAGTTAGAACAGGCACTCAGAGACATTTCCAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTGCTCGTTACAAACATAATACtctctc                                                                                     <  1:4096364/150‑1 (MQ=255)
                                                             gCACCGCCTAATGCCCCTGCACCCGCGCCGATTGCGGTGTTGCGGTCCCGTTTAGACCAGTTAGAACAGGCACTCAGAGACATTTCCAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTGCTCGTTACAAACATAATACTCTCTCCt                                                                                   >  1:1203205/1‑150 (MQ=255)
                                                                ccGCCTAATGCCCCTGCACCCGCGCCGATTGCGGTGTTGCGGTCCCGTTTAGACCAGTTAGAACAGGCACTCAGAGACATTTCCAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTGCTCGTTACAAACATAATACTCTCTCCTGaa                                                                                <  1:2663296/150‑1 (MQ=255)
                                                                  gCCTAATGCCCCTGCACCCGCGCCGATTGCGGTGTTGCGGTCCCGTTTAGACCAGTTAGAACAGGCACTCAGAGACATTTCCAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTGCTCGTTACAAACATAATACTCTCTCCTGAAtt                                                                              <  2:2425939/150‑1 (MQ=255)
                                                                    cTAATGCCCCTGCACCCGCGCCGATTGCGGTGTTGCGGTCCCGTTTAGACCAGTTAGAACAGGCACTCAGAGACATTTCCAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTGCTCGTTACAAACATAATACTCTCTCCTGAATTTa                                                                            >  1:3794171/1‑150 (MQ=255)
                                                                    cTAATGCCCCTGCACCCGCGCCGATTGCGGTGTTGCGGTCCCGTTTAGACCAGTTAGAACAGGCACTCAGAGACATTTCCAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTGCTCGTTACAAACATAATACTCTCTCCTGAATGTa                                                                            >  2:1049793/1‑150 (MQ=255)
                                                                     tAATGCCCCTGCACCCGCGCCGATTGCGGTGTTGCGGTCCCGTTTAGACCAGTTAGAACAGGCACTCAGAGACATTTCCAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTGCTCGTTACAAACATAATACTCTCTCCTGAATTTAt                                                                           >  1:3611495/1‑150 (MQ=255)
                                                                         gCCCCTGCACCCGCGCCGATTGCGGTGTTGCGGTCCCGTTTAGACCAGTTAGAACAGGCACTCAGAGACATTTCCAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTGCTCGTTACAAACATAATACTCTCTCCTGAATTTATGAtt                                                                       <  1:3265745/150‑1 (MQ=255)
                                                                           cccTGCACCCGCGCCGATTGCGGTGTTGCGGTCCCGTTTAGACCAGTTAGAACAGGCACTCAGAGACATTTCCAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTGCTCGTTACAAACATAATACTCTCTCCTGAATTTATGATTCa                                                                     <  2:2739182/150‑1 (MQ=255)
                                                                             cTGCACCCGCGCCGATTGCGGTGGTGCGGTCCCGTTTAGACCAGTTAGAACAGGCACTCAGAGACATTTCCAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTGCTCGTTACAAACATAATACTCTCTCCTGAATTTATGATTCACg                                                                   <  2:82452/150‑1 (MQ=255)
                                                                              tGCACCCGCGCCGATTGCGGTGTTGCGGTCCCGTTTAGACCAGTTAGAACAGGCACTCAGAGACATTTCCAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTGCTCGTTACAAACATAATACTCTCTCCTGAATTTATGATTCACgg                                                                  >  2:3314264/1‑150 (MQ=255)
                                                                               gCACCCGCGCCGATTGCGGTGTTGCGGTCCCGTTTAGACCAGTTAGAACAGGCACTCAGAGACATTGCCAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTGCTCGTTACAAACATAATACTCTCTCCTGAATTTATGATTCACGGa                                                                 <  2:30275/150‑1 (MQ=255)
                                                                                 aCCCGCGCCGATTGCGGTGTTGCGGTCCCGTTTAGACCAGTTAGAACAGGCACTCAGAGACATTTCCAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTGCTCGTTACAAACATAATACTCTCTCCTGAATTTATGATTCACGGAAg                                                               <  1:2211147/150‑1 (MQ=255)
                                                                                 aCCCGCGCCGATTGCGGTGTTGCGGTCCCGTTTAGACCAGTTAGAACAGGCACTCAGAGACATTTCCAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTGCTCGTTACAAACATAATACTCTCTCCTGAATTTATGATTCACGGAAg                                                               <  2:2712179/150‑1 (MQ=255)
                                                                                   ccGCGCCGATTGCGGTGTTGCGGTCCCGTTTAGACCAGTTAGAACAGGCACTCAGAGACATTTCCAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTGCTCGTTACAAACATAATACTCTCTCCTGAATTTATGATTCACGGAAGTa                                                             <  1:3946219/150‑1 (MQ=255)
                                                                                         cGATTGCGGTGTTGCGGTCCCGTTTAGACCAGTTAGAACAGGCACTCAGAGACATTTCCAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTGCTCGTTACAAACATAATACTCTCTCCTGAATTTATGATTCACGGAAGTAAGctct                                                       >  1:112381/1‑150 (MQ=255)
                                                                                            ttGCGGTGTTGCGGTCCCGTTTAGACCAGTTAGAACAGGCACTCAGAGACATTTCCAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTGCTCGTTACAAACATAATACTCTCTCCTGAATTTATGATTCACGGAAGTAAGCTCTcta                                                    <  2:3887121/150‑1 (MQ=255)
                                                                                             tGCGGTGTTGCGGTCCCGTTTAGACCAGTTAGAACAGGCACTCAGAGACATTTCCAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTGCTCGTTACAAACATAATACTCTCTCCTGAATTTATGATTCACGGAAGTAAGCTCTctat                                                   >  1:3632379/1‑150 (MQ=255)
                                                                                              gCGGTGTTGCGGTCCCGTTTAGACCAGTTAGAACAGGCACTCAGAGACATTTCCAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTGCTCGTTACAAACATAATACTCTCTCCTGAATTTATGATTCACGGAAGTAAGCTCTctata                                                  <  1:1891556/150‑1 (MQ=255)
                                                                                                     tGCGGTCCCGTTTAGACCAGTTAGAACAGGCACTCAGAGACATTTCCAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTGCTCGTTACAAACATAATACTCTCTCCTGAATTTATGATTCACGGAAGTAAGCTCTctataactataa                                           <  1:2023796/150‑1 (MQ=255)
                                                                                                      gCGGTCCCGTTTAGACCAGTTAGAACAGGCACTCAGAGACATTTCCAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTGCTCGTTACAAACATAATACTCTCTCCTGAATTTATGATTCACGGAAGTAAGCTCTCTATAACTataat                                          <  2:1607497/150‑1 (MQ=255)
                                                                                                      gCGGTCCCGTTTAGACCAGTTAGAACAGGCACTCAGAGACATTTCCAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTGCTCGTTACAAACATAATACTCTCTCCTGAATTTATGATTCACGGAAGTAAGCTCTCTATAACTataat                                          <  1:3725382/150‑1 (MQ=255)
                                                                                                       cGGTCCCGTTTAGACCAGTTAGAACAGGCACTCAGAGACATTTCCAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTGCTCGTTACAAACATAATACTCTCTCCTGAATTTATGATTCACGGAAGTAAGCTCTCTATAACTataata                                         <  2:3708852/150‑1 (MQ=255)
                                                                                                       cGGTCCCGTTTAGACCAGTTAGAACAGGCACTCAGAGACATTTCCAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTGCTCGTTACAAACATAATACTCTCTCCTGAATTTATGATTCACGGAAGTAAGCTCTCTATAACTataata                                         <  1:327306/150‑1 (MQ=255)
                                                                                                         gTCCCGTTTAGACCAGTTAGAACAGGCACTCAGAGACATTTCCAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTGCTCGTTACAAACATAATACTCTCTCCTGAATTTATGATTCACGGAAGTAAGCTCTCTAAAACTATAAAAGc                                       >  2:181736/1‑150 (MQ=255)
                                                                                                           cccGTTTAGACCAGTTAGAACAGGCACTCAGAGACATTTCCAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTGCTCGTTACAAACATAATACTCTCTCCTGAATTTATGATTCACGGAAGTAAGCTCTCTATAACTATAATAGCTa                                     <  2:2258152/150‑1 (MQ=255)
                                                                                                            ccGTTTAGACCAGTTAGAACAGGCACTCAGAGACATTTCCAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTGCTCGTTACAAACATAATACTCTCTCCTGAATTTATGATTCACGGAAGTAAGCTCTCTATAACTATAATAGCTaa                                    >  1:2222315/1‑150 (MQ=255)
                                                                                                             cGTTTAGACCAGTTAGAACAGGCACTCAGAGACATTTCCAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTGCTCGTTACAAACATAATACTCTCTCCTGAATTTATGATTCACGGAAGTAAGCTCTCTATAACTATAATAGCTAAg                                   >  1:1651880/1‑150 (MQ=255)
                                                                                                             cGTTTAGACCAGTTAGAACAGGCACTCAGAGACATTTCCAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTGCTCGTTACAAACATAATACTCTCTCCTGAATTTATGATTCACGGAAGTAAGCTCTCTATAACTATAATAGCTAAg                                   >  2:3987616/1‑150 (MQ=255)
                                                                                                                  aGACCAGTTAGAACAGGCACTCAGAGACATTTCCAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTGCTCGTTACAAACATAATACTCTCTCCTGAATTTATGATTCACGGAAGTAAGCTCTCTATAACTATAATAGCTAAGAATaa                              <  2:4052558/150‑1 (MQ=255)
                                                                                                                  aGACCAGTTAGAACAGGCACTCAGAGACATTTCCAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTGCTCGTTACAAACATAATACTCTCTCCTGAATTTATGATTCACGGAAGTAAGCTCTCTATAACTATAATAGCTAAGAATaa                              <  1:3169317/150‑1 (MQ=255)
                                                                                                                   gACCAGTTAGAACAGGCACTCAGAGACATTTCCAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTGCTCGTTACAAACATAATACTCTCTCCTGAATTTATGATTCACGGAAGTAAGCTCTCTATAACTATAATAGCTAAGAATAAg                             <  1:2643753/150‑1 (MQ=255)
                                                                                                                     ccAGTTAGAACAGGCACTCAGAGACATTTCCAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTGCTCGTTACAAACATAATACTCTCTCCTGAATTTATGATTCACGGAAGTAAGCTCTCTATAACTATAATAGCTAAGAATAAGTc                           <  2:604805/150‑1 (MQ=255)
                                                                                                                          tAGAACAGGCACTCAGAGACATTTCCAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTGCTCGTTACAAACATAATACTCTCTCCTGAATTTATGATTCACGGAAGTAAGCTCTCTATAACTATAATAGCTAAGAATAAGTCTGGTg                      >  1:2056311/1‑150 (MQ=255)
                                                                                                                          tAGAACAGGCACTCAGAGACATTTCCAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTGCTCGTTACAAACATAATACTCTCTCCTGAATTTATGATTCACGGAAGTAAGCTCTCTATAACTATAATAGCTAAGAATAAGTCTGGTg                      <  2:3003349/150‑1 (MQ=255)
                                                                                                                           aGAACAGGCACTCAGAGACATTTCCAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTGCTCGTTACAAACATAATACTCTCTCCTGAATTTATGATTCACGGAAGTAAGCTCTCTATAACTATAATAGCTAAGAATAAGTCTGGTGa                     >  1:786621/1‑150 (MQ=255)
                                                                                                                           aGAACAGGCACTCAGAGACATTTCCAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTGCTCGTTACAAACATAATACTCTCTCCTGAATTTATGATTCACGGAAGTAAGCTCTCTATAACTATAATAGCTAAGAATAAGTCTGGTGa                     >  2:2364682/1‑150 (MQ=255)
                                                                                                                             aaCAGGCACTCAGAGACATTTCCAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTGCTCGTTACAAACATAATACTCTCTCCTGAATTTATGATTCACGGAAGTAAGCTCTCTATAACTATAATAGCTAAGAATAAGTCTGGTGaaa                   <  2:3188801/150‑1 (MQ=255)
                                                                                                                               cAGGCACTCAGAGACATTTCCAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTGCTCGTTACAAACATAATACTCTCTCCTGAATTTATGATTCACGGAAGTAAGCTCTCTATAACTATAATAGCTAAGAATAAGTCTGGTGAAatt                 >  1:1579932/1‑150 (MQ=255)
                                                                                                                                aGGCACTCAGAGACATTTCCAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTGCTCGTTACAAACATAATACTCTCTCCTGAATTTATGATTCACGGAAGTAAGCTCTCTATAACTATAATAGCTAAGAATAAGTCTGGTGAAatta                <  1:4051757/150‑1 (MQ=255)
                                                                                                                                 ggCACTCAGAGACATTTCCAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTGCTCGTTACAAGCGTAATACTCTCTCCTGAATTTATGATTCACGGAAGTAAGCTCTCTATAACTATAATAGCTAAGAATAAGTCTGGTGAAattat               >  2:1745529/1‑150 (MQ=255)
                                                                                                                                  gCACTCAGAGACATTTCCAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTGCTCGTTACAAACATAATACTCTCTCCTGAATTTATGATTCACGGAAGTAAGCTCTCTATAACTATAATAGCTAAGAATAAGTCTGGTGAAattatt              >  2:481715/1‑150 (MQ=255)
                                                                                                                                  gCACTCAGAGACATTTCCAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTGCTCGTTACAAACATAATACTCTCTCCTGAATTTATGATTCACGGAAGTAAGCTCTCTATAACTATAATAGCTAAGAATAAGTCTGGTGAAattatt              <  2:3663064/150‑1 (MQ=255)
                                                                                                                                   cACTCAGAGACATTTCCAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTGCTCGTTACAAACATAATACTCTCTCCTGAATTTATGATTCACGGAAGTAAGCTCTCTATAACTATAATAGCTAAGAATAAGTCTGGTGAAATTATTc             >  1:1361963/1‑150 (MQ=255)
                                                                                                                                    aCTCAGAGACATTTCCAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTGCTCGTTACAAACATAATACTCTCTCCTGAATTTATGATTCACGGAAGTAAGCTCTCTATAACTATAATAGCTAAGAATAAGTCTGGTGAAATTATTcc            >  2:1779669/1‑150 (MQ=255)
                                                                                                                                    aCTCAGAGACATTTCCAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTGCTCGTTACAAACATAATACTCTCTCCTGAATTTATGATTCACGGAAGTAAGCTCTCTATAACTATAATAGCTAAGAATAAGTCTGGTGAAATTATTcc            <  1:366523/150‑1 (MQ=255)
                                                                                                                                    aCTCAGAGACATTTCCAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTGCTCGTTACAAACATAATACTCTCTCCTGAATTTATGATTCACGGAAGTAAGCTCTCTATAACTATAATAGCTAAGAATAAGTCTGGTGAAATTATTcc            >  2:1977768/1‑150 (MQ=255)
                                                                                                                                            aCATTTCCAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTGCTCGTTACAAACATAATACTCTCTCCTGAATTTATGATTCACGGAAGTAAGCTCTCTATAACTATAATAGCTAAGAATAAGTCTGGTGAAATTATTCCGTGAAATc    >  2:3114161/1‑150 (MQ=255)
                                                                                                                                             cATTTCCAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTGCTCGTTACAAACATAATACTCTCTCCTGAATTTATGATTCACGGAAGTAAGCTCTCTATAACTATAATAGCTAAGAATAAGTCTGGTGAAATTATTCCGTGAAATCt   <  1:2955138/150‑1 (MQ=255)
                                                                                                                                             cATTTCCAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTGCTCGTTACAAACATAATACTCTCTCCTGAATTTATGATTCACGGAAGTAAGCTCTCTATAACTATAATAGCTAAGAATAAGTCTGGTGAAATTATTCCGTGAAATCt   <  2:6630/150‑1 (MQ=255)
                                                                                                                                              aTTTCCAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTGCTCGTTACAAACATAATACTCTCTCCTGAATTTATGATTCACGGAAGTAAGCTCTCTATAACTATAATAGCTAAGAATAAGTCTGGTGAAATTATTCCGTGAAATCTg  <  1:2548969/150‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                  
AATAACACCTCCGACAGCTGCACCACCTAATGTACCCAACGTACTGCCATCGGTCAGTACTGCACCGCCTAATGCCCCTGCACCCGCGCCGATTGCGGTGTTGCGGTCCCGTTTAGACCAGTTAGAACAGGCACTCAGAGACATTGCCAAAGTAATTGCCAGAACAGCCGCGGTCATTTTTTTGCTCGTTACAAACATAATACTCTCTCCTGAATTTATGATTCACGGAAGTAAGCTCTCTATAACTATAATAGCTAAGAATAAGTCTGGTGAAATTATTCCGTGAAATCTG  >  CP009273/1337388‑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: ‑