Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,897,768 T→A 7.0% intergenic (+99/‑51) PP_2551 → / → PP_2552 LysR family transcriptional regulator/DOPA decarboxylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,897,7680TA 7.0% 200.9 / 3.4 68intergenic (+99/‑51)PP_2551/PP_2552LysR family transcriptional regulator/DOPA decarboxylase
Reads supporting (aligned to +/- strand):  ref base T (30/33);  new base A (2/3);  total (32/36)
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.98e-01

GCGCCAGTCGCAGTTCGGGGTGCAGGGCGGGCATGCGGGTTCGCCATAGCAGCTATGCGGTAAGCGAGGTTATTCGGCTGGGGATAGGTGCCTAGACTGGGGCATTGTGTTGATTGTGCGGCTTCTTCGCGGCTGTAGGCGCGGGTTTACCCGCGAAAGGGCCAGCACAGGCAATGGATAACCCCCAGGAGAGCTAGCGTGACCCCCGAACAATTCCGCCAGTACGGCCACCAACTGATCGACCTGATTGCCGACTACCGCCAGACCGTGGGCGAACGCCCGGTCATGGCCCAG  >  NC_002947/2897621‑2897914
                                                                                                                                                   |                                                                                                                                                  
gcgcCAGTCGCAGTTCGGGGTGCAGGGCGGGCATGCGGGTTCGCCATAGCAGCTATGCGGTAAGCGAGGTTATTCGGCTGGGGATAGGTGCCTAGACTGGGGCCTTGTGTTGATTGTTCGGCTTCTTCGCGGCTGTAGGCGCGGGTTTAc                                                                                                                                                  >  2:78421/1‑150 (MQ=255)
gcgcCAGTCGCAGTTCGGGGTGCAGGGCGGGCATGCGGGTTCGCCATAGCAGCTATGCGGTAAGCGAGGTTATTCGGCTGGGGATAGGTGCCTAGACTGGGGCATTGTGTTGATTGTGCGGCTTCTTCGCGGCTGTAGGCGCGGGTTTAc                                                                                                                                                  >  2:788514/1‑150 (MQ=255)
   ccAGTCGCAGTTCGGGGTGCAGGGCGGGCATGCGGGTTCGCCATAGCAGCTATGCGGTAAGCGAGGTTATTCGGCTGGGGATAGGTGCCTAGACTGGGGCATTGTGTTGATTGTGCGGCTTCTTCGCGGCTGTAGGCGCGGGTTTACCcg                                                                                                                                               >  1:254548/1‑150 (MQ=255)
      gTCGCAGTTCGGGGTGCAGGGCGGGCATGCGGGTTCGCCATAGCAGCTATGCGGTAAGCGAGGTTATTCGGCTGGGGATAGGTGCCTAGACTGGGGCATTGTGTTGATTGTGCGGCTTCTTCGCGGCTGTAGGCGCGGGTTTACCcgcgc                                                                                                                                            >  2:446328/1‑149 (MQ=255)
      gTCGCAGTTCGGGGTGCAGGGCGGGCATGCGGGTGCGCCATAGCAGCTATGCGGTAAGCGAGGTTATTCGGCTGGGGATAGGTGCCTAGACTGGGGCATTGTGTTGATTGTGCGGCTTCTTCGCGGCTGTAGGCGCGGGTTTACCCGCGa                                                                                                                                            >  2:910247/1‑150 (MQ=255)
         gCAGTTCGGGGTGCAGGGCGGGCATGCGGGTTCGCCATAGCAGCTATGCGGTAAGCGAGGTTATTCGGCTGGGGATAGGTGCCTAGACTGGGGCATTGTGTTGATTGTGCGGCTTCTTCGCGGCTGTAGGCGCGGGTTTACCCGCGAAAg                                                                                                                                         >  2:535156/1‑150 (MQ=255)
          cAGTTCGGGGTGCAGGGCGGGCATGCGGGTTCGCCATAGCAGCTATGCGGTAAGCGAGGTTATTCGGCTGGGGATAGGTGCCTAGACTGGGGCATTGTGTTGATTGTGCGGCTTCTTCGCGGCTGTAGGCGCGGGTTTACCCGCGAAAgg                                                                                                                                        <  1:690619/150‑1 (MQ=255)
            gTTCGGGGTGCAGGGCGGGCATGCGGGTTCGCCATAGCAGCTATGCGGTAAGCGAGGTTATTCGGCTGGGGATAGGTGCCTAGACTGGGGCATTGTGTTGATTGTGCGGCTTCTTCGCGGCTGTAGGCGCGGGTTTACCCGCGAAAGGGc                                                                                                                                      <  2:254548/150‑1 (MQ=255)
             ttCGGGGTGCAGGGCGGGCATGCGGGTTCGCCATAGCAGCTATGCGGTAAGCGAGGTTATTCGGCTGGGGATAGGTGCCTAGACTGGGGCATTGTGTTGATTGTGCGGCTTCTTCGCGGCTGTAGGCGCGGGTTTACCCGCGAAAGGGcc                                                                                                                                     <  1:535156/150‑1 (MQ=255)
                 gggTGCAGGGCGGGCATGCGGGTTCGCCATAGCAGCTATGCGGTAAGCGAGGTTATTCGGCTGGGGATAGGTGCCTAGACTGGGGCATTGTGTTGATTGTGCGGCTTCTTCGCGGCTGTAGGCGCGGGTTTACCCGCGAAAg                                                                                                                                         <  1:710300/142‑1 (MQ=255)
                 gggTGCAGGGCGGGCATGCGGGTTCGCCATAGCAGCTATGCGGTAAGCGAGGTTATTCGGCTGGGGATAGGTGCCTAGACTGGGGCATTGTGTTGATTGTGCGGCTTCTTCGCGGCTGTAGGCGCGGGTTTACCCGCGAAAg                                                                                                                                         >  2:710300/1‑142 (MQ=255)
                   gTGCAGGGCGGGCATGCGGGTTCGCCATAGCAGCTATGCGGTAAGCGAGGTTATTCGGCTGGGGATAGGTGCCTAGACTGGGGCATTGTGTTGATTGTGCGGCTTCTTCGCGGCTGTAGGCGCGGGTTTACCCGCGAAAGGGCCAGcaca                                                                                                                               <  2:1095953/150‑1 (MQ=255)
                    tGCAGGGCGGGCATGCGGGTTCGCCATAGCAGCTATGCGGTAAGCGAGGTTATTCGGCTGGGGATAGGTGCCTAGACTGGGGCATTGTGTTGATTGTGCGGCTTCTTCGCGGCTGTAGGCGCGGGTTTACCCGCGAAAGGGCCAGCACAg                                                                                                                              <  1:1052131/150‑1 (MQ=255)
                      cAGGGCGGGCATGCGGGTTCGCCATAGCAGCTATGCGGTAAGCGAGGTTATTCGGCTGGGGATAGGTGCCTAGACTGGGGCATTGTGTTGATTGTGCGGCTTCTTCGCGGCTGTAGGCGCGGGTTTACCCGCGAAAGGGCCAGCACAGGc                                                                                                                            >  2:319656/1‑150 (MQ=255)
                       aGGGCGGGCATGCGGGTTCGCCATAGCAGCTATGCGGTAAGCGAGGTTATTCGGCTGGGGATAGGTGCCTAGACTGGGGCATTGTGTTGATTGTGCGGCTTCTTCGCGGCTGTAGGCGCGGGTTTACCCGCGAAAGGGCCAGCACAGGCa                                                                                                                           >  2:1026841/1‑150 (MQ=255)
                            gggcATGCGGGTTCGCCATAGCAGCTATGCGGTAAGCGAGGTTATTCGGCTGGGGATAGGTGCCTAGACTGGGGCATTGTGTTGATTGTGCGGCTTCTTCGCGGCTGTAGGCGCGGGTTTACCCGCGAAAGGGCCAGCACAGGCAATGGa                                                                                                                      <  1:935505/150‑1 (MQ=255)
                                       ttCGCCATAGCAGCTATGCGGTAAGCGAGGTTATTCGGCTGGGGATAGGTGCCTAGACTGGGGCATTGTGTTGATTGTGCGGCTTCTTCGCGGCTGTAGGCGCGGGTTTACCCGCGAAAGGGCCAGCACAg                                                                                                                              <  1:680999/131‑1 (MQ=255)
                                       ttCGCCATAGCAGCTATGCGGTAAGCGAGGTTATTCGGCTGGGGATAGGTGCCTAGACTGGGGCATTGTGTTGATTGTGCGGCTTCTTCGCGGCTGTAGGCGCGGGTTTACCCGCGAAAGGGCCAGCACAg                                                                                                                              >  2:680999/1‑131 (MQ=255)
                                       ttCGCCATAGCAGCTATGCGGTAAGCGAGGTTATTCGGCTGGGGATAGGTGCCTAGACTGGGGCATTGTGTTGATTGTGCGGCTTCTTCGCGGCTGTAGGCGCGGGTTTACCCGCGAAAGGGCCAGCACAGGCAATGGATAACCCCCAgg                                                                                                           <  1:319656/150‑1 (MQ=255)
                                          gCCATAGCAGCTATGCGGTAAGCGAGGTTATTCGGCTGGGGATAGGTGCCTAGACTGGGGGATTGTGTTGATTGTGCGGCTTCTTCGCGGCTGTAGGCGCGGGTTTACCCGCGAAAGGGCCAGCACAGGCAATGGATAACCCCCAGgaga                                                                                                        <  2:126495/150‑1 (MQ=255)
                                          gCCATAGCAGCTATGCGGTAAGCGAGGTTATTCGGCTGGGGATAGGTGCCTAGACTGGGGCATTGTGTTGATTGTGCGGCTTCTTCGCGGCTGTAGGCGCGGGTTTACCCGCGAAAGGGCCAGCACAGGCAATGGATAACCCCCAGgaga                                                                                                        <  2:132979/150‑1 (MQ=255)
                                             aTAGCAGCTATGCGGTAAGCGAGGTTATTCGGCTGGGGATAGGTGCCTAGACTGGGGCATTGTGTTGATTGTGCGGCTTCTTCGCGGCTGTAGGCGCGGGTTTACCCGCGAAAGGGCCAGCACAGGCAATGGATAACCCCCAGGAGAGCt                                                                                                     <  2:947899/150‑1 (MQ=255)
                                              tAGCAGCTATGCGGTAAGCCATGGTATTCGGCTGGGGATAGGTGCCTAGACTGGGGCATTGTGTTGCTTGTGCGGCTTCTTCGCGGCTGTAGGCGCGGGTTTACCCGCGAAAGGGCCAGCACAGGCAATGGATAACCCCCAGGAGAGCTa                                                                                                    <  1:899095/150‑1 (MQ=255)
                                                   gcTATGCGGTAAGCGAGGTTATTCGGCTGGGGATAGGTGCCTAGACTGGGGCATTGTGTTGATTGTGCGGCTTCTTCGCGGCTGTAGGCGCGGGTTTACCCGCGAAAGGGCCAGCACAGGCAATGGATAACCCCCAGGAGAGCTAGCGTg                                                                                               >  2:579860/1‑150 (MQ=255)
                                                   gcTATGCGGTAAGCGAGGTTATTCGGCTGGGGATAGGTGCCTAGACTGGGGCATTGTGTTGATTGTGCGGCTTCTTCGCGGCTGTAGGCGCGGGTTTACCCGCGAAAGGGCCAGCACAGGCAATGGATAACCCCCAGGAGAGCTAGCGTg                                                                                               <  1:255792/150‑1 (MQ=255)
                                                    cTATGCGGTAAGCGAGGTTATTCGGCTGGGGATAGGTGCCTAGACTGGGGCATTGTGTTGATTGTGCGGCTTCTTCGCGGCTGTAGGCGCGGGTTTACCCGCGAAAGGGCCAGCACAGGCAATGGATAACCCCCAGGAGAGCTAGCGTGa                                                                                              <  1:788514/150‑1 (MQ=255)
                                                          ggTAAGCGAGGTTATTCGGCTGGGGATAGGTGCCTAGACTGGGGCATTGTGTTGATTGTGCGGCTTCTTCGCGGCTGTAGGCGCGGGTTTACCCGCGAAAGGGCCAGCACAGGCAATGGATAACCCCCAGGAGAGCTAGCGTGACCCCCg                                                                                        <  2:850049/150‑1 (MQ=255)
                                                                    gTTATTCGGCTGGGGATAGGTGCCTAGACTGGGGCATTGTGTTGATTGTGCGGCTTCTTCGCGGCTGTAGGCGCGGGTTTACCCGCGAAAGGGCCAGCACAGGCAATGGATAACCCCCAGGAGAGCTAGCGTGACCCCCGAACAATTCCg                                                                              >  1:736774/1‑150 (MQ=255)
                                                                     ttATTCGGCTGGGGATAGGTGCCTAGACTGGGGCATTGTGTTGATTGTGCGGCTTCTTCGCGGCTGTAGGCGCGGGTTTACCCGCGAAAGGGCCAGCACAGGCAATGGATAACCCCCAGGAGAGCTAGCGTGACCCCCGAACAATTCCGc                                                                             <  1:1053817/150‑1 (MQ=255)
                                                                                  gATNGGTGCCTAGACTGGGGCATTGTGTTGATTGTGCGGCTTCTTCGCGGCTGTAGGCGCGGGTTTACCCGCGAAAGGGCCAGCACAGGCAATGGATAACCCCCAGGAGAGCTAGCGTGACCCCCGAACAATTCCGCCAGTACGGccacc                                                                >  2:702979/1‑150 (MQ=255)
                                                                                       gTGCCTAGACTGGGGCATTGTGTTGATTGTGCGGCTTCTTCGCGGCTGTAGGCGCGGGTTTACCCGCGAAAGGGCCAGCACAGGCAATGGATAACCCCCAGGAGAGCTAGCGTGACCCCCGAACAATTCCGCCAGTACGGCCACCAACTg                                                           >  1:589562/1‑150 (MQ=255)
                                                                                        tGCCTAGACTGGGGCATTGTGTTGATTGTGCGGCTTCTTCGCGGCTGTAGGCGCGGGTTTACCCGCGAAAGGGCCAGCACAGGCAATGGATAACCCCCAGGAGAGCTAGCGTGACCCCCGAACAATTCCGCCAGTACGGCCACCAACTGa                                                          <  2:300789/150‑1 (MQ=255)
                                                                                         gCCTAGACTGGGGCATTGTGTTGATTGTGCGGCTTCTTCGCGGCTGTAGGCGCGGGTTTACCCGCGAAAGGGCCAGCGCAGGCAATGGATAACCCCCAGGAGAGCTAGCGTGAc                                                                                             >  1:816602/1‑114 (MQ=255)
                                                                                         gCCTAGACTGGGGCATTGTGTTGATTGTGCGGCTTCTTCGCGGCTGTAGGCGCGGGTTTACCCGCGAAAGGGCCAGCGCAGGCAATGGATAACCCCCAGGAGAGCTAGCGTGAc                                                                                             <  2:816602/114‑1 (MQ=255)
                                                                                          ccTAGACTGGGGCATTGTGTTGATTGTGCGGCTTCTTCGCGGCTGTAGGCGCGGGTTTACCCGCGAAAGGGCCAGCACAGGCAATGGATAACCCCCAGGAGAGCTAGCGTGACCCCCGAACAATTCCGCCAGTACGGCCACCAACTGATc                                                        <  1:1026841/150‑1 (MQ=255)
                                                                                                 tGGGGCATTGTGTTGATTGTGCGGCTTCTTCGCGGCTGTAGGCGCGGGTTTACCCGCGAAAGGGCCAGCACAGGCAATGGATAACCCCCAGGAGAGCTAGCGTGACCCCCGAACAATTCCGCCAGTACGGCCACCAACTGATCGACCTGa                                                 >  1:930950/1‑150 (MQ=255)
                                                                                                   gggCATTGTGTTGATTGTGCGGCTTCTTCGCGGCTGTAGGCGCGGGTTTACCCGCGAAAGGGCCAGCACAGGCAATGGATAACCCCCAGGAGAGCTAGCGTGACCCCCGAACAATTCCGCCAGTACGGCCACCAACTGATCGACCTGAtt                                               <  1:910247/150‑1 (MQ=255)
                                                                                                    ggCATTGTGTTGATTGTGCGGCTTCTTCGCGGCTGTAGGCGCGGGTTTACCCGCGAAAGGGCCAGCACAGGCAATGGATAACCCCCAGGAGAGCTAGCGTGACCCCCGCAGAATTCCGCCAGTACGGCCACCAACTGATCGACCTGCTTg                                              >  1:396934/1‑150 (MQ=255)
                                                                                                     gCATTGTGTTGATTGTGCGGCTTCTTCGCGGCTGTAGGCGCGGGTTTACCCGCGAAAGGGCCAGCACAGGCAATGGATAACCCCCAGGAGAGCTAGCGTGACCCCCGAACAATTCCGCCAGTACGGCCACCAACTGATCGACCTGATTGc                                             >  1:766426/1‑150 (MQ=255)
                                                                                                     gCATTGTGTTGATTGTGCGGCTTCTTCGCGGCTGTAGGCGCGGGTTTACCCGCGAAAGGGCCAGCACAGGCAATGGATAACCCCCAGGAGAGCTAGCGTGACCCCCGAACAATTCCGCCAGTACGGCCACCAACTGATCGACCTGATTGc                                             <  2:914181/150‑1 (MQ=255)
                                                                                                      cATTGTGTTGATTGTGCGGCTTCTTCGCGGCTGTAGGCGCGGGTTTACCCGCGAAAGGGCCAGCACAGGCAATGGATAACCCCCAGGAGAGCTAGCGTGACCCCCGAACAATTCCGCCAGTACGGCCACCAACTGATCGACCTGATTGcc                                            >  2:894692/1‑150 (MQ=255)
                                                                                                               gATTGTGCGGCTTCTTCGCGGCTGTAGGCGCGGGTTTACCCGCGAAAGGGCAAAAGCAGGCAATGGATAACCCCCAGGAGAGCTAGCGTGACCCCCGAACAATTCCGCCAGTACGGCCACCAACTGATCGACCTGATTGCCGACTACCGc                                   >  1:357088/1‑150 (MQ=255)
                                                                                                               ccgatctcGGCTTCTTCGCGGCTGTAGGCGCGGGTTAACCCGCGAAAGGGCCAGCACAGGCAATGGATAACCCCCAGGAGCGCTAGCGTGACCCCCGAACAATTCCGCCAGTTCGGCCACCAACTGATCGACCTGATTGCCGACTACCGc                                   <  1:790526/143‑1 (MQ=255)
                                                                                                                   gtgCGGCTTCTTCGCGGCTGTAGGCGCGGGTTTACCCGCGAAAGGGCCAGCACAGGCAATGGATAACCCCCAGGAGAGCTAGCGTGACCCCCGAACAATTCCGCCAGTACGGCCACCAACTGATCGACCTGATTGCCGACTACCGCCAGa                               >  1:1027963/1‑150 (MQ=255)
                                                                                                                       ggCTTCTTCGCGGCTGTAGGCGCGGGTTTACCCGCGAAAGGGCCAGCCCAGGCAATGGATAACCCCCAGGAGAGCTAGCGTGACCCCCGAACaa                                                                                   >  1:476092/1‑94 (MQ=255)
                                                                                                                       ggCTTCTTCGCGGCTGTAGGCGCGGGTTTACCCGCGAAAGGGCCAGCACAGGCAATGGATAACCCCCAGGAGAGCTAGCGTGACCCCCGAACaa                                                                                   <  2:476092/94‑1 (MQ=255)
                                                                                                                       ggCTTCTTCGCGGCTGTAGGCGCGGGTTTACCCGCGAAAGGGCCAGCACAGGCAATGGATAACCCCCAGGAGAGCTAGCGTGACCCCCGAACAATGCCGCCAGTACGGCCACCAACTGAGCGACCTGATTGCCGACTACCGCCAGATCgg                           >  2:790526/1‑149 (MQ=255)
                                                                                                                        gCTTCTTCGCGGCTGTAGGCGCGGGTTTACCCGCGAAAGGGCCAGCACAGGCAATGGATAACCCCCAGGAGAGCTAGCGTGACCCCCGAACAATTCCGCCAGTACGGCCACCAACTGATCGACCTGATTGCCGACTACCGCCAGACCggg                          <  1:783893/150‑3 (MQ=255)
                                                                                                                         cttcttCGCGGCTGTAGGCGCGGGTTTACCCGCGAAAGGGCCAGCACAGGCAATGGATAACCCCCAGGAGAGCTAGCGTGACCCCCGAACAATTCCGCCAGTACGGCCACCAACTGATCGACCTGATTGCCGACTACCGCCAGACCGTgg                         >  1:992703/1‑150 (MQ=255)
                                                                                                                             ttcGCGGCTGTAGGCGCGGGTTTACCCGCGAAAGGGCCAGCACAGGCAATGGATAACCCCCAGGAGAGCTAGCGTGACCCCCGAACAATTCCGCCAGTACGGCCACCAACTGATCGACCTGATTGCCGACTACCGCCAGACCGTGGGCGa                     >  2:251080/1‑150 (MQ=255)
                                                                                                                               cgcgGCTGTAGGCGCGGGTTTACCCGCGAAAGGGCCAGCACAGGCAATGGATAACCCCCAGGAGAGCTAGCGTGACCCCCGAACAATTCCGCCAGTACGGCCACCAACTGATCGACCTGATTGCCGACTACCGCCAGACCGTGGGCGAAc                   <  2:930950/150‑1 (MQ=255)
                                                                                                                                 cgGCTGTAGGCGCGGGTTTACCCGCGAAAGGGCCAGCACAGGCAATGGATAACCCCCAGGAGAGCTAGCGTGACCCCCGAACAATTCCGCCAGTACGGCCGCCAACTGATCGACCTGATTGCCGACTACCGCCAGACCGTGGGCGAACGc                 >  2:574171/1‑150 (MQ=255)
                                                                                                                                   gCTGTAGGCGCGGGTTTACCCGCGAAAGGGCCAGCACAGGCAATGGATAACCCCCAGGAGAGCTAGCGTGACCCCCGAACAATTCCGCCAGTACGGCCACCAACTGATCGACCTGATTGCCGACTACCGCCAGACCGTGGGCGAACGccc               <  1:251080/150‑1 (MQ=255)
                                                                                                                                    gatctggCGCGGGTTTACCCGCGAAAGGGCCAGCACAGGCAATGGATAACCCCCAGGAGAGCTAGCGTGACCCCCGAACAATTCCGCCAGTACGGCCACCCCCTGATCGACCTGATTGCCGACTACCGCCAGACCGTGGGCGAACGCCCg              <  1:256704/145‑1 (MQ=255)
                                                                                                                                     tGTAGGCGCGGGTTTACCCGCGAAAGGGCCAGCACAGGCAATGGATAACCCCCAGGAGAGCTAGCGTGACCCCCGAACAATTCCGCCAGTACGGCCACCAACTGATCGACCTGATTGCCGACTACCGCCAGACCGTGGGCGAACGCCCgg             >  1:595916/1‑150 (MQ=255)
                                                                                                                                        tttcgcgGGTAAACCCGCGAAAGGGCCAGCACAGGCAATGGATAACCCCCAGGAGAGCTAGCGTGACCCCCGAACAATTCCGCCAGTACGGCCACCAACTGATCGACCTGa                                                 <  2:623412/108‑1 (MQ=255)
                                                                                                                                        tttcgcgGGTAAACCCGCGAAAGGGCCAGCACAGGCAATGGATAACCCCCAGGAGAGCTAGCGTGACCCCCGAACAATCCCGCCAGTACGGCCACCAACTGATCGACCTGa                                                 >  1:623412/4‑111 (MQ=39)
                                                                                                                                         ttcgcgGGTAAACCCGCGAAAGGGCCAGCACAGGCAATGGATAACCCCCAGGAGAGCTAGCGTGACCCCCGAACAATTCCGCCAGTACGGCCACCAACTGATc                                                        >  2:917178/3‑103 (MQ=255)
                                                                                                                                         ttcgcgGGTAAACCCGCGAAAGGGCAAGCACAGGCAATGGATAACCCCCAGGAGAGCTAGCGTGACCCCCGAACAATTCCGCCAGTACGGCCACCAACTGATc                                                        <  1:917178/101‑1 (MQ=255)
                                                                                                                                         ggCGCGGGTTTACCCGCGAAAGGGCCAGCACAGGCAATGGATAACCCCCAGGAGAGCTAGCGTGACCCCCGAACAATTCCGCCAGTACGGCCACCACCTGATCGACCTGATTGCCGACTCCCGCTAGACCGTGGGCGAACGCCCGGTCAt         <  1:579860/150‑1 (MQ=255)
                                                                                                                                          gcgcGGGTTTACCCGCGAAAGGGCCAGCACAGGCAATGGATAACCCCCAGGAGAGCTAGCGTGACCCCCGAACAATTCCGCCAGTACGGCCACCAACTGATCGACCTGATTGCCGACTACCGCCCGACCGTGGTCGAACGCCCggagatc        >  2:256704/1‑145 (MQ=255)
                                                                                                                                          gcgcGGGTTTACCCGCGAAAGGGCCAGCACAGGCAATGGATAACCCCCAGGAGAGCTAGCGTGACCCCCGAACAATTCCGCCAGTACGGCCACCAACTGATCGACCTGATTGCCGACTACCGCCAGACCGTGGGCGAACGCCCGGTCATg        <  2:595916/150‑1 (MQ=255)
                                                                                                                                             cgGGTTTACCCGCGAAAGGGCCAGCACAGGCAATGGATAACCCCCAGGAGAGCTAGCGTGACCCCCGAACAATTCCGCCAGTACGGCCACCAACTGATCGACCTGATTGCCGACTACCGCCAGACCGTGGGCGAACGCCCGGTCATGGcc     <  2:589562/150‑1 (MQ=255)
                                                                                                                                                gTTTACCCGCGAAAGGGCCAGCACAGGCAATGGATAACCCCCAGGAGAGCTAGCGTGACCCCCGAACAATTCCGCCAGTACGGCCACCaa                                                              <  2:1034260/90‑1 (MQ=255)
                                                                                                                                                gTTTACCCGCGAAAGGGCCAGCACAGGCAATGGATAACCCCCAGGAGAGCTAGCGTGACCCCCGAACAATTCCGCCAGTACGGCCACCaa                                                              <  1:114525/90‑1 (MQ=255)
                                                                                                                                                gTTTACCCGCGAAAGGGCCAGCACAGGCAATGGATAACCCCCAGGAGAGCTAGCGTGACCCCCGAACAATTCCGCCAGTACGGCCACCaa                                                              >  2:114525/1‑90 (MQ=255)
                                                                                                                                                gTTTACCCGCGAAAGGGCCAGCACAGGCAATGGATAACCCCCAGGAGAGCTAGCGTGACCCCCGAACAATTCCGCCAGTACGGCCACCaa                                                              >  1:1034260/1‑90 (MQ=255)
                                                                                                                                                gTTTACCCGCGAAAGGGCCAGCACAGGCAATGGATAACCCCCAGGAGAGCTAGCGTGACCCCCGAACAATTCCGCCAGTACGGCCACCAACTGATCGACCTGATTGCCGACTACCGCCAGACCGTGGGCGAACGCCCGGTCATGGCCCAg  <  1:154517/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
GCGCCAGTCGCAGTTCGGGGTGCAGGGCGGGCATGCGGGTTCGCCATAGCAGCTATGCGGTAAGCGAGGTTATTCGGCTGGGGATAGGTGCCTAGACTGGGGCATTGTGTTGATTGTGCGGCTTCTTCGCGGCTGTAGGCGCGGGTTTACCCGCGAAAGGGCCAGCACAGGCAATGGATAACCCCCAGGAGAGCTAGCGTGACCCCCGAACAATTCCGCCAGTACGGCCACCAACTGATCGACCTGATTGCCGACTACCGCCAGACCGTGGGCGAACGCCCGGTCATGGCCCAG  >  NC_002947/2897621‑2897914

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: