Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,338,448 G→A 7.3% intergenic (+24/+118) PP_3812 → / ← PP_5632 lipase/allo‑isoleucine metabolism protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,338,4480GA 7.3% 187.0 / 3.3 66intergenic (+24/+118)PP_3812/PP_5632lipase/allo‑isoleucine metabolism protein
Reads supporting (aligned to +/- strand):  ref base G (31/29);  new base A (2/3);  total (33/33)
Fisher's exact test for biased strand distribution p-value = 6.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.60e-01

TAACCTCGGTGCGCTACAACGGAATGATTCATGACTATGGGCTGCTCAACCCGTTGAGCCAAGTGCCGGCAGTAAAGGCGGCGCTGCGGCAGGCGGGCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAAT  >  NC_002947/4338300‑4338589
                                                                                                                                                    |                                                                                                                                             
tAACCTCGGTGCGCTACAACGGAATGATTCATGACTATGGGCTGCTCAACCCGTTGAGCCAAGTGCCGGCAGTAAAGGCGGCGCTGCGGCAGGCGGGCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTgg                                                                                                                                              >  2:977305/1‑150 (MQ=255)
    cTCGGTGCGCTACAACGGAATGATTCATGACTATGGGCTGCTCAACCCGTTGAGCCAAGTGCCGGCAGTAAAGGCGGCGCTGCGGCAGGCGGGCAACGAATTTCCCGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCgg                                                                                                                                          <  2:558568/150‑1 (MQ=255)
    cTCGGTGCGCTACAACGGAATGATTCATGACTATGGGCTGCTCAACCCGTTGAGCCAAGTGCCGGCAGTAAAGGCGGCGCTGCGGCAGGCGGGCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCgg                                                                                                                                          >  1:147022/1‑150 (MQ=255)
    cTCGGTGCGCTACAACGGAATGATTCATGACTATGGGCTGCTCAACCCGTTGAGCCAAGTGCCGGCAGTAAAGGCGGCGCTGCGGCAGGCGGGCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCCAGGCgg                                                                                                                                          <  2:767166/150‑1 (MQ=255)
      cGGTGCGCTACAACGGAATGATTCATGACTATGGGCTGCTCAACCCGTTGAGCCAAGTGCCGGCAGTAAAGGCGGCGCTGCGGCAGGCGGGCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGt                                                                                                                                        <  2:1303843/150‑1 (MQ=255)
      cGGTGCCCTACAACGGAATGATTCATGACTATGGGCTGCTCAACCCGTTGAGCCAAGTGCCGGCAGTAAAGGCGGCGCTGCGGCAGGCAGGCAACGAACTTAAAGTGCATCTGTAGTGAGGCGATGAGCAATAGACCCGCCTGTGCGGGt                                                                                                                                        <  2:15861/150‑1 (MQ=255)
             cTACAACGGAATGATTCATGACTATGGGCTGCTCAACCCGTTGAGCCAAGTGCCGGCAGTAAAGGCGGCGCTGCGGCAGGCGGGCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGc                                                                                                                                 <  2:851384/150‑1 (MQ=255)
              tACAACGGAATGATTCATGACTATGGGCTGCTCAACCCGTTGAGCCAAGTGCCGGCAGTAAAGGCGGCGCTGCGGCAGGCGGGCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCt                                                                                                                                >  2:774896/1‑150 (MQ=255)
                cAACGGAATGATTCATGACTATGGGCTGCTCAACCCGTTGAGCCAAGTGCCGGCAGTAAAGGCGGCGCTGCGGCAGGCGGGCCACGAACTTCAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGGCTTTTGCTTc                                                                                                                              >  1:207084/1‑150 (MQ=255)
                   cGGAATGATTCATGACTATGGGCTGCTCAACCCGTTGAGCCAAGTGCCGGCAGTAAAGGCGGCGCTGCGGCAGGCGGGCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTAt                                                                                                                           <  2:1319966/150‑1 (MQ=255)
                    ggAATGATTCATGACTATGGGCTGCTCAACCCGTTGAGCCAAGTGCCGGCAGTAAAGGCGGCGCTGCGGCAGGCGGGCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTAtt                                                                                                                          >  1:896655/1‑150 (MQ=255)
                       aTGATTCATGACTATGGGCTGCTCAACCCGTTGAGCCAAGTGCCGGCAGTAAAGGCGGCGCTGCGGCAGGCGGGCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGc                                                                                                                       <  2:104340/150‑1 (MQ=255)
                           ttCATGACTATGGGCTGCTCAACCCGTTGAGCCAAGTGCCGGCAGTAAAGGCGGCGCTGCGGCAGGCGGGCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGaa                                                                                                                   >  1:534705/1‑150 (MQ=255)
                           ttCATGACTATGGGCTGCTCAACCCGTTGAGCCAAGTGCCGGCAGTAAAGGCGGCGCTGCGGCAGGCGGGCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGaa                                                                                                                   >  2:319575/1‑150 (MQ=255)
                                 aCTATGGGCTGCTCAACCCGTTGAGCCAAGTGCCGGCAGTAAAGGCGGCGCTGCGGCAGGCGGGCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGcc                                                                                                             >  2:808110/1‑150 (MQ=255)
                                   tATGGGCTGCTCAACCCGTTGAGCCAAGTGCCGGCAGTAAAGGCGGCGCTGCGGCAGGCGGGCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCg                                                                                                           <  1:977305/150‑1 (MQ=255)
                                      gggCTGCTCAACCCGTTGAGCCAAGTGCCGGCAGTAAAGGCGGCGCTGCGGCAGGCGGGCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAg                                                                                                        <  2:238254/150‑1 (MQ=255)
                                       ggCTGCTCAACCCGTTGAGCCAAGTGCCGGCAGTAAAGGCGGCGCTGCGGCAGGCGGGCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCt                                                                                                                                <  1:204635/125‑1 (MQ=255)
                                       ggCTGCTCAACCCGTTGAGCCAAGTGCCGGCAGTAAAGGCGGCGCTGCGGCAGGCGGGCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCt                                                                                                                                >  2:204635/1‑125 (MQ=255)
                                           gctCAACGGGTTGAGCCAAGTGCCGGCAGTAAAGATGGCTCTGCGGCAGGCGGGCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTc                                                                                                   >  1:337687/1‑150 (MQ=255)
                                               aacgggTTGAGCCAAGTGCCGGCAGTAAAGGCGGCGCTGCGGCAGGCGGGCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGc                                                                                                                       <  1:564482/121‑1 (MQ=255)
                                               aacgggTTGAGCCAAGTGCCGGCAGTAAAGGCGGCGCTGCGGCAGGCGGGCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGc                                                                                                                       >  2:564482/6‑126 (MQ=255)
                                                  ccGTTGAGCCAAGTGCCGGCAGTAAAGGCGGCGCTGCGGCAGGCGGGCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAAt                                                                                            <  2:869066/150‑1 (MQ=255)
                                                                    gCAGTAAAGGCGGCGCTGCGGCAGGCGGGCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGc                                                                                                                       >  1:911388/1‑105 (MQ=255)
                                                                    gCAGTAAAGGCGGCGCTGCGGCAGGCGGGCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGc                                                                                                                       <  2:911388/105‑1 (MQ=255)
                                                                        tAAAGGCGGCCATCCGGCAGGCGGGCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTcc                                                                      >  2:916054/1‑150 (MQ=255)
                                                                              cggcgCTGCGGCAGGCGGGCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCt                                                                <  1:808110/150‑1 (MQ=255)
                                                                                     gCGGCAGGCGGGCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCgg                                                         <  2:896655/150‑1 (MQ=255)
                                                                                     ccGGCAGGCGGGCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCgg                                                         <  1:5108/149‑1 (MQ=255)
                                                                                       ggCAGGCGGGCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGt                                                       >  2:800629/1‑150 (MQ=255)
                                                                                        gCAGGCGGGCAACGAACTTAAAGTGCATATGCAGTGAGGGGCTAAAAACCCCACGCTTTTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGGGCGGCAGGTCCGGCCCTTTCGCGGGtt                                                      <  2:37088/150‑2 (MQ=255)
                                                                                               ggCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCcgc                                               <  1:774896/150‑1 (MQ=255)
                                                                                                   aCGAACTTAAAGTGCATCTGCAGTGAGGCGATGCGCCAACGCCCACAGTTAACGGGGGTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGaaa                                           <  1:319575/150‑1 (MQ=255)
                                                                                                   aCGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGaaa                                           >  2:1101922/1‑150 (MQ=255)
                                                                                                     gAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAgg                                         <  1:1101922/150‑1 (MQ=255)
                                                                                                       aCTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGc                                       >  2:1232283/1‑150 (MQ=255)
                                                                                                         ttAAAGTGCATCTGCAGTGAGGCGATGAGCAAAAGACCCGCCCAGGCGGGTCtat                                                                                                                                    <  1:1059927/55‑3 (MQ=18)
                                                                                                         ttAAAGTGCATCTGCAGTGAGGCGATGAGCAAAAGACCCGCCCAGGCGGGTCtat                                                                                                                                    >  2:1059927/1‑53 (MQ=18)
                                                                                                             aGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAgg                                                                         <  1:576002/110‑1 (MQ=255)
                                                                                                             aGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAgg                                                                         >  2:576002/1‑110 (MQ=255)
                                                                                                              gTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGc                                                                                                              >  2:1145085/1‑72 (MQ=255)
                                                                                                              gTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGc                                                                                                              <  1:1145085/72‑1 (MQ=255)
                                                                                                              gTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCaaa                                >  2:1178916/1‑150 (MQ=255)
                                                                                                                 cATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATcc                             >  2:622759/1‑150 (MQ=255)
                                                                                                                        tgTGAGGCGATGAGCAATAGACCCGCCTGGGCGGTTCTTTTGCTTCTATTTGCGGAAAGAGCCCGCAGGTTTCAATTTATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCCGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAg                      <  2:628399/149‑1 (MQ=255)
                                                                                                                        aGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAg                      >  2:1081996/1‑150 (MQ=255)
                                                                                                                        aGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTATTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAg                      <  1:785363/150‑1 (MQ=255)
                                                                                                                          ggAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTCTGGTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTACACCCGCGAAAGGGCCGGCAAATCCTGCTCAggg                    <  1:727595/149‑1 (MQ=255)
                                                                                                                             ggCGATGAGCAATAGACGCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGTGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTc                 >  2:729464/1‑150 (MQ=255)
                                                                                                                                gATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCAGCAGGCTTCAATGAATCTCTCATTGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAAACCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCaaa              >  1:173392/1‑150 (MQ=255)
                                                                                                                                 aTGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAAc             <  1:622759/150‑1 (MQ=255)
                                                                                                                                  gaaGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCACAGGCTTCAATGAATCTGCCa                                                                                      >  2:479999/3‑76 (MQ=255)
                                                                                                                                  gaaGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCACAGGCTTCAATGAATCTGCCa                                                                                      <  1:479999/74‑1 (MQ=255)
                                                                                                                                   gAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCCTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAAccc           >  2:841157/1‑150 (MQ=255)
                                                                                                                                    ggCAATAGACCCGCCTGGGCGGGTCTTTTGGTTCTATTTGCGGGAACAGCGCGCCGGTTTCAATGAGTCGGCCAGCGGCGCGGCAGGTCCGGCCCTTTCGTGGGGAAACCCGCGAAAGGGCCGGCAAATCCTGCGCAGGGCTCAAACCCt          <  2:207084/149‑1 (MQ=255)
                                                                                                                                     gCAAAAGACCCGCCCAGGCGGGTCTATTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGc                                                                               <  2:336055/80‑1 (MQ=25)
                                                                                                                                     gCAAAAGACCCGCCCAGGCGGGTCTATTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGc                                                                               >  1:336055/1‑80 (MQ=25)
                                                                                                                                      cAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTc        >  1:402183/1‑150 (MQ=255)
                                                                                                                                       aaTAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTcc       >  1:930865/1‑150 (MQ=255)
                                                                                                                                       aaTAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTcc       >  1:788043/1‑150 (MQ=255)
                                                                                                                                        aTAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCTGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAGATCCTGCTCAGGGCTCAAACCCTTCca      >  1:1241365/1‑150 (MQ=255)
                                                                                                                                            aCCCGCCTGGGCGGGTTTTTTGGTTCGATTTGCGGAAACAGCCCGCAGGCTTCAATGAACCTTCCAGATGTGCGGCAGGTCCCGCCCTTTCGCGGGTAAACCTGCGAAAGAGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAAt  <  2:623692/150‑1 (MQ=255)
                                                                                                                                            aCCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAAt  >  1:628601/1‑150 (MQ=255)
                                                                                                                                              ccGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCttt                                                              >  2:172403/1‑88 (MQ=255)
                                                                                                                                              ccGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCttt                                                              <  1:172403/88‑1 (MQ=255)
                                                                                                                                                gCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGc                                                                   <  1:1338038/81‑1 (MQ=255)
                                                                                                                                                gCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGc                                                                   >  2:1338038/1‑81 (MQ=255)
                                                                                                                                                    |                                                                                                                                             
TAACCTCGGTGCGCTACAACGGAATGATTCATGACTATGGGCTGCTCAACCCGTTGAGCCAAGTGCCGGCAGTAAAGGCGGCGCTGCGGCAGGCGGGCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAAT  >  NC_002947/4338300‑4338589

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