Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,281,581 G→T 20.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,5810GT20.1% 804.5 / 135.1 282D48Y (GAT→TAT) nrfEheme lyase (NrfEFG) for insertion of heme into c552, subunit NrfE
Reads supporting (aligned to +/- strand):  ref base G (116/109);  new base T (24/33);  total (140/142)
Fisher's exact test for biased strand distribution p-value = 2.36e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

Only 100 of 287 total aligned reads displayed.

GGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCGATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCCTGCTGTCGTGGGAACTCAAGCTGGCAGCGGTGTGGGGCGGCCATGAAGGTTCGCTGCTGCTTTGGGTGCTGCTGCTTTCCGCCTGGAGCGCGCTGTTTGCCTGGCATTAT  >  CP009273/4281432‑4281730
                                                                                                                                                     |                                                                                                                                                     
       cGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCGATTtctc                                                                                                                                                <  2:57696/150‑1 (MQ=255)
         tgtTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCGATTTCTCgg                                                                                                                                              <  1:899306/150‑1 (MQ=255)
         tgtTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCGATTTCTCgg                                                                                                                                              >  1:725618/1‑150 (MQ=255)
           ttGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCGATTTCTCGGTc                                                                                                                                            <  1:3449526/150‑1 (MQ=255)
           ttGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCGATTTCTCGGTc                                                                                                                                            >  1:2551652/1‑150 (MQ=255)
           ttGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCGATTTCTCGGTc                                                                                                                                            <  2:3124422/150‑1 (MQ=255)
             gACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCGATTTCTCGGTCAt                                                                                                                                          <  2:3682850/150‑1 (MQ=255)
             gACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCGATTTCTCGGTCAt                                                                                                                                          <  2:2486030/150‑1 (MQ=255)
              aCCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCGGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCGATTTCTCGGTCAtt                                                                                                                                         <  2:4378572/150‑1 (MQ=255)
                 ccGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCGATTTCTCGGTCATTTAt                                                                                                                                      >  1:4559923/1‑150 (MQ=255)
                   gTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCGATTTCTCGGTCATTTATGt                                                                                                                                    <  1:26134/150‑1 (MQ=255)
                   gTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCGATTTCTCGGTCATTTATGt                                                                                                                                    <  1:26011/150‑1 (MQ=255)
                   gTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCGATTTCTCGGTCATTTATGt                                                                                                                                    <  1:2710596/150‑1 (MQ=255)
                    ttGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGTCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCGATTTCTCGGTCATTTATGTc                                                                                                                                   <  2:3047862/150‑1 (MQ=255)
                     tGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCGATTTCTCGGTCATTTATGTcg                                                                                                                                  >  1:2479270/1‑150 (MQ=255)
                          gCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCTATTTCTCGGTCATTTATGTCGCGCaa                                                                                                                             <  1:69118/150‑1 (MQ=255)
                          gCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCGATTTCTCGGTCATTTATGTCGCGCaa                                                                                                                             <  2:3588014/150‑1 (MQ=255)
                           ccTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCGATTTCTCGGTCATTTATGTCGCGCAAc                                                                                                                            >  1:2266023/1‑150 (MQ=255)
                              tCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCGATTTCTCGGTCATTTATGTCGCGCAACATa                                                                                                                         >  2:680024/1‑150 (MQ=255)
                               cgcgGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCGATTTCTCGGTCATTTATGTCGCGCAACATAg                                                                                                                        <  2:104480/150‑1 (MQ=255)
                               cgcgGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCGATTTCTCGGTCATTTATGTCGCGCAACATAg                                                                                                                        >  1:2298025/1‑150 (MQ=255)
                               cgcgGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCGATTTCTCGGTCATTTATGTCGCGCAACATAg                                                                                                                        >  1:2114108/1‑150 (MQ=255)
                                  gggAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCGATTTCTCGGTCATTTATGTCGCGCAACATAGCta                                                                                                                     <  1:665168/150‑1 (MQ=255)
                                   ggAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCTATTTCTCGGTCATTTATGTCGCGCAACATAGCtat                                                                                                                    <  2:798378/150‑1 (MQ=255)
                                   ggAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCGATTTCTCGGTCATTTATGTCGCGCAACATAGCtat                                                                                                                    >  1:4046327/1‑150 (MQ=255)
                                     aGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCGATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAg                                                                                                                  >  2:622191/1‑150 (MQ=255)
                                      gTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCGATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGc                                                                                                                 <  1:4630239/150‑1 (MQ=255)
                                             tGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCGATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCCTGCTGt                                                                                                          <  2:2070047/150‑1 (MQ=255)
                                             tGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCGATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCCTGCTGt                                                                                                          <  1:3627026/150‑1 (MQ=255)
                                               cgcTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCGATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCCTGCTGTCg                                                                                                        <  1:1298195/150‑1 (MQ=255)
                                               cgcTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCGATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCCTGCTGTCg                                                                                                        >  1:1083385/1‑150 (MQ=255)
                                                gcTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCTATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCCTGCTGTCGt                                                                                                       <  2:3017606/150‑1 (MQ=255)
                                                 cgggCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCGATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCCTGCTGTCGTg                                                                                                      <  2:1228191/148‑1 (MQ=255)
                                                  tGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCGATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCCTGCTGTCGTgg                                                                                                     >  2:942855/1‑150 (MQ=255)
                                                     ccTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCGATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCCTGCTGTCGTGGGaa                                                                                                  <  1:213603/150‑1 (MQ=255)
                                                       tGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCGATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCCTGCTGTCGTGGGAACt                                                                                                <  1:2996236/150‑1 (MQ=255)
                                                         ccATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCTATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCCTGCTGTCGTGGGAACTCa                                                                                              >  2:3945726/1‑150 (MQ=255)
                                                             gatgCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCGATGTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCCTGCTGTCGTGGGAACTCAAGCt                                                                                          >  2:1515320/1‑150 (MQ=255)
                                                              atgCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCGATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCCTGCTGTCGTGGGAACTCAAGCTg                                                                                         <  2:4106704/150‑1 (MQ=255)
                                                               tgCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCGATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCCTGCTGTCGTGGGAACTCAAGCTgg                                                                                        >  1:3733625/1‑150 (MQ=255)
                                                                    cTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCGATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCCTGCTGTCGTGGGAACTCAAGCTGGCAGCg                                                                                   >  1:3720344/1‑150 (MQ=255)
                                                                    cTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCGATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCCTGCTGTCGTGGGAACTCAAGCTGGCAGCg                                                                                   >  1:2004760/1‑150 (MQ=255)
                                                                    cTCACTTGCATCGGCATGCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTGTGGCGTACTGACGGATTGTTTTCTCATCAGCGCTTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCCTGCTGTCGTGGGAACTCAAGCTGGCAGCg                                                                                   >  2:2174612/1‑150 (MQ=255)
                                                                     tCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCGATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCCTGCTGTCGTGGGAACTCAATCTGGCAGCgg                                                                                  >  1:3297417/1‑150 (MQ=255)
                                                                      cACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCTATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCCTGCTGTCGTGGGAACTCAAGCTGGCAGCggg                                                                                 >  1:3879013/1‑149 (MQ=255)
                                                                      cACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCGATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCCTGCTGTCGTGGGAACTCAAGCTGGCAGCGgt                                                                                 >  2:2563724/1‑150 (MQ=255)
                                                                       tcTTGCATCGGAATTCTGGCGCAGTTCTCGCTCCTGCTGCTCGACTTTGGCGTACTGACGTATTGTTTTCTCATCAGCTATTTCTAGGTCATTTATGTCGCGCAACATTGCTATATCCTGCTGTCGTGGGAACTCAAGCTGTCAGCGgtg                                                                                <  1:2420074/149‑1 (MQ=255)
                                                                        cTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCGATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCCTGCTGTCGTGGGAACTCAAGCTGGCAGCGgtgt                                                                               <  2:2527311/150‑1 (MQ=255)
                                                                         ttGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCTATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCCTGCTGTCGTGGGAACTCAAGCTGGCAGCGGtgtg                                                                              <  2:1406918/150‑1 (MQ=255)
                                                                           gCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCATACTGACGTATTGTTTTCTCATCAGCTATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCCTGCTGTCGTGGGAACTCAAGCTGGCAGCGGTGTggg                                                                            >  1:2931093/1‑150 (MQ=255)
                                                                             aTCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCGATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCCTGCTGTCGTGGGAACTCAAGCTGGCAGCGGTGTGGggc                                                                          >  1:4259026/1‑150 (MQ=255)
                                                                              tCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCTATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCCTGCTGTCGTGGGAACTCAAGCTGGCAGCGGTGTGGggcg                                                                         >  1:236426/1‑150 (MQ=255)
                                                                               cGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCGATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCCTGCTGTCGTGGGAACTCAAGCTGGCAGCGGTGTGGggcgg                                                                        <  1:3991603/150‑1 (MQ=255)
                                                                                ggCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCGATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCCTGCTGTCGTGGGAACTCAAGCTGGCAGCGGTGTGGggcggc                                                                       >  2:1398142/1‑150 (MQ=255)
                                                                                 gCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCGATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCCTGCTGTCGTGGGACCTCAAGCTGGCAGCGGTGTGGGGCGGcc                                                                      >  2:2430164/1‑150 (MQ=255)
                                                                                   aTTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCGATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCCTGCTGTCGTGGGAACTCAAGCTGGCAGCGGTGTGGGGCGGCCAt                                                                    <  2:1748732/150‑1 (MQ=255)
                                                                                    ttCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCCTCAGCGATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCCTGCTGTCGTGGGAACTCAAGCTGGCAGCGGTGTGGGGCGGCCATg                                                                   <  1:2875419/150‑1 (MQ=255)
                                                                                        ggCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCGATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCCTGCTGTCGTGGGAACTCAAGCTGGCAGCGGTGTGGGGCGGCCATGAAgg                                                               >  1:3635740/1‑150 (MQ=255)
                                                                                         gcgcAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCTATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCCTGCTGTCGTGGGAACTCAAGCTGGCAGCGGTGTGGGGCGGCCATGAAGGt                                                              <  2:3616779/150‑1 (MQ=255)
                                                                                         gcgcAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCGATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCCTGCTGTCGTGGGAACTCAAGCTGGCAGCGGTGTGGGGCGGCCATGAAGGt                                                              <  2:201982/150‑1 (MQ=255)
                                                                                           gcAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCTATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCCTGCTGTCGTGGGAACTCAAGCTGGCAGCGGTGTGGGGCGGCCATGAAGGTTc                                                            <  1:1849745/150‑1 (MQ=255)
                                                                                              gTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCGATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCCTGCTGTCGTGGGAACTCAAGCTGGCAGCGGTGTGGGGCGGCCATGAAGGTTCgct                                                         >  1:3960935/1‑150 (MQ=255)
                                                                                                 cgcgCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCGATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCCTGCTGTCGTGGGAACTCAAGCTGGCAGCGGTGTGGGGCGGCCATGAAGGTTCgctgct                                                      >  1:2939657/1‑150 (MQ=255)
                                                                                                 cgcgCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCGATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCCTGCTGTCGTGGGAACTCAAGCTGGCAGCGGTGTGGGGCGGCCATGAAGGTTCgctgct                                                      <  2:2105470/150‑1 (MQ=255)
                                                                                                   cgcTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCGATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCCTGCTGTCGTGGGAACTCAAGCTGGCAGCGGTGTGGGGCGGCCATGAAGGTTCgctgctgc                                                    >  2:4507037/1‑150 (MQ=255)
                                                                                                   cgcTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCGATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCCTGCTGTCGTGGGAACTCAAGCTGGCAGCGGTGTGGGGCGGCCATGAAGGTTCgctgctgc                                                    >  2:3250012/1‑150 (MQ=255)
                                                                                                   cgcTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCGATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCCTGCTGTCGTGGGAACTCAAGCTGGCAGCGGTGTGGGGCGGCCATGAAGGTTCgctgctgc                                                    <  1:1917254/150‑1 (MQ=255)
                                                                                                     cTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCTATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCCTGCTGTCGTGGGAACTCAAGCTGGCAGCGGTGTGGGGCGGCCATGAAGGTTCGCTGCTGCtt                                                  >  1:1655267/1‑150 (MQ=255)
                                                                                                     cTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCTATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCCTGCTGTCGTGGGAACTCAAGCTGGCAGCGGTGTGGGGCGGCCATGAAGGTTCGCTGCTGCtt                                                  >  1:3230182/1‑150 (MQ=255)
                                                                                                     cTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCGATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCCTGCTGTCGTGGGAACTCAAGCTGGCAGCGGTGTGGGGCGGCCATGAAGGTTCGCTGCTGCtt                                                  <  1:988824/150‑1 (MQ=255)
                                                                                                       ccTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCGATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCCTGCTGTCGTGGGAACTCAAGCTGGCAGCGGTGTGGGGCGGCCATGAAGGTTCGCTGCTGCTTTg                                                <  2:3304753/150‑1 (MQ=255)
                                                                                                           ctgctCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCGATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCCTGCTGTCGTGGGAACTCAAGCTGGCAGCGGTGTGGGGCGGCCATGAAGGTTCGCTGCTGCTTTGGGtg                                            <  2:2757289/150‑1 (MQ=255)
                                                                                                           ctgctCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCGATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCCTGCTGTCGTGGGAACTCAAGCTGGCAGCGGTGTGGGGCGGCCATGAAGGTTCGCTACTGCTTTGGGtg                                            >  1:1441136/1‑150 (MQ=255)
                                                                                                            tgctCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCGATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCCTGCTGTCGTGGGAACTCAAGCTGGCAGCGGTGTGGGGCGGCCATGAAGGTTCGCTGCTGCTTTGGGtgc                                           >  2:4559245/1‑150 (MQ=255)
                                                                                                                cGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCGATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCCTGCTGTCGTGGGAACTCAAGCTGGCAGCGGTGTGGGGCGGCCATGAAGGTTCGCTGCTGCTTTGGGtgctgct                                       >  1:10631/1‑150 (MQ=255)
                                                                                                                cGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCGATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCCTGCTGTCGTGGGAACTCAAGCTGGCAGCGGTGTGGGGCGGCCATGAAGGTTCGCTGCTGCTTTGGGtgctgct                                       >  2:3065322/1‑150 (MQ=255)
                                                                                                                cGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCGATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCCTGCTGTCGTGGGAACTCAAGCTGGCAGCGGTGTGGGGCGGCCATGAAGGTTCGCTGCTGCTTTGGGtgctgct                                       <  1:4486936/150‑1 (MQ=255)
                                                                                                                cGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCGATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCCTGCTGTCGTGGGAACTCAAGCTGGCAGCGGTGTGGGGCGGCCATGAAGGTTCGCTGCTGCTTTGGGtgctgct                                       >  1:3583894/1‑150 (MQ=255)
                                                                                                                cGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCGATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCCTGCTGTCGTGGGAACTCAAGCTGGCAGCGGTGTGGGGCGGCCATGAAGGTTCGCTGCTGCTTTGGGtgctgct                                       <  2:1086760/150‑1 (MQ=255)
                                                                                                                 gCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCGATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCCTGCTGTCGTGGGAACTCAAGCTGGCAGCGGTGTGGGGCGGCCATGAAGGTTCGCTGCTGCTTTGGGtgctgctg                                      <  1:3703119/150‑1 (MQ=255)
                                                                                                                 gCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCGATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCCTGCTGTCGTGGGAACTCAAGCTGGCAGCGGTGTGGGGCGGCCATGAAGGTTCGCTGCTGCTTTGGGtgctgctg                                      <  2:3706279/150‑1 (MQ=255)
                                                                                                                   cTTTGGCGTACTGACGTATTGTTTTCTCATCAGCGATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCCTGCTGTCGTGGGAACTCAAGCTGGCAGCGGTGTGGGGCGGCCATGAAGGTTCGCTGCTGCTTTGGGTgctgctgct                                    >  2:106501/1‑150 (MQ=255)
                                                                                                                    tttGGCGTACTGACGTATTGTTTTCTCATCAGCGATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCCTGCTGTCGTGGGAACTCAAGCTGGCAGCGGTGTGGGGCGGCCATGAAGGTTCGCTGCTGCTTTGGGTGCTGCTGCtt                                   <  1:2362601/150‑1 (MQ=255)
                                                                                                                      tGGCGTACTGACGTATTGTTTTCTCATCAGCGATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCCTGCTGTCGTGGGAACTCAAGCTGGCAGCGGTGTGGGGCGGCCATGAAGGTTCGCTGCTGCTTTGGGTGCTGCTGCTTTc                                 <  1:3510592/150‑1 (MQ=255)
                                                                                                                      tGGCGTACTGACGTATTGTTTTCTCATCAGCGATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCCTGCTGTCGTGGGAACTCAAGCTGGCAGCGGTGTGGGGCGGCCATGAAGGTTCGCTGCTGCTTTGGGTGCTGCTGCTTTc                                 <  1:3411775/150‑1 (MQ=255)
                                                                                                                         cGTACTGTCGTCTTGTGGTGACCTCCGCGATTTCTAGGTCTTTTCTGACGCGCTTCCTAGCTATCCCCCGCTGTCGTGGGAACTCATTCTGGGAGGGGTGTGGGGCGGCCATGAAGGTTCGCTGCTGCTTTGGGTGCTGCTGCTTTCCGc                              <  2:2585190/150‑1 (MQ=255)
                                                                                                                          gTACTGACGTATTGTTTTCTCATCAGCGATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCCTGCTGTCGTGGGAACTCAAGCTGGCAGCGGTGTGGGGCGGCCATGAAGGTTCGCTGCTGCTTTGGGTGCTGCTGCTTTCCGcc                             <  2:417645/150‑1 (MQ=255)
                                                                                                                             cTGACGTATTGTTTTCTCATCAGCGATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCCTGCTGTCGTGGGAACTCAAGCTGGCAGCGGTGTGGGGCGGCCATGAAGGTTCGCTGCTGCTTTGGGTGCTGCTGCTTTCCGCCTgg                          >  2:232655/1‑150 (MQ=255)
                                                                                                                                  gTATTGTTTTCTCATCAGCGATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCCTGCTGTCGTGGGAACTCAAGCTGGCAGCGGTGTGGGGCGGCCATGAAGGTTCGCTGCTGCTTTGGGTGCTGCTGCTTTCCGCCTGGAgcgc                     >  2:1320320/1‑150 (MQ=255)
                                                                                                                                     ttGTTTTCTCATCAGCTATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCCTGCTGTCGTGGGAACTCAAGCTGGCAGCGGTGTGGGGCGGCCATGAAGGTTCGCTGCTGCTTTGGGTGCTGCTGCTTTCCGCCTGGAGCGCGCt                  <  2:3372008/150‑1 (MQ=255)
                                                                                                                                     ttGTTTTCTCATCAGCTATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCCTGCTGTCGTGGGAACTCAAGCTGGCAGCGGTGTGGGGCGGCCATGAAGGTTCGCTGCTGCTTTGGGTGCTGCTGCTCTCCGCCTGGAGCGCGCt                  <  2:3718484/150‑1 (MQ=255)
                                                                                                                                      tGTTTTCTCATCAGCTATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCCTGCTGTCGTGGGAACTCAAGCTGGCAGCGGTGTGGGGCGGCCATGAAGGTTCGCTGCTGCTTTGGGTGCTGCTGCTTTCCGCCTGGAGCGCGCTg                 <  1:4379950/150‑1 (MQ=255)
                                                                                                                                      tGTTTTCTCATCAGCTAGTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCCTGCTGTCGTGGGAACTCAAGCTGGCAGCGGTGTGGGGCGGCCATGAAGGTTCGCTGCTGCTTTGGGTGCTGCTGCTTTCCGCCTGGAGCGCGCTg                 <  2:2628659/150‑1 (MQ=255)
                                                                                                                                        ttttCTCATCAGCGATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCCTGCTGTCGTGGGAACTCAAGCTGGCAGCGGTGTGGGGCGGCCATGAAGGTTCGCTGCTGCTTTGGGTGCTGCTGCTTTCCGCCTGGAGCGCGCTGtt               <  2:1854838/150‑1 (MQ=255)
                                                                                                                                           tctcATCAGCTATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCCTGCTGTCGTGGGAACTCAAGCTGGCAGCGGTGTGGGGCGGCCATGAAGGTTCGCTGCTGCTTTGGGTGCTGCTGCTTTCCGCCTGGAGCGCGCTGTTTGc            <  2:322260/150‑1 (MQ=255)
                                                                                                                                               atcaGCGATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCCTGCTGTCGTGGGAACTCAAGCTGGCAGCGGTGTGGGGCGGCCATGAAGGTTCGCTGCTGCTTTGGGTGCTGCTGCTTTCCGCCTGGAGCGCGCTGTTTGCCTgg        >  1:3091437/1‑150 (MQ=255)
                                                                                                                                               atcaGCGATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCCTGCTGTCGTGGGAACTCAAGCTGGCAGCGGTGTGGGGCGGCCATGAAGGTTCGCTGCTGCTTTGGGTGCTGCTGCTTTCCGCCTGGAGCGCGCTGTTTGCCTgg        >  1:774370/1‑150 (MQ=255)
                                                                                                                                                 caGCGATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCCTGCTGTCGTGGGAACTCAAGCTGGCAGCGGTGTGGGGCGGCCATGAAGGTTCGCTGCTGCTTTGGGTGCTGCTGCTTTCCGCCTGGAGCGCGCTGTTTGCCTGGCa      >  1:1784599/1‑150 (MQ=255)
                                                                                                                                                     gATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCCTGCTGTCGTGGGAACTCAAGCTGGCAGCGGTGTGGGGCGGCCATGAAGGTTCGCTGCTGCTTTGGGTGCTGCTGCTTTCCGCCTGGAGCGCGCTGTTTGCCTGGCATTAt  >  2:1735374/1‑150 (MQ=255)
                                                                                                                                                     gATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCCTGCTGTCGTGGGAACTCAAGCTGGCAGCGGTGTGGGGCGGCCATGAAGGTTCGCTGCTGCTTTGGGTGCTGCTGCTTTCCGCCTGGAGCGCGCTGTTTGCCTGGCATTAt  >  1:515656/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCGATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCCTGCTGTCGTGGGAACTCAAGCTGGCAGCGGTGTGGGGCGGCCATGAAGGTTCGCTGCTGCTTTGGGTGCTGCTGCTTTCCGCCTGGAGCGCGCTGTTTGCCTGGCATTAT  >  CP009273/4281432‑4281730

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: