Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,960,435 C→A 5.7% intergenic (‑7/‑322) mnxG ← / → PP_3491 manganese‑oxidizing multicopper oxidase/enoly‑coenzyme A hydratase/isomerase family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,960,4350CA 5.7% 225.0 / 2.9 70intergenic (‑7/‑322)mnxG/PP_3491manganese‑oxidizing multicopper oxidase/enoly‑coenzyme A hydratase/isomerase family protein
Reads supporting (aligned to +/- strand):  ref base C (30/36);  new base A (2/2);  total (32/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.12e-01

ATCAGTGGCTGGTCCAGTGCCACCACGTTGGCGACCAGGGTGCGCTGGCACTGCACAGCCGCCTCGCTGGTTGCTTCCCAGCCAAGCAAGGCCAGTACAAGTGATGACAGGATGGGGCCTTTGAGCGGCGTCGACATGCTGCACCTCCGGGACACAGGATCTCTGCGGGGAGCTGGAGCAAGGTGTGTGCCAGAAGAAAAGTTCAATATTAATCAATAACTTGTACGATCAACCGCAACGGATTGGGGAGGATCCCCCAGCAGGTTCCCCCGAATCGGGGGACG  >  NC_002947/3960292‑3960575
                                                                                                                                               |                                                                                                                                            
aTCAGTGGCTGGTCCAGTGCCACCACGTTGGCGACCAGGGTGCGCTGGCACTGCACAGCCGCCTCGCTGGTTGCTTCCCAGCCAAGCAAGGCCAGTACAAGTGATGACAGGATGGGGCCTTTGAGCGGCGTCGACATGCTGCACCTCCgg                                                                                                                                        >  2:72334/1‑150 (MQ=255)
    gTGGCTGGTCCAGTGCCACCACGTTGGCGACCAGGGTGCGCTGGCACTGCACAGCCGCCTCGCTGGTTGCTTCCCAGCCAAGCAAGGCCAGTACAAGTGATGACAGGATGGGGCCTTTGAGCGGCGTCGACATGCTGCACCTCCGGGaca                                                                                                                                    >  1:1073037/1‑150 (MQ=255)
      ggCTGTTCCAGTGCCACCACGTTGGCGACCAGGGTTCGCTGGCACTGCACAGCCGCCTCGCTGGTTGCTTCCCAGCCAAGCAAGGCCAGTACAAGTGATGACAGGATGCGTCCTTTGAGCGGCGTCGACATGCTGCACCTCCGGGAcaca                                                                                                                                  <  1:540732/150‑1 (MQ=255)
       gCTGGTCCAGTGCCACCACGTTGGCGACCAGGGTGCGCTGGCACTGCACAGCCGCCTCGCTGGTTGCTTCCCAGCCAAGCAAGGCCAGTACAAGTGATGACAGGATGGGGCCTTTGAGCGGCGTCGACATGCTGCACCTCCGGGACACAg                                                                                                                                 >  1:1339473/1‑150 (MQ=255)
        cTGGTCCAGTGCCACCACGTTGGCGACCAGGGTGCGCTGGCACTGCACAGCCGCCTCGCTGGTTGCTTCCCAGCCAAGCAAGGCCAGTACAAGTGATGACAGGATGGGGCCTTTGAGCGGCGTCGACATGCTGCACCTCCGGGACACAgg                                                                                                                                <  1:72334/150‑1 (MQ=255)
           gTCCAGTGCCACCACGTTGGCGACCAGGGTGCGCTGGCACTGCACAGCCGCCTCGCTGGTTGCTTCCCAGCCAAGCAAGGCCAGTACAAGTGATGACAGGATGGGGCCTTTGAGCGGCGTCGACATGCTGCACCTCCGGGACACAGGAtc                                                                                                                             >  2:1163444/1‑150 (MQ=255)
            tCCAGTGCCACCACGTTGGCGACCAGGGTGCGCTGGCACTGCACAGCCGCCTCGCTGGTTGCTTCCCAGCCAAGCAAGGCCAGTACAAGTGATGACAGGATGGGGCCTTTGAGCGGCGTCGACATGCTGCACCTCCGGGACACAGGAtct                                                                                                                            >  2:910323/1‑150 (MQ=255)
                gTGCCACCACGTTGGCGACCAGGGTGCGCTGGCACTGCACAGCCGCCTCGCTGGTTGCTTCCCAGCCAAGCAAGGCCAGTACAAGTGATGACAGGATGGGGCCTTTGAGCGGCGTCGACATGCTGCACCTCCGGGACACAGGATCTCTGc                                                                                                                        <  1:860090/150‑1 (MQ=255)
                 tGCCACCACGTTGGCGACCAGGGTGCGCTGGCACTGCACAGCCGCCTCGATGGTTGCTTCCCAGCCAAGCAAGGCCAGTACAAGTGATGACAGGATGGGGCGTTTGAGCGGCGTCGACATGCTGCACCTCCGGGACACAGGATCTCTGCg                                                                                                                       >  2:952913/1‑150 (MQ=255)
                  gCCACCACGTTGGCGACCAGGGTGCGCTGGCACTGCACAGCCGCCTCGCTGGTTGCTTCCCAGCCAAGCAAGGCCAGTACAAGTGATGACAGGATGGGGCCTTTGAGCGGCGTCGACATGCTGCACCTCCGGGACACAGGATCTCTGCgg                                                                                                                      <  1:477839/150‑1 (MQ=255)
                     accacGTTGGCGACCAGGGTGCGCTGGCACTGCACAGCCGCCTCGCTGGTTGCTTCCCAGCCAAGCAAGGCCAGTACAAGTGATGACAGGATGGGGCCTTTGAGCGGCGTCGACATGCTGCACCTCCAGGACACACGATCTCTTCggtca                                                                                                                   <  2:1022635/150‑4 (MQ=255)
                      ccacGTTGGCGACCAGGGTGCGCTGGCACTGCACAGCCGCCTCGCTGGTTGCTTCCCAGCCCAGCAAGGCAAATGATAGTTGGGACAGGATGGGGCCTTTGAGCGGCGTCGACATGCTGCACCTCCGGGACACAGGATCTCTGCGGGGAg                                                                                                                  >  1:883399/1‑150 (MQ=255)
                      ccacGTTGGCGACCAGGGTGCGCTGGCACTGCACAGCCGCCTCGCTGGTTGCTTCCCAGCCAAGCAAGGCCAGTACAAGTGATGACAGGATGGGGCCTTTGAGCGGCGTCGACATGCTGCACCTCCGGGACACAGGATCTCTGCGGGGAg                                                                                                                  >  1:925860/1‑150 (MQ=255)
                         cGTTGGCGACCAGGGTGCGCTGGCACTGCACAGCCGCCTCGCTGGTTGCTTCCCAGCCAAGCAAGGCCAGTACAAGTGATGACAGGATGGGGCCTTTGAGCGGCGTCGACATGCTGCACCTCCGGGACACAGGATCTCTGCGGGGAGCTg                                                                                                               <  1:407273/150‑1 (MQ=255)
                           ttGGCGACCAGGGTGCGCTGGCACTGCACAGCCGCCTCGCTGGTTGCTTCCCAGCCAAGCAAGGCCAGTACAAGTGATGACAGGATGGGGCCTTTGAGCGGCGTCGACATGCTGCACCTcc                                                                                                                                          <  1:102880/121‑1 (MQ=255)
                           ttGGCGACCAGGGTGCGCTGGCACTGCACAGCCGCCTCGCTGGTTGCTTCCCAGCCAAGCAAGGCCAGTACAAGTGATGACAGGATGGGGCCTTTGAGCGGCGTCGACATGCTGCACCTcc                                                                                                                                          >  2:102880/1‑121 (MQ=255)
                                gACCAGGGTGCGCTGGCACTGCACAGCCGCCTCGCTGGTTGCTTCCCAGCCGAGCAACGCCATTACAGGTGATGACAGGATGGGGCCTTTGAGCGGCGTCGACATGCTGCACCTCCGGGACACAGGATCTCTGCGGGGAGCTGGAGCAAg                                                                                                        <  1:1236824/150‑1 (MQ=255)
                                gACCAGGGTGCGCTGGCACTGCACAGCCGCCTCGCTGGTTGCTTCCCAGCCAAGCAAGGCCAGTACAAGTGATGACAGGATGGGGCCTTTGAGCGGCGTCGACATGCTGCACCTCCGGGACACAGGATCTCTGCGGGGAGCTGGAGCAAg                                                                                                        <  1:1114140/150‑1 (MQ=255)
                                   cAGGGTGCGCTGGCACTGCACAGCCGCCTCGCTGGTTGCTTCCCAGCCAAGCAAGGCCAGTACAAGTGATGACAGGATGGGGCCTTTGAGCGGCGTCGACATGCTGCACCTCCGGGACACAGGATCTCTGCGGGGAGCTGGAGCAAGgtg                                                                                                     >  1:274719/1‑150 (MQ=255)
                                      ggTGCGCTGGCACTGCACAGCCGCCTCGATGGTTGCTTCCCAGCCAAGCAAGGCCAGTACAAGTGATGACAGGATGGGGCCTTTGAGCGGCGTCGACATGCTGCACCTCCGGGACACAGGATCTCTGCGGGGAGCTGGAGCAAGgtgtgt                                                                                                  <  1:952913/150‑1 (MQ=255)
                                            cTGGCACTGCACAGCCGCCTCGCTGGTTGCTTCCCAGCCAAGCAAGGCCAGTACAAGTGATGACAGGATGGGGCCTTTGAGCGGCGTCGACATGCTGCACCTCCGGGACACAGGATCTCTGCGGGGAGCTGGAGCAAGGTGTGTGCCaga                                                                                            <  1:459962/150‑1 (MQ=255)
                                                 aCTGCACAGCCGCCTCGCTGGTTGCTTCCCAGCCAAGCAAGGCCAGTACAAGTGATGACAGGATGGGGCCTTTGAGCGGCGTCGACATGCTGCACCTCCGGGACACAGGATCTCTGCGGGGAGCTGGAGCAAGGTGTGTGCCAGAAGaaa                                                                                       >  2:680196/1‑150 (MQ=255)
                                                  cTGCACAGCCGCCTCGCTGGTTGCTTCCCAGCCAAGCAAGGCCAGTACAAGTGATGACAGGATGGGGCCTTTGAGCGGCGTCGACATGCTGCACCTCCGGGa                                                                                                                                      >  2:499613/1‑102 (MQ=255)
                                                  cTGCACAGCCGCCTCGCTGGTTGCTTCCCAGCCAAGCAAGGCCAGTACAAGTGATGACAGGATGGGGCCTTTGAGCGGCGTCGACATGCTGCACCTCCGGGa                                                                                                                                      <  1:499613/102‑1 (MQ=255)
                                                  cTGCACAGCCGCCTCGCTGGTTGCTTCCCAGCCAAGCAAGGCCAGTACAAGTGATGACAGGATGGGGCCTTTGAGCGGCGTCGACATGCTGCACCTCCGGGACACAGGATCTCTGCGGGGAGCTGGAGCAAGGTGTGTGCCAGAAGaaaa                                                                                      <  2:1339473/150‑1 (MQ=255)
                                                    gCACAGCCGCCTCGCTGGTTGCTTCCCAGCCAAGCAAGGCCAGTACAAGTGATGACAGGATGGGGCCTTTGAGCGGCGTCGACATGCTGCACCTCCGGGACACAGGATCTCTGCGGGGAGCTGGAGCAAGGTGTGTGCCAGAAGAAAAGt                                                                                    >  1:876479/1‑150 (MQ=255)
                                                        aGCCGCCTCGCTGGTTGCTTCCCAGCCAAGCAAGGCCAGTACAAGTGCTGACAGGATGGGGCCTTTGAACGGCGTCGACATGCTGCACCTCCGGGACACAGGATCTCTGCGGGGAGCTGGAGCAAGGTGCGTGCTAGATGAAAAGTTaaa                                                                                >  1:217932/1‑150 (MQ=255)
                                                         gccgccTCGCTGGTTGCTTCCCAGCCAAGCAAGGCCAGTACAAGTGATGACAGGATGGGGCCTTTGAGCGGCGTCGACATGCTGCACCTCCGGGACACAGGATCTCTGCGGGGAGCTGGAGCAAGGTGTGTGCCAGAAGAAAAGTTCAat                                                                               <  2:1073037/150‑1 (MQ=255)
                                                             ccTCGCTGGTTGCTTCCCAGCCAAGCAAGGCCAGTACAAGTGATGACAGGATGGGGCCTTTGAGCGGCGTCGACATGCTGCACCTCCGGGACACAGGATCTCTGCGGGGAGCTGGAGCAAGGTGTGTGCCAGAAGAAAAGTTCAATATTa                                                                           >  2:226652/1‑150 (MQ=255)
                                                             ccTCGCTGGTTGCTTCCCAGCCAAGCAAGGCCAGTACAAGTGATGACAGGATGGGGCCTTTGAGCGGCGTCGACATGCTGCACCTCCGGGACACAGGATCTCTGCGGGGAGCTGGAGCAAGGTGTGTGCCAGAAGAAAAGTTCAATATTa                                                                           <  2:852740/150‑1 (MQ=255)
                                                                     gTTGCTTCCCAGCCAAGCAAGGCCAGTACAAGTGATGACAGGATGGGGCCTTTGAGCGGCGTCGACATGCTGCACCTCCGGGACACAGGATCTCTGCGGGGAGCTGGAGCAAGGTGTGCGCCAGAAGAAAAGTTCAATATTAATCAATaa                                                                   <  2:876479/150‑1 (MQ=255)
                                                                      ttGCTTCCCAGCCAAGCAAGGCCAGTACAAGTGATGACAGGATGGGGCCTTTGAGCGGCGTCGACATGCTGCACCTCCGGGACACAGGATCTCTGCGGGGAGCTGGAGCAAGGTGTGTGCCAGAAGAAAAGTTCAATATTAATCAATAAc                                                                  <  1:910323/150‑1 (MQ=255)
                                                                        gCTTCCCAGCCAAGCAAGGCCAGTACAAGTCATGACAGGATGGGGCCTTTGAGCGGCGTCGACATGCTGCACCTCCGGGACACAGGATCTCTGCGGGGAGCTGGAGCAAGGTGTGTGCCAGAAGAAAAGTTCAATATTAATCAATAACtt                                                                <  1:1072063/150‑1 (MQ=255)
                                                                            cccAGCCAAGCAAGGCCAGTACAAGTGATGACAGGATGGGGCCTTTGAGCGGCGTCGACATGCTGCACCTCCGGGACACAGGATCTCTGCGGGGAGCTGGAGCAAGGTGTGTGCCAGAAGAAAAGTTCAATATTAATCAATAACTTGTAc                                                            <  1:680196/150‑1 (MQ=255)
                                                                              cAGCCAAGCAAGGCCAGTACAAGTGATGACAGGATGGGGCCTTTGAGCGGCGTCGACATGCTGCACCTCCGGGACACAGGATCTCTGCGGGGAGCTGGAGCAAGGTGTGTGCCAGAAGAAAAGTTCAATATTAAt                                                                         >  2:119759/1‑135 (MQ=255)
                                                                              cAGCCAAGCAAGGCCAGTACAAGTGATGACAGGATGGGGCCTTTGAGCGGCGTCGACATGCTGCACCTCCGGGACACAGGATCTCTGCGGGGAGCTGGAGCAAGGTGTGTGCCAGAAGAAAAGTTCAATATTAAt                                                                         <  1:119759/135‑1 (MQ=255)
                                                                               aGCCAAGCAAGGCCAGTACAAGTGATGACAGGATGGGGCCTTTGAGCGGCGTCGACATGCTGCACCTCCGGGACACAGGATCTCTGCGGGGAGCTGGAGCAAGGTGTGTGCCAGAAGAAAAGTTCAATATTAATCAATAACTTGTACGAt                                                         <  2:274719/150‑1 (MQ=255)
                                                                                gCCAAGCAAGGCCAGTACAAGTGATGACAGGATGGGGCCTTTGAGCGGCGTCGACATGCTGCACCTCCGGGACACAGGATCTCTGCGGGGAGCTGGAGCAAGGTGTGTGCCAGAAGAAAAGTTCAATATTAATCAATAACTTGTACGATc                                                        >  1:152060/1‑150 (MQ=255)
                                                                                        agGCCAGTACAAGTGATGACAGGATGGGGCCTTTGAGCGGCGTCGACATGCTGCACCTCCGGGACACAGGATCTCTGCGGGGAGCTGGAGCAAGGTGTGTGCCAGAAGAAAAGTTCAATATTAATCAATAACTTGTACGATCAACCGCaa                                                >  2:731600/1‑150 (MQ=255)
                                                                                        agGCCAGTACAAGTGATGACAGGATGGGGCCTTTGAGCGGCGTCGACATGCTGCACCTCCGGGACACAGGATCTCTGCGGGGAGCTGGAGCAAGGTGTGTGCCAGAAGAAAAGTTCAATATTAATCAATAACTTGTACCATCAACCGCaa                                                <  2:925860/150‑1 (MQ=255)
                                                                                              gTACAAGTGATGACAGGATGGGGCCTTTGAGCGGCGTCGACATGCTGCACCTCCGGGACACAGGATCTCTGCTGAGAGCTGGAGCAAGGTGTGTGCCAGAAGAAAAGTTCAATATTAATCAATAACTTGTACGATCAACCGCAACGGAtt                                          <  1:1163444/150‑1 (MQ=255)
                                                                                               tACAAGTGATGACAGGATGGGGCCTTTGAGCGGCGTCGACATGCTGCACCTCCGGGACACAGGATCTCTGCGGGGAGCTGGAGCAAGGTGTGTGCCAGAAGAAAAGTTCAATATTAATCAATAACTTGTACGATCAACCGCAACGGATTg                                         <  1:226652/150‑1 (MQ=255)
                                                                                                 cAAGTGATGACAGGATGGGGCCTTTGAGCGGCGTCGACATGCTGCACCTCCGGGACACAGGATCTCTGCGGGGAGCTGGAGCAAGGTGTGTGCCAGAAGAAAAGTTCAATATTAATCAATAACTTGTACGATCAACCGCAACGGATTggg                                       <  1:731600/150‑1 (MQ=255)
                                                                                                   aGTGATGACAGGATGGGGCCTTTGAGCGGCGTCGACATGCTGCACCTCCGGGACACAGGATCTCTGCGGGGAGCTGGAGCAAGGTGTGTGCCAGAAGAAAAGTTCAATATTAATCAATAACTTGTACGATCAACCGCAACGGATTGGgga                                     <  2:334851/150‑1 (MQ=255)
                                                                                                    gTGATGACAGGATGGGGCCTTTGAGCGGCGCCGACATGCTGCACCTCCGGGACACAGGATCTCTGCGGGGAGCTGGAGTAAGGTGTGTGCCAGAAGAAAAGTTCAATATTAATCAATAACTTGTACGATCAACCGCAACGGATTGGggag                                    <  2:217932/150‑1 (MQ=255)
                                                                                                       atgaCAGGATGGGGCCTTTGAGCGGCGTCGACATGCTGCACCTCCGGGACACAGGATCTCTGCGGGGAGCTGGAGCAAGGTGTGTGCCAGAAGAAAAGTTCAATATTAATCAATAACTTGTACGATCAACCGCAACGGATTGGGGAGGAt                                 >  1:453586/1‑150 (MQ=255)
                                                                                                         gaCAGGATGGGGCCTTTGAGCGGCGTCGACATGCTGCACCTCCGGGACACAGGATCTCTGCGGGGAGCTGGAGCAAGGTGTGTGCCAGAAGAAAAGTTCaa                                                                                >  1:938199/1‑101 (MQ=255)
                                                                                                         gaCAGGATGGGGCCTTTGAGCGGCGTCGACATGCTGCACCTCCGGGACACAGGATCTCTGCGGGGAGCTGGAGCAAGGTGTGTGCCAGAAGAAAAGTTCaa                                                                                <  2:938199/101‑1 (MQ=255)
                                                                                                          aCAGGATGGGGCCTTTGAGCGGCGTCGACATGCTGCACCTCCGGGACACAGGATCTCTGCGGGGAGCTGGAGCAAGGTGTGTGCCAGAAGAAAAGTTCAATATTAATCAATAACTTGTACGATCAAcc                                                    <  1:1254593/128‑1 (MQ=255)
                                                                                                          aCAGGATGGGGCCTTTGAGCGGCGTCGACATGCTGCACCTCCGGGACACAGGATCTCTGCGGGGAGCTGGAGCAAGGTGTGTGCCAGAAGAAAAGTTCAATATTAATCAATAACTTGTACGATCAAcc                                                    >  2:1254593/1‑128 (MQ=255)
                                                                                                                 ggggCCTTTGAGCGGCGTCGACATGCTGCACCTCCGGGACACAGGATCTCTGCGGGGAGCTGGAGCAAGGTGTGTGCCAGAAGAAAAGTTCAATATTAATCAATAACTTGTACGATCAACCGCAACGGATTGGGGAGGATCCCCcagcag                       <  1:1377725/150‑1 (MQ=255)
                                                                                                                   ggCCTTTGAGCGGCGTCGACATGCTGCACCTCCGGGACACAGGATCTCTGCGGGGAGCTGGAGCAAGGTGTGTGCCAGAAGAAAAGTTCAATATTAATCAATAACTTGTACGATCAACCGCAACGGATTGGGGAGGATCCCCCAGCAGGt                     >  1:622346/1‑150 (MQ=255)
                                                                                                                    gCCTTTGAGCGGCGTCGACATGCTGCACCTCCGGGACACAGGATCTCTGCGGGGAGCTGGAGCAAGGTGTGTGCCAGAAGaaaa                                                                                      <  2:1409971/84‑1 (MQ=255)
                                                                                                                    gCCTTTGAGCGGCGTCGACATGCTGCACCTCCGGGACACAGGATCTCTGCGGGGAGCTGGAGCAAGGTGTGTGCCAGAAGaaaa                                                                                      >  1:1409971/1‑84 (MQ=255)
                                                                                                                    gCCTTTGAGCGGCGTCGACATGCTGCACCTCCGGGACACAGGATCTCTGCGGGGAGCTGGAGCAAGGTGTGTGCCAGAAGAAAAGTTCAATATTAATCAATAACTTGTACGATCaa                                                      >  2:1286612/1‑116 (MQ=255)
                                                                                                                    gCCTTTGAGCGGCGTCGACATGCTGCACCTCCGGGACACAGGATCTCTGCGGGGAGCTGGAGCAAGGTGTGTGCCAGAAGAAAAGTTCAATATTAATCAATAACTTGTACGATCaa                                                      <  1:1286612/116‑1 (MQ=255)
                                                                                                                      cTTTGAGCGGCGTCGACATGCTGCACCTCCGGGACACAGGATCTCTGCGGGGAGCTGGAGCAAGGTGTGTGCCAGAAGAAAAGTTCAATATTAATCAATAACTTGTACGATCAACCGCAACGGATTGGGGAGGATCCCCCAGCAGGTTcc                  <  2:883399/150‑1 (MQ=255)
                                                                                                                      cTTTGAGCGGCGTCGACATGCTGCACCTCCGGGACACAGGATCTCTGCGGGGAGCTGGAGCAAGGTGTGTGCCAGAAGAAAAGTTCAATATTAATCAATAACTTGTACGATCAACCGCAACGGATTGGGGAGGATCCCCCAGCAGGTTcc                  >  1:1023640/1‑150 (MQ=255)
                                                                                                                                 gTCGACATGCTGCACCTCCGGGACACAGGATCTCTGCGGGGAGCTGGAGCAAGGTGTGTGCCAGAAGAAAAGTTCaa                                                                                <  1:92874/77‑1 (MQ=255)
                                                                                                                                 gTCGACATGCTGCACCTCCGGGACACAGGATCTCTGCGGGGAGCTGGAGCAAGGTGTGTGCCAGAAGAAAAGTTCaa                                                                                >  2:92874/1‑77 (MQ=255)
                                                                                                                                    gACATGCTGCACCTCCGGGACACAGGATCTATACGTGAAGCGGGAGCAAGGTGTGATCACGACGAAAAGTTCCCCCACCACCGCAAGAGTGGACAGGCAACCGCAACGGATTGGGGAGGATCCCCCAGCAAATTCCACCGGATTGGGGGa    >  1:325244/1‑150 (MQ=255)
                                                                                                                                    gACATGCTGCAACTCCGGGACACAGGATCTCTGCGGGGAGCTGGAGCAAGGTGTGTGCCAGAAGAAAAGTTCAATATTAATCAATAACTTGTACGATCAACCGCAACGGATTGGGGAGGAt                                 >  2:1345233/1‑121 (MQ=255)
                                                                                                                                    gACATGCTGCAACTCCGGGACACAGGATCTCTGCGGGGAGCTGGAGCAAGGTGTGTGCCAGAAGAAAAGTTCAATATTAATCAATAACTTGTACGATCAACCGCAACGGATTGGGGAGGAt                                 >  2:1344967/1‑121 (MQ=255)
                                                                                                                                    gACATGCTGCAACTCCGGGACACAGGATCTCTGCGGGGAGCTGGAGCAAGGTGTGTGCCAGAAGAAAAGTTCAATATTAATCAATAACTTGTACGATCAACCGCAACGGATTGGGGAGGAt                                 <  1:1344967/121‑1 (MQ=255)
                                                                                                                                    gACATGCTGCAACTCCGGGACACAGGATCTCTGCGGGGAGCTGGAGCAAGGTGTGTGCCAGAAGAAAAGTTCAATATTAATCAATAACTTGTACGATCAACCGCAACGGATTGGGGAGGAt                                 <  1:1345233/121‑1 (MQ=255)
                                                                                                                                      cATGCTGCACCTCCGGGACACAGGATCTCTGCGGGGAGCTGGAGCAAGGTGTGTGCCAGAAGAAAAGTTCAATATTAATCAATAACTTGTACGATCAACCGCAACGGATTGGGGAGGATCCCCCAGCAGGTTCCCCCGAATCGGGGGACg  >  2:1112658/1‑150 (MQ=255)
                                                                                                                                         gctgcACCTCCGGGACACAGGATCTCTGCGGGGAGCTGGAGCAAGGTGTGTGCCAGAAGaaaa                                                                                      >  2:587664/1‑63 (MQ=255)
                                                                                                                                         gctgcACCTCCGGGACACAGGATCTCTGCGGGGAGCTGGAGCAAGGTGTGTGCCAGAAGaaaa                                                                                      >  2:1301777/1‑63 (MQ=255)
                                                                                                                                         gctgcACCTCCGGGACACAGGATCTCTGCGGGGAGCTGGAGCAAGGTGTGTGCCAGAAGaaaa                                                                                      <  1:1301777/63‑1 (MQ=255)
                                                                                                                                         gctgcACCTCCGGGACACAGGATCTCTGCGGGGAGCTGGAGCAAGGTGTGTGCCAGAAGaaaa                                                                                      <  1:587664/63‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                            
ATCAGTGGCTGGTCCAGTGCCACCACGTTGGCGACCAGGGTGCGCTGGCACTGCACAGCCGCCTCGCTGGTTGCTTCCCAGCCAAGCAAGGCCAGTACAAGTGATGACAGGATGGGGCCTTTGAGCGGCGTCGACATGCTGCACCTCCGGGACACAGGATCTCTGCGGGGAGCTGGAGCAAGGTGTGTGCCAGAAGAAAAGTTCAATATTAATCAATAACTTGTACGATCAACCGCAACGGATTGGGGAGGATCCCCCAGCAGGTTCCCCCGAATCGGGGGACG  >  NC_002947/3960292‑3960575

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