Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,844,258 C→A 5.8% T354K (ACG→AAG)  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,2580CA 5.8% 275.3 / 3.1 86T354K (ACG→AAG) PP_4259iron‑sulfur cluster‑binding protein
Reads supporting (aligned to +/- strand):  ref base C (39/41);  new base A (3/2);  total (43/43)
Fisher's exact test for biased strand distribution p-value = 6.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.77e-01

CGACAGCATCATGGACAAGATGAACTACCCCAAAGGGCTGATCAGCTACACCACCGAGCACAACCTGTCCGGGCAGAAGACCCACATGCTGCGCCCCCGTCTGATCGGCTATGCCGTGGTACTGCTGGTGATGATG‑‑ATCGGGTTGGCGACGGCCTTCGCCACCCGTTCGCTGGTCGGTTTCGACGTCAGCAAAGACCGCGTGTTGTACCGCGAAAACGCCCAGGGCCGGATCGAGAACGTGTACAGCCTGAAGGTAATGAACAAGGATCAGCGCGATCACCTCTACGTGCTGGACGC  >  NC_002947/4844109‑4844405
                                                                                                                                                       |                                                                                                                                                    
cGACAGCATCATGGACAAGATGAACTACCCCAAAGGGCTGATCAGCTACACCACCGAGCACAACCTGTCCGGGCAGAAGACCCACATGCTGCGCCCCCGTCTGATCGGCTATGCCGTGGTACTGCTGGTGATGATG‑‑ATCGGGTTGGCGAc                                                                                                                                                      <  2:1738037/150‑1 (MQ=255)
 gACAGCATCATGGACAAGATGAACTACCCCAAAGGGCTGATCAGCTACACCACCGAGCACAACCTGTCCGGGCAGAAGACCCACATGCTGCGCCCCCGTCTGATCGGCTATGCCGTGGTACTGCTGGTGATGATG‑‑ATCGGGTTGGCGACg                                                                                                                                                     <  2:980370/150‑1 (MQ=255)
      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                                                                                                                        >  2:1760538/1‑150 (MQ=255)
                             ccAAAGGGCTGATCAGCTACACCACCGAGCACAACCTGTCCGGGCAGAAGACCCACATGCTGCGCCCCCGTCTGATCGGCTATGCCGTGGTACTGCTGGTGATGATG‑‑ATCGGGTTGGCGACGGCCTTCGCCACCCGTTCGCTGGTCGGtt                                                                                                                        >  1:499210/1‑150 (MQ=255)
                                    gCTGATCAGCTACACCACCGAGCACAACCTGTCCGGGCAGAAGACCCACATGCTGCGCCCCCGTCTGATCGGCTATGCCGTGGTACTGCTGGTGATGATG‑‑ATCGGGTTGGCGACGGCCTTCGCCACCCGTTCGCTGGTCGGTTTCGACGt                                                                                                                 <  2:1048174/150‑1 (MQ=255)
                                        aTCAGCTACACCACCGAGCACAACCTGTCCGGGCAGAAGACCCACATGCTGCGCCCCCGTCTGATCGGCTATGCCGTGGTACTGCTGGTGATGATG‑‑ATCGGGTTGGCGACGGCCTTCGCCACCCGTTc                                                                                                                                   <  2:1454393/128‑1 (MQ=255)
                                        aTCAGCTACACCACCGAGCACAACCTGTCCGGGCAGAAGACCCACATGCTGCGCCCCCGTCTGATCGGCTATGCCGTGGTACTGCTGGTGATGATG‑‑ATCGGGTTGGCGACGGCCTTCGCCACCCGTTc                                                                                                                                   >  1:1454393/1‑128 (MQ=255)
                                        aTCAGCTACACCACCGAGCACAACCTGTCCGGGCAGAAGACCCACATGCTGCGCCCCCGTCTGATCGGCTATGCCGTGGTACTGCTGGTGATGATG‑‑ATCGGGGTGGCGAAGGCCGTCGCCAacccga                                                                                                                                    <  1:1620370/127‑4 (MQ=255)
                                        aTCAGCTACACCACCGAGCACAACCTGTCCGGGCAGAAGACCCACATGCTGCGCCCCCGTCTGATCGGCTATGCCGTGGTACTGCTGGTGATGATG‑‑ATCGGGGTGGCGAAGGCCGTCGCCAacccga                                                                                                                                    >  2:1620370/1‑124 (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                                                                                                                                               >  1:1343595/1‑106 (MQ=255)
                                                   caccGAGCACAACCTGTCCGGGCAGAAGACCCACATGCTGCGCCCCCGTCTGATCGGCTATGCCGTGGTACTGCTGGTGATGATG‑‑ATCGGGTTGGCGACGGCCTTc                                                                                                                                               <  2:1343595/106‑1 (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‑‑ATCGGGTTGGCGAc                                                                                                                                                      >  2:892517/1‑86 (MQ=255)
                                                                cTGTCCGGGCAGAAGACCCACATGCTGCGCCCCCGTCTGATCGGCTATGCCGTGGTACTGCTGGTGATGATG‑‑ATCGGGTTGGCGAc                                                                                                                                                      <  1:892517/86‑1 (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                                                                                                    <  2:704535/118‑1 (MQ=255)
                                                                                 ccACATGCTGCGCCCCCGTCTGATCGGCTATGCCGTGGTACTGCTGGTGATGATG‑‑ATCGGGTTGGCGACGGCCTTCGCCACCCGTTCGCTGGTCGGTTTCGACGTCAGCAAAGACcgc                                                                                                    >  1:704535/1‑118 (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)
                                                                                 ccACATGCTGCGCCCCCGTCTGATCGGCTATGCCGTGGTACTGCTGGTGATGATG‑‑ATCGGGTTGGCGACGGCCTTCGCCACCCGTTCGCTGGTCGGTTTCGACGTCAGCAAAGACCGCGTGTTGTACCGCGAAAACGCCCAGGGCCGGAt                                                                    >  1:557776/1‑150 (MQ=255)
                                                                                   acaTGCTGCGCCCCCGTCTGATCGGCTATGCCGTGGTACTGCTGGTGATGATG‑‑ATCGGGTTGGCGACGGCCTTCGCCACCCGTTCGCTGGTCGGTTTCGACGTCAg                                                                                                              <  1:956652/106‑1 (MQ=255)
                                                                                   acaTGCTGCGCCCCCGTCTGATCGGCTATGCCGTGGTACTGCTGGTGATGATG‑‑ATCGGGTTGGCGACGGCCTTCGCCACCCGTTCGCTGGTCGGTTTCGACGTCAg                                                                                                              >  2:956652/1‑106 (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                                                                                                 <  2:537415/100‑1 (MQ=255)
                                                                                                      gATCGGCTATGCCGTGGTACTGCTGGTGATGATG‑‑ATCGGGTTGGCGACGGCCTTCGCCACCCGTTCGCTGGTCGGTTTCGACGTCAGCAAAGACCGCgtg                                                                                                 >  1:537415/1‑100 (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                                              <  1:842592/150‑1 (MQ=255)
                                                                                                       aTCGGCTATGCCGTGGTACTGCTGGTGATGATG‑‑ATCGGGTTGGCGACGGCCTTCGCCACCCGTTCGCTGGTCGGTTTCGACGTCAGCAAAGACCGCGTGTTGTACCGCGAAAACGCCCAGGGCCGGATCGAGAACGTGTACAGCCTGAAg                                              <  2:499210/150‑1 (MQ=255)
                                                                                                       aTCGGCTATGCCGTGGTACTGCTGGTGATGATG‑‑ATCGGGTTGGCGACGGCCTTCGCCACCCGTTCGCTGGTCGGTTTCGACGTCAGCAAAGACCGCGTGTTGTACCGCGAAAACGCCCAGGGCCGGATCGAGAACGTGTACAGCCTGAAg                                              >  2:1620150/1‑150 (MQ=255)
                                                                                                         cGGCTATGCCGTGGTACTGCTGGTGATGATG‑‑ATCGGGTTGGCGACGGCCTTCGCCACCCGTTCGCTGGTCGGTTTCGACGTCAGCAAAGACCGCGTGTTGTACCGCGAAAACGCCCAGGGCCGGATCGAGAACGTGTACAGCCTGAAGGt                                            <  1:1620150/150‑1 (MQ=255)
                                                                                                             tATGCCGTGGTACTGCTGGTGATGATG‑‑ATCGGGTTGGCGACGGCCTTCGCCACCCGTTCGCTGGTCGGTTTCGACGTCAGCAAAGACCGCGTGTTGTACCGCGAAAACGCCCAGGGCCGGATCGAGAACGTGTACAGCCTGAAGGTAATg                                        <  2:557776/150‑1 (MQ=255)
                                                                                                               tGCCGTGGTACTGCTGGTGATGATG‑‑ATCGGGTTGGCGACGGCCTTCGCCACCCGTTCGCTGGTCGGTTTCGACGTCAGCAAAGACCGCGTGTTGTACCGCGAAAACGCCCAGGGCCGGATCGAGAACGTGTACAGCCTGAAGGTAATGaa                                      <  1:819647/150‑1 (MQ=255)
                                                                                                               tGCCGTGGTACTGCTGGTGATGATG‑‑ATCGGGTTGGCGACGGCCTTCGCCACCCGTTCGCTGGTCGGTTTCGACGTCAGCAAAGACCGCGTGTTGTACCGCGAAAACGCCCAGGGCCGGATCGAGAACGTGTACAGCCTGAAGGTAATGaa                                      <  1:1760538/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)
                                                                                                                                                       |                                                                                                                                                    
CGACAGCATCATGGACAAGATGAACTACCCCAAAGGGCTGATCAGCTACACCACCGAGCACAACCTGTCCGGGCAGAAGACCCACATGCTGCGCCCCCGTCTGATCGGCTATGCCGTGGTACTGCTGGTGATGATG‑‑ATCGGGTTGGCGACGGCCTTCGCCACCCGTTCGCTGGTCGGTTTCGACGTCAGCAAAGACCGCGTGTTGTACCGCGAAAACGCCCAGGGCCGGATCGAGAACGTGTACAGCCTGAAGGTAATGAACAAGGATCAGCGCGATCACCTCTACGTGCTGGACGC  >  NC_002947/4844109‑4844405

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