Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,659,951 G→A 6.2% intergenic (‑43/+32) PP_3223 ← / ← PP_3224 ABC transporter substrate‑binding protein/class II aldolase/adducin domain‑containing protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,659,9510GA 6.2% 182.6 / 2.5 64intergenic (‑43/+32)PP_3223/PP_3224ABC transporter substrate‑binding protein/class II aldolase/adducin domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base G (28/32);  new base A (2/2);  total (30/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.45e-01

GGCGATGTTGAGGGTCTTGCTGGTTGAGGTGTTGTCGGCCGGCGAGCAAGCGGTGAGGCCCAGGGCCAGGGTGCCGATCGCGAGCAGTTGGTGCAGTGGGGGCATGGGATTGATCCTTCAATGAATGGGGCTCTCCTGTTCCGGCCAGTTCGCCGGAACAGGCAAAAGGTGCCTTGGAGTCAGTCCTGGTAATCGAGGCCCTCGCGATCCAACTGCCGAATCAGCGCATCCCAATGCCGTTGGATACCTTCCCCCTGGCCATCGGAAAACCGCTTGGACTGCTCATGCACCTTGA  >  NC_002947/3659804‑3660098
                                                                                                                                                   |                                                                                                                                                   
ggCGATGTTGAGGGTCTTGCTGGTTGAGGTGTTGTCGGCCGGCGAGCAAGCGGTGAGGCCCAGGGCCAGGGTGCCGATCGCGAGCAGTTGGTGCAGTGGGGGCATGGGATTGATCCTTCAATGAATGGGGCTCTCCTGTTCCGGCCAGtt                                                                                                                                                   >  2:494139/1‑150 (MQ=255)
ggCGATGTTGAGGGTCTTGCTGGTTGAGGTGTTGTCGGCCGGCGAGCAAGCGGTGAGGCCCAGGGCCAGGGTGCCGATCGCGAGCAGTTGGTGCAGTGGGGGCATGGGATTGATCCTTCAATGAATGGGGCTCTCCTGTTCCGGCCAGtt                                                                                                                                                   <  1:552363/150‑1 (MQ=255)
   gATGTTGAGGGTCTTGCTGGTTGAGGTGTTGTCGGCCGGCGAGCAAGGGGTAAGGCCCAGGGCCAGGGTGCCGATCGCGAGCAGTTGGTGCAGTGGGGGCATGGGATTGATCCTTCAATGAATGGGGCTTTGTTGTTTGTGCCAGtttcc                                                                                                                                                >  2:1292169/1‑147 (MQ=255)
    aTGTTGAGGGTCGTGCTGGTGGAGGTGTGGTCGGCCGGCGAGCAAGCGGTGAGGCCCAGGGCCAGGGTGCCGATCGCGAGCAGTTGGTGCAGTGGGGGCATGGGATTGATCCTTCAATGAATGGGGCTCTCCTGTTCCGGCCAGTTCGcc                                                                                                                                               <  2:228833/150‑1 (MQ=255)
     tGTTGAGGGTCTTGCTGGTTGAGGTGTTGTCGGCCGGCGAGCAAGCGGTGAGGCCCAGGGCCAGGGTGCCGATCGCGAGCAGTTGGTGCAGTGGGGGCATGGGATTGATCCTTCAATGAATGGGGCTCTCCTGTTCCGGCCAGTTCGCCg                                                                                                                                              <  1:494139/150‑1 (MQ=255)
      gTTGAGGGTCTTGCTGGTTGAGGTGTTGTCGGCCGGCGAGCAAGCGGTGAGGCCCAGGGCCAGGGTGCCGATCGCGAGCAGTTGGTGCAGTGGGGGCATGGGATTGATCCTTCAATGAATGGGGCTCTCCTGTTCCGGCCAGTTCGCCgg                                                                                                                                             >  2:987175/1‑150 (MQ=255)
          aGGGTCTTGCTGGTTGAGGTGTTGTCGGCCGGCGAGCAAGCGGTGAGGCCCAGGGCCAGGGTGCCGATCGCGAGCAGTTGGTGCAGTGGGGGCATGGGATTGATCCTTCAATCAATGGGGCTCTCCTGTTCCGGCCAGTTCGCCGGAACa                                                                                                                                         <  1:1150404/150‑1 (MQ=255)
                   cTGGTTGAGGTGTTGTCGGCCGGCGAGCAAGCGGTGAGGCCCAGGGCCAGGGTGCCGATCGCGAGCAGTTGGTGCAGTGGGGGCATGGGATTGATCCTTCAATGAATGGGGCTCTCCTGTTCCGGCCAGTTCGCCGGAACAGGCAAAAgg                                                                                                                                >  1:396884/1‑150 (MQ=255)
                    tGGATGAGGTGCTTGCGGCCGGCGAGCAAGCGGTGAGGCCCAGGGCCAGGGTGCCGATCGCGAGCAGTTGGTGCAGTGGGGGCATGGGATTGATCCTTCAATGAATGGGGCTCTCCTGTTCCGGCCAGTTCGCCGGAACAGGCAAAAGGt                                                                                                                               <  2:325340/150‑1 (MQ=255)
                         gAGGTGTTGTCGGCCGGCGAGCAAGCGGTGAGGCCCAGGGCCAGGGTGCCGATCGCGAGCAGTTGGTGCAGTGGGGGCATGGGATTGATCCTTCAATGAATGGGGCTCTCCTGTTCCGGCCAGTTCGCCGGAACAGGCAAAAGGTGCCtt                                                                                                                          >  1:1418563/1‑150 (MQ=255)
                         gAGGTGTTGTCGGCCGGCGAGCAAGCGGTGAGGCCCAGGGCCAGGGTGCCGATCGCGAGCAGTTGGTGCAGTGGGGGCATGGGATTGATCCTTCAATGAATGGGGCTCTCCTGTTCCGGCCAGTTCGCCGGAACAGGCAAAAGGTGCCtt                                                                                                                          <  1:1292169/150‑1 (MQ=255)
                          aGGCGTTGTTGGCCGGCGGGCAGGCGGTGAGGCCCAGGGCCAGGGTGCCGATCGCGAGCAGTTGGTGCAGTGGGGGCATGAGACTGATTCTTCAATGAACGGGGCTTTCCTGTTCCGGCCAGTTTGCCGGAACAGGCAAACGGTGCCTTg                                                                                                                         <  2:10979/150‑1 (MQ=255)
                            gtgtTGTCGGCCGGCGAGCAAGCGGTGAGGCCCAGGGCCAGGGTGCCGATCGCGAGCAGTTGGTGCAGTGGGGGCATGGGATTGATCCTTCAATGAATGGGGCTCTCCTGTTCCGGCCAGTTCGCCGGAACAGGCAAAAGGTGCCTTGGa                                                                                                                       <  2:1418563/150‑1 (MQ=255)
                            gtgtTGTCGGCCGGCGAGCAAGCGGTGAGGCCCAGGGCCAGGGTGCCGATCGCGAGCAGTTGGTGCAGTGAGTGCATTGGATTGATCCTTCAATGAATGGGGCTCTCCTGTTCCGGCCAGTTCGCCGGAACAGGCAAAAGGTGCCTTGGa                                                                                                                       <  2:691806/150‑1 (MQ=255)
                             tgttgtCGGCCGGCGAGCAAGCGGTGAGGCCCAGGGCCAGGGTGCCGATCGCGAGCAGTTGGTGCAGTGGGGGCATGGGATTGATCCTTCAATGAATGGGGCTCTCCTGTTCCGGCGAACTGGCCGGAACAg                                                                                                                                        >  2:620571/1‑132 (MQ=255)
                             tgttgtCGGCCGGCGAGCAAGCGGTGAGGCCCAGGGCCAGGGTGCCGATCGCGAGCAGTTGGTGCAGTGGGGGCATGGGATTGATCCTTCAATGAATGGGGCTCTCCTGTTCCGGCGAACTGGCCGGAACAg                                                                                                                                        <  1:620571/132‑1 (MQ=255)
                                tgtCGGCCGGCGAGCAAGCGGTGAGGCCCAGGGCCAGGGTGCCGATCGCGAGCAGTTGGTGCAGTGGGGGCATGGGATTGATCCTTCAATGAATGGGGCTCTCCTGTTCCGGCCAGTTCGCCGGAACAGGCAAAAGGTGCCTTGGagtca                                                                                                                   <  2:396884/150‑1 (MQ=255)
                                     gccggcGAGCAAGCGGTGAGGCCCAGGGCCAGGGTGCCGATCGCGAGCAGTTGGTGCAGTGGGGGCATGGGATTGATCCTTCAATGAATGGGGCTCTCCTGTTCCGGCCAGTTCGCCGGAACAGGCAAAAGGTGCCTTGGAGTCAGTCCt                                                                                                              <  1:584709/150‑1 (MQ=255)
                                              cAAGCGGTGAGGCCCAGGGCCAGGGTGCCGATCGCGAGCAGTTGGTGCAGTGGGGGCATGGGATTGATCCTTCAATGAATGGGGCTCTCCTGTTCCGGCCAGTTCGCCGGAACAGGCAAAAGGTGCCTTGGAGTCAGTCCTGGTAATCGa                                                                                                     >  1:62778/1‑150 (MQ=255)
                                              cAAGCGGTGAGGCCCAGGGCCAGGGTGCCGATCGCGAGCAGTTGGTGCAGTGGGGGCATGGGATTGATCCTTCAATGAATGGGGCTCTCCTGTTCCGGCCAGTTCGCCGGAACAGGCAAAAGGTGCCTTGGAGTCAGTCCTGGTAATCGa                                                                                                     <  1:781703/150‑1 (MQ=255)
                                               aaGCGGTGAGGCCCAGGGCCAGGGTGCCGATCGCGAGCAGTTGGTGCAGTGGGGGCATGGGATTGATCCTTCCGTGATTGGGGCTCTCCTGTTCAGAGCAGTTGGGACGAGCATGGTAAAGGTGCCTTGGAGTCAGTCCTGGTAATCGAg                                                                                                    >  1:1385771/1‑150 (MQ=255)
                                               aaGCGGTGAGGCCCAGGGCCAGGGTGCCGATCGCGAGCAGTTGGTGCAGTGGGGGCATGGGATTGATCCTTCAATGAATGGGGCTCTCCTGTTCCGGCCAGTTCGCCGGAACAGGCAAAAGGTGCCTTGGa                                                                                                                       >  2:860487/1‑131 (MQ=255)
                                               aaGCGGTGAGGCCCAGGGCCAGGGTGCCGATCGCGAGCAGTTGGTGCAGTGGGGGCATGGGATTGATCCTTCAATGAATGGGGCTCTCCTGTTCCGGCCAGTTCGCCGGAACAGGCAAAAGGTGCCTTGGa                                                                                                                       <  1:860487/131‑1 (MQ=255)
                                               aaGCGGTGAGGCCCAGGGCCAGGGTGCCGATCGCGAGCAGTTGGTGCAGTGGGGGCATGGGATTGATCCTTCAATGAATGGGGCTCTCCTGTTCCGGCCAGTTCGCCGGAACAGGCAAAAGGTGCCTTGGAGTCAGTCCTGGTAATCGAg                                                                                                    >  2:935552/1‑150 (MQ=255)
                                               aaGCGGTGAGGCCCAGGGCCAGGGTGCCGATCGCGAGCAGTTGGTGCAGTGGGGGCATGGGATTGATCCTTCAATGAATGGGGCTCTCCTGTTCCGGCCAGTTCGCCGGAACAGGCAAAAGGTGCCTTGGAGTCAGTCCTGGTAATCGAg                                                                                                    <  1:987175/150‑1 (MQ=255)
                                                         gCCCAGGGCCAGGGTGCCGATCGCGAGCAGTTGGTGCAGTGGGGGCATGGGATTGATCCTTCAATGAATGGGGCTCTCCTGTTCCGGCCAGTTCGCCGGAACAGGCAAAAGGTGCCTTGGAGTCAGTCCTGGTAATCGAGGCCCTCGCGa                                                                                          >  1:523663/1‑150 (MQ=255)
                                                          cccAGGGCCAGGGTGCCGATCGCGAGCAGTTGGTGCAGTGGGGGCATGGGATTGATCCTTCAATGAATGGGGCTCTCCTGTTCCGGCCAGTTCGCCGGAACAGGCAAAAGGTGCCTTGGAGTCAGTCCTGGTAATCGAGGCCCTCGCGAt                                                                                         <  2:62778/150‑1 (MQ=255)
                                                           ccagggccagggTGCCGATCGCGAGCAGTTGGTGCAGTGGGGGCATGGGATTGATCCTTCAATGAATGGGGCTCTCCTGTTCCGGCGAACTGGCCGGAACAg                                                                                                                                        >  1:1274112/1‑102 (MQ=38)
                                                           ccagggccagggTGCCGATCGCGAGCAGTTGGTGCAGTGGGGGCATGGGATTGATCCTTCAATGAATGGGGCTCTCCTGTTCCGGCGAACTGGCCGGAACAg                                                                                                                                        <  2:1274112/102‑1 (MQ=38)
                                                           ccagggccagggTGCCGATCGCGAGCAGTTGGTGCAGTGGGGGCATGGGATTGATCCTTCAATGAATGGGGCTCTCCTGTTACGGCCAGTTCGCCGGAACAGGCAACACGTGACTTAGCATCAGTCCTGGTAATCGAGGCCCTCGCGATc                                                                                        <  2:591235/150‑1 (MQ=255)
                                                            cagggccagggTGCCGATCGCGAGCAGTTGGTGCAGTGGGGGCATGGGATTGATCCTTCAATGAATGGGGCTCTCCTGTTCCGGCCAGTTCGCCGGAACAGGCAAAAGGTGCCTTGGAGGCAGTCCTGGTAAGCTATGCCCTCGCGATcc                                                                                       >  1:1002284/1‑150 (MQ=255)
                                                             agggccagggTGCCGATCGCGAGCAGTTGGTGCAGTGGGGGCATGGGATTGATCCTTCAATGAATGGGGCTCTCCTGTTCCGGCCAGTTCGCCGGAACAGGCAAAAGGTGCCTTGGAGTCAGTcc                                                                                                               >  2:298288/1‑125 (MQ=255)
                                                             agggccagggTGCCGATCGCGAGCAGTTGGTGCAGTGGGGGCATGGGATTGATCCTTCAATGAATGGGGCTCTCCTGTTCCGGCCAGTTCGCCGGAACAGGCAAAAGGTGCCTTGGAGTCAGTcc                                                                                                               <  1:298288/125‑1 (MQ=255)
                                                                          cGATCGCGAGCAGTTGGTGCAGTGGGGGCATGGGATTGATCCTTCAATGAATGGGGCTCTCCTGTTCCGGCCAGTTCGCCGGAACAGGCAAAAGGTGCCTTGGAGTCAGTCCTGGTAATCGAGGCCCTCGCGATcg                                                                                       >  1:116370/1‑135 (MQ=255)
                                                                          cGATCGCGAGCAGTTGGTGCAGTGGGGGCATGGGATTGATCCTTCAATGAATGGGGCTCTCCTGTTCCGGCCAGTTCGCCGGAACAGGCAAAAGGTGCCTTGGAGTCAGTCCTGGTAATCGAGGCCCTCGCGATcg                                                                                       <  2:116370/136‑2 (MQ=255)
                                                                                cgAGCAGTTGGTGCAGTGGGGGCATGGGATTGATCCTTCAATGAATGGGGCTCTCCTGTTCCGGCCAGTTCGCCGGAACAGGCAAAAGGTGCCTTGGAGTCAGTCCTGGTAATCGAGGCCCTCGCGATCCAACTGCCGAATCAGCGCATc                                                                   >  1:914385/1‑150 (MQ=255)
                                                                                   gCAGTTGGTGCAGTGGGGGCATGGGATTGATCCTTCAATGAATGGGGCTCTCCTGTTCCGGCCAGTTCGCCGGAACAGGCAAAAGGTGCCTTGGAGTCAGTCCTGGTAATCGAGGCCCTCGCGATCCAACTGCCGAATCAGCGCATCCCa                                                                >  2:1322539/1‑150 (MQ=255)
                                                                                     acttGGTGCAGTGGGGGCATGGGATTGATCCTTCAATGAATGTGGCTCTCCTTTTCCGGCCAGTTCGCCGGAACAGGCAAAAGGTGCCTTTGAGTCAGTCCTGGTAATCGAGGCCCTCGCGATCCAACTGCCGAATCAGCGCATCCCAAt                                                              <  2:1385771/148‑1 (MQ=255)
                                                                                     aGTTGGTGCAGTGGGGGCATGGGATTGATCCTTCAATGAATGGGGCTCTCCTGTTCCGGCCAGTTCGCCGGAACAGGCAAAAGGTGCCTTGGAGTCAGTCCTGGTAATCGAGGCCCTCGCGATCCAACTGCCGAATCAGCGCATCCCAAt                                                              <  2:523663/150‑1 (MQ=255)
                                                                                       ttGGTGCAGTGGGGGCATGGGATTGATCCTTCAATGAATGGGGCTCTCCTGTTCCGGCCAGTTCGCCGGAACAGGCAAAAGGTGCCTTGGAGTCAGTCCTGGTAATCGAGGCCCTCGCGATCCAACTGCCGAATCAGCGCATCCAAATGc                                                            >  2:1071195/1‑150 (MQ=255)
                                                                                            gCAGTGGGGGCATGGGATTGATCCTTCAATGAATGGGGCTCTCCTGTTCCGGCCAGTTCGCCGGAACAGGCAAAAGGTGCCTTGGAGTCAGTCCTGGTAATCGAGGCCCTCGCGATCCAACTGCCGAATCAGCGCATCCCAATGCCGTTg                                                       <  2:914385/150‑1 (MQ=255)
                                                                                              aGTGGGGTCATGGGATTGATCCTTCCATGAATGGGGCTCTCCTGTTCCGGCCAGTTCGCCGGAACAGGCAAAAGGTGCCTTGGAGCGACACTAAGAGCTCGATTGTCTCGCGATCCAACTGCCGAATCAGCGCATTCCAATGCCGTTGGa                                                     <  2:703202/150‑1 (MQ=255)
                                                                                                tGGGGGCATGGGATTGATCCTTCAATGAATGGGGCTCTCCTGTTCCGGCCAGTTCGCCGGAACAGGCAAAAGGTGCCTTGGAGTCAGTCCTgg                                                                                                            <  2:942110/93‑1 (MQ=255)
                                                                                                tGGGGGCATGGGATTGATCCTTCAATGAATGGGGCTCTCCTGTTCCGGCCAGTTCGCCGGAACAGGCAAAAGGTGCCTTGGAGTCAGTCCTgg                                                                                                            >  1:942110/1‑93 (MQ=255)
                                                                                                         gggATTGATCCTTCAATGAATGGGGCTCTCCTGTTCCGGCCAGTTCGCCGGAACAGGCAAAAGGTGCCTTGGAGTCAGTCCTgg                                                                                                            >  1:422367/1‑84 (MQ=255)
                                                                                                         gggATTGATCCTTCAATGAATGGGGCTCTCCTGTTCCGGCCAGTTCGCCGGAACAGGCAAAAGGTGCCTTGGAGTCAGTCCTgg                                                                                                            <  2:422367/84‑1 (MQ=255)
                                                                                                          ggATTGATCCTTCATTGAATGGGTCTCTCCTGTTCCGGCCAGTTCCCCGTAACAGGCAAAAGGTGCCTTGGAGTCAGTCCTAGTAATCGCGGCCCTCGCGATCCAACTGCCGAATCAGCGCATCCCAATGCCGTTGGATACCTTCCCCCt                                         <  2:129628/150‑1 (MQ=255)
                                                                                                                 cccTTCAATGAATGGGGCTCTCCTGTTCCGGCCAGTTCGCCGGAACAGGCAAAAGGTGCCTTGGAGTCAGTCCTGGTAATCGAGGCCCTCGCGATCCAACTGCCGAATCAGCGCATCCCAATGCCGTTGGATACCTTCCCCCTGGCCATc                                  <  1:333798/149‑1 (MQ=255)
                                                                                                                   cTTCAATGAATGGGGCTCTCCTGTTCCGGCCAGTTCGCCGGAACAGGCAAAAGGTGCCTTGGAGTCAGTCCTGGTAATCGAGGCCCTCGCGATCCAACTGCCGAATCAGCGCATCCCAATGCCGTTGGATACCTTCCCCCTGGCCATCgg                                >  1:118719/1‑150 (MQ=255)
                                                                                                                    ttCAATGAATGGGGCTCTCCTGTTCCGGCCAGTTCGCCGGAACAGGCAAAAGGTGCCTTGGAGTCAGTCCTGGTAATCGAGGCCCTCGCGATCCAACTGCCGAATCAGCGCATCCCAATGcc                                                           >  1:568748/1‑122 (MQ=255)
                                                                                                                    ttCAATGAATGGGGCTCTCCTGTTCCGGCCAGTTCGCCGGAACAGGCAAAAGGTGCCTTGGAGTCAGTCCTGGTAATCGAGGCCCTCGCGATCCAACTGCCGAATCAGCGCATCCCAATGcc                                                           <  2:568748/122‑1 (MQ=255)
                                                                                                                    ttCAATGAATGGGGCTCTCCTGTTCCGGCCAGTTCGCCGGAACAGGCAAAAGGTGCCTTGGAGTCAGTCCTGGTAATCGAGGCCCTCGCGATCCAACTGCCGAATCAGCGCATCCCAATGCCGTTGGATACCTTCCCCCTGGCCATCGGa                               >  1:729680/1‑150 (MQ=255)
                                                                                                                     tCAATGAATGGGGCTCTCCTGTTCCGGCCAGTTCGCCGGAACAGGCAAAAGGTGCCTTGGAGTCAGTCCTGGTAATCGAGGCCCTCGCGATCCAACTGCCGAATCAGCGCATCCCAATGCCGTTGGATACCTTCCCCCTGGCCATCGGaa                              >  2:196984/1‑150 (MQ=255)
                                                                                                                              ggggCTCTCCTGTTCCGGCCAGTTCGCCGGAACAGGCAAAAGGTGCCTTGGAGTCAGTCCTGGTAATCGAGGCCCTCGCGATCCAACTGCCGAATCAGCGCATCCCAATGCCGTTGGATACCTTCCCCCTGGCCATCGGAAAACCGCTTg                     <  1:935552/150‑1 (MQ=255)
                                                                                                                                ggCTCTCCTGTTCCGGCCAGTTCGCCGGAACAGGCAAAAGGTGCCTTGGAGTCAGTCCt                                                                                                              <  2:253400/59‑1 (MQ=255)
                                                                                                                                ggCTCTCCTGTTCCGGCCAGTTCGCCGGAACAGGCAAAAGGTGCCTTGGAGTCAGTCCt                                                                                                              >  1:253400/1‑59 (MQ=255)
                                                                                                                                ggCTCTCCTGTTCCGGCCAGTTCGCCGGAACAGGCAAAAGGTGCCTTGGAGTCAGTCCTGGTAATCGAGGCCCTCGCGATCCAACTGCCGAATCAGCGCATCCCAATGCCGTTGGATACCTTCCCCCTGGCCATCGGAAAACCGCTTGGa                   <  1:376506/150‑1 (MQ=255)
                                                                                                                                     tcCTGTTCCGGCCAGTTCGCCGGAACAGGCAAAAGGTGCCTTGGAGTCAGTCCTGGTAATCGAGGCCCTCGCGATCCAACTGCCGAATCAgc                                                                        <  1:230439/92‑1 (MQ=255)
                                                                                                                                     tcCTGTTCCGGCCAGTTCGCCGGAACAGGCAAAAGGTGCCTTGGAGTCAGTCCTGGTAATCGAGGCCCTCGCGATCCAACTGCCGAATCAgc                                                                        >  2:230439/1‑92 (MQ=255)
                                                                                                                                     tcCTGTTCCGGCCAGTTCGCCGGAACAGGCAAAAGGTGCCTTGGAGTCAGTCCTGGTAATCGAGGCCCTCGCGATCCAACTGCCGAATCAGCGCATCCCAATGCCGTTGGATACCTTCCCCCTGGCCATCGGAAAACCGCTTGGACTGCt              >  1:852956/1‑150 (MQ=255)
                                                                                                                                          ttCCGGCCAGTTCGCCGGAACAGGCAAAAGGTGCCTTGGAGTCAGTCCTGGTAATCGAGGCCCTCGCGATCCAACTGCCGAATCAGCGCATCCCTATGCCGTTGGATACCTTCCCCCTGGCCATCGGAAAACCGCTTGGACTGCTCATGc         >  1:729059/1‑150 (MQ=255)
                                                                                                                                               gCCAGTTCGCCGGAACAGGCAAAAGGTGCCTTGGAGTCAGTCCTGGTAATCGAGGCCCTCGCGATCCAACTGCCGAATCAGCGCATCCCAATGCCGTTGGATACCTTCCCCCTGGCCATCGGAAAACCGCTTGGACTGCTCATGCACCtt    >  2:1089428/1‑150 (MQ=255)
                                                                                                                                               gCCAGTTCGCCGGAACAGGCAAAAGGTGCCTTGGAGTCAGTCCTGGTAATCGAGGCCCTCGCGATCCAACTGCCGAATCAGCGCATCCCAATGCCGTTGGATACCTTCCCCCTGGCCATCGGAAAACCGCTTGGACTGCTCATGCACCtt    <  2:852956/150‑1 (MQ=255)
                                                                                                                                                 cAGTTCGCCGGAACAGGCAAAAGGTGCCTTGGAGTCAGTCCTGGTAATCGAGGCCCTCGCGATCCAACTGCCGAATCAGCGCATCCCAATGCCGTTGGATACCTTCCCCCTGGCCATCGGAAAACCGCTTGGACTGCTCATGCACCTTGa  >  2:348149/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
GGCGATGTTGAGGGTCTTGCTGGTTGAGGTGTTGTCGGCCGGCGAGCAAGCGGTGAGGCCCAGGGCCAGGGTGCCGATCGCGAGCAGTTGGTGCAGTGGGGGCATGGGATTGATCCTTCAATGAATGGGGCTCTCCTGTTCCGGCCAGTTCGCCGGAACAGGCAAAAGGTGCCTTGGAGTCAGTCCTGGTAATCGAGGCCCTCGCGATCCAACTGCCGAATCAGCGCATCCCAATGCCGTTGGATACCTTCCCCCTGGCCATCGGAAAACCGCTTGGACTGCTCATGCACCTTGA  >  NC_002947/3659804‑3660098

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