Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,844,263 T→G 5.8% F356V (TTC→GTC)  PP_4259 → iron‑sulfur cluster‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,844,2630TG 5.8% 248.5 / 4.9 86F356V (TTC→GTC) PP_4259iron‑sulfur cluster‑binding protein
Reads supporting (aligned to +/- strand):  ref base T (42/39);  new base G (3/2);  total (45/41)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.87e-01

CATCATGGACAAGATGAACTACCCCAAAGGGCTGATCAGCTACACCACCGAGCACAACCTGTCCGGGCAGAAGACCCACATGCTGCGCCCCCGTCTGATCGGCTATGCCGTGGTACTGCTGGTGATGATG‑‑ATCGGGTTGGCGACGGCCTTCGCCACCCGTTCGCTGGTCGGTTTCGACGTCAGCAAAGACCGCGTGTTGTACCGCGAAAACGCCCAGGGCCGGATCGAGAACGTGTACAGCCTGAAGGTAATGAACAAGGATCAGCGCGATCACCTCTACGTGCTGGACGCCGCC  >  NC_002947/4844115‑4844409
                                                                                                                                                      |                                                                                                                                                   
catcatGGACAAGATGAACTACCCCAAAGGGCTGATCAGCTACACCACCGAGCACAACCTGTCCGGGCAGAAGACCCACATGCTGCGCCCCCGTCTGATCGGCTATGCCGTGGTACTGCTGGTGATGATG‑‑ATCGGGTTGGCGACGGCCtt                                                                                                                                                    <  1:414793/150‑1 (MQ=255)
        aCAAGATGAACTACCCCAAAGGGCTGATCAGCTACACCACCGAGCACAACCTGTCCGGGCAGAAGACCCACATGCTGCGCCCCCGTCTGATCGGCTATGCCGTGGTACTGCTGGTGATGATG‑‑ATCGGGTTGGCGACGGCCTTCGCCAccc                                                                                                                                           <  1:1298596/150‑1 (MQ=255)
        aCAAGATGAACTACCCCAAAGGGCTGATCAGCTACACCACCGAGCACAACCTGTCCGGGCAGAAGACCCACATGCTGCGCCCCCGTCTGATCGGCTATGCCGTGGTACTGACGGTTTGGATG‑‑ATCGGGTTGGCGACGGCCTTCGCCAccc                                                                                                                                           <  1:460087/150‑1 (MQ=255)
            gATGAACTACCCCAAAGGGCTGATCAGCTACACCACCGAGCACAACCTGTCCGGGCAGAAGACCCACATGCTGCGCCAACGTATGATTGCCTCACGGGCGGTACTTCGGGGGATGATG‑‑ATCGGGTTGGCGACGGCCTTCGCCACCCGTTc                                                                                                                                       >  1:717565/1‑150 (MQ=255)
              tGAACTACCCCAAAGGGCTGATCAGCTACACCACCGAGCACAACCTGTCCGGGCAGAAGACCCACATGCTGCGCCCCCGTCTGATCGGCTATGCCGTGGTACTGCTGGTGATGATG‑‑ATCGGGTTGGCGACGGCCTTCGCCACCCGTTCGc                                                                                                                                     >  2:369994/1‑150 (MQ=255)
               gAACTACCCCAAAGGGCTGATCAGCTACACCACCGAGCACAACCTGTCCGGGCAGAAGACCCACATGCTGCGCCCCCGTCTGATCGGCTATGCCGTGGTACTGCTGGTGATGATG‑‑ATCGGGTTGGCGACGGCCTTCGCCACCCGTTCGCt                                                                                                                                    <  1:9364/150‑1 (MQ=255)
                    aCCCCAAAGGGCTGATCAGCTACACCACCGAGCACAACCTGTCCGGGCAGAAGACCCACATGCTGCGCCCCCGTCTGATCGGCTATGCCGTGGTACTGCTGGTGATGATG‑‑ATCGGGTTGGCGACGGCCTTCGCCACCCGTTCGCTGGTCg                                                                                                                               <  1:175264/150‑1 (MQ=255)
                       ccAAAGGGCTGATCAGCTACACCACCGAGCACAACCTGTCCGGGCAGAAGACCCACATGCTGCGCCCCCGTCTGATCGGCTATGCCGTGGTACTGCTGGTGATGATG‑‑ATCGGGTTGGCGACGGCCTTCGCCACCCGTTCGCTGGTCGGtt                                                                                                                            >  1:499210/1‑150 (MQ=255)
                       ccAAAGGGCTGATCAGCTACACCACCGAGCACAACCTGTCCGGGCAGAAGACCCACATGCTGCGCCCCCGTCTGATCGGCTATGCCGTGGTACTGCTGGTGATGATG‑‑ATCGGGTTGGCGACGGCCTTCGCCACCCGTTCGCTGGTCGGtt                                                                                                                            >  2:1760538/1‑150 (MQ=255)
                              gCTGATCAGCTACACCACCGAGCACAACCTGTCCGGGCAGAAGACCCACATGCTGCGCCCCCGTCTGATCGGCTATGCCGTGGTACTGCTGGTGATGATG‑‑ATCGGGTTGGCGACGGCCTTCGCCACCCGTTCGCTGGTCGGTTTCGACGt                                                                                                                     <  2:1048174/150‑1 (MQ=255)
                                  aTCAGCTACACCACCGAGCACAACCTGTCCGGGCAGAAGACCCACATGCTGCGCCCCCGTCTGATCGGCTATGCCGTGGTACTGCTGGTGATGATG‑‑ATCGGGTTGGCGACGGCCTTCGCCACCCGTTc                                                                                                                                       >  1:1454393/1‑128 (MQ=255)
                                  aTCAGCTACACCACCGAGCACAACCTGTCCGGGCAGAAGACCCACATGCTGCGCCCCCGTCTGATCGGCTATGCCGTGGTACTGCTGGTGATGATG‑‑ATCGGGTTGGCGACGGCCTTCGCCACCCGTTc                                                                                                                                       <  2:1454393/128‑1 (MQ=255)
                                  aTCAGCTACACCACCGAGCACAACCTGTCCGGGCAGAAGACCCACATGCTGCGCCCCCGTCTGATCGGCTATGCCGTGGTACTGCTGGTGATGATG‑‑ATCGGGGTGGCGAAGGCCGTCGCCAacccga                                                                                                                                        >  2:1620370/1‑124 (MQ=255)
                                  aTCAGCTACACCACCGAGCACAACCTGTCCGGGCAGAAGACCCACATGCTGCGCCCCCGTCTGATCGGCTATGCCGTGGTACTGCTGGTGATGATG‑‑ATCGGGGTGGCGAAGGCCGTCGCCAacccga                                                                                                                                        <  1:1620370/127‑4 (MQ=255)
                                  aTCAGCTACACCACCGAGCACAACCTGTCCGGGCAGAAGACCCACATGCTGCGCCCCCGTCTGATCGGCTATGCCGTGGTACTGCTGGTGACGATGCCATTGGGAT‑‑CTATGGCCTTCGCCACCCGTTCGCTGGTCGGTTTCGACGTCAGc                                                                                                                 >  2:542927/1‑150 (MQ=255)
                                   tCAGCTACACCACCGAGCACAACCTGTCCGGGCAGAAGACCCACATGCTGCGCCCCCGTCTGATCGGCTATGCCGTGGTACTGCTGGTGATGATG‑‑ATCGGGTTGGCGACGGCCTTCGCCACCCGTTCGCTGGTCGGTTTCGACGTCAGCa                                                                                                                >  1:605636/1‑150 (MQ=255)
                                    cAGCTACACCACCGAGCACAACCTGTCCGGGCAGAAGACCCACATGCTGCGCCCCCGTCTGATCGGCTATGCCGTGGTACTGCTGGTGATGATG‑‑ATCGGGTTGGCGACGGCCTTCGCCACCCGTTCGCTGGTCGGTTTCGACGTCAGCaa                                                                                                               <  1:3568/150‑1 (MQ=255)
                                     aGCTACACCACCGAGCACAACCTGTCCGGGCAGAAGACCCACATGCTGCGCCCCCGTCTGATCGGCTATGCCGTGGTACTGCTGGTGATGATG‑‑ATCGGGTTGGCGACGGCCTTCGCCACCCGTTCGCTGGTCGGTTTCGACGTCAGCaaa                                                                                                              >  2:1513968/1‑150 (MQ=255)
                                      gCTACACCACCGAGCACAACCTGTCCGGGCAGAAGACCCACATGCTGCGCCCCCGTCTGATCGGCTATGCCGTGGTACTGCTGGTGATGATG‑‑ATTGGGTTGGCGACGGCCTTCGc                                                                                                                                                 >  2:1783924/1‑115 (MQ=255)
                                      gCTACACCACCGAGCACAACCTGTCCGGGCAGAAGACCCACATGCTGCGCCCCCGTCTGATCGGCTATGCCGTGGTACTGCTGGTGATGATG‑‑ATTGGGTTGGCGACGGCCTTCGc                                                                                                                                                 <  1:1783924/115‑1 (MQ=255)
                                      gCTACACCACCGAGCACAACCTGTCCGGGCAGAAGACCCACATGCTGCGCCCCCGTCTGATCGGCTATGCCGTGGTACTGCTGGTGATGATG‑‑ATCGGGTTGGCGACGGCCTTCGCCACCCGTTCGCTTGTCGGTTTCGACGTCAGCAAAg                                                                                                             >  2:881271/1‑150 (MQ=255)
                                          caccacCGAGCACAACCTGTCCGGGCAGAAGACCCACATGCTGCGCCCCCGTCTGATCGGCTATGCCGTGGTACTGCTGGTGATGATG‑‑ATCGGGTTGGCGACGGCCTTCGCCACCCGTTCGCTGGTCGGTTTCGACGTCAGCAAAGACcg                                                                                                         <  1:430067/150‑1 (MQ=255)
                                          caccacCGAGCACAACCTGTCCGGGCAGAAGACCCACATGCTGCGCCCCCGTCTGATCGGCTATGCCGTGGTACTGCTGGTGATGATG‑‑ATCGGGTTGGCGACGGCCTTCGCCACCCGTTCGCTGGTCGGTTTCGACGTCAGCAAAGACcg                                                                                                         <  2:605636/150‑1 (MQ=255)
                                             caccGAGCACAACCTGTCCGGGCAGAAGACCCACATGCTGCGCCCCCGTCTGATCGGCTATGCCGTGGTACTGCTGGTGATGATG‑‑ATCGGGTTGGCGACGGCCTTc                                                                                                                                                   <  2:1343595/106‑1 (MQ=255)
                                             caccGAGCACAACCTGTCCGGGCAGAAGACCCACATGCTGCGCCCCCGTCTGATCGGCTATGCCGTGGTACTGCTGGTGATGATG‑‑ATCGGGTTGGCGACGGCCTTc                                                                                                                                                   >  1:1343595/1‑106 (MQ=255)
                                             caccGAGCACAACCTGTCCGGGCAGAAGACCCACATGCTGCGCCCCCGTCTGATCGGCTATGCCGTGGTACTGCTGGTGATGATG‑‑ATCGGGTTGGCGACGGCCTTCGCCACCCGTTCGCTGGTCGGTTTCGACGTCAGCAAAGACCGCgt                                                                                                      >  1:1383613/1‑150 (MQ=255)
                                             caccGAGCACAACCTGTCCGGGCAGAAGACCCACATGCTGCGCACCCGTCTGATCGGCTATGCCGTGGTACTGCTGGTGATGATG‑‑ATCGGGTTGGCGACGGCCTTCGCCACCCGTTCGCTGGTCGGTTTCGACGTCAGCAAAGACCGCgt                                                                                                      >  2:1403423/1‑150 (MQ=255)
                                               ccacccacaACCCGTCCGAGCACATGACCCACATGCTTCCCCCCCGTATGATCAGCTATCCCGTGGTACTGCTGGTTATGATG‑‑ATCGGGTTTACGACGGCCTTCTCCACCCGTTCGCTCCTCTGTTTCGACGTCAGCAAAGACCGCgtgt                                                                                                    <  2:793274/145‑1 (MQ=255)
                                               ccGAGCACAACCTGTCCGGGCAGAAGACCCACATGCTGCGCCCCCGTCTGATCGGCTATGCCGTGGTACTGCTGGTGATGATG‑‑ATCGGGTTGGCGACGGCCTTCGCCACCCGTTCGCTGGTCGGTTTCGACGTCAGCAAAGACCGCgtgt                                                                                                    >  2:571216/1‑150 (MQ=255)
                                                       aaCCTGTCCGGGCAGAAGACCCACATGCTGCGCCCCCGTCTGATCGGCTATGCCGTGGTACTGCTGGTGATGATG‑‑ATCGGGTTGGCGACGGCCTTCGCCACCCGTTCGCTGGTCGGTTTCGACGTCAGCAAAGACCGCGTGTTGTACcgc                                                                                            >  1:1464512/1‑150 (MQ=255)
                                                          cTGTCCGGGCAGAAGACCCACATGCTGCGCCCCCGTCTGATCGGCTATGCCGTGGTACTGCTGGTGATGATG‑‑ATCGGGTTGGCGACGGCCTTCGCCACCCGTTCGCTGGTCGGTTTCGACGTCAGCAAAGACCGCGtgtt                                                                                                   <  1:1252207/140‑1 (MQ=255)
                                                          cTGTCCGGGCAGAAGACCCACATGCTGCGCCCCCGTCTGATCGGCTATGCCGTGGTACTGCTGGTGATGATG‑‑ATCGGGTTGGCGACGGCCTTCGCCACCCGTTCGCTGGTCGGTTTCGACGTCAGCAAAGACCGCGtgtt                                                                                                   >  2:1252207/1‑140 (MQ=255)
                                                          cTGTCCGGGCAGAAGACCCACATGCTGCGCCCCCGTCTGATCGGCTATGCCGTGGTACTGCTGGTGATGATG‑‑ATCGGGTTGGCGACGGCCTTCGCCACCCGTTCGCTGGTCGGTTTCGACGTCAGCAAAGACCGCGTGTTGTACCGCGaa                                                                                         >  2:1694199/1‑150 (MQ=255)
                                                                gggCAGAAGACCCACATGCTGCGCCCCCGTCTGATCGGCTATGCCGTGGTACTGCTGGTGATGATG‑‑ATCGGGTTGGCGACGGCCTTCGCCACCCGTTCGCTGGTCGGTTTCGAAGTACGCAAACTCCGAGTGTTGTACCCCGAAAACGcc                                                                                   >  1:336317/1‑150 (MQ=255)
                                                                 ggCAGAAGACCCACATGCTGCGCCCCCGTCTGATCGGCTATGCCGTGGTACTGCTGGTGATGATG‑‑ATCGGGTTGGCGACGGCCTTCGCCACCCGTTCGCTGGTCGGTTTCGACGTCAGCAAAGACCGCGTGTTGTACCGCGAAAACGccc                                                                                  >  1:480882/1‑150 (MQ=255)
                                                                     gaagaCCCATATGCTGCGCCCCCCTCTGATCGGCTATGCCGTGGTACTGCTGGTGATGATG‑‑ATCGGGTTGGCGACGGCCTTCGCCACCCGTTCGCTGGTCGGTTTCGACGTCAGCAAAGACCGCGTGTTGTACCGCGAAAACGCCCAggg                                                                              <  2:112764/150‑1 (MQ=255)
                                                                         aCCCACATGCTGCGCCCCCGTCTGATCGGCTATGCCGTGGTACTGCTGGTGATGATG‑‑ATCGGGTTGGCGACGGCCTTCGCCACCCGTTCGCTGGTCGGTTTCGACGTCAGCAAAGACCGCGTGTTGTACCGCGAAAACGCCCAGGGCCgg                                                                          >  2:819647/1‑150 (MQ=255)
                                                                           ccACATGCTGCGCCCCCGTCTGATCGGCTATGCCGTGGTACTGCTGGTGATGATG‑‑ATCGGGTTGGCGACGGCCTTCGCCACCCGTTCGCTGGTCGGTTTCGACGTCAGCAAAGACcgc                                                                                                        >  1:704535/1‑118 (MQ=255)
                                                                           ccACATGCTGCGCCCCCGTCTGATCGGCTATGCCGTGGTACTGCTGGTGATGATG‑‑ATCGGGTTGGCGACGGCCTTCGCCACCCGTTCGCTGGTCGGTTTCGACGTCAGCAAAGACcgc                                                                                                        <  2:704535/118‑1 (MQ=255)
                                                                           ccACATGCTGCGCCCCCGTCTGATCGGCTATGCCGTGGTACTGCTGGTGATGATG‑‑ATCGGGTTGGCGACGGCCTTCGCCACCCGTTCGCTGGTCGGTTTCGACGTCAGCAAAGACCGCGTGTTGTACCGCGAAAACGCCCAGGGCCGGAt                                                                        >  1:557776/1‑150 (MQ=255)
                                                                           ccACATGCTGCGCCCCCGTCTGATCGGCTATGCCGTGGTACTGCTGGTGATGATG‑‑ATCGGGTTGGCGACGGCCTTCGCCACCCGTTCGCTGGTCGGTTTCGACGTCAGCAAAGACCGCGTGTTGTACCGCGAAAACGCCCAGGGCCGGAt                                                                        >  1:866518/1‑150 (MQ=255)
                                                                           ccACATGCTGCGCCCCCGTCTGATCGGCTATGCCGTGGTACTGCTGGTGATGATG‑‑ATCGGGTTGGCGACGGCCTTCGCCACCCGTTCGCTGGTCGGTTTCGACGTCAGCAAAGACCGCGTGTTGTACCGCGAAAACGCCCAGGGCCGGAt                                                                        <  2:1021473/150‑1 (MQ=255)
                                                                           ccACATGCTGCGCCCCCGTCTGATCGGCTATGCCGTGGTACTGCTGGTGATGATG‑‑ATCGGGTTGGCGACGGCCTTCGCCACCCGTTCGCTGGTCGGTTTCGACGTCAGCAAAGACCGCGTGTTGTACCGCGAAAACGCCCAGGGCCGGAt                                                                        >  1:1021473/1‑150 (MQ=255)
                                                                             acaTGCTGCGCCCCCGTCTGATCGGCTATGCCGTGGTACTGCTGGTGATGATG‑‑ATCGGGTTGGCGACGGCCTTCGCCACCCGTTCGCTGGTCGGTTTCGACGTCAg                                                                                                                  >  2:956652/1‑106 (MQ=255)
                                                                             acaTGCTGCGCCCCCGTCTGATCGGCTATGCCGTGGTACTGCTGGTGATGATG‑‑ATCGGGTTGGCGACGGCCTTCGCCACCCGTTCGCTGGTCGGTTTCGACGTCAg                                                                                                                  <  1:956652/106‑1 (MQ=255)
                                                                               aTGCTGCGCCCCCGTCTGATCGGCTATGCCGTGGTACTGCTGGTGATGATG‑‑ATCGGGTTGGCGACGGCCTTCGCCACCCGTTCGCTGGTCGGTTTCGACGTCAGCAAAGACCGCGTGTTGTACCGCGAAAACGCCCAGGGCCGGATCgag                                                                    <  2:1464512/150‑1 (MQ=255)
                                                                                       cccccGTCTGATCGGCTATGCCGTGGTACTGCTGGTGATGATG‑‑ATCGGGTTGGCGACGGCCTTCGCCACCCGTTCGCTGGTCGGTTTCGACGTCAGCAAAGACCGCGTGTTGTACCGCGAAAACGCCCAGGGCCGGATCGAGAACGTGTa                                                            <  1:369994/150‑1 (MQ=255)
                                                                                         cccGTCTGATCGGCTATGCCGTGGTACTGCTGGTGATGATG‑‑ATCTGGTTGGCGACGGCCTTCGCCACCCGTTCGCTGGTCGGTTTCGACGTCAGCAAAGACCGCGTGTTGTACCGCGAAAACGCCCAGGGCCGGATCGAGAACGTGTACa                                                          >  1:1252610/1‑150 (MQ=255)
                                                                                          ccGTCTGATCGGCTATGCCGTGGTACTGCTGGTGATGATG‑‑ATCGGGTTGGCGACGGCCTTCGCCACCCGTTCGCTGGTCGGTTTCGACGTCAGCAAAGACCGCGTGTTGTACCGCGAAAACGCCCAGGGCCGGATc                                                                       <  1:1076005/136‑1 (MQ=255)
                                                                                          ccGTCTGATCGGCTATGCCGTGGTACTGCTGGTGATGATG‑‑ATCGGGTTGGCGACGGCCTTCGCCACCCGTTCGCTGGTCGGTTTCGACGTCAGCAAAGACCGCGTGTTGTACCGCGAAAACGCCCAGGGCCGGATc                                                                       >  2:1076005/1‑136 (MQ=255)
                                                                                          ccGTCTGATCGGCTATGCCGTGGTACTGCTGGTGATGATG‑‑ATCGGGTTGGCGACGGCCTTCGCCACCCGTTCGCTGGTCGGTTTCGACGTCAGCAAAGACCGCGTGTTGTACCGCGAAAACGCCCAGGGCCGGATCGAGAACGTGTACAg                                                         >  2:842592/1‑150 (MQ=255)
                                                                                             tCTGATCGGCTATGCCGTGGTACTGCTGGTGATGATG‑‑ATCGGGTTGGCGACGGCCTTCGCCACCCGTTCGCTGGTCGGTTTCGAc                                                                                                                       <  1:1546326/85‑1 (MQ=255)
                                                                                             tCTGATCGGCTATGCCGTGGTACTGCTGGTGATGATG‑‑ATCGGGTTGGCGACGGCCTTCGCCACCCGTTCGCTGGTCGGTTTCGAc                                                                                                                       >  2:1546326/1‑85 (MQ=255)
                                                                                                gATCGGCTATGCCGTGGTACTGCTGGTGATGATG‑‑ATCGGGTTGGCGACGGCCTTCGCCACCCGTTCGCTGGTCGGTTTCGACGTCAGCAAAGACCGCgtg                                                                                                     >  1:537415/1‑100 (MQ=255)
                                                                                                gATCGGCTATGCCGTGGTACTGCTGGTGATGATG‑‑ATCGGGTTGGCGACGGCCTTCGCCACCCGTTCGCTGGTCGGTTTCGACGTCAGCAAAGACCGCgtg                                                                                                     <  2:537415/100‑1 (MQ=255)
                                                                                                gATCGGCTATGCCGTGGTACTGCTGGTGATGATG‑‑ATCGGGTTGGCGACGGCCTTCGCCACCCGTTCGCTGGTCGGTTTCGACGTCAGCAAAGACCGCGTGTTGTACCGCGAAAACGCCCAGGGCCGGATCGAGAACGTGTACAGCCTGaa                                                   >  2:230905/1‑150 (MQ=255)
                                                                                                 aTCTGCTATGCCGTGGTACTGCTGGTGATGATG‑‑ATCGGGTTGGCGACGGCCTTCGCCACCCGTTCGCTGGTCGGTTTCGACGTCAGCAAAGACCGCGTGTTGTACCGCGAAAACGCCCAGGGCCGGATCGAGAACGTGTACAGCCTGAAg                                                  >  1:144258/1‑150 (MQ=255)
                                                                                                 aTCGGCTATGCTGTGGTACTGCTGGTGATGATG‑‑ATCGGGTTGGCGACGGCCTTCGCCACCCGTTCGCTGGTCGGTTTCGACGTCAGCAAAGACCGCGTGTTGTACCGCGAAAACGCCCAGGGCCGGATCGAGAACGTGTACAGCCTGAAg                                                  <  1:1561533/150‑1 (MQ=255)
                                                                                                 aTCGGCTATGCCGTGGTACTGCTGGTGATGATG‑‑ATCGGGTTGGCGACGGCCTTCGCCACCCGTTCGCTGGTCGGTTTCGACGTCAGCAAAGACCGCGTGTTGTACCGCGAAAACGCCCAGGGCCGGATCGAGAACGTGTACAGCCTGAAg                                                  <  2:499210/150‑1 (MQ=255)
                                                                                                 aTCGGCTATGCCGTGGTACTGCTGGTGATGATG‑‑ATCGGGTTGGCGACGGCCTTCGCCACCCGTTCGCTGGTCGGTTTCGACGTCAGCAAAGACCGCGTGTTGTACCGCGAAAACGCCCAGGGCCGGATCGAGAACGTGTACAGCCTGAAg                                                  >  2:1620150/1‑150 (MQ=255)
                                                                                                 aTCGGCTATGCCGTGGTACTGCTGGTGATGATG‑‑ATCGGGTTGGCGACGGCCTTCGCCACCCGTTCGCTGGTCGGTTTCGACGTCAGCAAAGACCGCGTGTTGTACCGCGAAAACGCCCAGGGCCGGATCGAGAACGTGTACAGCCTGAAg                                                  <  1:842592/150‑1 (MQ=255)
                                                                                                   cGGCTATGCCGTGGTACTGCTGGTGATGATG‑‑ATCGGGTTGGCGACGGCCTTCGCCACCCGTTCGCTGGTCGGTTTCGACGTCAGCAAAGACCGCGTGTTGTACCGCGAAAACGCCCAGGGCCGGATCGAGAACGTGTACAGCCTGAAGGt                                                <  1:1620150/150‑1 (MQ=255)
                                                                                                       tATGCCGTGGTACTGCTGGTGATGATG‑‑ATCGGGTTGGCGACGGCCTTCGCCACCCGTTCGCTGGTCGGTTTCGACGTCAGCAAAGACCGCGTGTTGTACCGCGAAAACGCCCAGGGCCGGATCGAGAACGTGTACAGCCTGAAGGTAATg                                            <  2:557776/150‑1 (MQ=255)
                                                                                                         tGCCGTGGTACTGCTGGTGATGATG‑‑ATCGGGTTGGCGACGGCCTTCGCCACCCGTTCGCTGGTCGGTTTCGACGTCAGCAAAGACCGCGTGTTGTACCGCGAAAACGCCCAGGGCCGGATCGAGAACGTGTACAGCCTGAAGGTAATGaa                                          <  1:1760538/150‑1 (MQ=255)
                                                                                                         tGCCGTGGTACTGCTGGTGATGATG‑‑ATCGGGTTGGCGACGGCCTTCGCCACCCGTTCGCTGGTCGGTTTCGACGTCAGCAAAGACCGCGTGTTGTACCGCGAAAACGCCCAGGGCCGGATCGAGAACGTGTACAGCCTGAAGGTAATGaa                                          <  1:819647/150‑1 (MQ=255)
                                                                                                          gCCGTGGTACTGCTGGTGATGATG‑‑ATCGGGTTGGCGACGGCCTTCGCCACCCGTTCGCTGGTCGGTTTCGACGTCAGCAAAGACCGCGTGTTGTACCGCGAAAACGCCCAGGGCCGGATCGAGAACGTGTACAGCCTGAAGGTAATGAAc                                         <  2:669092/150‑1 (MQ=255)
                                                                                                           ccGTGGTACTGCTGGTGATGATG‑‑ATCGGGTTGGCGACGGCCTTCGCCACCCGTTCGCTGGTCGGTTTCGACGTCAGCAAAGACCGCGTGTTGTACCGCGAAAACGCCCAGGGCCGGATCGAGAACGTGTACAGCCTGAAGGTAATGAACa                                        >  2:367289/1‑150 (MQ=255)
                                                                                                            cGTGGTACTGCTGGTGATGATG‑‑ATCGGGTTGGCGACGGCCTTCGCCACCCGTTCGCTGGTCGGTTTCGACGTCAGCAAAGACCGCGTGTTGTACCGCGAAAACGCCCAGGGCCGGATATTGGGCTTTTACAGCCTGAAGGTAATGAACaa                                       <  1:542927/150‑1 (MQ=255)
                                                                                                              tGGTACTGCTGGTGATGATG‑‑ATCGGGTTGGCGACGGCCTTCGCCACCCGTTCGCTGGTCGGTTTCGACGTCAGCAAAGACCGCGTGTTGTACCGCGAAAACGCCCAGGGCCGGATCGAGAACGTGTACAGCCTGAAGGTAATGAACAAgg                                     >  2:1794243/1‑150 (MQ=255)
                                                                                                                  aCTGCTGGTGATGATG‑‑ATCGGGTTGGCGACGGCCTTCGCCACCCGTTCGCTGGTCGGTTTCGACGTCAGCAAAGACCGCGTGTTGTACCGCGAAAACGCCCAGGGCCGGATCGAGAACGTGTACAGCCTGAAGGTAATGAACAAGGATCa                                 <  2:1013135/150‑1 (MQ=255)
                                                                                                                         gTGATGATG‑‑ATCGGGTTGGCGACGGCCTTCGCCACCCGTTCGCTGGTCGGTTTCGACGTCAGCAAAGACCGCGtgtt                                                                                                   >  2:1567756/1‑77 (MQ=255)
                                                                                                                         gTGATGATG‑‑ATCGGGTTGGCGACGGCCTTCGCCACCCGTTCGCTGGTCGGTTTCGACGTCAGCAAAGACCGCGtgtt                                                                                                   <  1:1567756/77‑1 (MQ=255)
                                                                                                                         gTGATGATG‑‑ATCGGGTTGGCGACGGCCTTCGCCACCCGTTCGCTGGTCGGTTTCGACGTCAGCAAAGACCGCGTGTTGTACCGCGAAAACGCCCAGGGCCGGATCGAGAACGTGTACAGCCTGAAGGTAATGAACAAGGATCAGCGCGAt                          <  2:144258/150‑1 (MQ=255)
                                                                                                                         gTGATGATG‑‑ATCGGGTTGGCGACGGCCTTCGCCACCCGTTCGCTGGTCGGTTTCGACGTCAGCAAAGACCGCGTGTTGTACCGCGAAAACGCCCAGGGCCGAGTCGGGAACGTGTACAGCCTGAAGGTAATGAACAAGGATCAGCGCGAt                          >  1:1667892/1‑150 (MQ=255)
                                                                                                                           gatgatg‑‑atCGGGGTGGCGACGGCCTTCGCCACCCGTTCGCTAGTCGGTTTCGACGTCAGCAAAGACCGCGTGTTGTACCGCGAAAACGCCCAGGGCCGGATCGAGAACGTGTACAGCCTGAAGGTAATGAACAAGGATCAGCGCGATCa                        >  2:1350954/1‑150 (MQ=255)
                                                                                                                              gatg‑‑atCGGGTTGGCGACGGCCTTCGCCACCCGTTCGCTGGTCGGTTTCGACGTCAGCAAAGACCGCGTGTTGTACCGCGAAAACGCCCAGGGCCGGATCGAGAACGTGTACAGCCTGAAGGTAATGAACAAGGATCAGCGCGATCACct                     >  2:1186844/1‑150 (MQ=255)
                                                                                                                                tg‑‑atCGGGTTGGCGACGGCCTTCGCCACCCGTTCGCTGGTCGGTTTCGACGTCAGCAAAGACCGCGTGTTGTACCGCGAAAACGcc                                                                                   >  1:1060801/1‑86 (MQ=255)
                                                                                                                                tg‑‑atCGGGTTGGCGACGGCCTTCGCCACCCGTTCGCTGGTCGGTTTCGACGTCAGCAAAGACCGCGTGTTGTACCGCGAAAACGcc                                                                                   <  2:1060801/86‑1 (MQ=255)
                                                                                                                                        gggtGGCGAAGGCCGTCGCCAACCCGATCGCTGGTCGGTTTCGACGTCAGCAAAGACCGCGTGTTGTACCGCGAAAACGCCCAGGGCCGGATCGAGAACGTGTACAGCCTGAAGGTAATGAACAAGGATCAGCGCGATCACCTCTAc                 <  2:1428155/144‑1 (MQ=255)
                                                                                                                                        gggtGGCGAAGGCCGTCGCCAACCCGATCGCTGGTCGGTTTCGACGTCAGCAAAGACCGCGTGTTGTACCGCGAAAACGCCCAGGGCCGGATCGAGAACGTGTACAGCCTGAAGGTAATGAACAAGGATCAGCGCGATCACCTCTAc                 >  1:1428155/4‑147 (MQ=255)
                                                                                                                                        gggtGGCGAAGGCCGTCGCCAACCGTTCGCTGGTCGGTTTCGACGTCAGCAAAGACCGCGTGTTGTACCGCGAAAACGCCCAGGGCCGGATCGAGAACGTGTACAGCCTGAAGGTAATGAACAAGGATCAGCGCGATCACCTCTACGTGc             >  2:126491/4‑150 (MQ=255)
                                                                                                                                        ggTTGGCGACGGCCTTCGCCACCCGTTCGCTGGTCGGTTTCGACGTCAGCAAAGACCGCGTGTTGTACCGCGAAAACGCCCAGGGCCGGATCGAGAACGTGTACAGCCTGAAGGTAATGAACAAGGATCAGCGCGATCACCTCTACGTGc             <  2:480882/150‑1 (MQ=255)
                                                                                                                                          ttGGCTACGGCCTCCGCCACCCGTTCGATTGTCGGTTTCGACGTCAGCAAATACCCCGGCTTGTACCGCGAAAACGCGTAGGGCCGTATCGAGAACGTGTACAGCCTGAAGGTAATGAACAAGGATCAGCGCGATCACCTCTACGTGCTg           <  2:336317/150‑1 (MQ=255)
                                                                                                                                          ttGGCGACGGCCTTCGCCACCCGTTCGCTGGTCGGTTTCGACGTCAGCAAAGACCGCGTGTTGTACCGCGAAAACGCCCAGGGCCGGATCGAGAACGTGTACAGCCTGAAGGTAATGAACAAGGATCAGCGCGATCACCTCTACGTGCTg           <  1:367289/150‑1 (MQ=255)
                                                                                                                                               gACGGCCTTCGCCACCCGTTCGCTGGTCGGTTTCGACGTCAGCAAAGACCGCGTGTTGTACCGCGAAAACGCCCAGGGCCGGATCGAGAACGTGTACAGCCTGAAGGTAATGAACAAGGATCAGCGCGATCACCTCTACGTGCTGGAcgc      <  2:717565/150‑1 (MQ=255)
                                                                                                                                                  ggCCTTCGCCACCCGTTCGCTGGTCGGTTTCGACGTCAGCAAAGACCGCGTGTTGTACCGCGAAAACGCCCAGGGCCGGATCGAGAACGTGTACAGCCTGAAGGTAATGAACAAGGATCAGCGCGATCACCTCTACGTGCTGGAcgccgc   >  1:962148/1‑150 (MQ=255)
                                                                                                                                                   gCCTTCGCCACCCGTTCGCTGGTCGGTTTCGACGTCAGCAAAGACCGCGTGTTGTACCGCGAAAACGCCCAGGGCCGGATCGAGAACGTGTACAGCCTGAAGGTAATGAACAAGGATCAGCGCGATCACCTCTACGTGCTGGACgccgcc  >  1:1619597/1‑150 (MQ=255)
                                                                                                                                                      |                                                                                                                                                   
CATCATGGACAAGATGAACTACCCCAAAGGGCTGATCAGCTACACCACCGAGCACAACCTGTCCGGGCAGAAGACCCACATGCTGCGCCCCCGTCTGATCGGCTATGCCGTGGTACTGCTGGTGATGATG‑‑ATCGGGTTGGCGACGGCCTTCGCCACCCGTTCGCTGGTCGGTTTCGACGTCAGCAAAGACCGCGTGTTGTACCGCGAAAACGCCCAGGGCCGGATCGAGAACGTGTACAGCCTGAAGGTAATGAACAAGGATCAGCGCGATCACCTCTACGTGCTGGACGCCGCC  >  NC_002947/4844115‑4844409

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: