Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,836,179 G→A 25.4% Q72* (CAG→TAG)  ynjH ← DUF1496 family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,836,1790GA25.4% 527.6 / 185.9 268Q72* (CAG→TAG) ynjHDUF1496 family protein
Reads supporting (aligned to +/- strand):  ref base G (87/113);  new base A (33/35);  total (120/148)
Fisher's exact test for biased strand distribution p-value = 4.84e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.93e-01

Only 100 of 269 total aligned reads displayed.

CATGAACATCAGGCGCTGGTCTGGTGCTCACCTGAAGAGGCGCTGCAATATCCGCTGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTA‑TGGCTTTACGCGCCGCCAGACCAGCGGATTAGTGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGTATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAATTTTGATCCTGATAAACGCAGCACTGAGTACAGGGCTGCGCTGACTGCCCGCCTGAACTGAACACTTCTGGCGGTACGTTTACCTCCACGTCCGGACGAT  >  CP009273/1836032‑1836328
                                                                                                                                                    |                                                                                                                                                     
  tGAACATCAGGCGCTGGTCAGGTGCTCACCTGAAGAGGCGCTGCAATATCCGCTGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTA‑TTGCTTTACCAGCACCCAGACCCGCGGAGTAGTGCAAAGGGTTTTGTAATCACGCTAGCAt                                                                                                                                                   >  2:4557306/1‑150 (MQ=255)
   gAACATCAGGCGCTGGTCTGGTGCTCACCTGAAGAGGCGCTGCAATATCCGCTGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTA‑TGGCTTTACGCGCCGCCAGACCAGCGGATTAGTGCTAAGGGTTTTGTCATCACGCTAGCAtt                                                                                                                                                  <  2:779517/150‑1 (MQ=255)
     aCATCAGGCGCTGGTCTGGTGCTCACCTGAAGAGGCGCTGCAATATCCGCTGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTA‑TGGCTTTACGCGCCGCCAGACCAGCGGATTAGTGCTAAGGGTTTTGTCATCACGCTGGCATTgc                                                                                                                                                <  1:3429413/150‑1 (MQ=255)
      cATCAGGCGCTGGTCTGGTGCTCACCTGAAGAGGCGCTGCAATATCCGCTGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTA‑TGGCTTTACGCGCCGCCAGACCAGCGGATTAGTGCTAAGGGTTTTGTCATCACGCTGGCATTgca                                                                                                                                               <  2:1785416/150‑1 (MQ=255)
      cATCAGGCGCTGGTCTGGTGCTCACCTGAAGAGGCGCTGCAATATCCGCTGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTA‑TGGCTTTACGCGCCGCCAGACCAGCGGATTAGTGCTAAGGGTTTTGTCATCACGCTAGCATTgca                                                                                                                                               <  2:5158155/150‑1 (MQ=255)
         cAGGCGCTGGTCTGGTGCTCACCTGAAGAGGCGCTGCAATATCCGCTGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTA‑TGGCTTTACGCGCCGCCAGACCAGCGGATTAGTGCTAAGGGTTTTGTCATCACGCTGGCATTgcagca                                                                                                                                            <  2:4018679/150‑1 (MQ=255)
          aGGCGCTGGTCTGGTGCTCACCTGAAGAGGCGCTGCAATATCCGCTGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTA‑TTGGTTTAAGCGCCGCCCGACCAGCGGATTAGTGCTAAGGGTTTTGTCATCACGCTGGCCTTGcagcag                                                                                                                                           >  1:5045467/1‑150 (MQ=255)
          aGGCGCTGGTCTGGTGCTCACCTGAAGAGGCGCTGCAATATCCGCTGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTA‑TGGCTTTACGCGCCGCCAGACCAGCGGATTATTGCGAAGGGTTTTGTCATCACACAGGCATTGcagcag                                                                                                                                           >  1:951862/1‑150 (MQ=255)
          aGGCGCTGGTCTGGTGCTCACCTGAAGAGGCGCTGCAATATCCGCTGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTA‑TGGCTTTACGCGCCGCCAGACCAGCGGATTAGTTCTAAGGGTTTTGTCATCACACTGGCATTGcagcag                                                                                                                                           >  1:3848858/1‑150 (MQ=255)
             cgcTGGTCTGGTGCTCACCTGAAGAGGCGCTGCAATATCCGCTGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTA‑TGGCTTTACGCGCCGCCAGACCAGCGGATTAGTGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGAAt                                                                                                                                        <  2:852500/150‑1 (MQ=255)
             cgcTGGTCTGGTGCTCACCTGAAGAGGCGCTGCAATATCCGCTGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTA‑TGGCTTTACGCGCCGCCAGACCAGCGGATTAGTGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGAAt                                                                                                                                        <  2:3751901/150‑1 (MQ=255)
              gcTGGTCTGGTGCTCACCTGAAGAGGCGCTGCAATATCCGCTGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTA‑TGGCTTTACGCGCCGCCAGACCAGCGGATTAGTGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGAAtt                                                                                                                                       <  1:4660491/150‑1 (MQ=255)
              gcTGGTCTGGTGCTCACCTGAAGAGGCGCTGCAATATCCGCTGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTA‑TGGCTTTACGCGCCGCCAGACCAGCGGATTAGTGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGAAtt                                                                                                                                       >  2:3899780/1‑150 (MQ=255)
                tggtctggtGCTCACCTGAAGAGGCGCTGCAATATCCGCTGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTA‑TGGCTTTACGCGCCGCCAGACCAGCGGATTAGTGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGAATTcc                                                                                                                                     <  2:4225174/150‑1 (MQ=255)
                 ggtctggtGCTCACCTGAAGAGGCGCTGCAATATCCGCTGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTA‑TTGCTTTACGCGCCGCCCGACCAGCGGGTTAGTGCTAAGGGTTTTTTCATCACGCTGGCATTTcagcagaaataca                                                                                                                                    >  1:3208871/1‑143 (MQ=255)
                 ggtctggtGCTCACCTGAAGAGGCGCTGCAATATCCGCTGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTA‑TGGCTTTACGCGCCGCCAGACCAGCGGATTAGTGCTAAGGGTTTTGTCATCACGCTGGCAATGCAGCAGAATTCCt                                                                                                                                    >  1:2202293/1‑150 (MQ=255)
                 ggtctggtGCTCACCTGAAGAGGCGCTGCAATATCCGCTGGCCACTGCTGACATTCCATTATTAGAGGCGTTTA‑TGGCTTTACGCGCCGCCAGACCAGCGGATTAGTGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGTATTCCt                                                                                                                                    <  1:1806516/150‑1 (MQ=255)
                     tggtGCTCACCTGAAGAGGCGCTGCAATATCCGCTGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTATTGGGTTTA‑GCGCCGCCACACCAGCGGATTATTGCTAAGGGGTTTTTCATAACGCTGGCCTTGCAGCAGAAttacatagg                                                                                                                                >  1:1900387/1‑143 (MQ=255)
                     tggtGCTCACCTGAAGAGGCGCTGCAATATCCGCTGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTA‑TGGCTTTACGCGCCGCCAGACCAGCGGATTAGTGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGAATTCCTTCgg                                                                                                                                <  1:1083175/150‑1 (MQ=255)
                       gtGCTCACCTGAAGAGGCGCTGCAATATCCGCTGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTA‑TGGCTTTACGCGCCGCCAGACCAGCGGATTAGTGCTAAGGGTTTTGTCATCACGCTAGCATTGCAGCAGAATTCCTTCGGCt                                                                                                                              <  2:2495778/150‑1 (MQ=255)
                       gtGCTCACCTGAAGAGGCGCTGCAATATCCGCTGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTA‑TGGCTTTACGCGCCGCCAGACCAGCGGATTAGTGCTAAGGGGTTTGTCATCACGCTAGCATTGCAGCAGAATTCCTTCGGCt                                                                                                                              >  2:2768995/1‑150 (MQ=255)
                          cTCACCTGAAGAGGCGCTGCAATATCCGCTGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTT‑TTGCTTTAAGCGCCGCCAGACCAGCGGATTAGTGCTAAGGGTTTTTTAATCACGCTAGCATTGCAGCAGAAATGCTTCGGGTTTa                                                                                                                           >  1:1465724/1‑150 (MQ=255)
                          cTCACCTGAAGAGGCGCTGCAATATCCGCTGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTA‑TGGCTTTACGCGCCGCCAGACCAGCGGATTAGTGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGAATTCCTTCGGCTTTa                                                                                                                           <  1:4998282/150‑1 (MQ=255)
                          cTCACCTGAAGAGGCGCTGCAATATCCGCTGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTA‑TGGCTTTACGCGCCGCCAGACCAGCGGATTAGTGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGAATTCCTTCGGCTTTa                                                                                                                           <  2:689075/150‑1 (MQ=255)
                                  aaGAGGCGCTGCAATATCCGCTGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTA‑TGGCTTTACGCGCCGCCAGACCAGCGGATTAGTGCTAAGGGTTTTGTAATCACGCTAGCATTGCAGCAGAATTACTTAGGCTTTAATTACCGc                                                                                                                   >  2:1312774/1‑150 (MQ=255)
                                   agagGCGCTGCAATATCCGCTGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTA‑TGGCTTTACGCGCCGCCAGACCAGCGGATTAGTGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGAATTCCTTCGGCGTTAATTACCGcc                                                                                                                  >  1:3055716/1‑150 (MQ=255)
                                       gcgcTGCAATATCCGCTGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTA‑TGGCTTTACGCGCCGCCAGACCAGCGGATTAGTGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGAATTCCTTCGGCTTTAATTACCGCCCCtt                                                                                                              <  2:2623592/150‑1 (MQ=255)
                                        cgcTGCAATATCCGCTGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTA‑TGGCTTTACGCGCCGCCAGACCAGCGGATTAGTGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGAATTCCTTCGGCTTTAATTACCGCCCCTTc                                                                                                             <  2:3346499/150‑1 (MQ=255)
                                        cgcTGCAATATCCGCTGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTA‑TGGCTTTACGCGCCGCCAGACCAGCGGATTAGTGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGAATTCCTTCGGCTTTAATTACCGCCCCTTc                                                                                                             >  2:4759703/1‑150 (MQ=255)
                                         gcTGCAATATCCGCTGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTA‑TGGCTTTACGCGCCGCCAGACCAGCGGATTAGTGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGAATTCCTTCGGCTTTAATTACCGCCCCTTCa                                                                                                            <  2:2289093/150‑1 (MQ=255)
                                         gcTGCAATATCCGCTGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTA‑TGGCTTTACGCGCCGCCAGACCAGCGGATTAGTGCTAAGGGTTTTGTCATCACGCTAGCATTGCAGCAGAATTCCTTCGGCTTTAATTACCGCCCCTTCa                                                                                                            <  2:2544618/150‑1 (MQ=255)
                                          cTGCAATATCCGCTGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTA‑TGGCTTTACGCGCCGCCAGACCAGCGGATTAGTGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGAATTCCTTCGGCTTTAATTACCGCCCCTTCAg                                                                                                           <  2:746611/150‑1 (MQ=255)
                                           tGCAATATCCGCTGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTA‑TGGCTTTACGCGCCGCCAGACCAGCGGATTAGTGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGAATTCCTTCGGCTTTAATTACCGCCCCTTCAGa                                                                                                          <  2:1823683/150‑1 (MQ=255)
                                            tcAATATACTCTGGCCCCTGCTGACAGTCCATTATTAGAGGCGTTTA‑TGTCTTTACGCTCCGCCAGACCAGCGGATTAGTGCTAAGGGTTTTGTCATCACGCTAGCATTGCAGCAGAATTCCTTCGGCTTTAATTACCGCCCCTTCAGaa                                                                                                         <  2:4497992/149‑1 (MQ=255)
                                            gCAATATCCGCTGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTA‑TGGCTTTACGCGCCGCCAGACCAGCGGATTAGTGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGAATTCCTTCGGCTTTAATTACCGCCCCATCAGaa                                                                                                         >  1:4340256/1‑150 (MQ=255)
                                             cAATATCCGCTGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTA‑TGGCTTTACGCGCCGCCAGACCAGCGGATTAGTGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGAATTCCTTCGGCTTTAATTACCGCCCCTTCAGaat                                                                                                        <  2:914263/150‑1 (MQ=255)
                                                   ccGCTGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTA‑TGGCTTTACGCGCCGCCAGACCAGCGGATTAGTGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGAATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAAtttt                                                                                                  <  1:2349673/150‑1 (MQ=255)
                                                     gCTGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTA‑TGGCTTTACGCGCCGCCAGACCAGCGGATTAGTGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGAATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAATTTTGa                                                                                                <  2:1381537/150‑1 (MQ=255)
                                                      cTGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTA‑TGGCTTTACGCGCCGCCAGACCAGCGGATTAGTGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGAATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAATTTTGAt                                                                                               <  2:2604575/150‑1 (MQ=255)
                                                       tGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTA‑TGGCTTTACGCGCCGCCAGACCAGCGGATTAGTGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGAATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAATTTTGATc                                                                                              <  1:4685401/150‑1 (MQ=255)
                                                       tGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTA‑TGGCTTTACGCGCCGCCAGACCAGCGGATTAGTGCTAAGGGTTTTGTCATCACGCTAGCATTGCAGCAGAATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAATTTTGATc                                                                                              <  2:2758751/150‑1 (MQ=255)
                                                         ggcccTGCTGACATTCCATTATTAGAGGCGTTTA‑TGGCTTTACGCGCCGCCAGACCAGCGGATTAGTGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGAATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAATTTTGATCCt                                                                                            <  1:4359948/148‑1 (MQ=255)
                                                         gCCCCTGCTGACATTCCATTATTAGAGGCGTTTA‑TGGCTTTACGCGCCGCCAGACCAGCGGATTAGTGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGAATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAATTTTGATCCt                                                                                            >  2:323907/1‑150 (MQ=255)
                                                          ccccTGCTGACATTCCATTATTAGAGGCGTTTA‑TGGCTTTACGCGCCGCCAGACCAGCGGATTAGTGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGAATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAATTTTGATCCTg                                                                                           >  2:5290452/1‑150 (MQ=255)
                                                             ctgctgACATTCCATTATTAGAGGCGTTTA‑TGGCTTTACGCGCCGCCAGACCAGCGGATTAGTGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGAATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAATTTTGATCCTGATa                                                                                        <  2:3984401/150‑1 (MQ=255)
                                                             ctgctgACATTCCATTATTAGAGGCGTTTA‑TGGCTTTACGCGCCCCCAGACCACCGGATTAGTGCAAAGGGTTTTGTCATCACGCAAGCATTGCAGCAGAATTCCTTCGGCTTTAATTACCGCCCCTTCAGAAAAATTTTGATCCTGATa                                                                                        >  2:5173576/1‑150 (MQ=255)
                                                              tgctgACATTCCATTATTAGAGGCGTTTA‑TGGCTTTACGCGCCGCCAGACCAGCGGATTAGTGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGAATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAATTTTGATCCTGATaa                                                                                       <  1:3502586/150‑1 (MQ=255)
                                                                 tgACATTCCATTATTAGAGGCGTTTA‑TGGCTTTACGCGCCGCCAGACCAGCGGATTAGTGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGAATTCCTTCGGCTTTAATTACCGCCCCATCAGAATAATTTTGATCCTGATAAACg                                                                                    >  2:2083227/1‑150 (MQ=255)
                                                                    cATTCCATTATTAGAGGCGTTTA‑TGGCTTTACGCGCCGCCAGACCAGCGGATTAGTGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGAATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAAGTTTGATCCTGATAAACgcag                                                                                 >  2:3843001/1‑150 (MQ=255)
                                                                        ccATTATTAGAGGCGTTTA‑TGGCTTTACGCGCCGCCAGACCAGCGGATTAGTGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGAATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAATTTTGATCCTGATAAACGCAGCACt                                                                             <  2:4684396/150‑1 (MQ=255)
                                                                        ccATTATTAGAGGCGTTTA‑TGGCTTTACGCGCCGCCAGACCAGCGGATTAGTGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGAATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAATTTTGATCCTGATAAACGCAGCACt                                                                             <  1:3708158/150‑1 (MQ=255)
                                                                         cATTATTAGAGGCGTTTA‑TGGCTTTACGCGCCGCCAGACCAGCGGATTAGTGCTAAGGGTTTTGTCATCACGCTGGCAATGCAGCAGAATTACTTCGGCATTAATTAACGCCCCATCAGAATAAATTTGATACAGATAAACGCAGCACTg                                                                            >  2:2419960/1‑150 (MQ=255)
                                                                          attattAGAGGCGTTTA‑TGGCTTTACGCGCCGCCAGACCAGCGGATTAGTGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGAATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAATTTTGATCCTGATAAACGCAGCACTGa                                                                           <  1:1159484/150‑1 (MQ=255)
                                                                            tattaGAGGCGTTTA‑TGGCTTTACGCGCCGCCAGACCAGCGGATTAGTGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGAATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAATTTTGATCCTGATAAACGCAGCACTGAGt                                                                         <  2:3542698/150‑1 (MQ=255)
                                                                              ttaGAGGCGTTTA‑TGGCTTTACGCGCCGCCAGACCAGCGGATTAGTGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGAATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAATTTTGATCCTGATAAACGCAGCACTGAGttc                                                                       <  1:72222/150‑3 (MQ=255)
                                                                               taGAGGCGTTTA‑TGGCTTTACGCGCCGCCAGACCAGCGGATTAGTGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGAATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAATTTTGATCCTGATAAACGCAGCACTGAGTACa                                                                      <  2:2911834/150‑1 (MQ=255)
                                                                                  agGCGTTTA‑TGGCTTTACGCGCCGCCAGACCAGCGGATTAGTGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGAATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAATTTTGATCCTGATAAACGCAGCACTGAGTACAggg                                                                   <  1:482149/150‑1 (MQ=255)
                                                                                     cGTTTA‑TGGCTTTACGCGCCGCCAGACCAGCGGATTAGTGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGAATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAATTTTGATCCTGATAAACGCAGCACTGAGTACAGGGCTg                                                                >  2:1879606/1‑150 (MQ=255)
                                                                                      gTTTA‑TGGCTTTACGCGCCGCCAGACCAGAGGATTAGTGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGAATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAATTTTGATCCTGATAAACGCAGCACTGAGTACAGGGCTgc                                                               >  1:4719884/1‑150 (MQ=255)
                                                                                       tttA‑TGGCTTTACGCGCCGCCAGACCAGCGGATTAGTGCTAAGGGTTTTGTCATCACGCTAGCATTGCAGCAGAATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAATTTTGATCCTGATAAACGCAGCACTGAGTACAGGGCTgcg                                                              <  1:5083349/150‑1 (MQ=255)
                                                                                         tA‑TGGCTTTACGCGCCGCCAGACCAGCGGATTAGTGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGAATTCCTTCGGCTTTAATTACCGCCCCCTCAGAATAATTTTGATCCTGATAAACGCAGCACTGAGTACAGGGCTGCGCt                                                            <  1:3883418/150‑1 (MQ=255)
                                                                                          a‑TGGCTTTACGCGCCGCCAGACCAGCGGATTAGTGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGAATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAATTTTGATCCTGATAAACGCGGCACTGAGTACAGGGCTGCGCTg                                                           >  2:309166/1‑150 (MQ=255)
                                                                                              gCTTTACGCGCCGCCAGACCAGCGGATTAGTGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGAATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAATTTTGATCCTGATAAACGCAGCACTGAGTACAGGGCTGCGCTGACt                                                        <  1:197368/150‑1 (MQ=255)
                                                                                               cTTTACGCGCCGCCAGACCAGCGGATTAGTGCTAAGGGTTTTGTCATCACGCTAGCATTGCAGCAGAATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAATTTTGATCCTGATAAACGCAGCACTGAGTACAGGGCTGCGCTGACTg                                                       >  2:3002969/1‑150 (MQ=255)
                                                                                                 ttACGCGCCGCCAGACCAGCGGATTAGTGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGAATTCCTTCGGCTTTAATTACCGCCCCATCAGAATAATTTTGATCCTGATAAACGAAGAACTGAGTACAGGGCTGCGCTGACTGcc                                                     >  1:1302527/1‑150 (MQ=255)
                                                                                                  tACGCGCCGCCAGACCAGCGGATTAGTGCTATGGGTTTTGTCATCACGCTGGCATTGCAGCAGAATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAATTTTGATCCTGATAAACGCAGCACTGAGTACAGGGCTGCGCTGACTGccc                                                    <  2:4992634/150‑1 (MQ=255)
                                                                                                  tACGCGCCGCCAGACCAGCGGATTAGTGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGAATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAATTTTGATCCTGATAAACGCAGCACTGAGTACAGGGCTGCGCTGACTGccc                                                    <  2:1221815/150‑1 (MQ=255)
                                                                                                  tACGCGCCGCCAGACCAGCGGATTAGTGCTAAGGGTTTTGTCATCACGCTAGCATTGCAGCAGAATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAATTTTGATCCTGATAAACGCAGCACTGAGTACAGGGCTGCGCTGACTGccc                                                    >  1:4918781/1‑150 (MQ=255)
                                                                                                  tACGCGCAGCCAGACAAGCGGATTAGTGCTAAGGGTTTTGTCATCACGCTGGAATTGCATCAGAATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAATTTTGATCCTGATAAACGCAGCACTGAGTACAGGGCTGCGCTGACTGccc                                                    <  1:3312095/150‑1 (MQ=255)
                                                                                                      cgccgcCAGACCTGCGTATTAGTGCTAAGGGTTTTGTCATCACGCTGGAATTGCAGCAGAATTCCTTCGGCTTTAATTACCGCCCATTCATAATAATTTTGATCCTTATAAACGCAGCACTGAGTACAGGGCTGCGCTGACTGCCCGCct                                                <  1:4100016/150‑1 (MQ=255)
                                                                                                            cAGACCAGCGGATTAGTGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGAATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAATTTTGATCCTGATAAACGCAGCACTGAGTACAGGGCTGCGCTGACTGCCCGCctgaactg                                          <  2:1918881/150‑1 (MQ=255)
                                                                                                               aCCCGCGGATTAGTGCTATGGTTTTTGTCATCACGCTGGCATTGCAGCAGAATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAATTTTGATCCTTATAAACGCAGCACTGAGTACAGGGCTGCGCTGACTGCCCGCCtgaactgaac                                       <  1:3626466/150‑1 (MQ=255)
                                                                                                                ccAGCGGATTAGTGCTAAGGGTTTTGTCATCACGCTAGCATTGCAGCAGAATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAATTTTGATCCTGATAAACGCAGCACTGAGTACAGGGCTGCGCTGACTGCCCGCCTGAACTGAaca                                      <  2:3442794/150‑1 (MQ=255)
                                                                                                                  aGCGGATTAGTGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGAATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAATTTTGATCCTGATAAACGCAGCACTGAGTACAGGGCTGCGCTGACTGCCCGCCTGAACTGAACACt                                    >  1:1029606/1‑150 (MQ=255)
                                                                                                                     ggATTAGTGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGAATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAATTTTGATCCTGATAAACGCAGCACTGAGTACAGGGCTGCGCTGACTGCCCGCCTGAACTGAACACTTCt                                 >  1:1291174/1‑150 (MQ=255)
                                                                                                                     ggATTAGTGCTAAGGGTTTTGTCATCACGCTAGCATTGCAGCAGAATTCCTTCGGCTTTAATTACCGCACCTTCAGAATAATTTTGATCCTGATAAACGCAGCACTGAGAACAGGGCTGCGCTGACTGCCCGCCTGAACTGAACACTTCt                                 >  1:1850079/1‑150 (MQ=255)
                                                                                                                      gATTATTGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGAATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAATTTTGATCCTGATAAACGCAGCACTGAGTACAGGGCTGCGCTGACTGCCCGCCTGAACTGAACACTTCTg                                <  2:2018137/150‑1 (MQ=255)
                                                                                                                      gATTAGTGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGAATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAATTTTGATCCTGATAAACGCAGCACTGAGTACAGGGCTGCGCTGACTGCCCGCCTGAACTGAACACTTCTg                                <  1:4888739/150‑1 (MQ=255)
                                                                                                                      gATTAGTGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGAATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAATTTTGATCCTGATAAACGCAGCACTGAGTACAGGGCTGCGCTGACTGCCCGCCTGAACTGAACACTTCTg                                <  1:815024/150‑1 (MQ=255)
                                                                                                                        ttAGTGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGAATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAATTTTGATCCTGATAAACGCAGCACTGAGTACAGGGCTGCGCTGACTGCCCGCCTGAACTGAACACTTCTGGc                              <  1:3899780/150‑1 (MQ=255)
                                                                                                                          aGTGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGAATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAATTTTGATCCTGATAAACGCAGCACTGAGTACAGGGCTGCGCTGACTGCCCGCCTGAACTGAACACTTCTGGCgg                            <  2:5045467/150‑1 (MQ=255)
                                                                                                                            tGCTAAGGGTTTTGTCATCACGCTAGCATTGCAGCAGAATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAATTTTGATCCTGATAAACGCAGCACTGAGTACAGGGCTGCGCTGACTGCCCGCCTGAACTGAACACTTCTGGCGGTa                          >  2:1675091/1‑150 (MQ=255)
                                                                                                                              cTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGAATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAATTTTGATCCTGATAAACGCAGCACTGAGTACAGGGCTGCGCTGACTGCCCGCCTGAACTGAACACTTCTGGCGGTACg                        >  1:2667848/1‑150 (MQ=255)
                                                                                                                                 aGGGTTTTGTCATCACGCTGGCATTGCAGCAGAATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAATTTTGATCCTGATAAACGCAGCACTGAGTACAGGGCTGCGCTGACTGCCCGCCTGAACTGAACACTTCTGGCGGTACGttt                     >  1:249711/1‑150 (MQ=255)
                                                                                                                                    gTTTTGTAATCACGCTGGCATTGCATCATAATTCCTTCGGCTTTAATTACCGCCCCTTCATAATAATTTTGATCCTGATAAACGCAGCACTTAGTACAGGGCTGCGCTGACTGCCCGCCTGAACTGAACACTTCTGGCGGTACGTTTAcc                  <  1:208728/150‑1 (MQ=255)
                                                                                                                                       ttGTCATCACGCTGGCATTGCAGCAGAATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAATTTTGATCCTGATAAACGCAGCACTGAGTACAGGGCTGCGCTGACTGCCCGCCTGAACTGAACACTTCTGGCGGTACGTTTACCTcc               <  2:4650423/150‑1 (MQ=255)
                                                                                                                                       ttGTCATCACGCTGGCATTGCAGCAGAATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAATTTTGATCCTGATAAACGCAGCACTGAGTACAGGGCTGCGCTGACTGCCCGCCTGAACTGAACACTTCTGGCGGTACGTTTACCTcc               <  2:1028882/150‑1 (MQ=255)
                                                                                                                                       ttGTCATCACGCTGGCATTGCAGCAGAATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAATTTTGATCCTGATAAACGCAGCACTGAGTACAGGGCTGCGCTGACTGCCCGCCTGAACTGAACACTTCTGGCGGTACGTTTACCTcc               <  2:4329607/150‑1 (MQ=255)
                                                                                                                                       ttGTCATCACGCTAGCATTGCAGCAGAATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAATTTTGATCCTGATAAACGCAGCACTGAGTACAGGGCTGCGCTGACTGCCCGCCTGAACTGAACACTTCTGGCGGTACGTTTACCTcc               <  1:2984571/150‑1 (MQ=255)
                                                                                                                                         gTCATCACGCTGGCATTGCAGCAGAATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAATTTTGATCCTGATAAACGCAGCACTGAGTACAGGGCTGCGCTGACTGCCCGCCTGAACTGAACACTTCTGGCGGTACGTTTACCTCCAc             <  2:195466/150‑1 (MQ=255)
                                                                                                                                         gTCATCACGCTGGCATTGCAGCAGAATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAATTTTGATCCTGATAAACGCAGCACTGAGTACAGGGCTGCGCTGACTGCCCGCCTGAACTGAACACTTCTGGCGGTACGTTTACCTCCAc             >  2:3689849/1‑150 (MQ=255)
                                                                                                                                         gTCATCACGCTGGCATTGCAGCAGAATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAATTTTGATCCTGATAAACGCAGCACTGAGTACAGGGCTGCGCTGACTGCCCGCCTGAACTGAACACTTCTGGCGGTACGTTTACCTCCAc             <  2:4658085/150‑1 (MQ=255)
                                                                                                                                           catcaCGCTGGCATTGCAGCAGAATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAATTTTGATCCTGATAAACGCAGCACTGAGTACAGGGCTGCGCTGACTGCCCGCCTGAACTGAACACTTCTGGCGGTACGTTTACCTCCACGt           >  2:439495/1‑150 (MQ=255)
                                                                                                                                           catcaCGCTGGCATTGCAGCAGAATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAATTTTGATCCTGATAAACGCAGCACTGAGTACAGGGCTGCGCTGACTGCCCGCCTGAACTGAACACTTCTGGCGGTACGTTTACCTCCACGt           >  2:4444135/1‑150 (MQ=255)
                                                                                                                                            atcaCGCTGGCATTGCAGCAGAATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAATTTTGATCCTGATAAACGCAGCACTGAGTACAGGGCTGCGCTGACTGCCCGCCTGAACTGAACACTTCTGGCGGTACGTTTACCTCCACGTc          >  2:1519318/1‑150 (MQ=255)
                                                                                                                                             tcaCGCTGGCATTGCAGCAGAATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAATTTTGATCCTGATAAACGCAGCACTGAGTACAGGGCTGCGCTGACTGCCCGCCTGAACTGAACACTTCTGGCGGTACGTTTACCTCCACGTcc         >  2:1168946/1‑150 (MQ=255)
                                                                                                                                             tcaCGCTGGCATTGCAGCAGAATTCCTTCGGCTTTAATTACCGCCCCTTCAGAAGAATTTTGATCCTGATAAACGCAGCACTGAGTACAGGGCTGCGCTGACTGCCCGCCTGAACTGAACACTTCTGGCGGTACGTTTACCTCCACGTcc         >  2:523647/1‑150 (MQ=255)
                                                                                                                                                 gCTGGCATTGCAGCAGAATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAATTTTGATCCTGATAAACGCAGCACTGAGTACAGGGCTGCGCTGACTGCCCGCCTGAACTGAACACTTCTGGCGGTACGTTTACCTCCACGTCCGGAc     <  1:4759703/150‑1 (MQ=255)
                                                                                                                                                  cTGGCATTGCAGCAGAATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAATTTTGATCCTGATAAACGCAGCACTGAGTACAGGGCTGCGCTGACTGCCCGCCTGAACTGAACACTTCTGGCGGTACGTTTACCTCCACGTCCGGACg    >  1:2849509/1‑150 (MQ=255)
                                                                                                                                                   tGGCATTGCAGCAGAATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAATTTTGATCCTGATAAACGCAGCACTGAGTACAGGGCTGCGCTGACTGCCCGCCTGAACTGAACACTTCTGGCGGTACGTTTACCTCCACGTCCGGACGa   >  1:2941299/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
CATGAACATCAGGCGCTGGTCTGGTGCTCACCTGAAGAGGCGCTGCAATATCCGCTGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTA‑TGGCTTTACGCGCCGCCAGACCAGCGGATTAGTGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGTATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAATTTTGATCCTGATAAACGCAGCACTGAGTACAGGGCTGCGCTGACTGCCCGCCTGAACTGAACACTTCTGGCGGTACGTTTACCTCCACGTCCGGACGAT  >  CP009273/1836032‑1836328

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