Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,281,581 G→T 11.1% D48Y (GAT→TAT)  nrfE → heme lyase (NrfEFG) for insertion of heme into c552, subunit NrfE

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,281,5810GT11.1% 947.4 / 58.8 274D48Y (GAT→TAT) nrfEheme lyase (NrfEFG) for insertion of heme into c552, subunit NrfE
Reads supporting (aligned to +/- strand):  ref base G (129/114);  new base T (14/17);  total (143/131)
Fisher's exact test for biased strand distribution p-value = 4.49e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 281 total aligned reads displayed.

GGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCGATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCCTGCTGTCGTGGGAACTCAAGCTGGCAGCGGTGTGGGGCGGCCATGAAGGTTCGCTGCTGCTTTGGGTGCTGCTGCTTTCCGCCTGGAGCGCGCTGTTTGCCTGGCATTA  >  CP009273/4281433‑4281729
                                                                                                                                                    |                                                                                                                                                    
ggTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCg                                                                                                                                                      >  1:3483635/1‑149 (MQ=255)
 gTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCGa                                                                                                                                                     >  2:3151369/1‑149 (MQ=255)
   cAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCGAtt                                                                                                                                                   <  1:3779144/149‑1 (MQ=255)
    aaCGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCTAttt                                                                                                                                                  >  1:458342/1‑149 (MQ=255)
      cGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCGATTtct                                                                                                                                                >  1:2353954/1‑149 (MQ=255)
       gtgtTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCGATTtctc                                                                                                                                               <  2:1616511/149‑1 (MQ=255)
       gtgtTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCGATTtctc                                                                                                                                               <  1:4151369/149‑1 (MQ=255)
        tgtTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCGATTTCTCg                                                                                                                                              <  1:4786562/149‑1 (MQ=255)
        tgtTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCGATTTCTCg                                                                                                                                              >  1:2838538/1‑149 (MQ=255)
            gACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCGATTTCTCGGTCa                                                                                                                                          <  1:106965/149‑1 (MQ=255)
            gACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCGATTTCTCGGTCa                                                                                                                                          >  2:1953469/1‑149 (MQ=255)
              ccccGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCGATTTCTCGGTCAtt                                                                                                                                        <  1:1263698/149‑1 (MQ=255)
              ccccGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCGATTTCTCGGTCAtt                                                                                                                                        >  1:1970589/1‑149 (MQ=255)
                ccGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCGATTTCTCGGTCATTTa                                                                                                                                      >  1:410601/1‑149 (MQ=255)
                 cGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCGATTTCTCGGTCATTTAt                                                                                                                                     <  2:15528/149‑1 (MQ=255)
                   ttGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCGATTTCTCGGTCATTTATGt                                                                                                                                   >  2:4334309/1‑149 (MQ=255)
                    tGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCGATTTCTCGGTCATTTATGTc                                                                                                                                  <  1:1528997/149‑1 (MQ=255)
                    tGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTGGGCGTACTGACGTATTGTTTTCTCATCAGCTATTTCTCGGTCATTTATGTc                                                                                                                                  <  2:797764/149‑1 (MQ=255)
                           cTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCGATTTCTCGGTCATTTATGTCGCGCAAc                                                                                                                           <  2:2240228/149‑1 (MQ=255)
                                 gggAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCGATTTCTCGGTCATTTATGTCGCGCAACATAGCt                                                                                                                     >  1:2103356/1‑149 (MQ=255)
                                 gggAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCGATTTCTCGGTCATTTATGTCGCGCAACATAGCt                                                                                                                     >  1:4502108/1‑149 (MQ=255)
                                  ggAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCGATTTCTCGGTCATTTATGTCGCGCAACATAGCta                                                                                                                    >  2:4752836/1‑149 (MQ=255)
                                    aGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCGATTTCTCGGTCATTTATGTCGCGCAACATAGCtata                                                                                                                  <  2:617443/149‑1 (MQ=255)
                                          gTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCGATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCctgc                                                                                                            >  1:2494436/1‑149 (MQ=255)
                                              cgcTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCTATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCCTGCTGTc                                                                                                        <  1:3966917/149‑1 (MQ=255)
                                                   gCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCGATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCCTGCTGTCGTggg                                                                                                   <  2:2516315/149‑1 (MQ=255)
                                                       gCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCGATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCCTGCTGTCGTGGGACCt                                                                                               >  1:1661196/1‑149 (MQ=255)
                                                        ccATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCGATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCCTGCTGTCGTGGGAACTc                                                                                              <  1:4013703/149‑1 (MQ=255)
                                                         cATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCGATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCCTGCTGTCGTGGGAACTCa                                                                                             <  1:1342983/149‑1 (MQ=255)
                                                            gctgCGTCTCTCTTGCATCGGCTTTCTGGCGGAGTTCGCGCTCCTGCTGCTCGCATTTGGCGTACTGACGTATTGTTTTCTCATCAGCGATTTCTCGGTCATGTTTGTCGCGGAACATAGCTATAGCCTGCTGTCGTGGGAACTCAAGc                                                                                          <  1:3621889/147‑1 (MQ=255)
                                                            gatgCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCGATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCCTGCTGTCGTGGGAACTCAAGc                                                                                          >  2:62759/1‑149 (MQ=255)
                                                              tgCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCGATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCCTGCTGTCGTGGGAACTCAAGCTg                                                                                        >  1:839177/1‑149 (MQ=255)
                                                               gCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCGATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCCTGCTGTCGTGGGAACTCAAGCTgg                                                                                       >  1:1920263/1‑149 (MQ=255)
                                                               gCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCGATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCCTGCTGTCGTGGGAACTCAAGCTgg                                                                                       >  1:4310524/1‑149 (MQ=255)
                                                                cGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCGATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCCTGCTGTCGTGGGAACTCAAGCTGGc                                                                                      <  1:2541617/149‑1 (MQ=255)
                                                                 gACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCGATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCCTGCTGTCGTGGGAACTCAAGCTGGCa                                                                                     <  1:3396884/149‑1 (MQ=255)
                                                                   cTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTGGGCGTACTGACGTATTGTTTTCTCATCAGCGATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCCTGCTGTCGGGGGAACTCACGCTGGCacc                                                                                   >  1:4627830/1‑147 (MQ=255)
                                                                    tCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCGATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCCTGCTGTCGTGGGAACTCAAGCTGGCAGCg                                                                                  >  2:3820654/1‑149 (MQ=255)
                                                                     cACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCGATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCCTGCTGTCGTGGGAACTCAAGCTGGCAGCgg                                                                                 <  2:410601/149‑1 (MQ=255)
                                                                      aCTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCGATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCCTGCTGTCGTGGGAACTCAAGCTGGCAGCGgt                                                                                <  2:2697398/149‑1 (MQ=255)
                                                                        ttGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCGATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCCTGCTGTCGTGGGAACTCCAGCTGGCAGCGgtgt                                                                              <  1:841990/149‑1 (MQ=255)
                                                                           cATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCGATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCCTGCTGTCGTGGGAACTCAAGCTGGCAGCGGTGTggg                                                                           <  1:906726/149‑1 (MQ=255)
                                                                            aTCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCGATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCCTGCTGTCGTGGGAACTCAAGCTGGCAGCGGTGTgggg                                                                          >  2:2927679/1‑149 (MQ=255)
                                                                             tCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCTATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCCTGCTGTCGTGGGAACTCAAGCTGGCAGCGGTGTGGggc                                                                         >  1:523721/1‑149 (MQ=255)
                                                                             tCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCGATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCCTGCTGTCGTGGGAACTCAAGCTGGCAGCGGTGTGGggc                                                                         >  1:3394529/1‑149 (MQ=255)
                                                                             tCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCGATTTCTCGGCCATTTATGTCGCGCAACATAGCTATAGCCTGCTGTCGTGGGAACTCAAGCTGGCAGCGGTGTGGggc                                                                         >  1:2461747/1‑149 (MQ=255)
                                                                              cGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCGATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCCTGCTGTCGTGGGAACTCAAGCTGGCAGCGGTGTGGggcg                                                                        >  2:878422/1‑149 (MQ=255)
                                                                               ggCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCGATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCCTGCTGTCGTGGGAACTCAAGCTGGCAGCGGTGTGGggcgg                                                                       >  1:2018801/1‑149 (MQ=255)
                                                                                 cATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCTATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCCTGCTGTCGTGGGAACTCAAGCTGGCAGCGGTGTGGGGCGGcc                                                                     <  2:3054686/149‑1 (MQ=255)
                                                                                 cATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCGATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCCTGCTGTCGTGGGAACTCAAGCTGGCAGCGGTGTGGGGCGGcc                                                                     <  2:2374352/149‑1 (MQ=255)
                                                                                   ttCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCTATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCCTGCTGTCGTGGGAACTCAAGCTGGCAGCGGTGTGGGGCGGCCAt                                                                   >  1:38548/1‑149 (MQ=255)
                                                                                   ttCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCGATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCCTGCTGTCGTGGGAACTCAAGCTGGCAGCGGTGTGGGGCGGCCAt                                                                   >  2:4167574/1‑149 (MQ=255)
                                                                                   ttCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCGATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCCTGCTGTCGTGGGAACTCAAGCTGGCAGCGGTGTGGGGCGGCCAt                                                                   >  2:2074889/1‑149 (MQ=255)
                                                                                    tCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCGATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCCTGCTGTCGTGGGAACTCAAGCTGGCAGCGGTGTGGGGCGGCCATg                                                                  <  1:3016207/149‑1 (MQ=255)
                                                                                      tGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCGATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCCTGCTGTCGTGGGAACTCAAGCTGGCAGCGGTGTGGGGCGGCCATGaa                                                                <  2:2213935/149‑1 (MQ=255)
                                                                                          gcAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCGATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCCTGCTGTCGTGGGAACTCAAGCTGGCAGCGGTGTGGGGCGGCCATGAAGGtt                                                            <  2:1049072/149‑1 (MQ=255)
                                                                                          gcAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCGATTTCTCGGTCATTTATGTCGCGCAACAGAGCTATAGCCTGCTGTCGTGGGAACTCAAGCTGGCAGCGGTGTGGGGCGGCCATGAAGGtt                                                            <  1:938241/149‑1 (MQ=255)
                                                                                            aGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCGATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCCTGCTGTCGTGGGAACTCAAGCTGGCAGCGGTGTGGGGCGGCCATGAAGGTTCg                                                          >  2:1181768/1‑149 (MQ=255)
                                                                                             gTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCGATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCCTGCTGTCGTGGGAACTCAAGCTGGCAGCGGTGTGGGGCGGCCATGAAGGTTCgc                                                         >  2:3412190/1‑149 (MQ=255)
                                                                                             gTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCGATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCCTGCTGTCGTGGGAACTCAAGCTGGCAGCGGTGTGGGGCGGCCATGAAGGTTCgc                                                         <  1:1569472/149‑1 (MQ=255)
                                                                                              gtCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCGATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCCTGCTGTCGTGGGAACTCAAGCTGGCAGCGGTGTGGGGCGGCCATGAAGGTTCgct                                                        <  1:4187040/148‑1 (MQ=255)
                                                                                                       ctgctgCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCGATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCCTGCTGTCGTGGGAACTCAAGCTGGCAGCGGTGTGGGGCGGCCATGAAGGTTCGCTGCTGCTTTg                                               >  2:1771095/1‑149 (MQ=255)
                                                                                                       ctgctgCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCGATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCCTGCTGTCGTGGGAACTCAAGCTGGCAGCGGTGTGGGGCGGCCATGAAGGTTCGCTGCTGCTTTg                                               <  2:2178333/149‑1 (MQ=255)
                                                                                                        tgctgcTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCGATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCCTGCTGTCGTGGGAACTCAAGCTGGCAGCGGTGTGGGGCGGCCATGAAGGTTCGCTGCTGCTTTgg                                              >  2:537144/1‑149 (MQ=255)
                                                                                                             ctCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCGATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCCTGCTGTCGTGGGAACTCAAGCTGGCAGCGGTGTGGGGCGGCCATGAAGGTTCGCTGCTGCTTTGGGtgct                                         >  1:713284/1‑149 (MQ=255)
                                                                                                                gCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCGATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCCTGCTGTCGTGGGAACTCAAGCTGGCAGCGGTGTGGGGCGGCCATGAAGGTTCGCTGCTGCTTTGGGtgctgct                                      >  2:3601992/1‑149 (MQ=255)
                                                                                                                   tttGGCGTACTGACGTATTGTTTTCTCATCAGCGATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCCTGCTGTCGTGGGAACTCAAGCTGGCAGCGGTGTGGGGCGGCCATGAAGGTTCGCTGCTGCTTTGGGTgctgctgct                                   <  2:955018/149‑1 (MQ=255)
                                                                                                                    ttGGCGTACTGACGTATTGTTTTCTCATCAGCGATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCCTGCTGTCGTGGGAACTCAAGCTGGCAGCGGTGTGGGGCGGCCATGAAGGTTCGCTGCTGCTTTGGGTGCTGCTGCtt                                  >  2:4079097/1‑149 (MQ=255)
                                                                                                                    ttGGCGTACTGACGTATTGTTTTCTCATCAGCGATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCCTGCTGTCGTGGGAACTCAAGCTGGCAGCGGTGTGGGGCGGCCATGAAGGTTCGCTGCTGCTTTGGGTGCTGCTGCtt                                  <  2:826458/149‑1 (MQ=255)
                                                                                                                      ggCGTACTGACGTATTGTTTTCTCATCAGCGATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCCTGCTGTCGTGGGAACTCAAGCTGGCAGCGGTGTGGGGCGGCCATGAAGGTTCGCTGCTGCTTTGGGTGCTGCTGCTTTc                                <  1:3409339/149‑1 (MQ=255)
                                                                                                                         gTACTGACGTATTGTTTTCTCATCAGCGATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCCTGCTGTCGTGGGAACTCAAGCTGGCAGCGGTGTGGGGCGGCCATGAAGGTTCGCTGCTGCTTTGGGTGCTGCTGCTTTCCGc                             >  1:2781194/1‑149 (MQ=255)
                                                                                                                          tACTGACGTATTGTTTTCTCATCAGCGATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCCTGCTGTCGTGGGAACTCAAGCTGGCAGCGGTGTGGGGCGGCCATGAAGGTTCGCTGCTGCTTTGGGTGCTGCTGCTTTCCGcc                            >  2:4001014/1‑149 (MQ=255)
                                                                                                                           aCTGACGTATTGTTTTCTCATCAGCGATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCCTGCTGTCGTGGGAACTCAAGCTGGCAGCGGTGTGGGGCGGCCATGAAGGTTCGCTGCTGCTTTGGGTGCTGCTGCTTTCCGCCt                           >  1:1138761/1‑149 (MQ=255)
                                                                                                                            cTGACGTATTGTTTTCTCATCAGCTATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCCTGCTGTCGTGGGAACTCAAGCTGGCAGCGGTGTGGGGCGGCCATGAAGGTTCGCTGCTGCTTTGGGTGCTGCTGCTTTCCGCCTg                          >  1:4675763/1‑149 (MQ=255)
                                                                                                                            cTGACGTATTGTTTTCTCATCAGCGATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCCTGCTGTCGTGGGAACTCAAGCTGGCAGCGGTGTGGGGCGGCCATGAAGGTTCGCTGCTGCTTTGGGTGCTGCTGCTTTCCGCCTg                          >  2:1519961/1‑149 (MQ=255)
                                                                                                                             tGACGTATTGTTTTCTCATCAGCGATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCCTGCTGTCGTGGGAACTCAAGCTGGCAGCGGTGTGGGGCGGCCATGAAGGTTCGCTGCTGCTTTGGGTGCTGCTGCTTTCCGCCTgg                         <  1:1790073/149‑1 (MQ=255)
                                                                                                                               aCGTATTGTTTTCTCATCAGCGATGTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCCTGCTGTCGTGGGAACTCAAGCTGGCAGCGGTGTGGGGCGGCCATGAAGGTTCGCTGCTGCTTTGGGTGCTGCTGCTTTCCGCCTGGAg                       >  2:555972/1‑149 (MQ=255)
                                                                                                                                  tATTGTTTTCTCATCAGCGATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCCTGCTGTCGTGGGAACTCAAGCTGGCAGCGGTGTGGGGCGGCCATGAAGGTTCGCTGCTGCTTTGGGTGCTGCTGCTTTCCGCCTGGAgcgc                    >  1:1038468/1‑149 (MQ=255)
                                                                                                                                  tATTGTTTTCTCATCAGCGATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCCTGCTGTCGTGGGAACTCAAGCTGGCAGCGGTGTGGGGCGGCCATGAAGGTTCGCTGCTGCTTTGGGTGCTGCTGCTTTCCGCCTGGAgcgc                    >  2:1176999/1‑149 (MQ=255)
                                                                                                                                  tATTGTTTTCTCATCAGCGATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCCTGCTGTCGTGGGAACTCAAGCTGGCAGCGGTGTGGGGCGGCCATGAAGGTTCGCTGCTGCTTTGGGTGCTGCTGCTTTCCGCCTGGAgcgc                    >  1:3110558/1‑149 (MQ=255)
                                                                                                                                      gTTTTCTCATCAGCGATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCCTGCTGTCGTGGGAACTCAAGCTGGCAGCGGTGTGGGGCGGCCATGAAGGTTCGCTGCTGCTTTGGGTGCTGCTGCTTTCCGCCTGGAGCGCGCTg                >  2:213913/1‑149 (MQ=255)
                                                                                                                                       ttttCTCATCAGCGATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCCTGCTGTCGTGGGAACTCAAGCTGGCAGCGGTGTGGGGCGGCCATGAAGGTTCGCTGCTGCTTTGGGTGCTGCTGCTTTCCGCCTGGAGCGCGCTGt               >  2:86268/1‑149 (MQ=255)
                                                                                                                                        tttCTCATCAGCGATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCCTGCTGTCGTGGGAACTCAAGCTGGCAGCGGTGTGGGGCGGCCATGAAGGTTCGCTGCTGCTTTGGGTGCTGCTGCTTTCCGCCTGGAGCGCGCTGtt              <  2:4418145/149‑1 (MQ=255)
                                                                                                                                         ttCTCATCAGCGATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCCTGCTGTCGTGGGAACTCAAGCTGGCAGCGGTGTGGGGCGGCCATGAAGGTTCGCTGCTGCTTTGGGTGCTGCTGCTTTCCGCCTGGAGCGCGCTGttt             >  2:4288270/1‑149 (MQ=255)
                                                                                                                                         ttCTCATCAGCGATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCCTGCTGTCGTGGGAACTCAAGCTGGCAGCGGTGTGGGGCGGCCATGAAGGTTCGCTGCTGCTTTGGGTGCTGCTGCTTTCCGCCTGGAGCGCGCTGttt             <  2:431083/149‑1 (MQ=255)
                                                                                                                                           ctcATCAGCGATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCCTGCTGTCGTGGGAACTCAAGCTGGCAGCGGTGTGGGGCGGCCATGAAGGTTCGCTGCTGCTTTGGGTGCTGCTGCTTTCCGCCTGGAGCGCGCTGTTTGc           <  2:1729596/149‑1 (MQ=255)
                                                                                                                                           ctcATCAGCGATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCCTGCTGTCGTGGGAACTCAAGCTGGCAGCGGTGTGGGGCGGCCATGAAGGTTCGCTGCTGCTTTGGGTGCTGCTGCTTTCCGCCTGGAGCGCGCTGTTTGc           <  1:3068553/149‑1 (MQ=255)
                                                                                                                                            tcatcaGCGATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCCTGCTGTCGTGGGAACTCAAGCTGGCAGCGGTGTGGGGCGGCCATGAAGGTTCGCTGCTGCTTTGGGTGCTGCTGCTTTCCGCCTGGAGCGCGCTGTTTGcc          <  2:4592253/149‑1 (MQ=255)
                                                                                                                                            tcatcaGCGATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCCTGCTGTCGTGGGAACTCAAGCTGGCAGCGGTGTGGGGCGGCCATGAAGGTTCGCTGCTGCTTTGGGTGCTGCTGCTTTCCGCCTGGAGCGCGCTGTTTGcc          <  2:696188/149‑1 (MQ=255)
                                                                                                                                             cctcaGCGATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCCTGCTGTCGTGGGAACTCAAGCTGGCAGCGGTGTGGGGCGGCCATGAAGGTTCGCTGCTGCTTTGGGTGCTGCTGCTTTCCGCCTGGAGCGCGCTGTTTGCCt         <  2:3665309/147‑1 (MQ=255)
                                                                                                                                              atcaGCGATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCCTGCTGTCGTGGGAACTCAAGCTGGCAGCGGTGTGGGGCGGCCATGAAGGTTCGCTGCTGCTTTGGGTGCTGCTGCTTTCCGCCTGGAGCGCGCTGTTTGCCTg        <  1:1760956/149‑1 (MQ=255)
                                                                                                                                              atcaGCGATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCCTGCTGTCGTGGGAACTCAAGCTGGCAGCGGTGTGGGGCGGCCATGAAGGTTCGCTGCTGCTTTGGGTGCTGCTGCTTTCCGCCTGGAGCGCGCTGTTTGCCTg        <  2:2771955/149‑1 (MQ=255)
                                                                                                                                               tcaGCTATTTCTCGGTCATTTAAGTCGCGCAACATAGCTATAGCCTGCTGTCGTGGGAACTCAAGCTGGCAGCGGTGTGGGGCGGCCATGAAGGTTCGCTGCTGCTTTGGGTGCTGCTGCTTTCCGCCTGGAGCGCGCTGTTTGCCTgg       >  2:2419885/1‑149 (MQ=255)
                                                                                                                                               tcaGCGATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCCTGCTGTCGTGGGAACTCAAGCTGGCAGCGGTGTGGGGCGGCCATGAAGGTTCGCTGCTGCTTTGGGTGCTGCTGCTTTCCGCCTGGAGCGCGCTGTTTGCCTgg       <  2:804156/149‑1 (MQ=255)
                                                                                                                                                caGCGATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCCTGCTGTCGTGGGAACTCAAGCTGGCAGCGGTGTGGGGCGGCCATGAAGGTTCGCTGCTGCTTTGGGTGCTGCTGCTTTCCGCCTGGAGCGCGCTGTTTGCCTGGc      >  1:3263045/1‑149 (MQ=255)
                                                                                                                                                caGCGATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCCTGCTGTCGTGGGAACTCAAGCTGGCAGCGGTGTGGGGCGGCCATGAAGGTTCGCTGCTGCTTTGGGTGCTGCTGCTTTCCGCCTGGAGCGCGCTGTTTGCCTGGc      >  1:4073989/1‑149 (MQ=255)
                                                                                                                                                  gCGATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCCTGCTGTCGTGGGAACTCAAGCTGGCAGCGGTGTGGGGCGGCCATGAAGGTTCGCTGCTGCTTTGGGTGCTGCTGCTTTCCGCCTGGAGCGCGCTGTTTGCCTGGCAt    >  2:3198852/1‑149 (MQ=255)
                                                                                                                                                  gCGATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCCTGCTGTCGTGGGAACTCAAGCTGGCAGCGGTGTGGGGCGGCCATGAAGGTTCGCTGCTGCTTTGGGTGCTGCTGCTTTCCGCCTGGAGCGCGCTGTTTGCCTGGCAt    <  1:3444188/149‑1 (MQ=255)
                                                                                                                                                  gCGATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCCTGCTGTCGTGGGAACTCAAGCTGGCAGCGGTGTGGGGCGGCCATGAAGGTTCGCTGCTGCTTTGGGTGCTGCTGCTTTCCGCCTGGAGCGCGCTGTTTGCCTGGCAt    >  2:1290639/1‑149 (MQ=255)
                                                                                                                                                    gATGTCTCGGTCAGTTATGTCGCGCAACATAGCTATAGCCTGCTGTCGTGGGAACTCAAGCTGGCAGCGGTGTGGGGCGGCCATGAAGGTTCGCTGCTGCTTTGGGTGCTGCTGCTTTCCGCCTGGAGCGCGCTGTTTGCCTGGCAttc  <  1:2961841/149‑2 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
GGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCGATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCCTGCTGTCGTGGGAACTCAAGCTGGCAGCGGTGTGGGGCGGCCATGAAGGTTCGCTGCTGCTTTGGGTGCTGCTGCTTTCCGCCTGGAGCGCGCTGTTTGCCTGGCATTA  >  CP009273/4281433‑4281729

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: