Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,163,250 T→A 6.1% R86W (AGG→TGG)  pssA ← CDP‑diacylglycerol‑‑serine O‑phosphatidyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,163,2500TA 6.1% 341.5 / 8.7 110R86W (AGG→TGG) pssACDP‑diacylglycerol‑‑serine O‑phosphatidyltransferase
Reads supporting (aligned to +/- strand):  ref base T (54/49);  new base A (4/3);  total (58/52)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

Only 100 of 110 total aligned reads displayed.

GCTGCCTTTGAGGTGCAGCACCCCGAACAGTTCGCGGGTCTGCACTGGCACGCCATGAATCGCGATGTCCAGCCCGCGCAACTGGCGCTGGGCCTGATACCAGGCGGCATTGCCCGGCTGACGCCCGGCCCCCAGCAGCCCGCGCCTGGCGCGGAACCAGTCGACCACGATGGTGATCTCCAGGCCTGGCCGCTCGGCCTTGGCCTGGTACAAGGCATCCAGGACTTCCTGGCCAGCCTCATCCTCCTGCAGGTACAACGCGACGATAACGATGCGCCGAGTCGCCGCACGGAT  >  NC_002947/4163104‑4163397
                                                                                                                                                   |                                                                                                                                                   
gctgcCTTTGAGGTGCAGCACCCCGAACAGTTCGCGGGTCTGCACTGGCACGCCATGAATCGCGATGTCCAGCCCGCGCAACTGGCGCTGGGCCTGATACCAGGCGGCATTGCCCGGCTGACGCCCGGCCCCCAGCAGCCCGCGCCTGgc                                                                                                                                                  >  2:1202355/1‑150 (MQ=255)
  tgcCTTTGAGGTGCAGCACCCCGAACAGTTCGCGGGTCTGCACTGGCACGCCATGAATCGCGATGTCCAGCCCGCGCAACTGGCGCTGGGCCTGATACCAGGCGGCATTGCCCGGCTGACGCCCGGCCCCCAGCAGCCCGCGCCTGgcgc                                                                                                                                                >  2:345823/1‑150 (MQ=255)
   gcCTTTGAGGTGCAGCACCCCGAACAGTTCGCGGGTCTGCACTGGCACGCCATGAATCGCGATGTCCAGCCCGCGCAACTGGCGCTGGGCCTGATACCAGGCGGCATTGCCCGGCTGACGCCCGGCCCCCAGCAGCCCGCGCCTGGcgcg                                                                                                                                               >  1:906956/1‑150 (MQ=255)
     cTTTGAGGTGCAGCACCCCGAACAGTGCGCGGGTCTGCACCGGCACGCCATGAATCGCGATGTCCAGCCCGCGCAACTGGCGCTGGGCCTGATACCAGGCGGCATTGCCCGGCTGACGCCCGGCCCCCAGCAGCCCGCGCCAGGCGCggg                                                                                                                                             <  2:544981/150‑2 (MQ=255)
      tttGAGGTGCAGCACCCCGAACAGTTCGCGGGTCTGGACTGGCACGCCATGAATCGCGATGTCCAGCCCGCGCAACTGGCGCTGGGCCTGATACCAGGCGGCATTGCCCGGCTGACGCCCGGCCCCCAGCAACCCGCGCCTGGCGCGGaa                                                                                                                                            >  1:782982/1‑150 (MQ=255)
      tttGAGGTGCAGCACCCCGAACAGTTCGCGGGTCTGCACTGGCACGCCATGAATCGCGATGTCCAGCCCGCGCAACTGGCGCTGGGCCTGATACCAGGCGGCATTGCCCGGCTGACGCCCGGCCCCCAGCAGCCCGCGCCTGGCGCGGaa                                                                                                                                            >  2:682714/1‑150 (MQ=255)
           ggTGGAGCACCCCGAATAGTTGGCGGGCATGCACTGGCACGCCATGAATCGCGGTGGCCAGCCCGCGCAACTGGCGCTGGGCCTGATACCAGGCGGCATTGCCCGGCTGACGCACGGCCCACAGCAGCCCGCGCCTGGCGCGGAAACAGt                                                                                                                                       >  2:375904/1‑150 (MQ=255)
            gTGCAGACCCGCGAACTGTTCGGGGTGCTGCACTGGCACGCCATGAATCGCGATGTCCAGCCCGCGCAACTGGCGCTGGGCCCGATACCAGGCGGCATTGCCCGGCTGACGCCCGGCCCCCAGCAGCCCGCGCCTGGCGCGGAACCAg                                                                                                                                        <  1:517397/148‑1 (MQ=255)
              gcagcaCCCCGAACAGTTCGCGGGTCTGCACTGGCACGCCATGAATCGCGATGTCCAGCCCGCGCAACTGGCGCTGGGCCTGATACCAGGCGGCATTGCCCGGCTGACGCCCGGCCCCCCGCAGCCCGCGCCTGGCGCGGAACCAGTCGa                                                                                                                                    >  2:864783/1‑150 (MQ=255)
                 gcaCGCCGCACAGCTCGCGGGTCTGCCCTGGCACGCCATGAATGGCGAGGGCCAGCCCGCGCAACTGGCGATGGGTCTGATACCAGCCGGCATTGCCCGGCTGACGCCCGGCCCCCAGCAGCCCGCGCCTGCAGCGCAACCACTCGACCa                                                                                                                                 >  2:89867/1‑150 (MQ=255)
                 gacccgcGAACTGTTCGGGGTGCTGCACTGGCACGCCATGAATCGCGATGTCCAGCCCGCGCAACTGGCGCTGGGCCTGATACCAGGCGGCATTGCCCGGCTGACGCCCGGCCCCCAGCAGCCCGCGCCTGGCGCGGAACCAGTCGACCa                                                                                                                                 >  2:110173/7‑150 (MQ=255)
                    ccccGAACAGTTCGCGGGTCTGCACTGGCACGCCATGAATCGCGATGTCCAGCCCGCGCAACTGGCGCTGGGCCTGATACCAGGCGGCATTGCCCGGCTGACGCCCGGCCCCCAGCAGCCCGCGCCTGGCGCGGAACCAGTCGACCACGa                                                                                                                              >  2:636144/1‑150 (MQ=255)
                     caccaaCAGTTCGCGGGCCTAAACTGGCACGACATGAATCTCGACGTGCAGCCCACGCAACTGGCGCTGGGCCTGATACCAGGCGGCATTGCCCGGCTGACGCCCGGCCCCCAGCAGCCCGCGCCTGGCGCGGAACCAGTCGACCACGAt                                                                                                                             >  2:524639/5‑150 (MQ=255)
                      ccGAACAGTTCGCGGGTCTGCACTGGCACGCCATGAATCGCGATGTCCAGCCCGCGCAACTGGCGCTGGGCCTGATACCAGGCGGCATTGCCCGGCTGACGCCCGGCCCCCAGCAGCCCGCGCCTGGCGCGGAACCAGTCGACCACGATg                                                                                                                            <  2:1101664/150‑1 (MQ=255)
                      ccGAACAGTTCGCGGGTCTGCACTGGCACGCCATGAATCGCGATGTCCAGCCCGCGCAACTGGCGCTGGGCCTGATACCAGGCGGCATTGCCCGGCTGACGCCCGGCCCCCAGCAGCCCGCGCCTGGCGCGGAACCAGTCGACCACGATg                                                                                                                            <  2:804733/150‑1 (MQ=255)
                      ccGAACAGTTCGCGGGTCTGCACTGGCACGCCATGAATCGCGATGTCCAGCCCGCGCAACTGGCGCTGGGCCTGATACCAGGCGGCATTGCCCGGCTGACGCCCGGCCCCCAGCAGCCCGCGCCTGGCGCGGAACCAGTCGACCACGATg                                                                                                                            <  1:636144/150‑1 (MQ=255)
                      ccGAACAGTTCGCGGGTCTGCACTGGCACGCCATGAATCGCGATGTCCAGCCCGCGCAACTGGCGCTGGGCCTGATACCAGGCGGCATTGCCCGGCTGACGCCCGGCCCCCAGCAGCCCGCGCCTGGCGCGGAACCAGTCGACCACGATg                                                                                                                            >  2:624073/1‑150 (MQ=255)
                      ccGAACAGTTCGCGGGTCTGCACTGGCACGCCATGAATCGCGATGTCCAGCCCGCGCAACTGGCGCTGGGCCTGATACCAGGCGGCATTGCCCGGCTGACGCCCGGCCCCCAGCAGCCCGCGCCTGGCGCGGAACCAGTCGACCACGATg                                                                                                                            <  1:302753/150‑1 (MQ=255)
                       cGAACAGTTCGCGGGTCTGCACTGGCACGCCATGAATCGCGATGGCCAGCCCGCGCAACTGGCGCTGGGCCTGATACCAGGCGGCATTGCCCGGCTGACGGCCGGCCCCCAGCAGCCCGCGCCTGGCGCGGAACCAGTCGACCACGATgg                                                                                                                           >  1:1247813/1‑150 (MQ=255)
                        gAACAGTTCGCGGGTCTGCACTGGCACGCCATGAATCGCGATGTCCAGCCCGCGCAACTGGCGCTGGGCCTGATACCAGGCGGCATTGCCCGGCTGACGCCCGGCCCCAGCAGCCCGTGCCTGTCGAGTAAGCAATCGGCCACAATGGTg                                                                                                                         >  2:1017015/1‑150 (MQ=255)
                           cAGTTCGCGGGTCTGCACTGGCACGCCATGAATCGCGATGTCCAGCCCGCGCAACTGGCGCTGGGCCTGATACCAGGCGGCATTTCCCGGCTGACGCCCGGCCCCCAGCAGCCCGCGCCTGGCGCGGAACCAGTCGACCACGATGGTGAt                                                                                                                       <  2:1025285/150‑1 (MQ=255)
                             gTTCGCGGGTCTGCACTGGCACGCCATGAATCGCGATGTCCAGCCCGCGCAACTGGCGCTGGGCCTGATACCAGGCGGCATTGCCCGGCTGACGCCCGGCCCCCAGCAGCCCGCGCCTGGCGCGGAACCAGTCGACCACGATGGTGAtct                                                                                                                     <  1:587327/150‑1 (MQ=255)
                              ttCGCGGGTCTGCACTGGCACGCCATGAATCGCGATGTCCAGCCCGCGCAACTGGCGCTGGGCCTGATACCAGGCGGCATTGCCCGGCTGACGCCCGGCCCCCAGCAGCCCGCGCCTGGCGCGGAACCAGTCGACCACGATgg                                                                                                                           >  1:645500/1‑143 (MQ=255)
                              ttCGCGGGTCTGCACTGGCACGCCATGAATCGCGATGTCCAGCCCGCGCAACTGGCGCTGGGCCTGATACCAGGCGGCATTGCCCGGCTGACGCCCGGCCCCCAGCAGCCCGCGCCTGGCGCGGAACCAGTCGACCACGATgg                                                                                                                           <  2:645500/143‑1 (MQ=255)
                              ttCGCGGGTCTGCACTGGCACGCCATGAATCGCGATGTCCAGCCCGCGCAACTGGCGCTGGGCCTGATACCAGGCGGCATTGCCCGGCTGACGCCCGGCCCCCAGCAGCCCGCGCCAGGCGCgggc                                                                                                                                            >  2:1362039/1‑124 (MQ=255)
                              ttCGCGGGTCTGCACTGGCACGCCATGAATCGCGATGTCCAGCCCGCGCAACTGGCGCTGGGCCTGATACCAGGCGGCATTGCCCGGCTGACGCCCGGCCCCCAGCAGCCCGCGCCAGGCGCgggc                                                                                                                                            <  1:1362039/126‑3 (MQ=255)
                              ttCGCGGGTCTGCACTGGCACGCCATGAATCGCGATGTCCAGCCCGCGCAACTGGCGCTGGGCCTGATACCAGGCGGCATTGCCCGGCTGACGCCCGGCCCCCAGCAGCCCGCGCCAGGCGCgg                                                                                                                                              >  1:450285/1‑124 (MQ=255)
                              ttCGCGGGTCTGCACTGGCACGCCATGAATCGCGATGTCCAGCCCGCGCAACTGGCGCTGGGCCTGATACCAGGCGGCATTGCCCGGCTGACGCCCGGCCCCCAGCAGCCCGCGCCAGGCGCgg                                                                                                                                              <  2:450285/124‑1 (MQ=255)
                              gtCGCGGGTGTGCACTGGCACGCCATGAATCGCGATGTCCAGCCCGCGCAACTGGCGCTGGGCCTGATACCAGGCGGCATTGCCCGGCTGACGCCCGGCCCCNAGCAGCCCGCGCCTGGCGNGGAACCAGTCGACCACGATGGTGAtctn                                                                                                                    <  1:89867/149‑2 (MQ=255)
                                cgcgGGTCTGCACTGGCACGCCATGAATCGCGATGTCCAGCCCGCGCAACTGGCGCTGGGCCGGATACCAGGCGGCATTGCCCGGCTGACGCCCGGCCCCCAGCAGCCCGCGCCTGGCGCGGAACCAGTCGACCACGATGGTGATCTCCa                                                                                                                  <  1:375904/150‑1 (MQ=255)
                                   gggTCTGCACTGGCACGCCATGAATCGCGATGTCCAGCCCGCGCAACTGGCGCTGGGCCTGATACCAGGCGGCATTGCCCGGCTGACGCCCGGCCCCCAGCAGCCCGCGCCTGGCGCGGAACCAGTCGACCACGATGGTGAtctc                                                                                                                    <  1:1378035/145‑1 (MQ=255)
                                   gggTCTGCACTGGCACGCCATGAATCGCGATGTCCAGCCCGCGCAACTGGCGCTGGGCCTGATACCAGGCGGCATTGCCCGGCTGACGCCCGGCCCCCAGCAGCCCGCGCCTGGCGCGGAACCAGTCGACCACGATGGTGAtctc                                                                                                                    >  2:1378035/1‑145 (MQ=255)
                                     gTCTGCACTGGCACGCCATGAATCGCGATGTCCAGCCCGCGCAACTGGCGCTGGGCCTGATACCAGGCGGCATTGCCCGGCTGACGCCCGGCCCCCAGCAGCCCTCGCCTGGCGCGGAACCAGTCGACCACGAt                                                                                                                             >  2:1038124/1‑134 (MQ=255)
                                     gTCTGCACTGGCACGCCATGAATCGCGATGTCCAGCCCGCGCAACTGGCGCTGGGCCTGATACCAGGCGGCATTGCCCGGCTGACGCCCGGCCCCCAGCAGCCCTCGCCTGGCGCGGAACCAGTCGACCACGAt                                                                                                                             <  1:1038124/134‑1 (MQ=255)
                                     gTCTGCACTGGCACGCCATGAATCGCGATGTCCAGCCCGCGCAACTGGCGCTGGGCCTGATACCAGGCGGCATTGCCCGGCTGACGCCCGGCCCCCAGCAGCCCGCGCCTGGCGCGGAACCAGTCGACCACGATGGTGATCTCCAg                                                                                                                 <  1:586275/146‑1 (MQ=255)
                                     gTCTGCACTGGCACGCCATGAATCGCGATGTCCAGCCCGCGCAACTGGCGCTGGGCCTGATACCAGGCGGCATTGCCCGGCTGACGCCCGGCCCCCAGCAGCCCGCGCCTGGCGCGGAACCAGTCGACCACGATGGTGATCTCCAg                                                                                                                 >  2:586275/1‑146 (MQ=255)
                                                gCACGCCATGAATCGCGATGTCCAGCCCGCGCAACTGGCGCTGGGCCTGATACCAGGCGGCATTGCCCGGCTGACGCCCGGCCCCCAGCAGCCCGCGCCTGGCGCGGAACCAGTCGACCACGATGGTGATCTCCAGGCCTGGCCGCTCgg                                                                                                   >  1:76670/1‑150 (MQ=255)
                                                  aCGCCATGAATCGCGATGTCCAGCCCGCGCAACTGGCGCTGGGCCTGATACCAGTCGGCATTGCCCGGCTGACGCCCGGCCCCCAGCAGCCCGCGCCTGGCGCGGAACCAGTCGACCACGATGGTGATCTCCAGGCCTGGCCGCTCGGcc                                                                                                 >  2:910043/1‑150 (MQ=255)
                                                  aCGCCATGAATCGCGATGTCCAGCCCGCGCAACTGGCGCTGGGCCTGATACCAGGCGGCATTGCCCGGCTGACGCCCGGCCCCCAGCAGCCCGCGCCTGGCGCGGa                                                                                                                                             <  2:442100/106‑1 (MQ=255)
                                                  aCGCCATGAATCGCGATGTCCAGCCCGCGCAACTGGCGCTGGGCCTGATACCAGGCGGCATTGCCCGGCTGACGCCCGGCCCCCAGCAGCCCGCGCCTGGCGCGGa                                                                                                                                             >  1:442100/1‑106 (MQ=255)
                                                    gCCATGAATCGCGATGTCCAGCCCGCGCAACTGGCGCTGGGCCTGATACCAGGCGGCATTGCCCGGCTGACGCCCGGCCCCCAGCAGCCCGCGCCTGGCGCGGAACCAGTCGACCACGATGGTGATCTCCAGGc                                                                                                               >  1:350199/1‑134 (MQ=255)
                                                      cATGAATCGCGATGTCCAGCCCGCGCAACTGGCGCTGGGCCTGATACCAGGCGGCATTGCCCGGCTGACGCCCGGCCACCAGCAGCCCGCGCCTGGCGCGGAACCAGTCGACCACGATGGTGATCTCCAGGCCTGGCCGCTCGGCCTTgg                                                                                             <  1:649068/150‑1 (MQ=255)
                                                         gAATCGCGATGTCCAGCCCGCGCAACTGGCGCTGGGCCTGATACCAGGCGGCATTGCCCGGCTGACGCCCGGCCCCCAGCAGCCCGCGCCTGGCGCGGAACCAGTCGACCACGATGGTGATCTCCAGGCCTGGCCGCTCGGCCTTGGCCt                                                                                          <  1:1332578/150‑1 (MQ=255)
                                                         gAATCGCGATGTCCAGCCCGCGCAACTGGCGCTGGGCCTGATACCAGGCGGCATTGCCCGGCTGACGCCCGGCCCCCAGCAGCCCGCGCCTGGCGCGGAACCAGTCGACCACGATGGTGATCTCCAGGCCTGGCCACTCGGCCTTGGCCt                                                                                          <  2:780488/150‑1 (MQ=255)
                                                              gcgATGTCCAGCCCGCGCAACTGGCGCTGGGCCTGATACCAGGCGGCATTGCCCGGCTGACGCCCGGCCCCCAGCAGCCCGCGCCTGGCGCGGAACCAGTCGACCACGATg                                                                                                                            <  1:766741/111‑1 (MQ=255)
                                                              gcgATGTCCAGCCCGCGCAACTGGCGCTGGGCCTGATACCAGGCGGCATTGCCCGGCTGACGCCCGGCCCCCAGCAGCCCGCGCCTGGCGCGGAACCAGTCGACCACGATg                                                                                                                            >  2:766741/1‑111 (MQ=255)
                                                                gATGTCCAGCCCGCGCAACTGGCGCTGGGCCTGACACCAGGCGGCATTGCCCGGCTGACGCCCGGCCCCAGCAGCCCGCGCCTGGCGCGGAACCAGTCGACCACGATGGTGATCTCCAGGCCTGGCCGCTCGGCCTTGGCCTGGTACAAg                                                                                  <  1:1017015/150‑1 (MQ=255)
                                                                    tCCAGCCCGCGCAACTGGCGCTGGGCCTGATACCAGGCGGCATTGCCCGGCTGACGCCCGGCCACCAGCAGCCCGCGCCAGGCGCGGAACCAGTCGACCACGATGGTGATCTCCAGGCCTGGCCGCTCGGCCTTGGCCTGGTACAAGGCa                                                                               >  1:38694/1‑150 (MQ=255)
                                                                            gcgcAACTGGCGCTGGGCCTGATACCAGGCGGCATTGCCCGGCTGACGCCCGGCCCCCAGCAGCCCGCGCCTGGCGCGGAACCAGTCGACCACGATGGTGATCTCCAGGCCTGGCCGCTCGGCCTTGGCCTGGTACAAGGCATCCa                                                                           <  1:75229/146‑1 (MQ=255)
                                                                            gcgcAACTGGCGCTGGGCCTGATACCAGGCGGCATTGCCCGGCTGACGCCCGGCCCCCAGCAGCCCGCGCCTGGCGCGGAACCAGTCGACCACGATGGTGATCTCCAGGCCTGGCCGCTCGGCCTTGGCCTGGTACAAGGCATCCa                                                                           >  2:75229/1‑146 (MQ=255)
                                                                            gcgcAACTGGCGCTGGGCCTGATACCAGGCGGCATTGCCCGGCTGACGCCCGGCCCCCAGCAGCCCGCGCCTGGCGCGGAACCAGTCGACCACGATGGTGATCTCCAGGCCTGGCCGCTCGGCCTTGGCCTGGTACAAGGCATCCAGGAc                                                                       >  1:633605/1‑150 (MQ=255)
                                                                                aaCTGGCGCTGGGCCTGATACCAGGCGGCATTGCCCGGCTGACGCCCGGCCCCCAGCAGCCCGCGCCTGGCGCGGAACCAGTCGACCACGATGGTGATCCAGAGGCCTGGCCGCTCGGCCTTGGCCTGGTACAAGGCATCCAGGACTTcc                                                                   >  2:398874/1‑150 (MQ=255)
                                                                                  cTGGCGCTGGGCCTGATACCAGGCGGCATTGCCCGGCTGACGCCCGGCCCCCAGCAGCCCGCGCCTGGCGCGGAACCAGTCGACCACGATGGTGATCTCCAGGCCTGGCCGCTCGGCCTTGGCCTGGTACAAGGCATCCAGGACTTCCTg                                                                 >  2:849176/1‑150 (MQ=255)
                                                                                   tGGCGCTGGGCCTGATACCAGGCGGCATTGCCCGGCTGACGCCCGGCCCCCAGCAGCCCGCGCCTGGCGCGGAACCAGTCGACCACGATGGTGATCTCCAGGCCTGGCCGCTCGGCCTTGGCCTGGTACAAGGCATCCAGGACTTCCTgg                                                                >  2:592519/1‑150 (MQ=255)
                                                                                   tGGCGCTGGGCCTGATACCAGGCGGCATTGCCCGGCTGACGCCCGGCCCCCAGCAGCCCGCGCCTGGCGCGGAACCAGTCGACCACGATGGTGATCTCCAGGCCTGGCCGCTCGGCCTTGGCCTGGTACAAGGCATCCAGGACTTCCTgg                                                                >  1:1364839/1‑150 (MQ=255)
                                                                                    ggCGCTGGGCCTGATACCAGGCGGCATTGCCCGGCTGACGCCCGGCCCCCAGCAGCCCGCGCCTGGCGCGGAACCAGTCGACCACGATGGTGATCTCCAGGCCTGGCCGCTCGGCCTTGGCCTGGTACAAGGCATCCAGGACTTCCTGGc                                                               <  2:963067/150‑1 (MQ=255)
                                                                                       gcTGGGCCTGATACCAGGCGGCATTGCCCGGCTGACGCCCGGCCCCCAGCAGCCCGCGCCTGGCGCGGAACCAGTCGACCACGATGGTGATCTCCAGGCCTGGc                                                                                                          >  2:1190047/1‑104 (MQ=255)
                                                                                       gcTGGGCCTGATACCAGGCGGCATTGCCCGGCTGACGCCCGGCCCCCAGCAGCCCGCGCCTGGCGCGGAACCAGTCGACCACGATGGTGATCTCCAGGCCTGGc                                                                                                          <  1:1190047/104‑1 (MQ=255)
                                                                                         tGGGCCTGATACCAGGCGGCATTGCCCGGCTGACGCCCGGCCCCCAGCAGCCCGCGCCTGGCGCGGAACCAGTCGACCACGATGGTGATCTCCAGGCCTGGCCGCTCGGCCTTGGCCTGGTACAAGGCATCCAGGACTTCCTGGCCAGcc                                                          <  2:1141666/150‑1 (MQ=255)
                                                                                           ggCCTGATACCAGGCGGCATTGCCCGGCTGACGCCCGGCCCCCAGCAGCCCGCGCCTGGCGCGGAACCAGTCGACCACGATg                                                                                                                            <  2:1235836/82‑1 (MQ=255)
                                                                                           ggCCTGATACCAGGCGGCATTGCCCGGCTGACGCCCGGCCCCCAGCAGCCCGCGCCTGGCGCGGAACCAGTCGACCACGATg                                                                                                                            >  1:1235836/1‑82 (MQ=255)
                                                                                           ggCCTGATACCAGGCGGCATTGCCCGGCTGACGCCCGGCCCCCAGCAGCCCGCGCCTGGCGCGGAACCAGTCGACCACGATGGTGATCTCCAg                                                                                                                 <  2:594156/93‑1 (MQ=255)
                                                                                           ggCCTGATACCAGGCGGCATTGCCCGGCTGACGCCCGGCCCCCAGCAGCCCGCGCCTGGCGCGGAACCAGTCGACCACGATGGTGATCTCCAg                                                                                                                 >  1:594156/1‑93 (MQ=255)
                                                                                           ggCCTGATACCAGGCGGCATTGCCCGGCTGACGCCCGGCCCCCAGCAGCCCGCGCCTGGCGCGGAACCAGTCGACCACGATGGTGATCTCCAGGCCTGGCCGCTCGGCCTTGGCCTGGTACaa                                                                                   >  1:1254290/1‑123 (MQ=255)
                                                                                           ggCCTGATACCAGGCGGCATTGCCCGGCTGACGCCCGGCCCCCAGCAGCCCGCGCCTGGCGCGGAACCAGTCGACCACGATGGTGATCTCCAGGCCTGGCCGCTCGGCCTTGGCCTGGTACAAGGCATCCAg                                                                          >  2:1324375/1‑132 (MQ=255)
                                                                                           ggCCTGATACCAGGCGGCATTGCCCGGCTGACGCCCGGCCCCCAGCAGCCCGCGCCTGGCGCGGAACCAGTCGACCACGATGGTGATCTCCAGGCCTGGCCGCTCGGCCTTGGCCTGGTACAAGGCATCCAg                                                                          <  1:1324375/132‑1 (MQ=255)
                                                                                              cttgtACCAGGCGGCATTGCCCGGCTGACGCCCGGCCCCCAGCAGCCCGCGCCTGGCGCGGAACCAGTCGACCACGATGGTGATCTCCAGGCCTGGCCGCTCGGCCTTGGCCTGGTACAAg                                                                                  <  1:499242/117‑1 (MQ=255)
                                                                                              cttgtACCAGGCGGCATTGCCCGGCTGACGCCCGGCCCCCAGCAGCCCGCGCCTGGCGCGGAACCAGTCGACCACGATGGTGATCTCCAGGCCTGGCCGCTCGGCCTTGGCCTGGTACAAg                                                                                  >  2:499242/5‑121 (MQ=255)
                                                                                                gATACCAGGCGGCATTGCCCGGCTGACGCCCGGCCCCCAGCAGCCCGCGCCTGGCGCGGAACCAGTCGACCACGATGGTGATCTCCAGGCCTGGCCGCTCGGCCTTGGCCTGGTACAAGGCATCCAGGa                                                                        <  2:349131/129‑1 (MQ=255)
                                                                                                 aTACCAGGCGGCATTGCCCGGCTGACGCCCGGCCCCCAGCAGCCCGCGCCTGGCGCGGAACCAGTCGACCACGATGGt                                                                                                                          <  1:1017059/78‑1 (MQ=255)
                                                                                                 aTACCAGGCGGCATTGCCCGGCTGACGCCCGGCCCCCAGCAGCCCGCGCCTGGCGCGGAACCAGTCGACCACGATGGt                                                                                                                          >  2:1017059/1‑78 (MQ=255)
                                                                                                   aCCAGGCGGCATTGCCCGGCTGACGCCCGGCCCCCAGCAGCCCGCGCCTGGCGCGGAACCAGTCGACCACGATGGTGATCTCCAGGCCTGGCCGCTCGGCCTTGGCCTGGTACAAGGCATCCAGGACTTCCTGGCCAGCCTCATcctcct                                                <  1:682714/150‑1 (MQ=255)
                                                                                                       ggcggcATTGCCCGGCTGACGCCCGGCCCCCAGCAGCCCGCGCCTGGCGCGGAACCAGTCGACCACGATGGTGATCTCCAGGCCTGGCCGCTCGGCCTTGGCCTg                                                                                         <  2:395874/105‑1 (MQ=255)
                                                                                                       ggcggcATTGCCCGGCTGACGCCCGGCCCCCAGCAGCCCGCGCCTGGCGCGGAACCAGTCGACCACGATGGTGATCTCCAGGCCTGGCCGCTCGGCCTTGGCCTg                                                                                         >  1:395874/1‑105 (MQ=255)
                                                                                                        gcggcATTGCCCGGCTGACGCCCGGCCCCCAGCAGCCCGCGCCTGGCGCGGAACCAGTCGACCACGATGGTGATCTCCAGGCCTGGCCGCTCGGCCTTGGCCTGGTACAAGGCATCCAGGACTTCCTGGCCAGCCTCATCCTCCTGCAgg                                           <  1:110173/150‑1 (MQ=255)
                                                                                                        gcggcATTGCCCGGCTGACGCCCGGCCCCCAGCAGCCCGCGCCTGGCGCGGAACCAGTCGACCACGATGGTGATCTCCAGGCCTGGCCGCTCGGCCTTGGCCTGGTACAAGGCATCCAGGACTTCCTGGCCAGCCTCATCCTCCTGCAgg                                           >  1:583724/1‑150 (MQ=255)
                                                                                                          ggcATTGCCCGGCTGACGCCCGGCCCCCAGCAGCCCGCGCCTGGCGCGGAACCAGTCGACCACGATgg                                                                                                                           >  1:791138/1‑68 (MQ=255)
                                                                                                              ttGCCCGGCTGACGCCCGGCCCCCAGCAGCCCGCGCCTGGCGCGGAACCAGTCGACCACGATGGTGATCTCCAGGCCTGGCCGCTCGGCCTTGGCCTGGTACAAGGCATCCAGGACTTCCTGGCCAGCCTCATCCTCCTGCAGGTACAAc                                     <  1:849176/150‑1 (MQ=255)
                                                                                                              ttGCCCGGCTGACGCCCGGCCCCCAGCAGCCCGCGCCTGGCGCGGAACCAGTCGACCACGATGGTGATCTCCAGGCCTGGCCGCTCGGCCTTGGCCTGGTACAAGGCATCCAGGACTTCCTGGCCAGCCTCATCCTCCTGCAGGTACAAc                                     <  2:1364839/150‑1 (MQ=255)
                                                                                                                  ccGGCTGACGCCCGGCCCCCAGCAGCCCGCGCCTGGCGCGGAACCAGTCGACCACGATGGTGATCTCCAGGCCTGGCCGCTCGGCCTTGGCCTGGTACAAGGCATCCAGGACTTCCTGGCCAGCCTCATCCTCCTTCa                                             >  1:294362/1‑138 (MQ=255)
                                                                                                                  ccGGCTGACGCCCGGCCCCCAGCAGCCCGCGCCTGGCGCGGAACCAGTCGACCACGATGGTGATCTCCAGGCCTGGCCGCTCGGCCTTGGCCTGGTACAAGGCATCCAGGACTTCCTGGCCAGCCTCATCCTCCTGc                                              >  1:185818/1‑137 (MQ=255)
                                                                                                                  ccGGCTGACGCCCGGCCCCCAGCAGCCCGCGCCTGGCGCGGAACCAGTCGACCACGATGGTGATCTCCAGGCCTGGCCGCTCGGCCTTGGCCTGGTACAAGGCATCCAGGACTTCCTGGCCAGCCTCATCCTCCTGc                                              <  2:185818/137‑1 (MQ=255)
                                                                                                                  ccGGCTGACGCCCGGCCCCCAGCAGCCCGCGCCTGGCGCGGAACCAGTCGACCACGATGGTGATCTCCAGGCCTGGCCGCTCGGCCTTGGCCTGGTACAAGGCATCCAGGACTTCCTGGCCAGCCTCATCCTCCTGCa                                             <  2:294362/138‑1 (MQ=255)
                                                                                                                         aCGCCCGGCCCCCAGCAGCCCGCGCCTGGCGCGGAACCAGTCGACCACGATGGTGATCTCCAGGCCTGGCCGCTCGGCCTTGGCCTGGTACAAGGCATCCAGGACTTCCTGGCCAGCCTCATCCTCCTGCAGGTACAACGCGACGATAAc                          <  1:592519/150‑1 (MQ=255)
                                                                                                                           gCCCGGCCCCCAGCAGCCCGCGCCTGGCGGGGAACCAGTCGACCACGATGGTGATCTCCAGGCCTGGCCGCTCGGCCTTGGCCTGGTACAAGGCATCCAGGACTTCCTGGCCAGCCTCATCCTCCTGCAGGTACAACGCGACGATAACGa                        >  1:124606/1‑150 (MQ=255)
                                                                                                                           gCCCGGCCCCCAGCAGCCCGCGCCTGGCGCGGAACCAGTCGACCACGATGGTGATCTCCAGGCCTGGCCGCTCGGCCTTGGCCTGGTACAAGGCATCCAGGCCTTCCTGGCCAGCCTCATCCTCCTGCAGGTACAACGCGACGATaa                           >  2:953204/1‑147 (MQ=255)
                                                                                                                           gCCCGGCCCCCAGCAGCCCGCGCCTGGCGCGGAACCAGTCGACCACGATGGTGATCTCCAGGCCTGGCCGCTCGGCCTTGGCCTGGTACAAGGCATCCAGGACGTCCTGGCCAGCCTCATCCTCCTGCAGGTACAACGCGACGATaa                           <  1:953204/147‑1 (MQ=255)
                                                                                                                              cGGCCCCCAGCAGCCCGCGCCTGGCGCGGAACCAGTCGACCACGATGGTGATCTCCAGGCCTGGCCGCTCGGCCTTGGCCTGGTACAAGGCATCCAGGACTTCCTGGCCAGCCTCATCCTCCTGCAGGTACAACGCGACGATAACGATgc                     >  1:393318/1‑150 (MQ=255)
                                                                                                                                 cccccAGCAGCCCGCGCCTGGCGCGGAACCAGTCGACCACGATGGTGATCTCCAGGCCTGGCCGCTCGGCCTTGGCCTGGTACAAGGCATCCa                                                                           <  1:548583/93‑1 (MQ=255)
                                                                                                                                 cccccAGCAGCCCGCGCCTGGCGCGGAACCAGTCGACCACGATGGTGATCTCCAGGCCTGGCCGCTCGGCCTTGGCCTGGTACAAGGCATCCa                                                                           >  2:548583/1‑93 (MQ=255)
                                                                                                                                   cccAGCAGCCCGCGCCTGGCGCGGAACCAGTCGACCACGATGGTGATCTCCAGGCCTGGCCGCTCGGCCTTGGCCTGGTACaa                                                                                   >  1:922289/1‑83 (MQ=255)
                                                                                                                                   cccAGCAGCCCGCGCCTGGCGCGGAACCAGTCGACCACGATGGTGATCTCCAGGCCTGGCCGCTCGGCCTTGGCCTGGTACaa                                                                                   <  2:922289/83‑1 (MQ=255)
                                                                                                                                   cccAGCAGCCCGCGCCTGGCGCGGAACCAGTCGACCACGATGGTGATCTCCAGGCCTGGCCGCTCGGCCTTGGCCTGGTACAAGGCATCCAGGACTTCCTGGCCAGCCTCATCCTCCTGCAGGTACAACGCGACGATAACGATGCGCCGa                <  2:906956/150‑1 (MQ=255)
                                                                                                                                    ccAGCAGCCCGCGCCTGGCGCGGAACCAGTCGACCACGATGGTGATCTCCAGGCCTGGCCGCTCGGCCTTGGCCTGGTACAAGGCATCCAGGACTTCCTGGCCAGCCTCATCCTCCTGCAgg                                           <  1:129341/122‑1 (MQ=255)
                                                                                                                                    ccAGCAGCCCGCGCCTGGCGCGGAACCAGTCGACCACGATGGTGATCTCCAGGCCTGGCCGCTCGGCCTTGGCCTGGTACAAGGCATCCAGGACTTCCTGGCCAGCCTCATCCTCCTGCAgg                                           >  2:129341/1‑122 (MQ=255)
                                                                                                                                          ttccGCGCCTGGCGCGGAACCAGTCGACCACGATGGTGATCTCCAGGCCTGGCCGCTCGGCCTTGGCCTGGTACAAGGCATCCAGGACTTCCTGGCCAGCCTCATCCTCCTGCAGGTACAACGCGACGATAACGATGCGCCGAGTcgccg         >  2:899553/3‑150 (MQ=255)
                                                                                                                                            ccGCGCCTGGCGCGGAACCAGTCGACCACGATGGTGATCTCCAGGCCTGGCCGCTCGGCCTTGGCCTGGTACAAGGCATCCAGGACTTCCTGGCCAGCCTCATCCTCCTGCAGGTACAACGCGACGATAACGATGCGCCGAGTCGCCGCa       >  1:146970/1‑150 (MQ=255)
                                                                                                                                             cgcgCCTGGCGCGGAACCAGTCGACCACGATGGTGATCTCCAGGCCTGGCCGCTCGGCCTTGGCCTGGTACAAGGCATCCAGGACTTCCTGGCCAGCCTCATCCTCCTGCAGGTACAACGCGACGATAACGATGCGCCGAGTcgccgccc      >  1:1208427/1‑148 (MQ=255)
                                                                                                                                                gcCTGGCGCGGAACCAGTCGACCACGATGGTGATCTCCAGGCCTGGCCGCTCGGCCTTGGCCTGGTACAAGGCATCCAGGACTTCCTGGCCAGCCTCATCCTCCTGCAGGTACAACGCGACGATAACGATGCGCCGAGTCGCCGCACGGa   >  2:246587/1‑150 (MQ=255)
                                                                                                                                                 ccTGGCGCGGAACCAGTCGACCACGATGGTGATCTCCAGGCCTGGCCGCTCGGCCTTGGCCTGGTACAAGGCATCCAGGACTTCCTGGCCAGCCTCATCCTCCTGCAGGTACAACGCGACGATAACGATGCGCCGAGTCGCCGCACGGAt  <  2:393318/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
GCTGCCTTTGAGGTGCAGCACCCCGAACAGTTCGCGGGTCTGCACTGGCACGCCATGAATCGCGATGTCCAGCCCGCGCAACTGGCGCTGGGCCTGATACCAGGCGGCATTGCCCGGCTGACGCCCGGCCCCCAGCAGCCCGCGCCTGGCGCGGAACCAGTCGACCACGATGGTGATCTCCAGGCCTGGCCGCTCGGCCTTGGCCTGGTACAAGGCATCCAGGACTTCCTGGCCAGCCTCATCCTCCTGCAGGTACAACGCGACGATAACGATGCGCCGAGTCGCCGCACGGAT  >  NC_002947/4163104‑4163397

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: