Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,338,447 T→C 7.0% intergenic (+23/+119) 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,4470TC 7.0% 196.0 / 2.8 67intergenic (+23/+119)PP_3812/PP_5632lipase/allo‑isoleucine metabolism protein
Reads supporting (aligned to +/- strand):  ref base T (32/30);  new base C (2/3);  total (34/33)
Fisher's exact test for biased strand distribution p-value = 6.73e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.00e-01

AGTAACCTCGGTGCGCTACAACGGAATGATTCATGACTATGGGCTGCTCAACCCGTTGAGCCAAGTGCCGGCAGTAAAGGCGGCGCTGCGGCAGGCGGGCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAAT  >  NC_002947/4338298‑4338589
                                                                                                                                                     |                                                                                                                                              
aGTAACCTCGGTGCGCTACAACGGAATGATTCATGACTATGGGCTGCTCAACCCGTTGAGCCAAGTGCCGGCAGTAAAGGCGGCGCTGCGGCAGGCGGGCAACGAACTTAAAGAGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCt                                                                                                                                                >  2:785363/1‑150 (MQ=255)
  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                                                                                                                                >  2:204635/1‑125 (MQ=255)
                                         ggCTGCTCAACCCGTTGAGCCAAGTGCCGGCAGTAAAGGCGGCGCTGCGGCAGGCGGGCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCt                                                                                                                                <  1:204635/125‑1 (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                                                                                                                       <  2:911388/105‑1 (MQ=255)
                                                                      gCAGTAAAGGCGGCGCTGCGGCAGGCGGGCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGc                                                                                                                       >  1:911388/1‑105 (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:788043/1‑150 (MQ=255)
                                                                                                                                         aaTAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTcc       >  1:930865/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)
                                                                                                                                                     |                                                                                                                                              
AGTAACCTCGGTGCGCTACAACGGAATGATTCATGACTATGGGCTGCTCAACCCGTTGAGCCAAGTGCCGGCAGTAAAGGCGGCGCTGCGGCAGGCGGGCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAAT  >  NC_002947/4338298‑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: