Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,277,897 T→A 5.1% L69Q (CTG→CAG)  nrfA → nitrite reductase, formate‑dependent, cytochrome

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,277,8970TA 5.1% 984.6 / 31.2 274L69Q (CTG→CAG) nrfAnitrite reductase, formate‑dependent, cytochrome
Reads supporting (aligned to +/- strand):  ref base T (147/111);  new base A (11/3);  total (158/116)
Fisher's exact test for biased strand distribution p-value = 1.64e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.03e-01

Only 100 of 277 total aligned reads displayed.

CTTTTTCACTTCTGTTCACGCTGAACAAACGGCTGCTCCCGCAAAACCTGTAACTGTGGAAGCGAAGAATGAAACCTTTGCCCCGCAGCATCCCGATCAATATCTCTCCTGGAAAGCCACCTCGGAACAGTCAGAGCGTGTTGACGCCCTGGCGGAAGATCCACGGCTGGTGATCCTGTGGGCGGGGTATCCCTTCTCGCGCGATTACAACAAGCCGCGTGGACATGCTTTTGCTGTGACCGATGTGCGTGAAACCCTGCGTACCGGTGCGCCGAAAAACGCTGAAGATGGTCCGCTAC  >  CP009273/4277748‑4278046
                                                                                                                                                     |                                                                                                                                                     
