Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,628,125 C→G 11.0% intergenic (+11/‑129) PP_3198 → / → PP_3199 2Fe‑2S ferredoxin/monooxygenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,628,1250CG11.0% 187.7 / 14.6 64intergenic (+11/‑129)PP_3198/PP_31992Fe‑2S ferredoxin/monooxygenase
Reads supporting (aligned to +/- strand):  ref base C (26/31);  new base G (3/4);  total (29/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.23e-01

CCACGGCGCGCGCTTCACGCCCGAAGACGGCACTTGCGTTCACGGCCCCTGCCGAGGCGAAGCCTTGGAGGCACTCGATCACAAGGAAGCGGGCGGCTGGCTGCTATTGCCGCTTGAGCAGCTCGATCGCTGACCTATCGCGCCGCGCGATAGGTCGCATCCCGCGACGCGATTGGTCATGGACCGGCGCAGGCGGGAAGCTGTGCATAGGCACAGGAAACGGCCGCCCCGCATGACGGGGCGCACAACAAGAACAAAACGTGAGACATGCCATGAATACTGT  >  NC_002947/3627982‑3628264
                                                                                                                                               |                                                                                                                                           
ccACGGCGCGCGCTTCACGCCCGAAGACGGCACTTGCGTTCACGGCCCCTGCCGAGGCGAAGCCTTGGAGGCACTCGATCACAAGGAAGCGGGCGGCTGGCTGCTATTGCCGCTTGAGCAGCTCGATCGCTGACCTATCGCGCCGCGCGa                                                                                                                                       >  2:1096886/1‑150 (MQ=255)
ccACGGCGCGCGCTTCACGCCCGAAGACGGCACTTGCGTTCACGGCCCCTGCCGAGGCGAAGCCTTGGAGGCACTCGATCACAAGGAAGCGGGCGGCTGGCTGCTATTGCCGCTTGAGCAGCTCGATCGCTGACCTATCGCGCCGCGCGa                                                                                                                                       <  1:139174/150‑1 (MQ=255)
  aCGGCGCGCGCTTCACGCCCGAAGACGGCACTTGCGTTCACGGCCCCTGCCGAGGCGAAGCCTTGGAGGCACTCGATCACAAGGAAGCGGGCGGCTGGCTGCTATTGCCGCTTGAGCAGCTCGATCGCTGACCTATCGCGCCGCGCGATa                                                                                                                                     <  1:277900/150‑1 (MQ=255)
   cGGCGCGCGCTTCACGCCCGAAGACGGCACTTGCGTTCACGGCCCCTGCCGAGGCGAAGCCTTGGAGGCACTCGATCACAAGGAAGCGGGCGGCTGGCTGCTATTGCCGCTTGAGCAGCTCGATCGCTGACCTATCGCGCGGCGCGATAg                                                                                                                                    <  2:146275/150‑1 (MQ=255)
   cGGCGCGCGCTTCACGCCCGAAGACGGCACTTGCGTTCACGGCCCCTGCCGAGGCGAAGCCTTGGAGGCACTCGATCACAAGGAAGCGGGCGGCTGGCTGCTATTGCCGCTTGAGCAGCTCGATCGCTGACCTATCGCGCGGCGCGATAg                                                                                                                                    <  2:726153/150‑1 (MQ=255)
     gcgcgcgcTTCACGCCCGAAGACGGCACTTGCGTTCACGGCCCCTGCCGAGGCGAAGCCTTGGAGGCACTCGATCACACGGATCCGTGCGGCTGGCTGCTATTGCCGCTTGAGCAGCTCGATCGCTGACCTATCGCGCCGCGCGATAGGt                                                                                                                                  >  1:74758/1‑150 (MQ=255)
           gcTTCACGCCCGAAGACGGCACTTGCGTTCACGGCCCCTGCCGAGGCGAAGCCTTGGAGGCACTCGATCACAAGGAAGCGGGCGGCTGGCTGCTATTGCCGCTTGAGCAGATCGATCGCTGACCTATCGCGCCGCGCGATAGGTCGCATc                                                                                                                            >  1:369095/1‑150 (MQ=255)
           gcTTCACGCCCGAAGACGGCACTTGCGTTCACGGCCCCTGCCGAGGCGAAGCCTTGGAGGCACTCGATCACAAGGAAGCGGGCGGCTGGCTGCTATTGCCGCGTGAGCAGCTCGATCGCTGACCTATCGCGCCGCGCGATAGGTCGCATc                                                                                                                            >  1:132276/1‑150 (MQ=255)
               cACGCCCGAAGACGGCACTTGCGTTCACGGCCCCTGCCGAGGCGAAGCCTTGGCGGGCCTCGATCACAAGGAAGCGGGCGGCTGGCTGCGATTGCGGCTTGAGCCGCTCGATTGCTTACCTATCGCGCCGCGCGATAGATCGCATCCcgc                                                                                                                        >  2:293877/1‑150 (MQ=255)
               cACGCCCGAAGACGGCACTTGCGTTCACGGCCCCTGCCGAGGCGAAGCCTTGGAGGCACTCGATCACAAGGAAGCGGGCGGCTGGCTGCTATTGCCGCTTGAGCAGCTCGATCGCTGACCTATCGCGCCGCGCGATAGGTCGCATCCcgc                                                                                                                        <  2:1048991/150‑1 (MQ=255)
                        aGACGGCACTTGCGTTCACGGCCCCTGCCGAGGCGAAGCCTTGGAGGCACTCGATCACAAGGAAGCGGGCGGCTGGCTGCTATTGCCGCTTGAGCAGCTCGATCGCTGACCTATCGCGCCGCGCGATAGGTCGCATCCCGCGACGCGAtt                                                                                                               <  2:830088/150‑1 (MQ=255)
                              cACTTGCGTTCACGGCCCCTGCCGAGGCGAAGCCTTGGAGGCACTCGATCACAAGGAAGCGGGCGGCTGGCTGCTATTGCCGCTTGAGCAGCTCGATCGCTGACCTATCGCGCCGCGCGATAGGTCGCATCCCGCGACGCGATTGGTCAt                                                                                                         <  2:369095/150‑1 (MQ=255)
                                 ttGCGTTCACGGCCCCTGCCGAGGCGAAGCCTTGGAGGCACTCGATCACAAGGAAGCGGGCGGCTGGCTGCTATTGCCGCTTGAGCAGCTCGATCGCTGACCTATCGCGCCGCGCGATAGGTCGCATccc                                                                                                                          <  2:455611/130‑1 (MQ=255)
                                 ttGCGTTCACGGCCCCTGCCGAGGCGAAGCCTTGGAGGCACTCGATCACAAGGAAGCGGGCGGCTGGCTGCTATTGCCGCTTGAGCAGCTCGATCGCTGACCTATCGCGCCGCGCGATAGGTCGCATccc                                                                                                                          >  1:455611/1‑130 (MQ=255)
                                    cGTTCACGGCCCCTGCCGAGGCGAAGCCTTGGAGGCACTCGATCACAAGGAAGCGGGCGGCTGGCTGCTATTGCCGCTTGAGCAGCTCGATCGCTGACCTATCGCGCCGCGCGATAGGTCGCATCCCGCGACGCGATTGGTCATGGACCg                                                                                                   <  2:503774/150‑1 (MQ=255)
                                     gTTCACGGCCCCTGCCGAGGCGAAGCCTTGGAGGCACTCGATCACAAGGAAGCGGGCGGCTGGCTGCTATTGCCGCTTGAGCAGCTCGATCGCTGACCTATCGCGCGGCGCGATa                                                                                                                                     <  2:1089599/115‑1 (MQ=255)
                                     gTTCACGGCCCCTGCCGAGGCGAAGCCTTGGAGGCACTCGATCACAAGGAAGCGGGCGGCTGGCTGCTATTGCCGCTTGAGCAGCTCGATCGCTGACCTATCGCGCGGCGCGATa                                                                                                                                     >  1:1089599/1‑115 (MQ=255)
                                      ttCACGGCCCCTGCCGAGGCGAAGCCTTGGAGGCACTCGATCACAAGGAAGCGGGCGGCTGGCTGCTATTGCCGCTTGAGCAGCTCGATCGCTGACCTATCGCGCGGCGCGATAGGTcag                                                                                                                               <  1:253096/120‑3 (MQ=255)
                                      ttCACGGCCCCTGCCGAGGCGAAGCCTTGGAGGCACTCGATCACAAGGAAGCGGGCGGCTGGCTGCTATTGCCGCTTGAGCAGCTCGATCGCTGACCTATCGCGCGGCGCGATAGGTcag                                                                                                                               >  2:253096/1‑118 (MQ=255)
                                         aCGGCCCCTGCCGAGGCGAAGCCTTGGAGGCACTCGATCACAAGGACGCGGGCGGCTGGCTGCTATTGCCGCTTGAGCAGCTCGATCGCTGACCTATCGCGCCGCGCGATAGGTCg                                                                                                                                >  1:536723/1‑116 (MQ=255)
                                         aCGGCCCCTGCCGAGGCGAAGCCTTGGAGGCACTCGATCACAAGGAAGCGGGCGGCTGGCTGCTATTGCCGCTTGAGCAGCTCGATCGCTGACCTATCGCGCCGCGCGATAGGTCg                                                                                                                                <  2:536723/116‑1 (MQ=255)
                                         aCGGCCCCTGCCGAGGCGAAGCCTTGGAGGCACTCGATCACAAGGAAGCGGGCGGCTGGCTGCTATTGCCGCTTGAGCAGCTCGATCGCTGACCTATCGCGCCGCGCGATAGGTCGCATCCCGCGACGCGATTGGTCATGGACCGGCGCa                                                                                              >  2:194267/1‑150 (MQ=255)
                                           ggCCCCTGCCGAGGCGAAGCCTTGGAGGCACTCGATCACAAGGAAGCGGGCGGCTGGCTGCTATTGCCGCTTGAGCAGCTCGATCGCTGACCTATCGCGCCGCGCGATAGGTCGCATCCCGCGACGCGATTGGTCATGGACCGGCGCAgg                                                                                            >  1:407286/1‑150 (MQ=255)
                                             ccccTGCCGAGGCGAAGTCTTGGCGGCTGTCGGTCCCCAGGTGTCAGTCAGTGGACGGCTGTTTCCGTTTTAGCAGCTCGATCGCTGACCTATCGCGCCGCGCGATCGGTCGCATCCCGCGACGCGATTGGTCATGGACCGGCGCAGGCg                                                                                          <  2:865419/150‑1 (MQ=255)
                                             ccccTGCCGAGGCGAAGCCTTGGAGGCACTCGATCACAAGGAAGCGGGCGGCTGGCTGCTATTGCCGCTTGAGCAGCTCGATCGCTGACCTATCGCGCCGCGCGATAGGTCGCATCCCGCGACGCGATTGGTCATGGACCGGCGCAGGCg                                                                                          >  1:397157/1‑150 (MQ=255)
                                              cccTGCCGAGGCGAAGCCTTGGAGGCACTCGATCACAAGGAAGCGGGCGGCTGGCTGCTATTGCCGCTTGAGCAGCTCGATCGCTGACCTATCGCGCCGCGCGATAGGTCGCATCCCGCGACGCGATTGGTCATGGACCGGCGCAGGCgg                                                                                         <  1:769381/150‑1 (MQ=255)
                                               ccTGCCGAGGCGAAGCCTTGGAGGCACTCGATCACAAGGAAGCGGGCGGCTGGCTGCTATTGCCGCTTGAGCAGCTCGATCGCTGACCTATCGCGCCGCGCGATAGGTCGCATCCCGCGACGCGATTGGTCATGGACCGGCGCAGGCggg                                                                                        >  1:689899/1‑150 (MQ=255)
                                                 tGCCGAGGCGAAGCCTTGGAGGCACTCGATCACAAGGAAGCGGGCGGCTGGCTGCTATTGCCGCTTGAGCAGCTCGATCGCTGACCTATCGCGCCGCGCGATAGGTCGCATCCCGCGACGCGATTGGTCATGGACCGGCGCAGGCGGGaa                                                                                      <  2:951321/150‑1 (MQ=255)
                                                             gCCTTGGAGGCACTAGATCACAAGGAAGCGGGCGGCTGGCTGTTATTGCCGCTTGAGCAGCTCGATCGCTGACCTATCGCGCCGCGCGATAGGTCGCATCCCGCGACGCGATTGGTCATGGACCGGCGCAGGCGGGAAGCTGTGCATAgg                                                                          <  2:1037720/150‑1 (MQ=255)
                                                                 tGGAGGCACTCGATCACAAGGAAGCGGGCGGCTGGCTGCTATTGCCGCTTGAGCAGCTCGATCGCTGACCTATCGCGCCGCGCGATAGGTCGCATCCCGCGACGCGATTGGTCATGGACCGGCGCAGGCGGGAAGCTGTGCATAGGcaca                                                                      <  2:397157/150‑1 (MQ=255)
                                                                 tGGAGGCACTCGATCACAAGGAAGCGGGCGGCTGGCTGCTATTGCCGCTTGAGCAGCTCGATCGCTGACCTATCGCGCCGCGCGATAGGTCGCATCCCGCGACGCGATTGGTCATGGACCGGCGCAGGCGGGAAGCTGTGCATAGGcaca                                                                      <  2:642856/150‑1 (MQ=255)
                                                                     ggCACTCGATCACAAGGAAGCGGGCGGCTGGCTGCTATTGCCGCTTGAGCAGCTCGATCGCTGACCTATCGCGCCGCGCGCTAGGTCGCATCCCGCGACGCGATTGGTCATGGACCGGCGCAGGCGGCAAGCTGTGCATAGGCACCGGaa                                                                  >  2:797/1‑150 (MQ=255)
                                                                      gCACTCGATCACAAGGAAGCGGGCGGCTGGCTGCTATTGCCGCTTGAGCAGCTCGATCGCTGACCTATCGCGCCGCGCGATAGGTc                                                                                                                                 <  2:4843/86‑1 (MQ=255)
                                                                      gCACTCGATCACAAGGAAGCGGGCGGCTGGCTGCTATTGCCGCTTGAGCAGCTCGATCGCTGACCTATCGCGCCGCGCGATAGGTc                                                                                                                                 >  1:4843/1‑86 (MQ=255)
                                                                         cTCGATCACAAGGAAGCGGGCGGCTGGCTGCTATTGCCGCTTGAGCAGCTCGATCGCTGACCTATCGCGCCGCGCGATAGGTCGCATCCCGCGACGCGATTGGTCATGGACCGGCGCAGGCGGGAAGCTGTGCATAGGCACAGGAAACgg                                                              <  2:689899/150‑1 (MQ=255)
                                                                         cTCGATCACAAGGAAGCGGGCGGCTGGCTGCTATTGCCGCTTGAGCAGCTCGATCGCTGACATATCGCGCCGCGCGATAGGTCGAATCCCGCGACGCGATTGGTCATGGACCGGCGCAGGCGGGAAGCTGTGCATAGGCACAGGAAACgg                                                              <  2:1064336/150‑1 (MQ=255)
                                                                                acaAGGAAGCGGGCGGCTGGCTGCTATTGCCGCTTGAGCAGCTCGATCGCTGACCTATCGCGCCGCGCGATAGGTCGCATCCCGCGACGCGATTGGTCATGGACCGGCGCAGGCGGGAAGCtgt                                                                                 <  1:747905/124‑1 (MQ=255)
                                                                                acaAGGAAGCGGGCGGCTGGCTGCTATTGCCGCTTGAGCAGCTCGATCGCTGACCTATCGCGCCGCGCGATAGGTCGCATCCCGCGACGCGATTGGTCATGGACCGGCGCAGGCGGGAAGCtgt                                                                                 >  2:747905/1‑124 (MQ=255)
                                                                                  aaGGAAGCGGGCGGCTGGCTGCTATTGCCGCTTGAGCAGCTCGATCGCTGACCTATCGCGCCGCGCGATAGGTCGCATCCCGCGACGCGATTGGTCATGGACCGGCGCAGGCGGGAAGCTGTGCATAGGCACAGGAAACGGCCGCCCCGc                                                     >  1:1187318/1‑150 (MQ=255)
                                                                                    ggAAGCGGGCGGCTGGCTTCTATTGCCGCTTGAGCAGCTCGATCGCTGACCTATCGCGCCGCGCGATATGTCGCATCCCGCGACGCGATTGGTCATGGACCGGCGCAGGCGGGAAGCTGTGCATAGGCACAg                                                                     <  2:260465/132‑1 (MQ=255)
                                                                                    ggAAGCGGGCGGCTGGCTGCTATTGCCGCTTGAGCAGCTCGATCGCTGACCTATCGCGCCGCGCGATATGTCGCATCCCGCGACGCGATTGGTCATGGACCGGCGCAGGCGGGAAGCTGTGCATAGGCACAg                                                                     >  1:260465/1‑132 (MQ=255)
                                                                                    ggAAGCGGGCGGCTGGCTGCTATTGCCGCTTGAGCAGCTCGATCGCTGACCTATCGCGCCGCGCGATAGGTCGCATCCCGCGACGCGATTGGTCATGGACCGGCGCAGGCGGGAAGCTGTGCATAGGCACAGGAAACGGCCGCCCCGCAt                                                   >  1:591906/1‑150 (MQ=255)
                                                                                    ggAAGCGGGCGGCTGGCTGCTATTGCCGCTTGAGCAGCTCGATCGCTGACCTATCGCGCCGCGCGATAGGTCGCATCCCGCGACGCGATTGGTCATGGACCGGCGCAGGCGGGAAGCTGTGCATAGGCACAGGAAACGGCCGCCCCGCAt                                                   >  2:949875/1‑150 (MQ=255)
                                                                                      aaGCGGGCGGCTGGCTGCTATTGCCGCTTGAGCAGCTCGATCGCTGACCTATCGCGCCGCGCGATAGGTCGCATCCCGCGACGCGATTGGTCATGGACCGGCGCAGGCGGGAAGCTGTGCATAGGCACAGGAAACGGCCGCCCCGCATGa                                                 <  2:869366/150‑1 (MQ=255)
                                                                                            gcggcTGGCTGCTATTGCCGTTTGAGAAGCTCGATCGCTGACCTATCGCGCCGCGCGATAGGTCGCATCCCGCGACGCGATTGGTCATGGACCGGCGCAGGCNGGAAGCTGTGCATAGGCACAGGAAACGGCCGCCCCGCATGACGGGgc                                           <  1:797/150‑1 (MQ=255)
                                                                                            gcggcTGGCTGCTATTGCCGCTTGAGCAGCTCGATCGCTGACCTATCGCGCCGCGCGATAGGTCGCATCCCGCGACGCGATTGGTCATGGACCGGCGCAGGCGGGAAGCTGTGCATAGGCACAGGAAACGGCCGCCCCGCATGACGGGgc                                           <  2:407286/150‑1 (MQ=255)
                                                                                                     tgctATTGCCGCTTGAGCAGCTCGATCGCTGACCTATCGCGCCGCGCGATAGGTCGCATCCCGCGACGCGATTGGTCATGGACCGGCGCAGGCGGGAAGCTGTGCATAGGCACAGGAAACGGCCGCCCCGCATGACGGGGCGCAcaacaa                                  <  2:166381/150‑1 (MQ=255)
                                                                                                        tATTGCCGCTTGAGCAGCTCGATCGCTGACCTATCGCGCCGCGCGATAGGTCGCATCCCGCGACGCGATTGGTCATGGACCGGCGCAGGCGGGAAGCTGTGCATAGGCACAGGAAACGGCCGCCCCGCATGACGGGGCGCACAACAAGa                                <  1:458907/149‑1 (MQ=255)
                                                                                                        tATTGCCGCTTGAGCAGCTCGATCGCTGACCTATCGCGCCGCGCGATAGGTCGCATCCCGCGACGCGATTGGTCATGGACCGGCGCAGGCGGGAAGCTGTGCATAGGCACAGGAAACGGCCGCCCCGCATGACGGGGCGCACAACAAGa                                >  2:458907/1‑149 (MQ=255)
                                                                                                               gCTTGAGCAGCTCGATCGCTGACCTATCGCGCCGCGCGATAGGTCGCATCCCGCGACGCGATTGGTCATg                                                                                                        >  1:1019450/1‑70 (MQ=255)
                                                                                                               gCTTGAGCAGCTCGATCGCTGACCTATCGCGCCGCGCGATAGGTCGCATCCCGCGACGCGATTGGTCATg                                                                                                        <  2:1019450/70‑1 (MQ=255)
                                                                                                                   gAGCAGCTCGATCGCTGACCTATCGCGCCGCGCGATAGGTCGCATCCCGCGACGCGATTGGTCATGGACCGGCGCAGGCGGGAAGCTGTGCATAGGCACAGGAAACGGCCGCCCCGCATGACGGGGCGCACAACAAGAACAAAACGTgag                    <  1:1096886/150‑1 (MQ=255)
                                                                                                                     gcagcTCGATCGCTGACCTATCGCGCCGCGCGATAGGTCGCATCCCGCGACGCGATTGGTCATGGACCGGCGCAGGCGGGAAGCTGTGCATAGGCACAGGAAACGGCCGCCCCGCATGACGGGGCGCACAACAAGAACAAAACGTGAGAc                  >  2:737299/1‑150 (MQ=255)
                                                                                                                        gcTCGATCGCTGACCTATCGCGCCGCGCGATAGGTCGCATCCCGCGACGCGATTGGTCATGGACCGGCGCAGGCGGGAAGCTGTGCATAGGCACAGGAAACGGCCGCCCCGCATGACGGGGCGCACAACAAGAACAAAACGTGAGACATg               >  2:1180493/1‑150 (MQ=255)
                                                                                                                          tCGATCGCTGACCTATCGCGCCGCGCGATAGGTCGCATCCCGCGACGCGATTGGTCATGGACCGGCGCAGGCGGGAAGCTGTGCATAGGCACAGGAAACGGCCGCCCCGCATGACGGGGCGCACAACAAGAACAAAACGTGAGACATGcc             >  1:240624/1‑150 (MQ=255)
                                                                                                                            gATCGCTGACCTATCGCGCCGCGCGATAGGTCGCATCCCGCGACGCGATTGGTCATGGACCGGCGCAGGCGGGaa                                                                                      >  1:1046218/1‑75 (MQ=255)
                                                                                                                            gATCGCTGACCTATCGCGCCGCGCGATAGGTCGCATCCCGCGACGCGATTGGTCATGGACCGGCGCAGGCGGGaa                                                                                      <  2:1046218/75‑1 (MQ=255)
                                                                                                                            gATCGCTGACCTATCGCGCCGCGCGATAGGTCGCATCCCGCGACGCGATTGGGCATGGACCGGCGCAGGCGGGAAGCTGTGCATAGGCACAGGAAACGGCCGCCCCGCATGACGGGGCGCACAACAAGAACAAAACGTGAGACATGCCAt           <  1:334781/150‑1 (MQ=255)
                                                                                                                                  tGACCTATCGCGCCGCGCGATAGGTCGCATCCCGCGACGCGATTGGTCATGGACCGGCGCAGGCGGGAAGCTGTGCATAGGCACAGGAAACGGCCGCCCCGCATGACGGGGCGCACAACAAGAACAAAACGTGAGACATGCCATGAATAc     <  2:240624/150‑1 (MQ=255)
                                                                                                                                  tGACCTATCGCGCCGCGCGATAGGTCGCATCCCGCGACGCGATTGGTCATGGACCGGCGCAGGCGGGAAGCTGTGCATAGGCACAGGAAACGGCCGCCCCGCATGACGGGGCGCACAACAAGAACAAAACGTGAGACATGCCATGAATAc     >  2:1021801/1‑150 (MQ=255)
                                                                                                                                    aCCTATCGCGCCGCGCGATAGGTCGCATCCCGCGACGCGATTGGTCATGGACCGGCGCAGGCGGGAAGCTGTGCATAGGCACAGGAAACGGCCGCCCCGCATGACGGGGCGCACAACAAGAACAAAACGTGAGACATGCCATGAATACtg   <  1:737299/150‑1 (MQ=255)
                                                                                                                                     ccTATCGCGCGGCGCGATAGGTCGCATCCCGCGACGCGATTGGTCATGGACCGGCGCAGGCGGGAAGCTGTGCATAGGCACAGGAAACGGCCGCCCCGCATGACGGGGCGCACAACAAGAACAAAACGTGAGACATGCCATGAATACtgt  >  1:728798/1‑150 (MQ=255)
                                                                                                                                       tATCGCGCCGCGCGATAGGTCGCATCCCGCGACGCGATTGGTCATGGACCGGCGCAGGCGGGAAGCTGTGCATAGGCACAGGAAACGGCCGCCCCGCATGACGGGGCGCACAACAAGa                                <  1:176188/118‑1 (MQ=255)
                                                                                                                                       tATCGCGCCGCGCGATAGGTCGCATCCCGCGACGCGATTGGTCATGGACCGGCGCAGGCGGGAAGCTGTGCATAGGCACAGGAAACGGCCGCCCCGCATGACGGGGCGCACAACAAGa                                >  2:176188/1‑118 (MQ=255)
                                                                                                                                               |                                                                                                                                           
CCACGGCGCGCGCTTCACGCCCGAAGACGGCACTTGCGTTCACGGCCCCTGCCGAGGCGAAGCCTTGGAGGCACTCGATCACAAGGAAGCGGGCGGCTGGCTGCTATTGCCGCTTGAGCAGCTCGATCGCTGACCTATCGCGCCGCGCGATAGGTCGCATCCCGCGACGCGATTGGTCATGGACCGGCGCAGGCGGGAAGCTGTGCATAGGCACAGGAAACGGCCGCCCCGCATGACGGGGCGCACAACAAGAACAAAACGTGAGACATGCCATGAATACTGT  >  NC_002947/3627982‑3628264

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