Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,184,175 G→T 14.0% Q88H (CAG→CAT PP_4564 → lipoprotein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,184,1750GT14.0% 177.4 / 14.0 71Q88H (CAG→CATPP_4564lipoprotein
Reads supporting (aligned to +/- strand):  ref base G (32/29);  new base T (4/6);  total (36/35)
Fisher's exact test for biased strand distribution p-value = 5.14e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.35e-01

GCTGAGCGACTACCTGGACATGAACCGCGAGCAAAGGGCGCTGCTGGATGATCGGTTGCGCGAGCACCTGGCCTGGCATTGCAAGACTCAACTGCCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGACCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGTGAGGCAATTGGTCGATTGGCCGAGGAAATCACCCCTTCGGCCACCGAACTGCTGCGCGGCATGAGCGACAACCAGGTAGCGCAGATGCGTGAGGCCTTCCGCGA  >  NC_002947/5184028‑5184321
                                                                                                                                                   |                                                                                                                                                  
gCTGAGCGACTACCTGGACATGAACCGCGAGCAAAGGGCGCTGCTGGATGATCGGTTGCGCGAGCACCTGGCCTGGCATTGCAAGACTCAACTGCCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACCTGGTCGTCGGCGACCATGt                                                                                                                                                  <  2:1261292/150‑1 (MQ=255)
gCTGAGCGACTACCTGGACATGAACCGCGAGCAAAGGGCGCTGCTGGATGATCGGTTGCGCGAGCACCTGGCCTGGCATTGCAAGACTCAACTGCCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGt                                                                                                                                                  >  2:1135583/1‑150 (MQ=255)
 cTGAGCGACTACCTGGACATGAACCGCGAGCAAAGGGCGCTGCTGGATGATCGGTTGCGCGAGCACCTGGCCTGGCATTGCAAGACTCAACTGCCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACCTGGTCGTCGGCGACCATGTc                                                                                                                                                 <  2:1023789/150‑1 (MQ=255)
 cTGAGCGACTACCTGGACATGAACCGCGAGCAAAGGGCGCTGCTGGATGATCGGTTGCGCGAGCACCTGGCCTGGCATTGCAAGACTCAACTGCCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTc                                                                                                                                                 <  2:1321262/150‑1 (MQ=255)
  tGAGCGACTACCTGGACATGAACCGCGAGCAAAGGGCGCTGCTGGATGATCGGTTGCGCGAGCACCTGGCCTGGCATTGCAAGACTCAACTGCCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCa                                                                                                                                                >  1:400312/1‑150 (MQ=255)
        aCTACCTGGACATGAACCGCGAGCAAAGGGCGCTGCTGGATGATCGGTTGCGCGAGCACCTGGCCTGGCATTGCAAGACTCAACTGCCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGAcc                                                                                                                                          <  1:584320/150‑1 (MQ=255)
            ccTGGACATGAACCGCGAGCAAAGGGCGCTGCTGGATGATCGGTTGCGCGAGCACCTGGCCTGGCATTGCAAGACTCAACTGCCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGACCAGGc                                                                                                                                      <  1:1209226/150‑1 (MQ=255)
             cTGGACATGAACCGCGAGCAAAGGGCGCTGCTGGGTTATCGGTTGCGCGAGCACCTGGCCTGGCATTGCAAGACTCAACTGCCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGACCAGGcc                                                                                                                                     >  1:114457/1‑150 (MQ=255)
             cTGGACATGAACCGCGAGCAAAGGGCGCTGCTGGATGATCGGTTGCGCGAGCACCTGGCCTGGCATTGCAAGACTCAACTGCCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACCTGGTCGTCGGCGACCATGTCAc                                                                                                                                               <  1:1326493/140‑1 (MQ=255)
             cTGGACATGAACCGCGAGCAAAGGGCGCTGCTGGATGATCGGTTGCGCGAGCACCTGGCCTGGCATTGCAAGACTCAACTGCCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACCTGGTCGTCGGCGACCATGTCAc                                                                                                                                               >  2:1326493/1‑140 (MQ=255)
             cTGGACATGAACCGCGAGCAAAGGGCGCTGCTGGATGATCGGTTGCGCGAGCACCTGGCCTGGCATTGCAAGACTCAACTGCCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGACCAGGcc                                                                                                                                     >  2:798567/1‑150 (MQ=255)
              tGGACATGAACCGCGAGCAAAGGGCGCTGCTGGATGATCGGTTGCGCGAGCACCTGGCCTGGCATTGCAAGACTCAACTGCCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGACCAGGccc                                                                                                                                    >  2:820394/1‑150 (MQ=255)
                  cATGAACCGCGAGCAAAGGGCGCTGCTGGATGATCGGTTGCGCGAGCACCTGGCCTGGCATTGCAAGACTCAACTGCCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACCTGGTCGTCGGCGACCATGTc                                                                                                                                                 <  2:864039/133‑1 (MQ=255)
                  cATGAACCGCGAGCAAAGGGCGCTGCTGGATGATCGGTTGCGCGAGCACCTGGCCTGGCATTGCAAGACTCAACTGCCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACCTGGTCGTCGGCGACCATGTc                                                                                                                                                 >  1:864039/1‑133 (MQ=255)
                  cATGAACCGCGAGCAAAGGGCGCTGCTGGATGATCGGTTGCGCGAGCACCTGGCCTGGCATTGCAAGACTCAACTGCCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGACCAGGCCCTgca                                                                                                                                <  1:487243/150‑1 (MQ=255)
                       aCCGCGAGCAAAGGGCGCTGCTGGATGATCGGTTGCGCGAGCACCTGGCCGGGCATTGCAAGACTCAACTGCCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACCTGGTCGCCGACGACCAGGTCACCGACCAGGCCCTGCAGCaaa                                                                                                                           <  2:18520/150‑1 (MQ=255)
                          gcgAGCAAAGGGCGCTGCTGGATGATCGGTTGCGCGAGCACCTGGCCTGGCATTGCAAGACTCAACTGCCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGACCAGGCCCTg                                                                                                                                  <  2:486848/140‑1 (MQ=255)
                          gcgAGCAAAGGGCGCTGCTGGATGATCGGTTGCGCGAGCACCTGGCCTGGCATTGCAAGACTCAACTGCCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGACCAGGCCCTg                                                                                                                                  >  1:486848/1‑140 (MQ=255)
                           cgAGCAAAGGGCGCTGCTGGATGATCGGTTGCGCGAGCACCTGCCGTGGCATGGCAAGACTCAACTGCCCGGCTACCTCGACTGGCGCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGACCAGGCCCTGCAGCAAAGCAc                                                                                                                       <  2:425810/150‑1 (MQ=255)
                                 aaGGGCGCTGCTGGATGATCGGTTGCGCGAGCACCTGGCCTGGCATTGCAAGACTCAACTGCCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGACCAGGCCCTGCAGCAAAGCACCCGCGa                                                                                                                 <  1:820394/150‑1 (MQ=255)
                                        ctgctgGATGATCGGTTGCGCGAGCACCTGGCCTGGCATTGCAAGACTCAACTGCCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGACCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGt                                                                                                          <  2:400312/150‑1 (MQ=255)
                                                tgatCGGTTGCGCGAGCACCTGGCCTGGCATTGCAAGACTCAACTGCCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGACCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGTGAGGCAAt                                                                                                  <  2:723597/150‑1 (MQ=255)
                                                 gatCGGTTGCGCGAGCACCTGGCCTGGCATTGCAAGACTCAACTGCCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGACCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGTGAGGCAAtt                                                                                                 >  1:625969/1‑150 (MQ=255)
                                                                   ctggcctggcATTGCAAGACTCAACTGCCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGACCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGTGAGGCAATTGGTCGATTGGCCGAGGaa                                                                               >  1:372518/1‑150 (MQ=255)
                                                                   ctggcctggcATTGCAAGACTCAACTGCCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGACCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGTGAGGCAATTGGTCGATTGGCCGAGGaa                                                                               >  1:1002669/1‑150 (MQ=255)
                                                                         tggcATTGCAAGACTCAACTGCCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGACCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGTGAGGCAATTg                                                                                                <  2:1257127/127‑1 (MQ=255)
                                                                         tggcATTGCAAGACTCAACTGCCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGACCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGTGAGGCAATTg                                                                                                >  1:1257127/1‑127 (MQ=255)
                                                                         tggcATTGCAAGACTCAACTGCCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGACCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGTGAGGCAATTGGTCGATTGGCCGAGGAAATCAcc                                                                         >  2:973913/1‑150 (MQ=255)
                                                                          ggcATTGCAAGACTCAACTGCCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGACCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGTGAGGCAATTGGTCGATTGGCCGAGGAAATCAccc                                                                        >  2:1256789/1‑150 (MQ=255)
                                                                             aTTGCAAGACTCAACTGCCCGGCTACCNCGACTGGCTCGANCGGCTANGTGACATGGTCGCCGACGACCAGGTCACCGACCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGTGAGGCAATTGGTCGATTGGCCGAGGAAATCACCCCtt                                                                     >  1:1021207/1‑150 (MQ=255)
                                                                               tGCAAGACTCAACTGCCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGACCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGTGAGGCAATTGGTCGATTGGCCGAGGAAATCACCCCTTCg                                                                   <  1:1256789/150‑1 (MQ=255)
                                                                                 cAAGACTCAACTGCCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGACCAGGCCCTGCAGCaaa                                                                                                                           <  1:1331253/92‑1 (MQ=255)
                                                                                 cAAGACTCAACTGCCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGACCAGGCCCTGCAGCaaa                                                                                                                           >  2:1331253/1‑92 (MQ=255)
                                                                                 cAAGACTCAACTGCCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGAAGACCAGGTCACCGACCAGGCCCTGCAGAACAGAACCCGCGAGCTCAGTGCGGCAATGGGTCGTGTGGGGGCGGAGATCACTCCTTCGGc                                                                 >  2:588983/1‑150 (MQ=255)
                                                                                    gACTCAACTGCCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAg                                                                                                                                                    >  2:296636/1‑64 (MQ=255)
                                                                                    gACTCAACTGCCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAg                                                                                                                                                    <  1:296636/64‑1 (MQ=255)
                                                                                    gACTCAACTGCCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGACCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGTGAGGCAATTGGTCGATTGGCCGAGGAAATCACCCCTTCGGCCAc                                                              <  2:932408/150‑1 (MQ=255)
                                                                                     aCTCAACTGCCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGACCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGTGAGGCAATTGGTCGATTGGCCGAGGAAATCACCCCTTCGGCCAcc                                                             >  1:1249814/1‑150 (MQ=255)
                                                                                         aaCTGCCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGACCAGGc                                                                                                                                      >  1:1301033/1‑73 (MQ=255)
                                                                                         aaCTGCCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGACCAGGc                                                                                                                                      <  2:1301033/73‑1 (MQ=255)
                                                                                         aaCTGCCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGACCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGTGAGGCAATTGGTCGATTGGCCGAGGAAATCACCCCTTCGGCCACCGAAc                                                         <  2:1384402/150‑1 (MQ=255)
                                                                                            tGCCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGACCAGGCCCTGCAGCAAAGCACCcgc                                                                                                                   >  1:1057794/1‑89 (MQ=255)
                                                                                            tGCCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGACCAGGCCCTGCAGCAAAGCACCcgc                                                                                                                   <  2:1057794/89‑1 (MQ=255)
                                                                                              cccGGCTACCTCGACGGTCCCGACCGGCTACCTGACATGGTCGCCGACGACCAGGTCACCGACCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGTGAGGCAATTGGTCGATTGGCCGAGGAAATCACCCCTTCGGCCACCGAActgctg                                                    >  2:414442/1‑150 (MQ=255)
                                                                                               ccGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGACCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGTGAGGCAATTGGTCGATTGGCCGAGGAAATCACCCCTTCGGCCACCGAACtgctgc                                                   <  1:973913/150‑1 (MQ=255)
                                                                                                      ccTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGACCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGTGAGGCAATTGGTCGATTGGCCGAGGAAATCACCCCTTCGGCCACCGAACTGCTGCGCGGCAt                                            >  2:601132/1‑150 (MQ=255)
                                                                                                       cTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGACCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGTGAGGCAATTGGTCGATTGGCCGAg                                                                                  <  1:1255633/111‑1 (MQ=255)
                                                                                                       cTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGACCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGTGAGGCAATTGGTCGATTGGCCGAg                                                                                  >  2:1255633/1‑111 (MQ=255)
                                                                                                       cTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGACCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGTGAGGCAATTGGTCGATTGGCCGAGGAAATCACCCCTTCGGCCACCGAACTGCTGCGCGGCATg                                           >  2:1247781/1‑150 (MQ=255)
                                                                                                                cTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGACCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGTGAGGCAATTGGTCGATTGGCCGAGGAAATCACCCCTTCGGCCACCGAACTGCTGCGCGGCATGAGCGACAAc                                  >  1:1226554/1‑150 (MQ=255)
                                                                                                                   gACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGACCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGTGAGGCAATTGGTCGATTGGCCGAGGAAATCACCCCTTCGGCCACCGaa                                                          <  2:1158086/123‑1 (MQ=255)
                                                                                                                   gACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGACCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGTGAGGCAATTGGTCGATTGGCCGAGGAAATCACCCCTTCGGCCACCGAACTGCTGCGCGGCATGAGCGACAACCAg                               <  2:372518/150‑1 (MQ=255)
                                                                                                                   gACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGACCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGTGAGGCAATGGGTCGATTGGCCGAGGAAATCACCCCTTCGGCCACCGaa                                                          >  1:1158086/1‑123 (MQ=255)
                                                                                                                    aCCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGACCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGTGAGGCAATTGGTCGATTGGCCGAGGAAATCACCCCTTCGGCCACCGAACTGCTGCGCGGCATGAGCGACAACCAgg                              <  1:503125/150‑1 (MQ=255)
                                                                                                                       ggCTACGTGACATGGTCGCCGACGACCAGGTCACCGACCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGTGAGGCAATTGGTCGATTGGCCGAGGAAATCACCCCTTCGGCCACCGAACTGCTGCGCGGCATGAGCGACAACCAGGTAg                           >  2:782140/1‑150 (MQ=255)
                                                                                                                           aCGTGACATGGTCGCCGACGACCAGGTCACCGACCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGTGAGGCAATTGGTCGATTGGCCGAGGAAATCACCCCTTCGGCCACCGAACTGCTGCGCGGCATGAGCGACAACCAGGTAGCGCa                       >  2:1082562/1‑150 (MQ=255)
                                                                                                                             gTGACCTGGTCGTCGGCGACCATGTCACCGACCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGTGAGGCAATTGGTc                                                                                             <  1:541372/78‑1 (MQ=37)
                                                                                                                             gTGACCTGGTCGTCGGCGACCATGTCACCGACCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGTGAGGCAATTGGTc                                                                                             >  2:541372/1‑78 (MQ=25)
                                                                                                                               gACCTGGTCGTCGGCGACCATGTCACCGACCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGTGAGGCAATTGGTCGATTGGCCGAGGAAATCACCCCTTCGGCCACCGAACTGCTGCGCGGCATGAGCGACAACCAGGTAGCGCAGATg                   >  2:149081/1‑150 (MQ=255)
                                                                                                                               gACATGGTCGCCGACGACCAGGTCACCGACCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGTGAGGCAAt                                                                                                  <  1:668999/71‑1 (MQ=255)
                                                                                                                               gACATGGTCGCCGACGACCAGGTCACCGACCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGTGAGGCAAt                                                                                                  >  2:668999/1‑71 (MQ=255)
                                                                                                                                 cATGGTCGCCGACGACCAGGTCACCGACCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGTGAGGCAATTGGTCGATTGGCCGAGGAAATCACCCCTTCGGCCACCGAACTGCTGCGCGGCATGAGCGACAACCAg                               <  1:779529/136‑1 (MQ=255)
                                                                                                                                 cATGGTCGCCGACGACCAGGTCACCGACCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGTGAGGCAATTGGTCGATTGGCCGAGGAAATCACCCCTTCGGCCACCGAACTGCTGCGCGGCATGAGCGACAACCAg                               >  2:779529/1‑136 (MQ=255)
                                                                                                                                  aTGGTCGCCGACGACCAGGTCACCGACCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGTGAGGCAATTGGTCGATTGGCCGAGGAAATCACCCCTTCGGCCACCGAACTGCTGCGCGGCATGAGCGACAACCAGGTAGCGCa                       >  2:1128503/1‑143 (MQ=255)
                                                                                                                                  aTGGTCGCCGACGACCAGGTCACCGACCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGTGAGGCAATTGGTCGATTGGCCGAGGAAATCACCCCTTCGGCCACCGAACTGCTGCGCGGCATGAGCGACAACCAGGTAGCGCa                       <  1:1128503/143‑1 (MQ=255)
                                                                                                                                  aTGGTCGCCGACGACCAGGTCACCGACCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGTGAGGCAATTGGTCGATTGGCCGAGGAAATCACCCCTTCGGCCACCGAACTGCTGCGCGGCATGAGCGACAACCAGGTAGCGCAGATGCGt                <  2:1249814/150‑1 (MQ=255)
                                                                                                                                  aTGGTCGCCGACGACCAGGTCACCGACCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGTGAGGCAATTGGTCGATTGGCCGAGGAAATCACCCCTTCGGCCACCGAACTGCTGCGCGGCATGAGCGACAACCAGGTAGCGCAGATGCGt                <  1:577703/150‑1 (MQ=255)
                                                                                                                                    ggTCGCCGACGACCAGGTCACCGACCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGTGAGGCAATTGGTCGATTGGCCGAGGAAATCACCCCTTCGGCCACCGAACTGCTGCGCGGCATGAGCGACAACCAGGTAGCGCAGATGCGTGa              <  1:414442/150‑1 (MQ=255)
                                                                                                                                          cgacgaCCAGGTCACCGACCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGTGAGGCAATTGGTCGATTGGCCGAGGAAATCACCCCTTCGGCCACCGAACTGCTGCGCGGCATGAGCGACACCCAGGTAGCGCAGATGCGTGAGGCCtt        >  1:536865/1‑150 (MQ=255)
                                                                                                                                          cgacgaCCAGGTCACCGACCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGTGAGGCAATTGGTCGATTGGCCGAGGAAATCACCCCTTCGGCCACCGAACTGCTGCGCGGCATGAGCGACAACCAGGTAGCGCAGATGCGTGAGGCCtt        >  1:626798/1‑150 (MQ=255)
                                                                                                                                            gcgacCATGTCACCGACCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGTGAGGCAATTGGTCGATTGGCCGAGGAAATCACCCCTTCGGCCACCGAACTGCTGCGCGGCATGAGCGACAACCAGGTAGCGCAGATGCGTGAGGCCTTcc      <  1:149081/149‑1 (MQ=255)
                                                                                                                                             cgacCAGGTCAACGCCCGGGCCCTGCAGCCCAATATCCGCCTGTTCAATGAGTCGATTTTGCGGGTTGCCGAGGAAATCACCCCTGCGGCCACCGAACTGCTGCGCGGCATGAGCGACAACCAGGTAGCGCAGATGCGTGAGGCCTTCcg     <  2:1021207/150‑1 (MQ=255)
                                                                                                                                                ccAGGTCACCGACCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGTGAGGCAATTGGTCGATTGGCCGAGGAAATCACCCCTTCGGCCACCGAACTGCTGCGCGGCATGAGCGACAACCAGGTAGCGCAGATGCGTGAGGCCTTCCGCGa  >  2:594331/1‑150 (MQ=255)
                                                                                                                                                   ggTCACCGACCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGTGAGGCAATTGGTCGATTGGCCGAGGAAATCACCCCTTCGGCCACCGAACTGCTGCGCGGCAt                                            >  2:190379/1‑105 (MQ=255)
                                                                                                                                                   ggTCACCGACCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGTGAGGCAATTGGTCGATTGGCCGAGGAAATCACCCCTTCGGCCACCGAACTGCTGCGCGGCAt                                            <  1:190379/105‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
GCTGAGCGACTACCTGGACATGAACCGCGAGCAAAGGGCGCTGCTGGATGATCGGTTGCGCGAGCACCTGGCCTGGCATTGCAAGACTCAACTGCCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGACCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGTGAGGCAATTGGTCGATTGGCCGAGGAAATCACCCCTTCGGCCACCGAACTGCTGCGCGGCATGAGCGACAACCAGGTAGCGCAGATGCGTGAGGCCTTCCGCGA  >  NC_002947/5184028‑5184321

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