ctttttcACTTCTGTTCACGCTGAACAAACGGCTGCTCCCGCAAAACCTGTAACTGTGGAAGCGAAGAATGAAACCTTTGCCCCGCAGCATCCCGATCAATATCTCTCCTGGAAAGCCACCTCGGAACAGTCAGAGCGTGTTGACGCCCt                                                                                                                                                       <  1:646602/150‑1 (MQ=255)
     tcACTTCTGTTCACGCTGAACAAACGGCTGCTCCCGCAAAACCTGTAACTGTGGAAGCGAAGAATGAAACCTTTGCCCCGCAGCATCCCGATCAATATCTCTCCTGGAAAGCCACCTCGGAACAGTCAGAGCGTGTTGACGCCCTGGCgg                                                                                                                                                  >  1:3394421/1‑150 (MQ=255)
      cACTTCTGTTCACGCTGAACAAACGGCTGCTCCCGCAAAACCTGTAACTGTGGAAGCGAAGAATGAAACCTTTGCCCCGCAGCATCCCGATCAATATCTCTCCTGGAAAGCCACCTCGGAACAGTCAGAGCGTGTTGACGCCCTGGCGGa                                                                                                                                                 <  1:3048038/150‑1 (MQ=255)
        cTTCTGTTCACGCTGAACAAACGGCTGCTCCCGCAAAACCTGTAACTGTGGAAGCGAAGAATGAAACCTTTGCCCCGCAGCATCCCGATCAATATCTCTCCTGGAAAGCCACCTCGGAACAGTCAGAGCGTGTTGACGCCCTGGCGGAAg                                                                                                                                               >  2:3395861/1‑150 (MQ=255)
         ttCTGTTCACGCTGAACAAACGGCTGCTCCCGCAAAACCTGTAACTGTGGAAGCGAAGAATGAAACCTTTGCCCCGCAGCATCCCGATCAATATCTCTCCTGGAAAGCCACCTCGGAACAGTCAGAGCGTGTTGACGCCCTGGCGGAAGa                                                                                                                                              <  2:831133/150‑1 (MQ=255)
               tCACGCTGAACAAACGGCTGCTCCCGCAAAACCTGTAACTGTGGAAGCGAAGAATGAAACCTTTGCCCCGCAGCATCCCGATCAATATCTCTCCTGGAAAGCCACCTCGGAACAGTCAGAGCGTGTTGACGCCCTGGCGGAAGATCCACg                                                                                                                                        >  2:3173097/1‑150 (MQ=255)
                cACGCTGAACAAACGGCTGCTCCCGCAAAACCTGTAACTGTGGAAGCGAAGAATGAAACCTTTGCCCCGCAGCATCCCGATCAATATCTCTCCTGGAAAGCCACCTCGGAACAGTCAGAGCGTGTTGACGCCCTGGCGGAAGATCCACgg                                                                                                                                       >  1:2861771/1‑150 (MQ=255)
                   gCTGAACAAACGGCTGCTCCCGCAAAACCTGTAACTGTGGAAGCGAAGAATGAAACCTTTGCCCCGCAGCATCCCGATCAATATCTCTCCTGGAAAGCCACCTCGGAACAGTCAGAGCGTGTTGACGCCCTGGCGGAAGATCCACGGCTg                                                                                                                                    >  1:179548/1‑150 (MQ=255)
                    cTGAACAAACGGCTGCTCCCGCAAAACCTGTAACTGTGGAAGCGAAGAATGAAACCTTTGCCCCGCAGCATCCCGATCAATATCTCTCCTGGAAAGCCACCTCGGAACAGTCAGAGCGTGTTGACGCCCTGGCGGAAGATCCACGGCTgg                                                                                                                                   <  1:134796/150‑1 (MQ=255)
                      gAACAAACGGCTGCTCCCGCAAAACCTGTAACTGTGGAAGCGAAGAATGAAACCTTTGCCCCGCAGCATCCCGATCAATATCTCTCCTGGAAAGCCACCTCGGAACAGTCAGAGCGTGTTGACGCCCTGGCGGAAGATCCACGGCTGGTg                                                                                                                                 >  2:189193/1‑150 (MQ=255)
                          aaaCGGCTGCTCCCGCAAAACCTGTAACTGTGGAAGCGAAGAATGAAACCTTTGCCCCGCAGCATCCCGATCAATATCTCTCCTGGAAAGCCACCTCGGAACAGTCAGAGCGTGTTGACGCCCTGGCGGAAGATCCACGGCTGGTGATcc                                                                                                                             >  2:1631103/1‑150 (MQ=255)
                           aaCGGCTGCTCCCGCAAAACCTGTAACTGTGGAAGCGAAGAATGAAACCTTTGCCCCGCAGCATCCCGATCAATATCTCTCCTGGAAAGCCACCTCGGAACAGTCAGAGCGTGTTGACGCCCTGGCGGAAGATCCACGGCTGGTGATCCt                                                                                                                            >  2:3056054/1‑150 (MQ=255)
                            aCGGCTGCTCCCGCAAAACCTGTAACTGTGGAAGCGAAGAATGAAACCTTTGCCCCGCAGCATCCCGATCAATATCTCTCCTGGAAAGCCACCTCGGAACAGTCAGAGCGTGTTGACGCCCTGGCGGAAGATCCACGGCTGGTGATCCtg                                                                                                                           <  1:244804/150‑1 (MQ=255)
                                ctgctCCCGCAAAACCTGTAACTGTGGAAGCGAAGAATGAAACCTTTGCCCCGCAGCATCCCGATCAATATCTCTCCTGGAAAGCCACCTCGGAACAGTCAGAGCGTGTTGACGCCCTGGCGGAAGATCCACGGCTGGTGATCCTGTggg                                                                                                                       <  1:1565615/150‑1 (MQ=255)
                                ctgctCCCGCAAAACCTGTAACTGTGGAAGCGAAGAATGAAACCTTTGCCCCGCAGCATCCCGATCAATATCTCTCCTGGAAAGCCACCTCGGAACAGTCAGAGCGTGTTGACGCCCTGGCGGAAGATCCACGGCTGGTGATCCTGTggg                                                                                                                       >  1:2790895/1‑150 (MQ=255)
                                  gctCCCGCAAAACCTGTAACTGTGGAAGCGAAGAATGAAACCTTTGCCCCGCAGCATCCCGATCAATATCTCTCCTGGAAAGCCACCTCGGAACAGTCAGAGCGTGTTGACGCCCTGGCGGAAGATCCACGGCTGGTGATCCTGTGGGCg                                                                                                                     <  1:103353/150‑1 (MQ=255)
                                  gctCCCGCAAAACCTGTAACTGTGGAAGCGAAGAATGAAACCTTTGCCCCGCAGCATCCCGATCAATATCTCTCCTGGAAAGCCACCTCGGAACAGTCAGAGCGTGTTGACGCCCAGGCGGAAGATCCACGGCTGGTGATCCTGTGGGCg                                                                                                                     >  2:2438033/1‑150 (MQ=255)
                                   ctCCCGCAAAACCTGTAACTGTGGAAGCGAAGAATGAAACCTTTGCCCCGCAGCATCCCGATCAATATCTCTCCTGGAAAGCCACCTCGGAACAGTCAGAGCGTGTTGACGCCCTGGCGGAAGATCCACGGCTGGTGATCCTGTGGGCgg                                                                                                                    >  2:2464220/1‑150 (MQ=255)
                                   ctCCCGCAAAACCTGTAACTGTGGAAGCGAAGAATGAAACCTTTGCCCCGCAGCATCCCGATCAATATCTCTCCTGGAAAGCCACCTCGGAACAGTCAGAGCGTGTTGACGCCCTGGCGGAAGATCCACGGCTGGTGATCCTGTGGGCgg                                                                                                                    <  1:1717379/150‑1 (MQ=255)
                                     cccGCAAAACCTGTAACTGTGGAAGCGAAGAATGAAACCTTTGCCCCGCAGCATCCCGATCAATATCTCTCCTGGAAAGCCACCTCGGAACAGTCAGAGCGTGTTGACGCCCTGGCGGAAGATCCACGGCTGGTGATCCTGTGGGCgggg                                                                                                                  <  1:1943610/150‑1 (MQ=255)
                                     cccGCAAAACCTGTAACTGTGGAAGCGAAGAATGAAACCTTTGCCCCGCAGCATCCCGATCAATATCTCTCCTGGAAAGCCACCTCGGAACAGTCAGAGCGTGTTGACGCCCCGGCGGAAGATCCACGGCTGGTGATCCTGTGGGCgggg                                                                                                                  <  2:572855/150‑1 (MQ=255)
                                       cGCAAAACCTGTAACTGTGGAAGCGAAGAATGAAACCTTTGCCCCGCAGCATCCCGATCAATATCTCTCCTGGAAAGCCACCTCGGAACAGTCAGAGCGTGTTGACGCCCTGGCGGAAGATCCACGGCTGGTGATCCTGTGGGCGGGGTa                                                                                                                <  1:3526851/150‑1 (MQ=255)
                                        gCAAAACCTGTAACTGTGGAAGCGAAGAATGAAACCTTTGCCCCGCAGCATCCCGATCAATATCTCTCCTGGAAAGCCACCTCGGAACAGTCAGAGCGTGTTGACGCCCTGGCGGAAGATCCACGGCTGGTGATCCTGTGGGCGGGGTAt                                                                                                               >  1:866211/1‑150 (MQ=255)
                                          aaaaCCTGTAACTGTGGAAGCGAAGAATGAAACCTTTGCCCCGCAGCATCCCGATCAATATCTCTCCTGGAAAGCCACCTCGGAACAGTCAGAGCGTGTTGACGCCCTGGCGGAAGATCCACGGCTGGTGATCCTGTGGGCGGGGTATcc                                                                                                             >  1:750149/1‑150 (MQ=255)
                                             aCCTGTAACTGTGGAAGCGAAGAATGAAACCTTTGCCCCGCAGCATCCCGATCAATATCTCTCCTGGAAAGCCACCTCGGAACAGTCAGAGCGTGTTGACGCCCTGGCGGAAGATCCACGGCTGGTGATCCTGTGGGCGGGGTATCCCtt                                                                                                          <  2:1335253/150‑1 (MQ=255)
                                               cTGTAACTGTGGAAGCGAAGAATGAAACCTTTGCCCCGCAGCATCCCGATCAATATCTCTCCTGGAAAGCCACCTCGGAACAGTCAGAGCGTGTTGACGCCCTGGCGGAAGATCCACGGCTGGTGATCCTGTGGGCGGGGTATCCCTtct                                                                                                        >  2:3668649/1‑150 (MQ=255)
                                               cTGTAACTGTGGAAGCGAAGAATGAAACCTTTGCCCCGCAGCATCCCGATCAATATCTCTCCTGGAAAGCCACCTCGGAACAGTCAGAGCGTGTTGACGCCCTGGCGGAAGATCCACGGCTGGTGATCCTGTGGGCGGGGTATCCCTtct                                                                                                        <  2:4269558/150‑1 (MQ=255)
                                               cTGTAACTGTGGAAGCGAAGAATGAAACCTTTGCCCCGCAGCATCCCGATCAATATCTCTCCTGGAAAGCCACCTCGGAACAGTCAGAGCGTGTTGACGCCCTGGCGGAAGATCCACGGCTGGTGATCCTGTGGGCGGGGTATCCCTtct                                                                                                        >  1:21705/1‑150 (MQ=255)
                                                tGTAACTGTGGAAGCGAAGAATGAAACCTTTGCCCCGCAGCATCCCGATCAATATCTCTCCTGGAAAGCCACCTCGGAACAGTCAGAGCGTGTTGACGCCCTGGCGGAAGATCCACGGCTGGTGATCCTGTGGGCGGGGTATCCCTtctc                                                                                                       <  1:1695244/150‑1 (MQ=255)
                                                  tAACTGTGGAAGCGAAGAATGAAACCTTTGCCCCGCAGCATCCCGATCAATATCTCTCCTGGAAAGCCACCTCGGAACAGTCAGAGCGTGTTGACGCCCTGGCGGAAGATCCACGGCTGGTGATCCTGTGGGCGGGGTATCCCTTCTcgc                                                                                                     >  2:3976299/1‑150 (MQ=255)
                                                  tAACTGTGGAAGCGAAGAATGAAACCTTTGCCCCGCAGCATCCCGATCAATATCTCTCCTGGAAAGCCACCTCGGAACAGTCAGAGCGTGTTGACGCCCTGGCGGAAGATCCACGGCTGGTGATCCTGTGGGCGGGGTATCCCTTCTcgc                                                                                                     >  1:4213445/1‑150 (MQ=255)
                                                     cTGTGGAAGCGAAGAATGAAACCTTTGCCCCGCAGCATCCCGATCAATATCTCTCCTGGAAAGCCACCTCGGAACAGTCAGAGCGTGTTGACGCCCTGGCGGAAGATCCACGGCTGGTGATCCTGTGGGCGGGGTATCCCTTCTcgcgcg                                                                                                  <  2:3809306/150‑1 (MQ=255)
                                                       gtgGAAGCGAAGAATGAAACCTTTGCCCCGCAGCATCCCGATCAATATCTCTCCTGGAAAGCCACCTCGGAACAGTCAGAGCGTGTTGACGCCCTGGCGGAAGATCCACGGCTGGTGATCCTGTGGGCGGGGTATCCCTTCTCGCGCGAt                                                                                                <  2:943293/150‑1 (MQ=255)
                                                        tgGAAGCGAAGAATGAAACCTTTGCCCCGCAGCATCCCGATCAATATCTCTCCTGGAAAGCCACCTCGGAACAGTCAGAGCGTGTTGACGCCCTGGCGGAAGATCCACGGCTGGTGATCCTGTGGGCGGGGTATCCCTTCTCGCGCGAtt                                                                                               <  2:2371591/150‑1 (MQ=255)
                                                           aaGCGAAGAATGAAACCTTTGCCCCGCAGCATCCCGATCAATATCTCTCCTGGAAAGCCACCTCGGAACAGTCAGAGCGTGTTGACGCCCTGGCGGAAGATCCACGGCTGGTGATCCTGTGGGCGGGGTATCCCTTCTCGCGCGATTaca                                                                                            <  2:4602681/150‑1 (MQ=255)
                                                            aGCGAAGAATGAAACCTTTGCCCCGCAGCATCCCGATCAATATCTCTCCTGGAAAGCCACCTCGGAACAGTCAGAGCGTGTTGACGCCCTGGCGGAAGATCCACGGCTGGTGATCCTGTGGGCGGGGTATCCCTTCTCGCGCGATTacaa                                                                                           >  1:4002934/1‑150 (MQ=255)
                                                             gCGAAGAATGAAACCTTTGCCCCGCAGCATCCCGATCAATATCTCTCCTGGAAAGCCACCTCGGAACAGTCAGAGCGTGTTGACGCCCTGGCGGAAGATCCACGGCTGGTGATCCTGTGGGCGGGGTATCCCTTCTCGCGCGATTacaac                                                                                          <  1:4338130/150‑1 (MQ=255)
                                                               gaagaaTGAAACCTTTGCCCCGCAGCATCCCGATCAATATCTCTCCTGGAAAGCCACCTCGGAACAGTCAGAGCGTGTTGACGCCCTGGCGGAAGATCCACGGCTGGTGATCCTGTGGGCGGGGTATCCCTTCTCGCGCGATTAcaacaa                                                                                        <  2:1843012/150‑1 (MQ=255)
                                                               gaagaaTGAAACCTTTGCCCCGCAGCATCCCGATCAATATCTCTCCTGGAAAGCCACCTCGGAACAGTCAGAGCGTGTTGACGCCCTGGCGGAAGATCCACGGCTGGTGATCCTGTGGGCGGGGTATCCCTTCTCGCGCGATTAcaacaa                                                                                        >  2:3307158/1‑150 (MQ=255)
                                                               gaagaaTGAAACCTTTGCCCCGCAGCATCCCGATCAATATCTCTCCTGGAAAGCCACCTCGGAACAGTCAGAGCGTGTTGACGCCCAGGCGGAAGATCCACGGCTGGTGATCCTGTGGGCGGGGTATCCCTTCTCGCGCGATTAcaacaa                                                                                        >  2:3103340/1‑150 (MQ=255)
                                                                aagaaTGAAACCTTTGCCCCGCAGCATCCCGATCAATATCTCTCCTGGAAAGCCACCTCGGAACAGTCAGAGCGTGTTGACGCCCTGGCGGAAGATCCACGGCTGGTGATCCTGTGGGCGGGGTATCCCTTCTCGCGCGATTACAACAAg                                                                                       >  1:1459694/1‑150 (MQ=255)
                                                                 agaaTGAAACCTTTGCCCCGCAGCATCCCGATCAATATCTCTCCTGGAAAGCCACCTCGGAACAGTCAGAGCGTGTTGACGCCCTGGCGGAAGATCCACGGCTGGTGATCCTGTGGGCGGGGTATCCCTTCTCGCGCGATTACAACAAGc                                                                                      <  1:2455329/150‑1 (MQ=255)
                                                                   aaTGAAACCTTTGCCCCGCAGCATCCCGATCAATATCTCTCCTGGAAAGCCACCTCGGAACAGTCAGAGCGTGTTGACGCCCTGGCGGAAGATCCACGGCTGGTGATCCTGTGGGCGGGGTATCCCTTCTCGCGCGATTACAACAAGCcg                                                                                    >  1:565860/1‑150 (MQ=255)
                                                                     tGAAACCTTTGCCCCGCAGCATCCCGATCAATATCTCTCCTGGAAAGCCACCTCGGAACAGTCAGAGCGTGTTGACGCCCTGGCGGAAGATCCACGGCTGGTGATCCTGTGGGCGGGGTATCCCTTCTCGCGCGATTACAACACGCcgcg                                                                                  >  2:2945460/1‑150 (MQ=255)
                                                                     tGAAACCTTTGCCCCGCAGCATCCCGATCAATATCTCTCCTGGAAAGCCACCTCGGAACAGTCAGAGCGTGTTGACGCCCTGGCGGAAGATCCACGGCTGGTGATCCTGTGGGCGGGGTATCCCTTCTCGCGCGATTACAACAAGCcgcg                                                                                  <  1:630528/150‑1 (MQ=255)
                                                                        aaCCTTTGCCCCGCAGCATCCCGATCAATATCTCTCCTGGAAAGCCACCTCGGAACAGTCAGAGCGTGTTGACGCCCAGGCGGAAGATCCACGGCTGGTGATCCTGTGGGCGGGGTATCCCTTCTCGCGCGATTACAACAAGCCGCGTgg                                                                               >  1:2093869/1‑150 (MQ=255)
                                                                         aCCTTTGCCCCGCAGCATCCCGATCAATATCTCTCCTGGAAAGCCACCTCGGAACAGTCAGAGCGTGTTGACGCCCTGGCGGAAGATCCACGGCTGGTGATCCTGTGGGCGGGGTATCCCTTCTCGCGCGATTACAACAAGCCGCGTGGa                                                                              >  2:4521516/1‑150 (MQ=255)
                                                                              tGCCCCGCAGCATCCCGATCAATATCTCTCCTGGAAAGCCACCTCGGAACAGTCAGAGCGTGTTGACGCCCTGGCGGAAGATCCACGGCTGGTGATCCTGTGGGCGGGGTATCCCTTCTCGCGCGATTACAACAAGCCGCGTGGACATGc                                                                         >  2:2848439/1‑150 (MQ=255)
                                                                                  ccGCAGCATCCCGATCAATATCTCTCCTGGAAAGCCACCTCGGAACAGTCAGAGCGTGTTGACGCCCTGGCGGAAGATCCACGGCTGGTGATCCTGTGGGCGGGGTATCCCTTCTCGCGCGATTACAACAAGCCGCGTGGACATGCtttt                                                                     <  2:3580535/150‑1 (MQ=255)
                                                                                    gcagcaTCCCGATCAATATCTCTCCTGGAAAGCCACCTCGGAACAGTCAGAGCGTGTTGACGCCCTGGCGGAAGATCCACGGCTGGTGATCCTGTGGGCGGGGTATCCCTTCTCGCGCGATTACAACAAGCCGCGTGGACATGCTTTTGc                                                                   >  1:3832299/1‑150 (MQ=255)
                                                                                       gcaTCCCGATCAATATCTCTCCTGGAAAGCCACCTCGGAACAGTCAGAGCGTGTTGACGCCCTGGCGGAAGATCCACGGCTGGTGATCCTGTGGGCGGGGTATCCCTTCTCGCGCGATTACAACAAGCCGCGTGGACATGCTTTTGCtgt                                                                >  1:4660744/1‑150 (MQ=255)
                                                                                        caTCCCGATCAATATCTCTCCTGGAAAGCCACCTCGGAACAGTCAGAGCGTGTTGACGCCCTGGCGGAAGATCCACGGCTGGTGATCCTGTGGGCGGGGTATCCCTTCTCGCGCGATTACAACAAGCCGCGTGGACATGCTTTTGCtgtg                                                               <  1:3193489/150‑1 (MQ=255)
                                                                                        caTCCCGATCAATATCTCTCCTGGAAAGCCACCTCGGAACAGTCAGAGCGTGTTGACGCCCTGGCGGAAGATCCACGGCTGGTGATCCTGTGGGCGGGGTATCCCTTCTCGCGCGATTACAACAAGCCGCGTGGACATGCTTTTGCtgtg                                                               <  1:1155833/150‑1 (MQ=255)
                                                                                         aTCCCGATCAATATCTCTCCTGGAAAGCCACCTCGGAACAGTCAGAGCGTGTTGACGCCCTGGCGGAAGATCCACGGCTGGTGATCCTGTGGGCGGGGTATCCCTTCTCGCGCGATTACAACAAGCCGCGTGGACATGCTTTTGCTGTGa                                                              >  1:2901836/1‑150 (MQ=255)
                                                                                         aTCCCGATCAATATCTCTCCTGGAAAGCCACCTCGGAACAGTCAGAGCGTGTTGACGCCCTGGCGGAAGATCCACGGCTGGTGATCCTGTGGGCGGGGTATCCCTTCTCGCGCGATTACAACAAGCCGCGTGGACATGCTTTTGCTGTGa                                                              <  2:3920440/150‑1 (MQ=255)
                                                                                           cccGATCAATATCTCTCCTGGAAAGCCACCTCGGAACAGTCAGAGCGTGTTGACGCCCTGGCGGAAGATCCACGGCTGGTGATCCTGTGGGCGGGGTATCCCTTCTCGCGCGATTACAACAAGCCGCGTGGACATGCTTTTGCTGTGAcc                                                            <  2:4780184/150‑1 (MQ=255)
                                                                                           cccGATCAATATCTCTCCTGGAAAGCCACCTCGGAACAGTCAGAGCGTGTTGACGCCCTGGCGGAAGATCCACGGCTGGTGATCCTGTGGGCGGGGTATCCCTTCTCGCGCGATTACAACAAGCCGCGTGGACATGCTTTTGCTGTGAcc                                                            >  2:2256050/1‑150 (MQ=255)
                                                                                            ccGATCAATATCTCTCCTGGAAAGCCACCTCGGAACAGTCAGAGCGTGTTGACGCCCTGGCGGAAGATCCACGGCTGGTGATCCTGTGGGCGGGGTATCCCTTCTCGCGCGATTACAACAAGCCGCGTGGACATGCTTTTGCTGTGACCg                                                           <  2:4373827/150‑1 (MQ=255)
                                                                                              gATCAATATCTCTCCTGGAAAGCCACCTCGGAACAGTCAGAGCGTGTTGACGCCCTGGCGGAAGATCCACGGCTGGTGATCCTGTGGGCGGGGTATCCCTTCTCGCGCGATTACAACAAGCCGCGTGGACATGCTTTTGCTGTGACCGAt                                                         <  1:4298647/150‑1 (MQ=255)
                                                                                              gATCAATATCTCTCCTGGAAAGCCACCTCGGAACAGTCAGAGCGTGTTGACGCCCTGGCGGAAGATCCACGGCTGGTGATCCTGTGGGCGGGGTATCCCTTCTCGCGCGATTACAACAAGCCGCGTGGACATGCTTTTGCTGTGACCGAt                                                         <  2:3751142/150‑1 (MQ=255)
                                                                                               aTCAATATCTCTCCTGGAAAGCCACCTCGGAACAGTCAGAGCGTGTTGACGCCCTGGCGGAAGATCCACGGCTGGTGATCCTGTGGGCGGGGTATCCCTTCTCGCGCGATTACAACAAGCCGCGTGGACATGCTTTTGCTGTGACCGAtg                                                        >  1:3697757/1‑150 (MQ=255)
                                                                                                tCAATATCTCTCCTGGAAAGCCACCTCGGAACAGTCAGAGCGTGTTGACGCCCTGGCGGAAGATCCACGGCTGGTGATCCTGTGGGCGGGGTATCCCTTCTCGCGCGATTACAACAAGCCGCGTGGACATGCTTTTGCTGTGACCGAtgt                                                       <  1:2364524/150‑1 (MQ=255)
                                                                                                tCAATATCTCTCCTGGAAAGCCACCTCGGAACAGTCAGAGCGTGTTGACGCCCAGGCGGAAGATCCACGGCTGGTGATCCTGTGGGCGGGGTATCCCTTCTCGCGCGATTACAACAAGCCGCGTGGACATGCTTTTGCTGTGACCGAtgt                                                       >  2:1527216/1‑150 (MQ=255)
                                                                                                 cAATATCTCTCCTGGAAAGCCACCTCGGAACAGTCAGAGCGTGTTGACGCCCTGGCGGAAGATCCACGGCTGGTGATCCTGTGGGCGGGGTATCCCTTCTCGCGCGATTACAACAAGCCGCGTGGACATGCTTTTGCTGTGACCGAtgtg                                                      <  1:1744781/150‑1 (MQ=255)
                                                                                                 cAATATCTCTCCTGGAAAGCCACCTCGGAACAGTCAGAGCGTGTTGACGCCCTGGCGGAAGATCCACGGCTGGTGATCCTGTGGGCGGGGTATCCCTTCTCGCGCGATTACAACAAGCCGCGTGGACATGCTTTTGCTGTGACCGAtgtg                                                      >  2:3861763/1‑150 (MQ=255)
                                                                                                  aaTATCTCTCCTGGAAAGCCACCTCGGAACAGTCAGAGCGTGTTGACGCCCTGGCGGAAGATCCACGGCTGGTGATCCTGTGGGCGGGGTATCCCTTCTCGCGCGATTACAACAAGCCGCGTGGACATGCTTTTGCTGTGACCGATGTGc                                                     <  1:1314246/150‑1 (MQ=255)
                                                                                                    tatCTCTCCTGGAAAGCCACCTCGGAACAGTCAGAGCGTGTTGACGCCCTGGCGGAAGATCCACGGCTGGTGATCCTGTGGGCGGGGTATCCCTTCTCGCGCGATTACAACAAGCCGCGTGGACATGCTTTTGCTGTGACCGATGTGCGt                                                   >  2:2207526/1‑150 (MQ=255)
                                                                                                     atCTCTCCTGGAAAGCCACCTCGGAACAGTCAGAGCGTGTTGACGCCCTGGCGGAAGATCCACGGCTGGTGATCCTGTGGGCGGGGTATCCCTTCTCGCGCGATTACAACAAGCCGCGTGGACATGCTTTTGCTGTGACCGATGTGCGTg                                                  <  2:1751779/150‑1 (MQ=255)
                                                                                                      tctctcCTGGAAAGCCACCTCGGAACAGTCAGAGCGTGTTGACGCCCTGGCGGAAGATCCACGGCTGGTGATCCTGTGGGCGGGGTATCCCTTCTCGCGCGATTACAACAAGCCGCGTGGACATGCTTTTGCTGTGACCGATGTGCGTGa                                                 >  2:1701303/1‑150 (MQ=255)
                                                                                                      tctctcCTGGAAAGCCACCTCGGAACAGTCAGAGCGTGTTGACGCCCTGGCGGAAGATCCACGGCTGGTGATCCTGTGGGCGGGGTATCCCTTCTCGCGCGATTACAACAAGCCGCGTGGACATGCTTTTGCTGTGACCGATGTGCGTGa                                                 >  2:3244351/1‑150 (MQ=255)
                                                                                                      tctctcCTGGAAAGCCACCTCGGAACAGTCAGAGCGTGTTGACGCCCTGGCGGAAGATCCACGGCTGGTGATCCTGTGGGCGGGGTATCCCTTCTCGCGCGATTACAACAAGCCGCGTGGACATGCTTTTGCTGTGACCGATGTGCGTGa                                                 >  2:752122/1‑150 (MQ=255)
                                                                                                        tctcCTGGAAAGCCACCTCGGAACAGTCAGAGCGTGTTGACGCCCTGGCGGAAGATCCACGGCTGGTGATCCTGTGGGCGGGGTATCCCTTCTCGCGCGATTACAACAAGCCGCGTGGACATGCTTTTGCTGTGACCGATGTGCGTGaaa                                               >  1:733717/1‑150 (MQ=255)
                                                                                                             tGGAAAGCCACCTCGGAACAGTCAGAGCGTGTTGACGCCCTGGCGGAAGATCCACGGCTGGTGATCCTGTGGGCGGGGTATCCCTTCTCGCGCGATTACAACAAGCCGCGTGGACATGCTTTTGCTGTGACCGATGTGCGTGAAACCCTg                                          >  1:4369641/1‑150 (MQ=255)
                                                                                                             tGGAAAGCCACCTCGGAACAGTCAGAGCGTGTTGACGCCCTGGCGGAAGATCCACGGCTGGTGATCCTGTGGGCGGGGTATCCCTTCTCGCGCGAGTACAACACGCCGCGTGGACATGCTTTTGCTGTGACCGATGTGCGTGAAACCCTg                                          <  2:3394421/150‑1 (MQ=255)
                                                                                                                 aaGCCACCTCGGAACAGTCAGAGCGTGTTGACGCCCTGGCGGAAGATCCACGGCTGGTGATCCTGTGGGCGGGGTATCCCTTCTCGCGCGATTACAACAAGCCGCGTGGACATGCTTTTGCTGTGACCGATGTGCGTGAAACCCTGCGTa                                      <  1:4289905/150‑1 (MQ=255)
                                                                                                                    ccACCTCGGAACAGTCAGAGCGTGTTGACGCCCTGGCGGAAGATCCACGGCTGGTGATCCTGTGGGCGGGGTATCCCTTCTCGCGCGATTACAACAAGCCGCGTGGACATGCTTTTGCTGTGACCGATGTGCGTGAAACCCTGCGTACCg                                   <  1:2700103/150‑1 (MQ=255)
                                                                                                                      aCCTCGGAACAGTCAGAGCGTGTTGACGCCCTGGCGGAAGATCCACGGCTGGTGATCCTGTGGGCGGGGTATCCCTTCTCGCGCGATTACAACAAGCCGCGTGGACATGCTTTTGCTGTGACCGATGTGCGTGAAACCCTGCGTACCGGt                                 <  1:1037517/150‑1 (MQ=255)
                                                                                                                      aCCTCGGAACAGTCAGAGCGTGTTGACGCCCTGGCGGAAGATCCACGGCTGGTGATCCTGTGGGCGGGGTATCCCTTCTCGCGCGATTACAACAAGCCGCGTGGACATGCTTTTGCTGTGACCGATGTGCGTGAAACCCTGCGTACCGGt                                 >  2:1825748/1‑150 (MQ=255)
                                                                                                                        cTCGGAACAGTCAGAGCGTGTTGACGCCCTGGCGGAAGATCCACGGCTGGTGATCCTGTGGGCGGGGTATCCCTTCTCGCGCGATTACAACAAGCCGCGTGGACATGCTTTTGCTGTGACCGATGTGCGTGAAACCCTGCGTACCGGTgc                               >  2:2621407/1‑150 (MQ=255)
                                                                                                                          cGGAACAGTCAGAGCGTGTTGACGCCCTGGCGGAAGATCCACGGCTGGTGATCCTGTGGGCGGGGTATCCCTTCTCGCGCGATTACAACAAGCCGCGTGGACATGCTTTTGCTGTGACCGATGTGCGTGAAACCCTGCGTACCGGTgcgc                             >  1:2071106/1‑150 (MQ=255)
                                                                                                                           ggAACAGTCAGAGCGTGTTGACGCCCAGGCGGAAGATCCACGGCTGGTGATCCTGTGGGCGGGGTATCCCTTCTCGCGCGATTACAACAAGCCGCGTGGACATGCTTTTGCTGTGACCGATGTGCGTGAAACCCTGCGTACCGGTGCGcc                            <  1:574663/150‑1 (MQ=255)
                                                                                                                            gAACAGTCAGAGCGTGTTGACGCCCTGGCGGAAGATCCACGGCTGGTGCTCCTGTGGGCGGGGTCTCACTTCTCGCGCGATTACAACAAGCCGCGTGGACATGCTTTTGCTGTGACCGATGTGCGTGAAACCCTGCGTAGCGGTGCGCCg                           >  2:2563568/1‑150 (MQ=255)
                                                                                                                                aGTCAGAGCGTGTTGACGCCCTGGCGGAAGATCCACGGCTGGTGATCCTGTGGGCGGGGTATCCCTTCTCGCGCGATTACAACAAGCCGCGTGGACATGCTTTTGCTGTGACCGATGTGCGTGAAACCCTGCGTACCGGTGCGCCGaaaa                       >  2:3148141/1‑150 (MQ=255)
                                                                                                                                aGTCAGAGCGTGTTGACGCCCTGGCGGAAGATCCACGGCTGGTGATCCTGTGGGCGGGGTATCCCTTCTCGCGCGATTACAACAAGCCGCGTGGACATGCTTTTGCTGTGACCGATGTGCGTGAAACCCTGCGTACCGGTGCGCCGaaaa                       >  1:3713820/1‑150 (MQ=255)
                                                                                                                                aGTCAGAGCGTGTTGACGCCCTGGCGGAAGATCCACGGCTGGTGATCCTGTGGGCGGGGTATCCCTTCTCGCGCGATTACAACAAGCCGCGTGGACATGCTTTTGCTGTGACCGATGTGCGTGAAACCCTGCGTACCGGTGCGCCGaaaa                       <  1:3787599/150‑1 (MQ=255)
                                                                                                                                   cAGAGCGTGTTGACGCCCTGGCGGAAGATCCACGGCTGGTGATCCTGTGGGCGGGGTATCCCTTCTCGCGCGATTACAACAAGCCGCGTGGACATGCTTTTGCTGTGACCGATGTGCGTGAAACCCTGCGTACCGGTGCGCCGACAAACg                    >  2:4164652/1‑150 (MQ=255)
                                                                                                                                    agagCGTGTTGACGCCCTGGCGGAAGATCCACGGCTGGTGATCCTGTGGGCGGGGTATCCCTTCTCGCGCGATTACAACAAGCCGCGTGGACATGCTTTTGCTGTGACCGATGTGCGTGAAACCCTGCGTACCGGTGCGCCGAAAAACGc                   >  2:1225253/1‑150 (MQ=255)
                                                                                                                                      agCGTGTTGACGCCCTGGCGGAAGATCCACGGCTGGTGATCCTGTGGGCGGGGTATCCCTTCTCGCGCGATTACAACAAGCCGCGTGGACATGCTTTTGCTGTGACCGATGTGCGTGAAACCCTGCGTACCGGTGCGCCGAAAAACGCTg                 <  1:3240741/150‑1 (MQ=255)
                                                                                                                                       gCGTGTTGACGCCCTGGCGGAAGATCCACGGCTGGTGATCCTGTGGGCGGGGTATCCCTTCTCGCGCGATTACAACAAGCCGCGTGGACATGCTTTTGCTGTGACCGATGTGCGTGAAACCCTGCGTACCGGTGCGCCGAAAAACGCTGa                <  1:317491/150‑1 (MQ=255)
                                                                                                                                       gCGTGTTGACGCCCTGGCGGAAGATCCACGGCTGGTGATCCTGTGGGCGGGGTATCCCTTCTCGCGCGATTACAACAAGCCGCGTGGACATGCTTTTGCTGTGACCGATGTGCGTGAAACCCTGCGTACCGGTGCGCCGAAAAACGCTGa                <  1:968344/150‑1 (MQ=255)
                                                                                                                                        cGTGTTGACGCCCTGGCGGAAGATCCACGGCTGGTGATCCTGTGGGCGGGGTATCCCTTCTCGCGCGATTACAACAAGCCGCGTGGACATGCTTTTGCTGTGACCGATGTGCGTGAAACCCTGCGTACCGGTGCGCCGAAAAACGCTGaa               <  1:3173077/150‑1 (MQ=255)
                                                                                                                                        cGTGTTGACGCCCTGGCGGAAGATCCACGGCTGGTGATCCTGTGGGCGGGGTATCCCTTCTCGCGCGATTACAACAAGCCGCGTGGACATGCTTTTGCTGTGACCGATGTGCGTGAAACCCTGCGTACCGGTGCGCCGAAAAACGCTGaa               >  1:2531326/1‑150 (MQ=255)
                                                                                                                                         gtgtTGACGCCCTGGCGGAAGATCCACGGCTGGTGATCCTGTGGGCGGGGTATCCCTTCTCGCGCGATTACAACAAGCCGCGTGGACATGCTTTTGCTGTGACCGATGTGCGTGAAACCCTGCGTACCGGTGCGCCGAAAAACGCTGAAg              >  2:1975268/1‑150 (MQ=255)
                                                                                                                                         gtgtTGACGCCCTGGCGGAAGATCCACGGCTGGTGATCCTGTGGGCGGGGTATCCCTTCTCGCGCGATTACAACAAGCCGCGTGGACATGCTTTTGCTGTGACCGATGTGCGTGAAACCCTGCGTACCGGTGCGCCGAAAAACGCTGAAg              <  1:4590382/150‑1 (MQ=255)
                                                                                                                                          tgtTGACGCCCTGGCGGAAGATCCACGGCTGGTGATCCTGTGGGCGGGGTATCCCTTCTCGCGCGATTACAACAAGCCGCGTGGACATGCTTTTGCTGTGACCGATGTGCGTGAAACCCTGCGTACCGGTGCGCCGAAAAACGCTGAAGa             >  2:4200934/1‑150 (MQ=255)
                                                                                                                                             tGACGCCCTGGCGGAAGATCCACGGCTGGTGATCCTGTGGGCGGGGTATCCCTTCTCGCGCGATTACAACAAGCCGCGTGGACATGCTTTTGCTGTGACCGATGTGCGTGAAACCCTGCGTACCGGTGCGCCGAAAAACGCTGAAGATgg          <  2:2781738/150‑1 (MQ=255)
                                                                                                                                                cGCCCTGGCGGAAGATCCACGGCTGGTGATCCTGTGGGCGGGGTATCCCTTCTCGCGCGATTACAACAAGCCGCGTGGACATGCTTTTGCTGTGACCGATGTGCGTGAAACCCTGCGTACCGGTGCGCCGAAAAACGCTGAAGATGGTcc       >  1:1064107/1‑150 (MQ=255)
                                                                                                                                                 gCCCTGGCGGAAGATCCACGGCTGGTGATCCTGTGGGCGGGGTATCCCTTCTCGCGCGATTACAACAAGCCGCGTGGACATGCTTTTGCTGTGACCGATGTGCGTGAAACCCTGCGTACCGGTGCGCCGAAAAACGCTGAAGATGGTCCg      <  1:991145/150‑1 (MQ=255)
                                                                                                                                                 gCCCTGGCGGAAGATCCACGGCTGGTGATCCTGTGGGCGGGGTATCCCTTCTCGCGCGATTACAACAAGCCGCGTGGACATGCTTTTGCTGTGACCGATGTGCGTGAAACCCTGCGTACCGGTGCGCCGAAAAACGCTGAAGATGGTCCg      <  1:2464595/150‑1 (MQ=255)
                                                                                                                                                     tGGCGGAAGATCCACGGCTGGTGATCCTGTGGGCGGGGTATCCCTTCTCGCGCGATTACAACAAGCCGCGTGGACATGCTTTTGCTGTGACCGATGTGCGTGAAACCCTGCGTACCGGTGCGCCGAAAAACGCTGAAGATGGTCCGCTAc  <  1:3307158/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CTTTTTCACTTCTGTTCACGCTGAACAAACGGCTGCTCCCGCAAAACCTGTAACTGTGGAAGCGAAGAATGAAACCTTTGCCCCGCAGCATCCCGATCAATATCTCTCCTGGAAAGCCACCTCGGAACAGTCAGAGCGTGTTGACGCCCTGGCGGAAGATCCACGGCTGGTGATCCTGTGGGCGGGGTATCCCTTCTCGCGCGATTACAACAAGCCGCGTGGACATGCTTTTGCTGTGACCGATGTGCGTGAAACCCTGCGTACCGGTGCGCCGAAAAACGCTGAAGATGGTCCGCTAC  >  CP009273/4277748‑4278046

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: