Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,571,658 C→A 6.8% A73A (GCG→GCT yddB ← putative TonB‑dependent outer membrane receptor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,571,6580CA 6.8% 1143.0 / 34.1 322A73A (GCG→GCTyddBputative TonB‑dependent outer membrane receptor
Reads supporting (aligned to +/- strand):  ref base C (143/157);  new base A (12/10);  total (155/167)
Fisher's exact test for biased strand distribution p-value = 6.60e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.89e-01

Only 100 of 322 total aligned reads displayed.

ACTCATTCGCTGGGTTCAGGTTATTAGTTGCACTAATACCGTCAATCAAATAGGCATTCTGGTAGGGCGACGCACCGTGAATAGAGATTTTCTCAGGGCGAATATCTCCCTGGTTCAACGAGGTACTTTGCGTTGAATCCATGCGTACCGCAGGGTTGGTTCTCAGCAGATCGCTGATATTGCCATTCCCGGTCGGTAAACGCTCGATGCTTTCACTGGTGTAATGCGTGTTGCCGTTTACGGGTACAGGGACAGGAGCATAGACCGTCATCTCTTCAAAAAAATGCATGTACAT  >  CP009273/1571510‑1571804
                                                                                                                                                    |                                                                                                                                                  
     ttCGCTGGGTTCAGGTTATTAGTTGCACTAATACCGTCAATCAAATAGGCATTCTGGTAGGGCGACGCACCGTGAATAGAGATTTTCTCAGGGCGAATATCTCCCTGGTTCAACGAGGTACTTTGCGTTGAATCCATGCGTACCGCAggg                                                                                                                                              <  2:2601323/150‑1 (MQ=255)
      tCGCTGGGTTCAGGTTATTAGTTGCACTAATACCGTCAATCAAATAGGCATTCTGGTAGGGCGACGCACCGTGAATAGAGATTTTCTCAGGGCGAATATCTCCCTGGTTCAACGAGGTACTTTGCGTTGAATCCATGCGTACCGCAGggt                                                                                                                                             >  1:69155/1‑150 (MQ=255)
      tCGCTGGGTTCAGGTTATTAGTTGCACTAATACCGTCAATCAAATAGGCATTCTGGTAGGGCGACGCACCGTGAATAGAGATTTTCTCAGGGCGAATATCTCCCTGGTTCAACGAGGTACTTTGCGTTGAATCCATGCGTACCGCAGggt                                                                                                                                             <  1:1234716/150‑1 (MQ=255)
      tCGCTGGGTTCAGGTTATTAGTTGCACTAATACCGTCAATCAAATAGGCATTCTGGTAGGGCGACGCACCGTGAATAGAGATTTTCTCAGGGCGAATATCTCCCTGGTTCAACGAGGTACTTTGCGTTGAATCCATGCGTACCGCAGggt                                                                                                                                             >  1:3592530/1‑150 (MQ=255)
      tCGCTGGGTTCAGGTTATTAGTTGCACTAATACCGTCAATCAAATAGGCATTCTGGTAGGGCGACGCACCGTGAATAGAGATTTTCTCAGGGCGAATATCTCCCTGGTTCAACGAGGTACTTTGCGTTGAATCCATGCGTACCGCAGggt                                                                                                                                             <  2:2988579/150‑1 (MQ=255)
             gatCAGGTAATTAGTTGCACTAATACCGTCAATCAACTAGGCATTCTGGTAGGGCGACGCACCGTGAATAGAGATTTTCTCAGGGCGAATATCTCCCTGGTTCAACGAGGTACTTTGCGTTGAATCCATGCGTACCGCAGGGTTGGTtct                                                                                                                                      <  2:5088008/148‑1 (MQ=255)
             gTTCAGGTTATTAGTTGCACTAATACCGTCAATCAAATAGGCATTCTGGTAGGGCGACGCACCGTGAATAGAGATTTTCTCAGGGCGAATATCTCCCTGGTTCAACGAGGTACTTTGCGTTGAATCCATGCGTACCGCAGGGTTGGTtct                                                                                                                                      <  2:705094/150‑1 (MQ=255)
               tCAGGTTATTAGTTGCACTAATACCGTCAATCAAATAGGCATTCTGGTAGGGCGACGCACCGTGAATAGAGATTTTCTCAGGGCGAATATCTCCCTGGTTCAACGAGGTACTTTGCGTTGAATCCATGCGTACCGCAGGGTTGGTTCTca                                                                                                                                    >  2:3791194/1‑150 (MQ=255)
                 aGGTTATTAGTTGCACTAATACCGTCAATCAAATAGGCATTCTGGTAGGGCGACGCACCGTGAATAGAGATTTTCTCAGGGCGAATATCTCCCTGGTTCAACGAGGTACTTTGCGTTGAATCCATGCGTACCGCAGGGTTGGTTCTcagc                                                                                                                                  <  2:4727623/150‑1 (MQ=255)
                 aGGTTATTAGTTGCACTAATACCGTCAATCAAATAGGCATTCTGGTAGGGCGACGCACCGTGAATAGAGATTTTCTCAGGGCGAATATCTCCCTGGTTCAACGAGGTACTTTGCGTTGAATCCATGCGTACCGCAGGGTTGGTTCTcagc                                                                                                                                  >  2:394783/1‑150 (MQ=255)
                  ggTTATTAGTTGCACTAATACCGTCAATCAAATAGGCATTCTGGTAGGGCGACGCACCGTGAATAGAGATTTTCTCAGGGCGAATATCTCCCTGGTTCAACGAGGTACTTTGCGTTGAATCCATGCGTACCGCAGGGTTGGTTCTcagca                                                                                                                                 <  2:2532063/150‑1 (MQ=255)
                    ttattaGTTGCACTAATACCGTCAATCNAATAGGCATTCTGGTAGGGCGACGCACCGTGAATAGAGATTTTCTCAGGGCGAATATCTCCCTGGTTCAACGAGGTACTTTGCGTTGAATCCATGCGTACCGCAGGGTTGGTTCTCAGCAGa                                                                                                                               >  1:3654708/1‑150 (MQ=255)
                     tattaGTTGCACTAATACCGTCAATCAAATAGGCATTCTGGTAGGGCGACGCACCGTGAATAGAGATTTTCTCAGGGCGAATATCTCCCTGGTTCAACGAGGTACTTTGCGTTGAATCCATGCGTACAGCAGGGTTGGTTCTCAGCAGAt                                                                                                                              >  2:772930/1‑150 (MQ=255)
                      attaGTTGCACTAATACCGTCAATCAAATAGGCATTCTGGTAGGGCGACGCACCGTGAATAGAGATTTTCTCAGGGCGAATATCTCCCTGGTTCAACGAGGTACTTTGCGTTGAATCCATGCGTACCGCAGGGTTGGTTCTCAGCAGATc                                                                                                                             <  1:1083524/150‑1 (MQ=255)
                         aGTTGCACTAATACCGTCAATCAAATAGGCATTCTGGTAGGGCGACGCACCGTGAATAGAGATTTTCTCAGGGCGAATATCTCCCTGGTTCAACGAGGTACTTTGCGTTGAATCCATGCGTACAGCAGGGTTGGTTCTCAGCAGATCGCt                                                                                                                          >  1:4778585/1‑150 (MQ=255)
                            tGCACTAATACCGTCAATCAAATAGGCATTCTGGTAGGGCGACGCACCGTGAATAGAGATTTTCTCAGGGCGAATATCTCCCTGGTTCAACGAGGTACTTTGCGTTGAATCCATGCGTACCGCAGGGTTGGTTCTCAGCAGATCGCTGat                                                                                                                       >  1:4050545/1‑150 (MQ=255)
                            tGCACTAATACCGTCAATCAAATAGGCATTCTGGTAGGGCGACGCACCGTGAATAGAGATTTTCTCAGGGCGAATATCTCCCTGGTTCAACGAGGTACTTTGCGTTGAATCCATGCGTACCGCAGGGTTGGTTCTCAGCAGATCGCTGat                                                                                                                       <  2:1793349/150‑1 (MQ=255)
                            tGCACTAATACCGTCAATCAAATAGGCATTCTGGTAGGGCGACGCACCGTGAATAGAGATTTTCTCAGGGCGAATATCTCCCTGGTTCAACGAGGTACTTTGCGTTGAATCCATGCGTACCGCAGGGTTGGTTCTCAGCAGATCGCTGat                                                                                                                       >  2:4063945/1‑150 (MQ=255)
                             gCACTAATACCGTCAATCAAATAGGCATTCTGGTAGGGCGACGCACCGTGAATAGAGATTTTCTCAGGGCGAATATCTCACTGGTTCAACGAGGTACTTTGCGTTGAATCCATGAGTACCGCAGGGTTGGTTCTCAGCAGATCGCTGata                                                                                                                      >  2:2938665/1‑150 (MQ=255)
                               aCTAATACCGTCAATCAAATAGGCATTCTGGTAGGGCGACGCACCGTGAATAGAGATTTTCTCAGGGCGAATATCTCCCTGGTTCAACGAGGTACTTTGCGTTGAATCCATGCGTACCGCAGGGTTGGTTCTCAGCAGATCGCTGATAtt                                                                                                                    <  1:139939/150‑1 (MQ=255)
                               aCTAATACCGTCAATCAAATAGGCATTCTGGTAGGGCGACGCACCGTGAATAGAGATTTTCTCAGGGCGAATATCTCCCTGGTTCAACGAGGTACTTTGCGTTGAATCCATGCGTACCGCAGGGTTGGTTCTCAGCAGATCGCTGATAtt                                                                                                                    <  2:3299956/150‑1 (MQ=255)
                                cTAATACCGTCAATCACATAGGCATTCTGGTAGGGCGACGCACCGTGAATAGAGATTTTCTCAGGGCGAATATCTCCCTGGTTCAACGAGGTACTTTGCGTTGAATCCATGCGTACAGCAGGGTTGGTTCTCAGCAGATCGCTGATATTg                                                                                                                   <  2:4513274/150‑1 (MQ=255)
                                 tAATACCGTCAATCAAATAGGCATTCTGGTAGGGCGACGCACCGTGAATAGAGATTTTCTCAGGGCGAATATCTCCCTGGTTCAACGAGGTACTTTGCGTTGAATCCATGCGTACAGCAGGGTTGGTTCTCAGCAGATCGCTGATATTGc                                                                                                                  <  1:756934/150‑1 (MQ=255)
                                  aaTACCGTCAATCAAATAGGCATTCTGGTAGGGCGACGCACCGTGAATAGAGATTTTCTCAGGGCGAATATCTCCCTGGTTCAACGAGGTACTTTGCGTTGAATCCATGCGTACCGCAGGGTTGGTTCTCAGCAGATCGCTGATATTGcc                                                                                                                 <  2:3114739/150‑1 (MQ=255)
                                    tACCGTCAATCAAATAGGCTTTCTGGTAGGGCGACGCACCGTGAATAGAGATTTTCTCAGGGCGAATATCTCCCTGGTTCAACGAGGTACTTTGCGTTGAATCCATGCGTACCGCAGGGTTGGTTCTCAGCAGATCGCTGATATTGCCAt                                                                                                               <  1:1571397/150‑1 (MQ=255)
                                     aCCGTCAATCAAATAGGCATTCTGGTAGGGCGACGCACCGTGAATAGAGATTTTCTCAGGGCGAATATCTCCCTGGTTCAACGAGGTACTTTGCGTTGAATCCATGCGTACCGCAGGGTTGGTTCTCAGCAGATCGCTGATATTGCCAtt                                                                                                              >  2:3620940/1‑150 (MQ=255)
                                     aCCGTCAATCAAATAGGCATTCTGGTAGGGCGACGCACCGTGAATAGAGATTTTCTCAGGGCGAATATCTCCCTGGTTCAACGAGGTACTTTGCGTTGAATCCATGCGTACCGCAGGGTTGGTTCTCAGCAGATCGCTGATATTGCCAtt                                                                                                              >  2:1558167/1‑150 (MQ=255)
                                     aCCGTCAATCAAATAGGCATTCTGGTAGGGCGACGCACCGTGAATAGAGATTTTCTCAGGGCGAATATCTCCCTGGTTCAACGAGGTACTTTGCGTTGAATCCATGCGTACCGCAGGGTTGGTTCTCAGCAGATCGCTGATATTGCCAtt                                                                                                              >  2:3672799/1‑150 (MQ=255)
                                       cGTCAATCAAATAGGCATTCTGGTAGGGCGACGCACCGTGAATAGAGATTTTCTCAGGGCGAATATCTCCCTGGTTCAACGAGGTACTTTGCGTTGAATCCATGCGTACCGCAGGGTTGGTTCTCAGCAGATCGCTGATATTGCCATTcc                                                                                                            >  1:3075102/1‑150 (MQ=255)
                                       cGTCAATCAAATAGGCATTCTGGTAGGGCGACGCACCGTGAATAGAGATTTTCTCAGGGCGAATATCTCCCTGGTTCAACGAGGTACTTTGCGTTGAATCCATGCGTACCGCAGGGTTGGTTCTCAGCAGATCGCTGATATTGCCATTcc                                                                                                            >  1:3481797/1‑150 (MQ=255)
                                        gTCAATCAAATAGGCATTCTGGTAGGGCGACGCACCGTGAATAGAGATTTTCTCAGGGCGAATATCTCCCTGGTTCAACGAGGTACTTTGCGTTGAATCCATGCGTACCGCAGGGTTGGTTCTCAGCAGATCGCTGATATTGCCATTccc                                                                                                           <  1:2882987/150‑1 (MQ=255)
                                          caatcaaATAGGCATTCTGGTAGGGCGACGCACCGTGAATAGAGATTTTCTCAGGGCGAATATCTCCCTGGTTCAACGAGGTACTTTGCGTTGAATCCATGCGTACCGCAGGGTTGGTTCTCAGCAGATCGCTGATATTGCCATTCCcgg                                                                                                         >  1:2258311/1‑150 (MQ=255)
                                            atcaaATAGGCATTCTGGTAGGGCGACGCACCGTGAATAGAGATTTTCTCAGGGCGAATATCTCCCTGGTTCAACGAGGTACTTTGCGTTGAATCCATGCGTACCGCAGGGTTGGTTCTCAGCAGATCGCTGATATTGCCATTCCcggtc                                                                                                       >  1:4409026/1‑150 (MQ=255)
                                             tcaaATAGGCATTCTGGTAGGGCGACGCACCGTGAATAGAGATTTTCTCAGGGCGAATATCTCCCTGGTTCAACGAGGTACTTTGCGTTGAATCCATGCGTACCGCAGGGTTGGTTCTCAGCAGATCGCTGATATTGCCATTCCcggtcg                                                                                                      >  1:2731109/1‑150 (MQ=255)
                                              caaATAGGCATTCTGGTAGGGCGACGCACCGTGAATAGAGATTTTCTCAGGGCGAATATCTCCCTGGTTCAACGAGGTACTTTGCGTTGAATCCATGCGTACCGCAGGGTTGGTTCTCAGCAGATCGCTGATATTGcc                                                                                                                 >  1:4632295/1‑138 (MQ=255)
                                               aaaTAGGCATTCTGGTAGGGCGACGCACCGTGAATAGAGATTTTCTCAGGGCGAATATCTCCCTGGTTCAACGAGGTACTTTGCGTTGAATCCATGCGTACCGCAGGGTTGGTTCTCAGCAGATCGCTGATATTGCCATTCCcggtcggt                                                                                                    <  2:2336233/150‑1 (MQ=255)
                                               aaaTAGGCATTCTGGTAGGGCGACGCACCGTGAATAGAGATTTTCTCAGGGCGAATATCTCCCTGGTTCAACGAGGTACTTTGCGTTGAATCCATGCGTACCGCAGGGTTGGTTCTCAGCAGATCGCTGATATTGCCATTCCcggtcggt                                                                                                    <  2:2836136/150‑1 (MQ=255)
                                                aaTAGGCATTCTGGTAGGGCGACGCACCGTGAATAGAGATTTTCTCAGGGCGAATATCTCCCTGGTTCAACGAGGTACTTTGCGTTGAATCCATGCGTACCGCAGGGTTGGTTCTCAGCAGATCGCTGATATTGCCATTCCCGGTCGGTa                                                                                                   <  1:1696114/150‑1 (MQ=255)
                                                   aGGCATTCTGGTAGGGCGACGCACCGTGAATAGAGATTTTCTCAGGGCGAATATCTCCCTGGTTCAACGAGGTACTTTGCGTTGAATCCATGCGTACCGCAGGGTTGGTTCTCAGCAGATCGCTGATATTGCCATTCCCGGTCGGTAAAc                                                                                                >  2:1199116/1‑150 (MQ=255)
                                                    ggCATTCTGGTAGGGCGACGCACCGTGAATAGAGATTTTCTCAGGGCGAATATCTCCCTGGTTCAACGAGGTACTTTGCGTTGAATCCATGCGTACCGCAGGGTTGGTTCTCAGCAGATCGCTGATATTGCCATTCCCGGTCGGTAAACg                                                                                               >  2:1517028/1‑150 (MQ=255)
                                                          cTGGTAGGGCGACGCACCGTGAATAGAGATTTTCTCAGGGCGAATATCTCCCTGGTTCAACGAGGTACTTTGCGTTGAATCCATGCGTACCGCAGGGTTGGTTCTCAGCAGATCGCTGATATTGCCATTCCCGGTCGGTAAACGCTCGAt                                                                                         <  2:2051495/150‑1 (MQ=255)
                                                          cTGGTAGGGCGACGCACCGTGAATAGAGATTTTCTCAGGGCGAATATCTCCCTGGTTCAACGAGGTACTTTGCGTTGAATCCATGCGTACCGCAGGGTTGGTTCTCAGCAGATCGCTGATATTGCCATTCCCGGTCGGTAAACGCTCGAt                                                                                         <  2:2884742/150‑1 (MQ=255)
                                                           tGGTAGGGCGACGCACCGTGAATAGAGATTTTCTCAGGGCGAATATCTCCCTGGTTCAACGAGGTACTTTGCGTTGAATCCATGCGTACCGCAGGGTTGGTTCTCAGCAGATCGCTGATATTGCCATTCCCGGTCGGTAAACGCTCGATg                                                                                        <  1:1759349/150‑1 (MQ=255)
                                                            ggTAGGGCGACGCACCGTGAATAGAGATTTTCTCAGGGCGAATATCTCCCTGGTTCAACGAGGTACTTTGCGTTGAATCCATGCGTACCGCAGGGTTGGTTCTCAGCAGATCGCTGATATTGCCATTCCCGGTCGGTAAACGCTCGATGc                                                                                       >  1:1629073/1‑150 (MQ=255)
                                                            ggTAGGGCGACGCACCGTGAATAGAGATTTTCTCAGGGCGAATATCTCCCTGGTTCAACGAGGTACTTTGCGTTGAATCCATGCGTACCGCAGGGTTGGTTCTCAGCAGATCGCTGATATTGCCATTCCCGGTCGGTAAACGCTCGATGc                                                                                       <  2:44807/150‑1 (MQ=255)
                                                             gTAGGGCGACGCACCGTGAATAGAGATTTTCTCAGGGCGAATATCTCCCTGGTTCAACGAGGTACTTTGCGTTGAATCCATGCGTACCGCAGGGTTGGTTCTCAGCAGATCGCTGATATTGCCATTCCCGGTCGGTAAACGCTCGATGCt                                                                                      <  2:1073401/150‑1 (MQ=255)
                                                              tAGGGCGACGCACCGTGAATAGAGATTTTCTCAGGGCGAATATCTCCCTGGTTCAACGAGGTACTTTGCGTTGAATCCATGCGTACCGCAGGGTTGGTTCTCAGCAGATCGCTGATATTGCCATTCCCGGTCGGTAAACGCTCGATGCtt                                                                                     >  2:1740590/1‑150 (MQ=255)
                                                               aGGGCGACGCACCGTGAATAGAGATTTTCTCAGGGCGAATATCTCCCTGGTTCAACGAGGTACTTTGCGTTGAATCCATGCGTACCGCAGGGTTGGTTCTCAGCAGATCGCTGATATTGCCATTCCCGGTCGGTAAACGCTCGATGCttt                                                                                    >  1:2772005/1‑150 (MQ=255)
                                                                gggCGACGCACCGTGAATAGAGATTTTCTCAGGGCGAATATCTCCCTGGTTCAACGAGGTACTTTGCGTTGAATCCATGCGTACCGCAGGGTTGGTTCTCAGCAGATCGCTGATATTGCCATTCCCGGTCGGTAAACGCTCGATGCTTTc                                                                                   <  2:92267/150‑1 (MQ=255)
                                                                gggCGACGCACCGTGAATAGAGATTTTCTCAGGGCGAATATCTCCCTGGTTCAACGAGGTACTTTGCGTTGAATCCATGCGTACCGCAGGGTTGGTTCTCAGCAGATCGCTGATATTGCCATTCCCGGTCGGTAAACGCTCGATGCTTTc                                                                                   <  2:1859717/150‑1 (MQ=255)
                                                                 ggCGACGCACCGTGAATAGAGATTTTCTCAGGGCGAATATCTCCCTGGTTCAACGAGGTACTTTGCGTTGAATCCATGCGTACCGCAGGGTTGGTTCTCAGCAGATCGCTGATATTGCCATTCCCGGTCGGTAAACGCTCGATGCTTTCa                                                                                  >  2:1165673/1‑150 (MQ=255)
                                                                    gACGCACCGTGAATAGAGATTTTCTCAGGGCGAATATCTCCCTGGTTCAACGAGGTACTTTGCGTTGAATCCATGCGTACCGCAGGGTTGGTTCTCAGCAGATCGCTGATATTGCCATTCCCGGTCGGTAAACGCTCGATGCTTTCACTg                                                                               >  1:1009800/1‑150 (MQ=255)
                                                                    gACGCACCGTGAATAGAGATTTTCTCAGGGCGAATATCTCCCTGGTTCAACGAGGTACTTTGCGTTGAATCCATGCGTACCGCAGGGTTGGTTCTCAGCAGATCGCTGATATTGCCATTCCCGGTCGGTAAACGCTCGATGCTTTCACTg                                                                               >  2:872585/1‑150 (MQ=255)
                                                                       gCACCGTGAATAGAGATTTTCTCAGGGCGAATATCTCCCTGGTTCAACGAGGTACTTTGCGTTGAATCCATGCGTACCGCAGGGTTGGTTCTCAGCAGATCGCTGATATTGCCATTCCCGGTCGGTAAACGCTCGATGCTTTCACTGgtg                                                                            >  2:3240509/1‑150 (MQ=255)
                                                                       gCACCGTGAATAGAGATTTTCTCAGGGCGAATATCTCCCTGGTTCAACGAGGTACTTTGCGTTGAATCCATGCGTACCGCAGGGTTGGTTCTCAGCAGATCGCTGATATTGCCATTCCCGGTCGGTAAACGCTCGATGCTTTCACTGgtg                                                                            <  2:452077/150‑1 (MQ=255)
                                                                         aCCGTGAATAGAGATTTTCTCAGGGCGAATATCTCCCTGGTTCAACGAGGTACTTTGCGTTGAATCCATGCGTACCGCAGGGTTGGTTCTCAGCAGATCGCTGATATTGCCATTCCCGGTCGGTAAACGCTCGATGCTTTCACTGGTGTa                                                                          >  2:2184705/1‑150 (MQ=255)
                                                                          ccGTGAATAGAGATTTTCTCAGGGCGAATATCTCCCTGGTTCAACGAGGTACTTTGCGTTGAATCCATGCGTACCGCAGGGTTGGTTCTCAGCAGATCGCTGATATTGCCATTCCCGGTCGGTAAACGCTCGATGCTTTCACTGGTGTaa                                                                         <  2:2719024/150‑1 (MQ=255)
                                                                            gTGAATAGAGATTTTCTCAGGGCGAATATCTCCCTGGTTCAACGAGGTACTTTGCGTTGAATCCATGCGTACCGCAGGGTTGGTTCTCAGCAGATCGCTGATATTGCCATTCCCGGTCGGTAAACGCTCGATGCTTTCACTGGTGTAATg                                                                       >  1:5102961/1‑150 (MQ=255)
                                                                             tGAATAGAGATTTTCTCAGGGCGAATATCTCCCTGGTTCAACGAGGTACTTTGCGTTGAATCCATGCGTACCGCAGGGTTGGTTCTCAGCAGATCGCTGATATTGCCATTCCCGGTCGGTAAACGCTCGATGCTTTCACTGGTGTAATGc                                                                      >  2:2475134/1‑150 (MQ=255)
                                                                              gAATAGAGATTTTCTCAGGGCGAATATCTCCCTGGTTCAACGAGGTACTTTGCGTTGAATCCATGCGTACCGCAGGGTTGGTTCTCAGCAGATCGCTGATATTGCCATTCCCGGTCGGTAAACGCTCGATGCTTTCACTGGTGTAATGCg                                                                     <  2:4705626/150‑1 (MQ=255)
                                                                               aaTAGAGATTTTCTCAGGGCGAATATCTCCCTGGTTCAACGAGGTACTTTGCGTTGAATCCATGCGTACCGCAGGGTTGGTTCTCAGCAGATCGCTGATATTGCCATTCCCGGTCGGTAAACGCTCGATGCTTTCACTGGTGTAATGCgt                                                                    <  2:3481797/150‑1 (MQ=255)
                                                                                aTAGAGATTTTCTCAGGGCGAATATCTCCCTGGTTCAACGAGGTACTTTGCGTTGAATCCATGCGTACAGCAGGGTTGGTTCTCAGCAGATCGCTGATATTGCCATTCCCGGTCGGTAAACGCTCGATGCTTTCACTGGTGTAATGCgtg                                                                   >  1:2153899/1‑150 (MQ=255)
                                                                                aTAGAGATTTTCTCAGGGCGAATATCTCCCTGGTTCAACGAGGTACTTTGCGTTGAATCCATGCGTACAGCAGGGTTGGTTCTCAGCAGATCGCTGATATTGCCATTCCCGGTCGGTAAACGCTCGATGCTTTCACTGGTGTAATGCgtg                                                                   <  2:4778585/150‑1 (MQ=255)
                                                                                     gaTTTTCTCAGGGCGAATATCTCCCTGGTTCAACGAGGTACTTTGTGTTGAATCCATGCGTACCGCAGGGTTGGTTCTCAGCAGATCGCTGATATTGCCATTCCCGGTCGGTAAACGCTCGATGCTTTCACTGGTGTAATGCGTGTTGcc                                                              <  1:2451274/150‑1 (MQ=255)
                                                                                        tttCTCAGGGCGAATATCTCCCTGGTTCAACGAGGTACTTTGCGTTGAATCCATGCGTACCGCAGGGTTGGTTCTCAGCAGATCGCTGATATTGCCATTCCCGGTCGGTAAACGCTCGATGCTTTCACTGGTGTAATGCGTGTTGCCGtt                                                           >  2:4340459/1‑150 (MQ=255)
                                                                                         ttCTCAGGGCGAATATCTCCCTGGTTCAACGAGGTACTTTGCGTTGAATCCATGCGTACCGCAGGGTTGGTTCTCAGCAGATCGCTGATATTGCCATTCCCGGTCGGTAAACGCTCGATGCTTTCACTGGTGTAATGCGTGTTGCCGttt                                                          <  1:3325733/150‑1 (MQ=255)
                                                                                          tctcAGGGCGAATATCTCCCTGGTTCAACGAGGTACTTTGCGTTGAATCCATGCGTACCGCAGGGTTGGTTCTCAGCAGATCGCTGATATTGCCATTCCCGGTAGGTAAACGCTCGATGCTTTCACTGGTGTAATGCGTGTTGCCGTTTa                                                         >  1:3234271/1‑150 (MQ=255)
                                                                                           ctcAGGGCGAATATCTCCCTGGTTCAACGAGGTACTTTGCGTTGAATCCATGCGTACCGCAGGGTTGGTTCTCAGCAGATCGCTGATATTGCCATTCCCGGTCGGTAAACGCTCGATGCTTTCACTGGTGTAATGCGTGTTGCCGTTTAc                                                        >  2:4667909/1‑150 (MQ=255)
                                                                                             cAGGGCGAATATCTCCCTGGTTCAACGAGGTACTTTGCGTTGAATCCATGCGTACCGCAGGGTTGGTTCTCAGCAGATCGCTGATATTGCCATTCCCGGTCGGTAAACGCTCGATGCTTTCACTGGTGTAATGCGTGTTGCCGTTTACgg                                                      >  2:4429868/1‑150 (MQ=255)
                                                                                              aGGGCGAATATCTCCCTGGTTCAACGAGGTACTTTGCGTTGAATCCATGCGTACCGCAGGGTTGGTTCTCAGCAGATCGCTGATATTGCCATTCCCGGTCGGTAAACGCTCGATGCTTTCACTGGTGTAATGCGTGTTGCCGTTTACggg                                                     <  1:991900/150‑1 (MQ=255)
                                                                                              aGGGCGAATATCTCCCTGGTTCAACGAGGTACTTTGCGTTGAATCCATGCGTACCGCAGGGTTGGTTCTCAGCAGATCGCTGATATTGCCATTCCCGGTCGGTAAACGCTCGATGCTTTCACTGGTGTAATGCGTGTTGCCGTTTACggg                                                     <  1:1855607/150‑1 (MQ=255)
                                                                                                ggCGAATATCTCCCTGGTTCAACGAGGTACTTTGCGTTGAATCCATGCGTACCGCAGGGTTGGTTCTCAGCAGATCGCTGATATTGCCATTCCCGGTCGGTAAACGCTCGATGCTTTCACTGGTGTAATGCGTGTTGCCGTTTACGGGTa                                                   >  2:4197499/1‑150 (MQ=255)
                                                                                                  cGAATATCTCCCTGGTTCAACGAGGTACTTTGCGTTGAATCCATGCGTACCGCAGGGTTGGTTCTCAGCAGATCGCTGATATTGCCATTCCCGGTCGGTAAACGCTCGATGCTTTCACTGGTGTAATGCGTGTTGCCGTTTACGGGTACa                                                 <  1:4571824/150‑1 (MQ=255)
                                                                                                    aaTATCTCCCTGGTTCAACGAGGTACTTTGCGTTGAATCCATGCGTACCGCAGGGTTGGTTCTCAGCAGATCGCTGATATTGCCATTCCCGGTCGGTAAACGCTCGATGCTTTCACTGGTGTAATGCGTGTTGCCGTTTACGGGTACAgg                                               <  1:2997848/150‑1 (MQ=255)
                                                                                                       atCTCCCTGGTTCAACGAGGTACTTTGCGTTGAATCCATGCGTACCGCAGGGTTGGTTCTCAGCAGATCGCTGATATTGCCATTCCCGGTCGGTAAACGCTCGATGCTTTCACTGGTGTAATGCGTGTTGCCGTTTACGGGTACAGGGAc                                            >  1:2101925/1‑150 (MQ=255)
                                                                                                         ctcCCTGGTTCAACGAGGTACTTTGCGTTGAATCCATGCGTACCGCAGGGTTGGTTCTCAGCAGATCGCTGATATTGCCATTCCCTGTCGGTAAACGCTCGATGCTTTCACTGGTGTAATGCGTGTTGCCGTTTACGGGTACAGGGACAg                                          <  2:1959250/150‑1 (MQ=255)
                                                                                                         ctcCCTGGTTCAACGAGGTACTTTGCGTTGAATCCATGCGTACCGCAGGGTTGGTTCTCAGCAGATCGCTGATATTGCCATTCCCGGTCGGTAAACGCTCGATGCTTTCACTGGTGTAATGCGTGTTGCCGTTTACGGGTACAGGGACAg                                          >  1:2602986/1‑150 (MQ=255)
                                                                                                          tcCCTGGTTCAACGAGGTACTTTGCGTTGAATCCATGCGTACCGCAGGGTTGGTTCTCAGCAGATCGCTGATATTGCCATTCCCGGTCGGTAAACGCTCGATGCTTTCACTGGTGTAATGCGTGTTGCCGTTTACGGGTACAGGGACAgg                                         >  1:1968776/1‑150 (MQ=255)
                                                                                                          tcCCTGGTTCAACGAGGTACTTTGCGTTGAATCCATGCGTACCGCAGGGTTGGTTCTCAGCAGATCGCTGATATTGCCATTCCCGGTCGGTAAACGCTCGATGCTTTCACTGGTGTAATGCGTGTTGCCGTTTACGGGTACAGGGACAgg                                         <  1:3514206/150‑1 (MQ=255)
                                                                                                            ccTGGTTCAACGAGGTACTTTGCGTTGAATCCATGCGTACCGCAGGGTTGGTTCTCAGCAGATCGCTGATATTGCCATTCCCGGTCGGTAAACGCTCGATGCTTTCACTGGTGTAATGCGTGTTGCCATTTACGGGTACAGGGACAGGAg                                       <  1:647952/150‑1 (MQ=255)
                                                                                                                gTTCAACGAGGTACTTTGCGTTGAATCCATGCGTACCGCAGGGTTGGTTCTCAGCAGATCGCTGATATTGCCATTCCCGGTCGGTAAACGCTCGATGCTTTCACTGGTGTAATGCGTGTTGCCGTTTACGGGTACAGGGACAGGAGCATa                                   <  1:2981695/150‑1 (MQ=255)
                                                                                                                gTTCAACGAGGTACTTTGCGTTGAATCCATGCGTACCGCAGGGTTGGTTCTCAGCAGATCGCTGATATTGCCATTCCCGGTCGGTAAACGCTCGATGCTTTCACTGGTGTAATGCGTGTTGCCGTTTACGGGTACAGGGACAGGAGCATa                                   >  1:4902063/1‑150 (MQ=255)
                                                                                                                gTTCAACGAGGTACTTTGCGTTGAATCCATGCGTACCGCAGGGTTGGTTCTCAGCAGATCGCTGATATTGCCATTCCCGGTCGGTAAACGCTCGATGCTTTCACTGGTGTAATGCGTGTTGCCGTTTACGGGTACAGGGACAGGAGCATa                                   >  1:5082083/1‑150 (MQ=255)
                                                                                                                 ttCAACGAGGTACTTTGCGTTGAATCCATGCGTACCGCAGGGTTGGTTCTCAGCAGATCGCTGATATTGCCATTCCCGGTCGGTAAACGCTCGATGCTTTCACTGGTGTAATGCGTGTTGCCGTTTACGGGTACAGGGACAGGAGCATAg                                  >  2:1749691/1‑150 (MQ=255)
                                                                                                                  tCAACGAGGTACTTTGCGTTGAATCCATGCGTACCGCAGGGTTGGTTCTCAGCAGATCGCTGATATTGCCATTCCCGGTCGGTAAACGCTCGATGCTTTCACTGGTGTAATGCGTGTTGCCGTTTACGGGTACAGGGACAGGAGCATAGa                                 >  1:2682510/1‑150 (MQ=255)
                                                                                                                    aaCGAGGTACTTTGCGTTGAATCCATGCGTACCGCAGGGTTGGTTCTCAGCAGATCGCTGATATTGCCATTCCCGGTCGGTAAACGCTCGATGCTTTCACTGGTGTAATGCGTGTTGCCGTTTACGGGTACAGGGACAGGAGCATAGAcc                               <  2:4092751/150‑1 (MQ=255)
                                                                                                                      cGAGGTACTTTGCGTTGAATCCATGCGTACCGCAGGGTTGGTTCTCAGCAGATCGCTGATATTGCCATTCCCGGTCGGTAAACGCTCGATGCTTTCACTGGTGTAATGCGTGTTGCCGTTTACGGGTACAGGGACAGGAGCATAGACCGt                             >  1:2430768/1‑150 (MQ=255)
                                                                                                                       gAGGTACTTTGCGTTGAATCCATGCGTACCGCAGGGTTGGTTCTCAGCAGATCGCTGATATTGCCATTCCCGGTCGGTAAACGCTCGATGCTTTCACTGGTGTAATGCGTGTTGCCGTTTACGGGTACAGGGACAGGAGCATAGACCGTc                            <  1:447505/150‑1 (MQ=255)
                                                                                                                        aGGTACTTTGCGTTGAATCCATGCGTACCGCAGGGTTGGTTCTCAGCAGATCGCTGATATTGCCATTCCCGGTCGGTAAACGCTCGATGCTTTCACTGGTGTAATGCGTGTTGCCGTTTACGGGTACAGGGACAGGAGCATAGACCGTCa                           <  1:2580676/150‑1 (MQ=255)
                                                                                                                            aCTTTGCGTTGAATCCATGCGTACCGCAGGGTTGGTTCTCAGCAGATCGCTGATATTGCCATTCCCGGTCGGTAAACGCTCGATGCTTTCACTGGTGTAATGCGTGTTGCCGTTTACGGGTACAGGGACAGGAGCATAGACCGTCAtctc                       <  1:1191413/150‑1 (MQ=255)
                                                                                                                             cTTTGCGTTGAATCCATGCGTACCGCAGGGTTGGTTCTCAGCAGATCGCTGATATTGCCATTCCCGGTCGGTAAACGCTCGATGCTTTCACTGGTGTAATGCGTGTTGCCGTTTACGGGTACAGGGACAGGAGCATAGACCGTCATctct                      <  1:423786/150‑1 (MQ=255)
                                                                                                                                tGCGTTGAATCCATGCGTACCGCAGGGTTGGTTCTCAGCAGATCGCTGATATTGCCATTCCCGGTCGGTAAACGCTCGATGCTTTCACTGGTGTAATGCGTGTTGCCGTTTACGGGTACAGGGACAGGAGCATAGACCGTCATCTCTTCa                   <  2:334800/150‑1 (MQ=255)
                                                                                                                                    ttGAATCCATGCGTACCGCAGGGTTGGTTCTCAGCAGATCGCTGATATTGCCATTCCCGGTCGGTAAACGCTCGATGCTTTCACTGGTGTAATGCGTGTTGCCGTTTACGGGTACAGGGACAGGAGCATAGACCGTCATCTCTTCaaaaa               >  2:138974/1‑150 (MQ=255)
                                                                                                                                      gAATCCATGCGTACCGCAGGGTTGGTTCTCAGCAGATCGCTGATATTGCCATTCCCGGTCGGTAAACGCTCGATGCTTTCACTGGTGTAATGCGTGTTGCCGTTTACGGGTACAGGGACAGGAGCATAGACCGTCATCTCTTCaaaaaaa             >  1:5038210/1‑150 (MQ=255)
                                                                                                                                        aTCCATGCGTACCGCAGGGTTGGTTCTCAGCAGATCGCTGATATTGCCATTCCCGGTCGGTAAACGCTCGATGCTTTCACTGGTGTAATGCGTGTTGCCGTTTACGGGTACAGGGACAGGAGCATAGACCGTCATCTCTTCAAAAAAATg           >  2:1989399/1‑150 (MQ=255)
                                                                                                                                        aTCCATGCGTACCGCAGGGTTGGTTCTCAGCAGATCGCTGATATTGCCATTCCCGGTCGGTAAACGCTCGATGCTTTCACTGGTGTAATGCGTGTTGCCGTTTACGGGTACAGGGACAGGAGCATAGACCGTCATCTCTTCAAAAAAATg           >  1:3808295/1‑150 (MQ=255)
                                                                                                                                          ccATGCGTACCGCAGGGTTGGTTCTCAGCAGATCGCTGATATTGCCATTCCCGGTCGGTAAACGCTCGATGCTTTCACTGGTGTAATGCGTGTTGCCGTTTACGGGTACAGGGACAGGAGCATAGACCGTCATCTCTTCAAAAAAATGCa         <  2:1502685/150‑1 (MQ=255)
                                                                                                                                          ccATGCGTACCGCAGGGGTGGTTCTCAGCAGATCGCTGATATTGCCATTCCCGGTCGGTAAACGCTCGATGCTTTCACTGGTGTAATGCGTGTTGCCGTTTACGGGTACAGGGACAGGAGCATAGACCGTCATCTCTTCGAAAAAATGCa         <  2:2623874/150‑1 (MQ=255)
                                                                                                                                              gCGTACCGCAGGGTTGGTTCTCAGCAGATCGCTGATATTGCCATTCCCGGTCGGTAAACGCTCGATGCTTTCACTGGTGTAATGCGTGTTGCCGTTTACGGGTACAGGGACAGGAGCATAGACCGTCATCTCTTCAAAAAAATGCATGTa     <  2:3654737/150‑1 (MQ=255)
                                                                                                                                                 tACCGCAGGGTTGGTTCTCAGCAGATCGCTGATATTGCCATTCCCGGTCGGTAAACGCTCGATGCTTTCACTGGTGTACTGCGTGTTGCCGTTTACGGGTACAGGGACAGGAGCATAGACCGTCATCTCTTCAAAAAAATGCATGTACAt  >  1:2020982/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
ACTCATTCGCTGGGTTCAGGTTATTAGTTGCACTAATACCGTCAATCAAATAGGCATTCTGGTAGGGCGACGCACCGTGAATAGAGATTTTCTCAGGGCGAATATCTCCCTGGTTCAACGAGGTACTTTGCGTTGAATCCATGCGTACCGCAGGGTTGGTTCTCAGCAGATCGCTGATATTGCCATTCCCGGTCGGTAAACGCTCGATGCTTTCACTGGTGTAATGCGTGTTGCCGTTTACGGGTACAGGGACAGGAGCATAGACCGTCATCTCTTCAAAAAAATGCATGTACAT  >  CP009273/1571510‑1571804

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: