Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,338,437 T→A 7.6% intergenic (+13/+129) 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,4370TA 7.6% 187.7 / 4.5 65intergenic (+13/+129)PP_3812/PP_5632lipase/allo‑isoleucine metabolism protein
Reads supporting (aligned to +/- strand):  ref base T (32/27);  new base A (2/3);  total (34/31)
Fisher's exact test for biased strand distribution p-value = 6.59e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCGTGACAGTAACCTCGGTGCGCTACAACGGAATGATTCATGACTATGGGCTGCTCAACCCGTTGAGCCAAGTGCCGGCAGTAAAGGCGGCGCTGCGGCAGGCGGGCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCA  >  NC_002947/4338291‑4338585
                                                                                                                                                  |                                                                                                                                                    
gcgTGACAGTAACCTCGGTGCGCTACAACGGAATGATTCATGACTATGGGCTGCTCAACCCGTTGAGCCAAGTGCCGGCAGTAAAGGCGGCGCTGCGGCAGGCGGGCAACGAACTTAAAGTGTAGCTGCAGTGAGGCGATGAGCAATAGa                                                                                                                                                   <  1:1108339/150‑1 (MQ=255)
  gTGACAGTAACCTCGGTGCGCTACAACGGAATGATTCATGACTATGGGCTGCTCAACCCGTTGAGCCAAGTGCCGGCAGTAAAGGCGGCGCTGCGGCAGGCGGGCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGAcc                                                                                                                                                 >  1:628399/1‑150 (MQ=255)
     aCAGTAACCTCGGTGCGCTACAACGGAATGATTCATGACTATGGGCTGCTCAACCCGTTGAGCCAAGTGCCGGCAGTAAAGGCGGCGCTGCGGCAGGCGGGCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGc                                                                                                                                              >  1:1303843/1‑150 (MQ=255)
       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                                                                                                                            <  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                                                                                                                   <  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:930865/1‑150 (MQ=255)
                                                                                                                                                aaTAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTcc   >  1:788043/1‑150 (MQ=255)
                                                                                                                                                 aTAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCTGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAGATCCTGCTCAGGGCTCAAACCCTTCca  >  1:1241365/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
GCGTGACAGTAACCTCGGTGCGCTACAACGGAATGATTCATGACTATGGGCTGCTCAACCCGTTGAGCCAAGTGCCGGCAGTAAAGGCGGCGCTGCGGCAGGCGGGCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCA  >  NC_002947/4338291‑4338585

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: