Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,902,115 C→A 6.1% intergenic (‑267/+403) yobF ← / ← yebO DUF2527 family heat‑induced protein/putative inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,902,1150CA 6.1% 1165.0 / 32.3 311intergenic (‑267/+403)yobF/yebODUF2527 family heat‑induced protein/putative inner membrane protein
Reads supporting (aligned to +/- strand):  ref base C (125/167);  new base A (6/13);  total (131/180)
Fisher's exact test for biased strand distribution p-value = 4.73e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.91e-01

Only 100 of 331 total aligned reads displayed.

TTTTGCACACTTATGCAACGGAATACGCGCCGATATATCATTGTGCTTAACCTTGCCAGTTCAGGCAGATACTTAAACTGGCGTATTTTCTAACATAGTTCAATAAGTTCAGTTGTTCCAGGCCAGTGGGAGAAGTTATTACATAGTGCGTGCAATATCACATTTTTTGCTATGCAATGAATAAAAAGTTATATCACTTTTTCTCATAAAACAGTCAGTTAACGGCTATTAATTACCCTAAAGAGAAGTCAATCCCCAAAGGGATTGTAAATTTAAAATAAGAAAAATTGATGAATGA  >  CP009273/1901967‑1902264
                                                                                                                                                    |                                                                                                                                                     
  ttGCACACTTATGCAACGGAATACGCGCCGATATATCATTGTGCTTAACCTTGCCAGTTCAGGCAGATACTTAAACTGGCGTATTTTCTAACATAGTTCAATAAGTTCAGTTGTTCCAGGCCAGTGGGAGAAGTTATTACATAgtgcgtg                                                                                                                                                    >  1:3459691/1‑150 (MQ=255)
  ttGCACACTTATGCAACGGAATACGCGCCGATATATCATTGTGCTTAACCTTGCCAGTTCAGGCAGATACTTAAACTGGCGTATTTTCTAACATAGTTCAATAAGTTCAGTTGTTCCAGGCCAGTGGGAGAAGTTATTACATAgtgcgtg                                                                                                                                                    >  1:3460209/1‑150 (MQ=255)
  ttGCACACTTATGCAACGGAATACGCGCCGATATATCATTGTGCTTAACCTTGCCAGTTCAGGCAGATACTTAAACTGGCGTATTTTCTAACATAGTTCAATAAGTTCAGTTGTTCCAGGCCAGTGGGAGAAGTTATTACATAgtgcgtg                                                                                                                                                    >  1:3460457/1‑150 (MQ=255)
   tGCACACTTATGCAACGGAATACGCGCCGATATATCATTGTGCTTAACCTTGCCAGTTCAGGCAGATACTTAAACTGGCGTATTTTCTAACATAGTTCAATAAGTTCAGTTGTTCCAGGCCAGTGGGAGAAGTTATTACATAgtgcgtgc                                                                                                                                                   >  1:2340044/1‑150 (MQ=255)
      acacTTATGCAACGGAATACGCGCCGATATATCATTGTGCTTAACCTTGCCAGTTCAGGCAGATACTTAAACTGGCGTATTTTCTAACATAGTTCAATAAGTTCAGTTGTTCCAGGCCAGTGGGAGAAGTTATTACATAGTGCGTGCAat                                                                                                                                                <  1:671909/150‑1 (MQ=255)
      acacTTATGCAACGGAATACGCGCCGATATATCATTGTGCTTAACCTTGCCAGTTCAGGCAGATACTTAAACTGGCGTATTTTCTAACATAGTTCAATAAGTTCAGTTGTTCCAGGCCAGTGGGAGAAGTTATTACATAGTGCGTGCAat                                                                                                                                                <  2:402502/150‑1 (MQ=255)
        acTTATGCAACGGAATACGCGCCGATATATCATTGTGCTTAACCTTGCCAGTTCAGGCAGATACTTAAACTGGCGTATTTTCTAACATAGTTCAATAAGTTCAGTTGTTCCAGGCCAGTGGGAGAAGTTATTACATAGTGCGTGCAatat                                                                                                                                              >  1:2531652/1‑150 (MQ=255)
        acTTATGCAACGGAATACGCGCCGATATATCATTGTGCTTAACCTTGCCAGTTCAGGCAGATACTTAAACTGGCGTATTTTCTAACATAGTTCAATAAGTTCAGTTGTTCCAGGCCAGTGGGAGAAGTTATTACATAGTGCGTGCAatat                                                                                                                                              >  1:2170840/1‑150 (MQ=255)
        acTTATGCAACGGAATACGCGCCGATATATCATTGTGCTTAACCTTGCCAGTTCAGGCAGATACTTAAACTGGCGTATTTTCTAACATAGTTCAATAAGTTCAGTTGTTCCAGGCCAGTGGGAGAAGTTATTACATAGTGCGTGCAatat                                                                                                                                              >  1:1191527/1‑150 (MQ=255)
           tATGCAACGGAATACGCGCCGATATATCATTGTGCTTAACCTTGCCAGTTCAGGCAGATACTTAAACTGGCGTATTTTCTAACATAGTTCAATAAGTTCAGTTGTTCCAGGCCAGTGGGAGAAGTTATTACATAGTGCGTGCAATATcac                                                                                                                                           <  1:639249/150‑1 (MQ=255)
            aTGCAACGGAATACGCGCCGATATATCATTGTGCTTAACCTTGCCAGTTCAGGCAGATACTTAAACTGGCGTATTTTCTAACATAGTTCAATAAGTTCAGTTGTTCCAGGCCAGTGGGAGAAGTTATTACATAGTGCGTGCAATATcaca                                                                                                                                          >  1:2545357/1‑150 (MQ=255)
             tGCAACGGAATACGCGCCGATATATCATTGTGCTTAACCTTGCCAGTTCAGGCAGATACTTAAACTGGCGTATTTTCTAACATAGTTCAATAAGTTCAGTTGTTCCAGGCCAGTGGGAGAAGTTATTACATAGTGCGTGCAATATCACAt                                                                                                                                         >  2:4277553/1‑150 (MQ=255)
             tGCAACGGAATACGCGCCGATATATCATTGTGCTTAACCTTGCCAGTTCAGGCAGATACTTAAACTGGCGTATTTTCTAACATAGTTCAATAAGTTCAGTTGTTCCAGGCCAGTGGGAGAAGTTATTACATAGTGCGTGCAATATCACAt                                                                                                                                         <  2:317163/150‑1 (MQ=255)
             tGCAACGGAATACGCGCCGATATATCATTGTGCTTAACCTTGCCAGTTCAGGCAGATACTTAAACTGGCGTATTTTCTAACATAGTTCAATAAGTTCAGTTGTTCCAGGCCAGTGGGAGAAGTTATTACATAGTGCGTGCAATATCACAt                                                                                                                                         <  2:2614726/150‑1 (MQ=255)
               cAACGGAATACGCGCCGATATATCATTGTGCTTAACCTTGCCAGTTCAGGCAGATACTTAAACTGGCGTATTTTCTAACATAGTTCAATAAGTTCAGTTGTTCCAGGCCAGTGGGAGAAGTTATTACATAGTGCGTGCAATATCACAttt                                                                                                                                       <  1:2400436/150‑1 (MQ=255)
               cAACGGAATACGCGCCGATATATCATTGTGCTTAACCTTGCCAGTTCAGGCAGATACTTAAACTGGCGTATTTTCTAACATAGTTCAATAAGTTCAGTTGTTCCAGGCCAGTGGGAGAAGTTATTACATAGTGCGTGCAATATCACAttt                                                                                                                                       <  2:2057902/150‑1 (MQ=255)
               cAACGGAATACGCGCCGATATATCATTGTGCTTAACCTTGCCAGTTCAGGCAGATACTTAAACTGGCGTATTTTCTAACATAGTTCAATAAGTTCAGTTGTTCCAGGCCAGTGGGAGAAGTTATTACATAGTGAGTGCAATATCACAttt                                                                                                                                       >  2:1114516/1‑150 (MQ=255)
                aaCGGAATACGCGCCGATATATCATTGTGCTTAACCTTGCCAGTTCAGGCAGATACTTAAACTGGCGTATTTTCTAACATAGTTCAATAAGTTCAGTTGTTCCAGGCCAGTGGGAGAAGTTATTACATAGTGCGTGCAATATCACAtttt                                                                                                                                      <  1:4020688/150‑1 (MQ=255)
                aaCGGAATACGCGCCGATATATCATTGTGCTTAACCTTGCCAGTTCAGGCAGATACTTAAACTGGCGTATTTTCTAACATAGTTCAATAAGTTCAGTTGTTCCAGGCCAGTGGGAGAAGTTATTACATAGTGCGTGCAATATCACAtttt                                                                                                                                      <  1:1772073/150‑1 (MQ=255)
                   gtaaTACGCGCCGATATATCATTGTGCTTAACCTTGCCAGTTCAGGCAGATACTTAAACTGGCGTATTTTCTAACATAGTTCAATAAGTTCAGTTGTTCCAGGCCAGTGGGAGAAGTTATTACATAGTGCGTGCAATATCACATTTTTTg                                                                                                                                   <  2:3557000/148‑1 (MQ=255)
                          gcgcCGATATATCATTGTGCTTAACCTTGCCAGTTCAGGCAGATACTTAAACTGGCGTATTTTCTAACATAGTTCAATAAGTTCAGTTGTTCCAGGCCAGTGGGAGAAGTTATTACATAGTGCGTGCAATATCACATTTTTTGCTATGCa                                                                                                                            >  2:3969730/1‑150 (MQ=255)
                             ccGATATATCATTGTGCTTAACCTTGCCAGTTCAGGCAGATACTTAAACTGGCGTATTTTCTAACATAGTTCAATAAGTTCAGTTGTTCCAGGCCAGTGGGAGAAGTTATTACATAGTGCGTGCAATATCACATTTTTTGCTATGCAATg                                                                                                                         <  2:4430508/150‑1 (MQ=255)
                               gATATATCATTGTGCTTAACCTTGCCAGTTCAGGCAGATACTTAAACTGGCGTATTTTCTAACATAGTTCAATAAGTTCAGTTGTTCCAGGCCAGTGGGAGAAGTTATTACATAGTGCGTGCAATATCACATTTTTTGCTATGCAATGaa                                                                                                                       >  1:4305712/1‑150 (MQ=255)
                                atatatCATTGTGCTTAACCTTGCCAGTTCAGGCAGATACTTAAACTGGCGTATTTTCTAACATAGTTCAATAAGTTCAGTTGTTCCAGGCCAGTGGGAGAAGTTATTACATAGTGCGTGCAATATCACATTTTTTGCTATGCAATGAAt                                                                                                                      >  2:1076528/1‑150 (MQ=255)
                                 tatatCATTGTGCTTAACCTTGCCAGTTCAGGCAGATACTTAAACTGGCGTATTTTCTAACATAGTTCAATAAGTTCAGTTGTTCCAGGCCAGTGGGAGAAGTTATTACATAGTGCGTGCAATATCACATTTTTTGCTATGCAATGAATa                                                                                                                     <  2:2767698/150‑1 (MQ=255)
                                  atatCATTGTGCTTAACCTTGCCAGTTCAGGCAGATACTTAAACTGGCGTATTTTCTAACATAGTTCAATAAGTTCAGTTGTACCAGGCCAGTGGGAGAAGTTATTACATAGTGCGTGCAATATCACATTTTTTGCTATGCAATGAATaa                                                                                                                    >  2:753692/1‑150 (MQ=255)
                                    atCATTGTGCTTAACCTTGCCAGTTCAGGCAGATACTTAAACTGGCGTATTTTCTAACATAGTTCAATAAGTTCAGTTGTTCCAGGCCAGTGGGAGAAGTTATTACATAGTGCGTGCAATATCACATTTTTTGCTATGCAATGAATaaaa                                                                                                                  <  1:2176134/150‑1 (MQ=255)
                                      cATTGTGCTTAACCTTGCCAGTTCAGGCAGATACTTAAACTGGCGTATTTTCTAACATAGTTCAATAAGTTCAGTTGTTCCAGGCCAGTGGGAGAAGTTATTACATAGTGCGTGCAATATCACATTTTTTGCTATGCAATGAATAAAAAg                                                                                                                <  2:2668824/150‑1 (MQ=255)
                                        ttGTGCTTAACCTTGCCAGTTCAGGCAGATACTTAAACTGGCGTATTTTCTAACATAGTTCAATAAGTTCAGTTGTTCCAGGCCAGTGGGAGAAGTTATTACATAGTGAGTGCAATATCACATTTTTTGCTATGCAATGAATAAAAAGtt                                                                                                              <  2:2284301/150‑1 (MQ=255)
                                           tgCTTAACCTTGCCAGGTCAGGCAGATACTTAAACTGGCGTATTTTCTAACATAGTTCAATAAGTTCAGTTGTTCCAGGCCAGTGGGAGAAGTTATTACATAGTGCGTGCAATATCACATTTTTTGCTATGCAATGAATAAAAAGTtata                                                                                                           <  2:606035/150‑1 (MQ=255)
                                            gCTTAACCTTGCCAGTTCAGGCAGATACTTAAACTGGCGTATTTTCTAACATAGTTCAATAAGTTCAGTTGTTCCAGGCCAGTGGGAGAAGTTATTACATAGTGCGTGCAATATCACATTTTTTGCTATGCAATGAATAAAAAGTTatat                                                                                                          <  1:508415/150‑1 (MQ=255)
                                             cTTAACCTTGCCAGTTCAGGCAGATACTTAAACTGGCGTATTTTCTAACATAGTTCAATAAGTTCAGTTGTTCCAGGCCAGTGGGAGAAGTTATTACATAGTGCGTGCAATATCACATTTTTTGCTATGCAATGAATAAAAAGTTATATc                                                                                                         >  2:1288914/1‑150 (MQ=255)
                                                aaCCTTGCCAGTTCAGGCAGATACTTAAACTGGCGTATTTTCTAACATAGTTCAATAAGTTCAGTTGTTCCAGGCCAGTGGGAGAAGTTATTACATAGTGCGTGCAATATCACATTTTTTGCTATGCAATGAATAAAAAGTTATATCACt                                                                                                      <  2:2341408/150‑1 (MQ=255)
                                                  ccTTGCCAGTTCAGGCAGATACTTAAACTGGCGTATTTTCTAACATAGTTCAATAAGTTCAGTTGTTCCAGGCCAGTGGGAGAAGTTATTACATAGTGCGTGCAATATCACATTTTTTGCTATGCAATGAATAAAAAGTTATATCACttt                                                                                                    <  1:4538403/150‑1 (MQ=255)
                                                       caaGTTCAGGCAGATACTTAAACTGGCGTATTTTCTAACATAGTTCAATAAGTTCAGTTGTTCCAGGCCAGTGGGAGAAGTTATTACATAGTGCGTGCAATATCACATTTTTTGCTATGCAATGAATAAAAAGTTATATCACTTTTtctc                                                                                               <  2:3011529/148‑1 (MQ=255)
                                                        cAGTTCAGGCAGATACTTAAACTGGCGTATTTTCTAACATAGTTCAATAAGTTCAGTTGTTCCAGGCCAGTGGGAGAAGTTATTACATAGTGCGTGCAATATCACATTTTTTGCTATGCAATGAATAAAAAGTTATATCACTTTTTCTCa                                                                                              <  2:1936254/150‑1 (MQ=255)
                                                              aGGCAGATACTTAAACTGGCGTATTTTCTAACATAGTTCAATAAGTTCAGTTGTTCCAGGCCAGTGGGAGAAGTTATTACATAGTGCGTGCAATATCACATTTTTTGCTATGCAATGAATAAAAAGTTATATCACTTTTTCTCATAAAAc                                                                                        <  1:3971917/150‑1 (MQ=255)
                                                              aGGCAGATACTTAAACTGGCGTATTTTCTAACATAGTTCAATAAGTTCAGTTGTTCCAGGCCAGGGGGAGAAGTTATTACATAGTGCGTGCAATATCACATTTTTTGCTATGCAATGAATAAAAAGTTATATCACTTTTTCTCTTAAAAc                                                                                        >  1:1263057/1‑150 (MQ=255)
                                                               ggCAGATACTTAAACTGGCGTATTTTCTAACATAGTTCAATAAGTTCAGTTGTTCCAGGCCAGTGGGAGAAGTTATTACATAGTGCGTGCAATATCACATTTTTTGCTATGCAATGAATAAAAAGTTATATCACTTTTTCTCATAAAAca                                                                                       >  2:947863/1‑150 (MQ=255)
                                                                 cAGATACTTAAACTGGCGTATTTTCTAACATAGTTCAATAAGTTCAGTTGTTCCAGGCCAGTGGGAGAAGTTATTACATAGTGCGTGCAATATCACATTTTTTGCTATGCAATGAATAAAAAGTTATATCACTTTTTCTCATAAAAcagt                                                                                     >  2:1832459/1‑150 (MQ=255)
                                                                          aaaCTGGCGTATTTTCTAACATAGTTCAATAAGTTCAGTTGTTCCAGGCCAGTGGGAGAAGTTATTACATAGTGCGTGCAATATCACATTTTTTGCTATGCAATGAATAAAAAGTTATATCACTTTTTCTCATAAAACAGTCAGTTAACg                                                                            >  1:1076774/1‑150 (MQ=255)
                                                                            aCTGGCGTATTTTCTAACATAGTTCAATAAGTTCAGTTGTTCCAGGCCAGTGGGAGAAGTTATTACATAGTGCGTGCAATATCACATTTTTTGCTATGCAATGAATAAAAAGTTATATCACTTTTTCTCATAAAACAGTCAGTTAACGGc                                                                          >  2:2428979/1‑150 (MQ=255)
                                                                             cTGGCGTATTTTCTAACATAGTTCAATAAGTTCAGTTGTTCCAGGCCAGTGGGAGAAGTTATTACATAGTGCGTGCAATATCACATTTTTTGCTATGCAATGAATAAAAAGTTATATCACTTTTTCTCATAAAACAGTCAGTTAACGGCt                                                                         >  1:1894458/1‑150 (MQ=255)
                                                                               ggCGTATTTTCTAACATAGTTCAATAAGTTCAGTTGTTCCAGGCCAGTGGGAGAAGTTATTACATAGTTCGTGCATTATCACATTTTTTGCTATGCAATGAATAAAAAGTTATATCACTTTTTCTCATAAAACAGTCAGTTAACGGCTat                                                                       <  1:3700464/150‑1 (MQ=255)
                                                                               ggCGTATTTTCTAACATAGTTCAATAAGTTCAGTTGTTCCAGGCCAGTGGGAGAAGTTATTACATAGTGCGTGCAATATCACATTTTTTGCTATGCAATGAATAAAAAGTTATATCACTTTTTCTCATAAAACAGTCAGTTAACGGCTat                                                                       <  1:1886686/150‑1 (MQ=255)
                                                                                 cGTATTTTCTAACATAGTTCAATAAGTTCAGTTGTTCCAGGCCAGTGGGAGAAGTTATTACATAGTGCGTGCAATATCACATTTTTTGCTATGCAATGAATAAAAAGTTATATCACTTTTTCTCATAAAACAGTCAGTTAACGGCTatta                                                                     >  1:4211476/1‑150 (MQ=255)
                                                                                 cGTATTTTCTAACATAGTTCAATAAGTTCAGTTGTTCCAGGCCAGTGGGAGAAGTTATTACATAGTGCGTGCAATATCACATTTTTTGCTATGCAATGAATAAAAAGTTATATCACTTTTTCTCATAAAACAGTCAGTTAACGGCTatta                                                                     >  2:81866/1‑150 (MQ=255)
                                                                                 cGTATTTTCTAACATAGTTCAATAAGTTCAGTTGTTCCAGGCCAGTGGGAGAAGTTATTACATAGTGCGTGCAATATCACATTTTTTGCTATGCAATGAATAAAAAGTTATATCACTTTTTCTCATAAAACAGTCAGTTAACGGCTatta                                                                     <  2:708668/150‑1 (MQ=255)
                                                                                  gTATTTTCTAACATAGTTCAATAAGTTCAGTTGTTCCAGGCCAGTGGGAGAAGTTATTACATAGTGCGTGCAATATCACATTTTTTGCTATGCAATGAATAAAAAGTTATATCACTTTTTCTCATAAAACAGTCAGTTAACGGCTattaa                                                                    >  2:1488893/1‑150 (MQ=255)
                                                                                       ttCTAACATAGTTCAATAAGTTCAGTTGTTCCAGGCCAGTGGGAGAAGTTATTACATAGTGCGTGCAATATCACATTTTTTGCTATGCAATGAATAAAAAGTTATATCACTTTTTCTCATAAAACAGTCAGTTAACGGCTATTAATTAcc                                                               <  1:4277952/150‑1 (MQ=255)
                                                                                         cTAACATAGTTCAATAAGTTCAGTTGTTCCAGGCCAGTGGGAGAAGTTATTACATAGTGCGTGCAATATCACATTTTTTGCTATGCAATGAATAAAAAGTTATATCACTTTTTCTCATAAAACAGTCAGTTAACGGCTATTAATTACCCt                                                             <  2:1198013/150‑1 (MQ=255)
                                                                                           aaCATAGTTCAATAAGTTCAGTTGTTCCAGGCCAGTGGGAGAAGTTATTACATAGTGCGTGCAATATCACATTTTTTGCTATGCAATGAATAAAAAGTTATATCACTTTTTCTCATAAAACAGTCAGTTAACGGCTATTAATTACCCTaa                                                           >  2:2881950/1‑150 (MQ=255)
                                                                                                 gTTCAATAAGTTCAGTTGTTCCAGGCCAGTGGGAGAAGTTATTACATAGTGCGTGCAATATCACATTTTTTGCTATGCAATGAATAAAAAGTTATATCACTTTTTCTCATAAAACAGTCAGTTAACGGCTATTAATTACCCTAAAGAGaa                                                     <  1:314237/150‑1 (MQ=255)
                                                                                                  ttCAATAAGTTCAGTTGTTCCAGGCCAGTGGGAGAAGTTATTACATAGTGCGTGCAATATCACATTTTTTGCTATGCAATGAATAAAAAGTTATATCACTTTTTCTCATAAAACAGTCAGTTAACGGCTATTAATTACCCTAAAGAGAAg                                                    <  1:4406664/150‑1 (MQ=255)
                                                                                                  ttCAATAAGTTCAGTTGTTCCAGGCCAGTGGGAGAAGTTATTACATAGTGCGTGCAATATCACATTTTTTGCTATGCAATGAATAAAAAGTTATATCACTTTTTCTCATAAAACAGTCAGTTAACGGCTATTAATTACCCTAAAGAGAAg                                                    <  1:3818916/150‑1 (MQ=255)
                                                                                                   tCAATAAGTTCAGTTGTTCCAGGCCAGTGGGAGAAGTTATTACATAGTGCGTGCAATATCACATTTTTTGCTATGCAATGAATAAAAAGTTATATCACTTTTTCTCATAAAACAGTCAGTTAACGGCTATTAATTACCCTAAAGAGAAGt                                                   <  1:2293485/150‑1 (MQ=255)
                                                                                                    cAATAAGTTCAGTTGTTCCAGGCCAGTGGGAGAAGTTATTACATAGTGCGTGCAATATCACATTTTTTGCTATGCAATGAATAAAAAGTTATATCACTTTTTCTCATAAAACAGTCAGTTAACGGCTATTAATTACCCTAAAGAGAAGTc                                                  <  1:1429048/150‑1 (MQ=255)
                                                                                                      aTAAGTTCAGTTGTTCCAGGCCAGTGGGAGAAGTTATTACATAGTGCGTGCAATATCACATTTTTTGCTATGCAATGAATAAAAAGTTATATCACTTTTTCTCATAAAACAGTCAGTTAACGGCTATTAATTACCCTAAAGAGAAGTCaa                                                >  2:3349392/1‑150 (MQ=255)
                                                                                                      aTAAGTTCAGTTGTTCCAGGCCAGTGGGAGAAGTTATTACATAGTGCGTGCAATATCACATTTTTTGCTATGCAATGAATAAAAAGTTATATCACTTTTTCTCATAAAACAGTCAGTTAACGGCTATTAATTACCCTAAAGAGAAGTCaa                                                >  2:3480768/1‑150 (MQ=255)
                                                                                                      aTAAGTTCAGTTGTTCCAGGCCAGTGGGAGAAGTTATTACATAGTGCGTGCAATATCACATTTTTTGCTATGCAATGAATAAAAAGTTATATCACTTTTTCTCATAAAACAGTCAGTTAACGGCTATTAATTACCCTAAAGAGAAGTCaa                                                <  1:2737000/150‑1 (MQ=255)
                                                                                                           ttCAGTTGTTCCAGGCCAGTGGGAGAAGTTATTACATAGTGCGTGCAATATCACATTTTTTGCTATGCAATGAATAAAAAGTTATATCACTTTTTCTCATAAAACAGTCAGTTAACGGCTATTAATTACCCTAAAGAGAAGTCAATcccc                                           <  1:102357/150‑1 (MQ=255)
                                                                                                           ttCAGTTGTTCCAGGCCAGTGGGAGAAGTTATTACATAGTGCGTGCAATATCACATTTTTTGCTATGCAATGAATAAAAAGTTATATCACTTTTTCTCATAAAACAGTCAGTTAACGGCTATTAATTACCCTAAAGAGAAGTCAATcccc                                           <  1:2375239/150‑1 (MQ=255)
                                                                                                             cAGTTGTTCCAGGCCAGTGGGAGAAGTTATTACATAGTGCGTGCAATATCACATTTTTTGCTATGCAATGAATAAAAAGTTATATCACTTTTTCTCATAAAACAGTCAGTTAACGGCTATTAATTACCCTAAAGAGAAGTCAATCCCCaa                                         >  1:515175/1‑150 (MQ=255)
                                                                                                               gttgttCCAGGCCAGTGGGAGAAGTTATTACATAGTGCGTGCAATATCACATTTTTTGCTATGCAATGAATAAAAAGTTATATCACTTTTTCTCATAAAACAGTCAGTTAACGGCTATTAATTACCCTAAAGAGAAGTCAATCCCCAAAg                                       <  2:1124657/150‑1 (MQ=255)
                                                                                                               gttgttCCAGGCCAGTGGGAGAAGTTATTACATAGTGCGTGCAATATCACATTTTTTGCTATGCAATGAATAAAAAGTTATATCACTTTTTCTCATAAAACAGTCAGTTAACGGCTATTAATTACCCTAAAGAGAAGTCAATCCCCAAAg                                       >  2:151732/1‑150 (MQ=255)
                                                                                                                 tgttCCAGGCCAGTGGGAGAAGTTATTACATAGTGCGTGCAATATCACATTTTTTGCTATGCAATGAATAAAAAGTTATATCACTTTTTCTCATAAAACAGTCAGTTAACGGCTATTAATTACCCTAAAGAGAAGTCAATCCCCAAAggg                                     >  2:876555/1‑150 (MQ=255)
                                                                                                                    tCCAGGCCAGTGGGAGAAGTTATTACATAGTGCGTGCAATATCACATTTTTTGCTATGCAATGAATAAAAAGTTATATCACTTTTTCTCATAAAACAGTCAGTTAACGGCTATTAATTACCCTAAAGAGAAGTCAATCCCCAAAGGGAtt                                  <  1:1537150/150‑1 (MQ=255)
                                                                                                                     ccAGGCCAGTGGGAGAAGTTATTACATAGTGCGTGCAATATCACATTTTTTGCTATGCAATGAATAAAAAGTTATATCACTTTTTCTCATAAAACAGTCAGTTAACGGCTATTAATTACCCTAAAGAGAAGTCAATCCCCAAAGGGATTg                                 <  2:3460457/150‑1 (MQ=255)
                                                                                                                     ccAGGCCAGTGGGAGAAGTTATTACATAGTGCGTGCAATATCACATTTTTTGCTATGCAATGAATAAAAAGTTATATCACTTTTTCTCATAAAACAGTCAGTTAACGGCTATTAATTACCCTAAAGAGAAGTCAATCCCCAAAGGGATTg                                 <  2:3460209/150‑1 (MQ=255)
                                                                                                                     ccAGGCCAGTGGGAGAAGTTATTACATAGTGCGTGCAATATCACATTTTTTGCTATGCAATGAATAAAAAGTTATATCACTTTTTCTCATAAAACAGTCAGTTAACGGCTATTAATTACCCTAAAGAGAAGTCAATCCCCAAAGGGATTg                                 <  2:3459691/150‑1 (MQ=255)
                                                                                                                      cAGGCCAGTGGGAGAAGTTATTACATAGTGCGTGCTATATCACATTTTTTGCTATGCAATGAATAAAAAGTTATATCACTTTTTCTCATAAAACAGTCAGTTAACGGCTATTAATTACCCTAAAGAGAAGTCAATCCCCAAAGGGATTGt                                <  1:3674569/150‑1 (MQ=255)
                                                                                                                       aGGCCAGTGGGAGAAGTTATTACATAGTGCGTGCAATATCACATTTTTTGCTATGCAATGAATAAAAAGTTATATCACTTTTTCTCATAAAACAGTCAGTTAACGGCTATTAATTACCCTAAAGAGAAGTCAATCCCCAAAGGGATTGTa                               >  2:3544173/1‑150 (MQ=255)
                                                                                                                        ggCCAGTGGGAGAAGTTATTACATAGTGCGTGCAATATCACATTTTTTGCTATGCCATGAATAAAAAGTTATATCACTTTTTCTCATAAAACAGTCAGTTAACCGCTATTAATTACCCTAAAGAGAAGTCAATCCCCAAAGGGCTTGTaa                              >  2:3896614/1‑150 (MQ=255)
                                                                                                                         gCCAGTGGGAGAAGTTATTACATAGTGCGTGCAATATCACATTTTTTGCTATGCAATGAATAAAAAGTTATATCACTTTTTCTCATAAAACAGTCAGTTAACGGCTATTAATTACCCTAAAGAGAAGTCAATCCCCAAAGGGATTGTaaa                             >  1:3530339/1‑150 (MQ=255)
                                                                                                                         gCCAGTGGGAGAAGTTATTACATAGTGCGTGCAATATCACATTTTTTGCTATGCAATGAATAAAAAGTTATATCACTTTTTCTCATAAAACAGTCAGTTAACGGCTATTAATTACCCTAAAGAGAAGTCAATCCCCAAAGGGATTGTaaa                             <  1:3971664/150‑1 (MQ=255)
                                                                                                                          ccAGTGGGAGAAGTTATTACATAGTGCGTGCAATATCACATTTTTTGCTATGCAATGAATAAAAAGTTATATCACTTTTTCTCATAAAACAGTCAGTTAACGGCTATTAATTACCCTAAAGAGAAGTCAATCCCCAAAGGGATTGTAAAt                            <  1:1729433/150‑1 (MQ=255)
                                                                                                                           cAGTGGGAGAAGTTATTACATAGTGCGTGCAATATCACATTTTTTGCTATGCAATGAATAAAAAGTTATATCACTTTTTCTCATAAAACAGTCAGTTAACGGCTATTAATTACCCTAAAGAGAAGTCAATCCCCAAAGGGATTGTAAAtt                           <  2:1586030/150‑1 (MQ=255)
                                                                                                                           cAGTGGGAGAAGTTATTACATAGTGCGTGCAATATCACATTTTTTGCTATGCAATGAATAAAAAGTTATATCACTTTTTCTCATAAAACAGTCAGTTAACGGCTATTAATTACCCTAAAGAGAAGTCAATCCCCAAAGGGATTGTAAAtt                           <  1:1926564/150‑1 (MQ=255)
                                                                                                                            aGTGGGAGAAGTTATTACATAGTGCGTGCAATATCACATTTTTTGCTATGACATGAATAAAAAGTTATATCACTTTTTCTCATAAAACATTCAGTTAACGGCTATTAATTACCCTAAAGAGAAGTCAATCCCCAAAGGGATTGTAAAttt                          >  1:2926971/1‑150 (MQ=255)
                                                                                                                                  agaAGTTATTACATAGTGCGTGCAATATCACATTTTTTGCTATGCAATGAATAAAAAGTTATATCACTTTTTCTCATAAAACAGTCAGTTAACGGCTATTAATTACCCTAAAGAGAAGTCAATCCCCAAAGGGATTGTAAATTTAAAATa                    >  2:3247588/1‑150 (MQ=255)
                                                                                                                                   gaAGTTATTACATAGTGCGTGCAATATCACATTTTTTGCTATGCAATGAATAAAAAGTTATATCACTTTTTCTCATAAAACAGTCAGTTAACGGCTATTAATTACCCTAAAGAGAAGTCAATCCCCAAAGGGATTGTAAATTTAAAATaa                   >  2:512206/1‑150 (MQ=255)
                                                                                                                                   gaAGTTATTACATAGTGCGTGCAATATCACATTTTTTGCTATGCAATGAATAAAAAGTTATATCACTTTTTCTCATAAAACAGTCAGTTAACGGCTATTAATTACCCTAAAGAGAAGTCAATCCCCAAAGGGATTGTAAATTTAAAATaa                   >  2:512247/1‑150 (MQ=255)
                                                                                                                                      gTTATTACATAGTGCGTGCAATATCACATTTTTTGCTATGCAATGAATAAAAAGTTATATCACTTTTTCTCATAAAACAGTCAGTTAACGGCTATTAATTACCCTAAAGAGAAGTCAATCCCCAAAGGGATTGTAAATTTAAAATAAGaa                >  1:3745757/1‑150 (MQ=255)
                                                                                                                                      gTTATTACATAGTGCGTGCAATATCACATTTTTTGCTATGCAATGAATAAAAAGTTATATCACTTTTTCTCATAAAACAGTCAGTTAACGGCTATTAATTACCCTAAAGAGAAGTCAATCCCCAAAGGGATTGTAAATTTAAAATAAGaa                <  1:2210250/150‑1 (MQ=255)
                                                                                                                                       ttattaCATAGTGCGTGCAATATCACATTTTTTGCTATGCAATGAATAAAAAGTTATATCACTTTTTCTCATAAAACAGTCAGTTAACGGCTATTAATTACCCTAAAGAGAAGTCAATCCCCAAAGGGATTGTAAATTTAAAATAAGaaa               >  2:1625816/1‑150 (MQ=255)
                                                                                                                                       ttattaCATAGTGCGTGCAATATCACATTTTTTGCTATGCAATGAATAAAAAGTTATATCACTTTTTCTCATAAAACAGTCAGTTAACGGCTATTAATTACCCTAAAGAGAAGTCAATCCCCAAAGGGATTGTAAATTTAAAATAAGaaa               >  2:4163132/1‑150 (MQ=255)
                                                                                                                                         attaCATAGTGCGTGCAATATCACATTTTTTGCTATGCAATGAATAAAAAGTTATATCACTTTTTCTCATAAAACAGTCAGTTAACGGCTATTAATTACCCTAAAGAGAAGTCAATCCCCAAAGGGATTGTAAATTTAAAATAAGaaaaa             <  2:1726427/150‑1 (MQ=255)
                                                                                                                                         attaCATAGTGCGTGCAATATCACATTTTTTGCTATGCAATGAATAAAAAGTTATATCACTTTTTCTCATAAAACAGTCAGTTAACGGCTATTAATTACCCTAAAGAGAAGTCAATCCCCAAAGGGATTGTAAATTTAAAATAAGaaaaa             <  2:2517098/150‑1 (MQ=255)
                                                                                                                                         attaCATAGTGCGTGCAATATCACATTTTTTGCTATGCAATGAATAAAAAGTTATATCACTTTTTCTCATAAAACAGTCAGTTAACGGCTATTAATTACCCTAAAGAGAAGTCAATCCCCAAAGGGATTGTAAATTTAAAATAAGaaaaa             >  1:3257895/1‑150 (MQ=255)
                                                                                                                                           taCATAGTGAGTGCAATATCACATTTTTTGCTATGCAATGAATAAAAAGTTATATCACTTTTTCTCATAAAACAGTCAGTTAACGGCTATTAATTACCCTAAAGAGAAGTCAATCCCCAAAGGGATTGTAAATTTAAAATAAGAAAAAtt           <  1:974916/150‑1 (MQ=255)
                                                                                                                                            aCATAGTGCGTGCAATATCACATTTTTTGCTATGCAATGAATAAAAAGTTATATCACTTTTTCTCATAAAACAGTCAGTTAACGGCTATTAATTACCCTAAAGAGAAGTCAATCCCCAAAGGGATTGTAAATTTAAAATAAGAAAAATtg          <  1:3897851/150‑1 (MQ=255)
                                                                                                                                            aCATAGTGCGTGCAATATCACATTTTTTGCTATGCAATGAATAAAAAGTTATATCACTTTTTCTCATAAAACAGTCAGTTAACGGCTATTAATTACCCTAAAGAGAAGTCAATCCCCAAAGGGATTGTAAATTTAAAATAAGAAAAATtg          <  2:1773341/150‑1 (MQ=255)
                                                                                                                                              aTAGTGCGTTCAATATCACATTTTTTGCTATGCAATGAATAAAAAGTTATATCACTTTTTCTCATAAAACAGTCAGTTAACGGCTATTAATTACCCTAAAGAGAAGTCAATCCCCAAAGGGATTGTAAATTTAAAATAAGAAAAATtgat        <  1:663504/150‑1 (MQ=255)
                                                                                                                                              aTAGTGCGTGCAATATCACATTTTTTGCTATGCAATGAATAAAAAGTTATATCACTTTTTCTCATAAAACAGTCAGTTAACGGCTATTAATTACCCTAAAGAGAAGTCAATCCCCAAAGGGATTGTAAATTTAAAATAAGAAAAATtgat        <  1:3666924/150‑1 (MQ=255)
                                                                                                                                               tAGTGCGTGCAATATCACATTTTTTGCTATGCAATGAATAAAAAGTTATATCACTTTTTCTCATAAAACAGTCAGTTAACGGCTATTAATTACCCTAAAGAGAAGTCAATCCCCAAAGGGATTGTAAATTTAAAATAAGAAAAATtgatg       >  1:1581087/1‑150 (MQ=255)
                                                                                                                                               tAGTGCGTGCAATATCACATTTTTTGCTATGCAATGAATAAAAAGTTATATCACTTTTTCTCATAAAACAGTCAGTTAACGGCTATTAATTACCCTAAAGAGAAGTCAATCCCCAAAGGGATTGTAAATTTAAAATAAGAAAAATtgatg       <  2:777875/150‑1 (MQ=255)
                                                                                                                                                aGTGCGTGCAATATCACATTTTTTGCTATGCAATGAATAAAAAGTTATATCACTTTTTCTCATAAAACAGTCAGTTAACGGCTATTAATTACCCTAAAGAGAAGTCAATCCCCAAAGGGATTGTAAATTTAAAATAAGAAAAATtgatga      >  1:1495180/1‑150 (MQ=255)
                                                                                                                                                aGTGCGTGCAATATCACATTTTTTGCTATGCAATGAATAAAAAGTTATATCACTTTTTCTCATAAAACAGTCAGTTAACGGCTATTAATTACCCTAAAGAGAAGTCAATCCCCAAAGGGATTGTAAATTTAAAATAAGAAAAATtgatga      >  1:886658/1‑150 (MQ=255)
                                                                                                                                                 gtgcgtgcAATATCACATTTTTTGCTATGCAATGAATAAAAAGTTATATCACTTTTTCTCATAAAACAGTCAGTTAACGGCTATTAATTACCCTAAAGAGAAGTCAATCCCCAAAGGGATTGTAAATTTAAAATAAGAAAAATTGatgaa     <  2:221102/150‑1 (MQ=255)
                                                                                                                                                    cgtgcAATATCACATTTTTTGCTATGCAATGAATAAAAAGTTATATCACTTTTTCTCATAAAACAGTCAGTTAACGGCTATTAATTACCCTAAAGAGAAGTCAATCCCCAAAGGGATTGTAAATTTAAAATAAGAAAAATTGatgaatga  <  2:4241062/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
TTTTGCACACTTATGCAACGGAATACGCGCCGATATATCATTGTGCTTAACCTTGCCAGTTCAGGCAGATACTTAAACTGGCGTATTTTCTAACATAGTTCAATAAGTTCAGTTGTTCCAGGCCAGTGGGAGAAGTTATTACATAGTGCGTGCAATATCACATTTTTTGCTATGCAATGAATAAAAAGTTATATCACTTTTTCTCATAAAACAGTCAGTTAACGGCTATTAATTACCCTAAAGAGAAGTCAATCCCCAAAGGGATTGTAAATTTAAAATAAGAAAAATTGATGAATGA  >  CP009273/1901967‑1902264

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: