Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,727,826 A→C 6.4% V113G (GTC→GGC)  plsB ← glycerol‑3‑phosphate acyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,727,8260AC 6.4% 186.0 / 4.1 63V113G (GTC→GGC) plsBglycerol‑3‑phosphate acyltransferase
Reads supporting (aligned to +/- strand):  ref base A (32/27);  new base C (2/2);  total (34/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.34e-01

GTACGGTCAGCAGCCGGCGCAGGCGCCCGGTTACCGCCCAGCTGTCGGCGAACAGCAGCTTCCAGGGGCTGGACTCGCTGGCCGGG‑GTCTGGCCCCAGAACACGCTGACAGGAATGATCTGCGCATCTTCTTCGGCATGCTGGCTGACGGCAGCCACCAGGCGCTCGAGCGTGGGCGGGGCGCCGCTTTTGTCCTGGCGGCCGAGCCAGTCCGGGTCGGGGGTCAGGTAGAAGAACGCGGCGGGTTCCTGCAGTGGGCCGACGGCTACCGGCAGCACCGGGCGCGGCAGCCC  >  NC_002947/1727680‑1727971
                                                                                                                                                   |                                                                                                                                                 
gTACGGTCAGCAGCCGGCGCAGGCGCCCGGTTACCGCCCAGCTGTCGGCGAACAGCAGCTTCCAGGGGCTGGACTCGCTGGCCGGG‑GTCTGGCCCCAGAACACGCTGACAGGAATGATCTGCGCATCTTCTTCGGCATGCTGGCTGACg                                                                                                                                                 >  2:258161/1‑149 (MQ=255)
 tACGGTCAGCAGCCGGCGCAGGCGCCCGGTTACCGCCCAGCTGTCGGCGAACAGCAGCTTCCAGGGGCTGGACTCGCTGGCCGGG‑GTCTGGCCCCAGAACACGCTGACAGGAATGATCTGCGCATCTTCTTCGGCATGCTGGCTGACgg                                                                                                                                                >  1:118032/1‑149 (MQ=255)
             ccGGCGCAGGCGCCCGGTTACCGCCCAGCTGTCGGCGAACAGCAGCTTCCAGGGGCTGGACTCGCTGGCCGGG‑GTCTGGCCCCAGAACACGCTGACAGGAATGATCTGCGCATCTTCTTCGGCATGCTGGCTGACg                                                                                                                                                 >  3:183701/1‑136 (MQ=255)
             ccGGCGCAGGCGCCCGGTTACCGCCCAGCTGTCGGCGAACAGCAGCTTCCAGGGGCTGGACTCGCTGGCCGGG‑GTCTGGCCCCAGAACACGCTGACAGGAATGATCTGCGCATCTTCTTCGGCATGCTGGCTGACg                                                                                                                                                 <  4:183701/136‑1 (MQ=255)
             ccGGCGCAGGCGCCCGGTTACCGCCCAGCTGTCGGCGAACAGCAGCTTCCAGGGGCTGGACTCGCTGGCCGGG‑GTCTGGCCCCAGAACACGCTGACAGGAATGATCTGCGCATCTTCTTCGGCATGCTGGCTGACg                                                                                                                                                 <  8:100455/136‑1 (MQ=255)
             ccGGCGCAGGCGCCCGGTTACCGCCCAGCTGTCGGCGAACAGCAGCTTCCAGGGGCTGGACTCGCTGGCCGGG‑GTCTGGCCCCAGAACACGCTGACAGGAATGATCTGCGCATCTTCTTCGGCATGCTGGCTGACg                                                                                                                                                 >  7:100455/1‑136 (MQ=255)
             ccGGCGCAGGCGCCCGGTTACCGCCCAGCTGTCGGCGAACAGCAGCTTCCAGGGGCTGGACTCGCTGGCCGGG‑GTCTGGCCCCAGAACACGCTGACAGGAATGATCTGCGCATCTTCTTCGGCATGCTGGCTGACGGCAGccac                                                                                                                                         >  8:22829/1‑144 (MQ=255)
             ccGGCGCAGGCGCCCGGTTACCGCCCAGCTGTCGGCGAACAGCAGCTTCCAGGGGCTGGACTCGCTGGCCGGG‑GTCTGGCCCCAGAACACGCTGACAGGAATGATCTGCGCATCTTCTTCGGCATGCTGGCTGACGGCAGccac                                                                                                                                         <  7:22829/144‑1 (MQ=255)
               ggCGCAGGCGCCCGGTTACCGCCCAGCTGTCGGCGAACAGCAGCTTCCAGGGGCTGGACTCGCTGGCCGGG‑GTCTGGCCCCAGAACACGCTGACAGGAATGATCTGCGCATCTTCTTCGGCATGCTGGCTGACGGCAGCCACCAGgc                                                                                                                                    >  7:205418/1‑147 (MQ=255)
               ggCGCAGGCGCCCGGTTACCGCCCAGCTGTCGGCGAACAGCAGCTTCCAGGGGCTGGACTCGCTGGCCGGG‑GTCTGGCCCCAGAACACGCTGACAGGAATGATCTGCGCATCTTCTTCGGCATGCTGGCTGACGGCAGCCACCAGgc                                                                                                                                    <  8:205418/147‑1 (MQ=255)
                gcgcAGGCGCCCGGTTACCGCCCAGCTGTCGGCGAACAGCAGCTTCCAGGGGCTGGACTCGCTGGCCGGG‑GTCTGGCCCCAGAACACGCTGACAGGAATGATCTGCGCATCTTCTTCGGCATGCTGGCTGACGGCAGCCACCAGGCGCt                                                                                                                                 <  7:134181/149‑1 (MQ=255)
                gcgcAGGCGCCCGGTTACCGCCCAGCTGTCGGCGAACAGCAGCTTCCAGGGGCTGGACTCGCTCGCCGGG‑GTCTGGCCCCAGAACACGCTGACAGGAATGATCTGCGCATCTTCTTCGGCATGCTGGCTGACGGCAGCCACCAGGCGCt                                                                                                                                 >  4:71286/1‑149 (MQ=255)
                      gcgcCCGGTTACCGCCCAGCTGTCGGCGAACAGCAGCTTCCAGGGGCTGGACTCGCTGGCCGGG‑GTCTGGCCCCAGAACACGCTGACAGGAATGATCTGCGCATCTTCTTCGGCATGCTGGCTGCCGTCAGccagcatg                                                                                                                                     <  6:18369/139‑3 (MQ=255)
                      gcgcCCGGTTACCGCCCAGCTGTCGGCGAACAGCAGCTTCCAGGGGCTGGACTCGCTGGCCGGG‑GTCTGGCCCCAGAACACGCTGACAGGAATGATCTGCGCATCTTCTTCGGCATGCTGGCTGCCGTCAGccagcatg                                                                                                                                     >  5:18369/1‑137 (MQ=255)
                      gcgcCCGGTTACCGCCCAGCTGTCGGCGAACAGCAGCTTCCAGGGGCTGGACTCGCTGGCCGGG‑GTCTGGCCCCAGAACACGCTGACAGGAATGATCTGCGCATCTTCTTCGGCATGCTGGCTGACGGCAGCCACCAGGCGCTCGAGCg                                                                                                                           >  8:46422/1‑149 (MQ=255)
                        gcCCGGTTACCGCCCAGCTGTCGGCGAACAGCAGCTTCCAGGGGCTGGACTCGCTGGCCGGG‑GTCTGGCCCCAGAACACGCTGACAGGAATGATCTGCGCATCTTCTTCGGCATGCTGGCTGACGGCAGCCACCAGGCGCTCGAGCGTg                                                                                                                         <  5:219113/149‑1 (MQ=255)
                         cccGGTTACCGCCCAGCTGTCGGCGAACAGCAGCTTCCAGGGGCTGGACTCGCTGGCCGGG‑GTCTGGCCCCAGAACACGCTGACAGGAATGATCTGCGCATCTTCTTCGGCATGCTGGCTGACGGCAGCCACCAGGCGCTCGAGCGTgg                                                                                                                        <  7:46422/149‑1 (MQ=255)
                                 ccGCCCAGCTGTCGGCGAACAGCAGCTTCCAGGGGCTGGACTCGCTGGCCGGG‑GTCTGGCCCCAGAACACGCTGACAGGAATGATCTGCGCATCTTCTTCGGCATGCTGGCTGCCGTCAGccagcat                                                                                                                                      <  3:255440/127‑2 (MQ=255)
                                 ccGCCCAGCTGTCGGCGAACAGCAGCTTCCAGGGGCTGGACTCGCTGGCCGGG‑GTCTGGCCCCAGAACACGCTGACAGGAATGATCTGCGCATCTTCTTCGGCATGCTGGCTGCCGTCAGccagcat                                                                                                                                      >  4:255440/1‑126 (MQ=255)
                                     ccAGCTGTCGGCGAACAGCAGCTTCCAGGGGCTGGACTCGCTGGCCGGG‑GTCTGGCCCCAGAACACGCTGACAGGAATGATCTGCGCATCTTCTTCGGCATGCTGGCTGACGGCAGCCACCAGGCGCTCGAGCGTGGGCGGGGcgccgc                                                                                                            >  7:245452/1‑149 (MQ=255)
                                     ccAGCTGTCGGCGAACAGCAGCTTCCAGGGGCTGGACTCGCTGGCCGGG‑GTCTGGCCCCAGAACACGCTGACAGGAATGATCTGCGCATCTTCTTCGGCATGCTGGCTGACGGCAGCCACCAGGCGCTCGAGCGTGGGCGGGGcgccgc                                                                                                            >  2:146067/1‑149 (MQ=255)
                                           gTCGGCGAACAGCAGCTTCCAGGGGCTGGACTCGCTGGCCGGGAGTCTGGCCCCAGAACACGCTGACAGGAATGATCTGCGCATCTTCTTCGGCATGCTGGCTGACGGCAGCCACCAGGCGCTCGAGCGTGGGCGGGGCGCCGCgttct                                                                                                       >  1:52953/1‑147 (MQ=255)
                                                   aCAGCAGCTTCCAGGGGCTGGACTCGCTGGCCGGG‑GTCTGGCCCCAGAACACGCTGACAGGAATGATCTGTGCATCTTCTTCGGCATGCTGGCTGACGGCAGCCACCAGGCGCTCGAGCGTGGGCGGGGCGCCGCTTTTGTCCTggcgg                                                                                              >  7:85401/1‑149 (MQ=255)
                                                    cagcagCTTCCAGGGGCTGGACTCGCTGGCCGGG‑GTCTGGCCCCAGAACACGCTGACAGGAATGATCTGCGCATCTTCTTCGGCATGCTGGCTGACGGCAGCCACCAGGCGCTCGAGCGTGGGCGGGGCGCCGCTTTTGTCCTggcggc                                                                                             <  7:68032/149‑1 (MQ=255)
                                                     agcagcTTCCAGGGGCTGGACTCGCTGGCCGGG‑GTCTGGCCCCAGAACACGCTGACAGGAATGATCTGCGCATCTTCTTCGGCATGCTGGCTGACGGCAGCCACCAGGCGCTCGAGCGTGGGCGGGGCGCCGCTTTTGTCCTGGCGGcc                                                                                            >  1:129080/1‑149 (MQ=255)
                                                       cagcTTCCAGGGGCTGGACTCGCTGGCCGGG‑GTCTGGCCCCAGAACACGCTGACAGGAATGATCTGCGCATCTTCTTCGGCATGCTGGCTGACGGCAGCCACCAGGCGCTCGAGCGTGGGCGGCGCGCCGCTTTTGTCCTGGCGGcc                                                                                            >  4:239021/1‑147 (MQ=255)
                                                         gcTTCCAGGGGCTGGACTCGCTGGCCGGG‑GTCTGGCCCCAGAACACGCTGACAGGAATGATCTGCGCATCTTCTTCGGCATGCTGGCTGACGGCAGCCACCAGGCGCTCGAGCGTGGGCGGGGCGCCGCTTTTGTCCTGGCGGCCGAGc                                                                                        <  8:138561/149‑1 (MQ=255)
                                                           ttCCAGGGGCTGGACTCGCTGGCCGGG‑GTCTGGCCCCAGAACACGCTGACAGGAATGATCTGCGCATCTTCTTCGGCATGCTGGCTGACGGCAGCCACCAGGCGCTCGAGCGTGGGCGGGGCGCCGCTTTTGTCCt                                                                                                   <  6:21709/136‑1 (MQ=255)
                                                           ttCCAGGGGCTGGACTCGCTGGCCGGG‑GTCTGGCCCCAGAACACGCTGACAGGAATGATCTGCGCATCTTCTTCGGCATGCTGGCTGACGGCAGCCACCAGGCGCTCGAGCGTGGGCGGGGCGCCGCTTTTGTCCt                                                                                                   >  5:21709/1‑136 (MQ=255)
                                                           ttCCAGGGGCTGGACTCGCTGGCCGGG‑GTCTGGCCCCAGAACACGCTGACAGGAATGATCTGCGCATCTTCTTCGGCATGCTGGCTGACGGCAGCCACCAGGCGCTCGAGCGTGGGCGGGGCGCCGCTTTTGTCCTGGCGGCCGAGc                                                                                        >  2:141487/1‑147 (MQ=255)
                                                               aGGGGCTGGACTCGCTGGCCGGG‑GTCTGGCCCCAGAACACGCTGACAGGAATGATCTGCGCATCTTCTTCGGCATGCTGGCTGACGGCAGCCACCAGGCGCTCGAGCGTGGGCGGGGCGCCGCTTTTGTCCTGGCGGCCGAGCCAGTcc                                                                                  <  4:127236/149‑1 (MQ=255)
                                                                ggggCTGGACTCGCTGGCCGGG‑GTCTGGCCCCAGAACACGCTGACAGGAATGATCTGCGCATCTTCTTCGGCATGCTGGCTGACGGCAGCCACCAGGCGCTCGAGCGTGGGCGGGGCGCCGCTTTTGTCCTGGCGGCCGAGCCAGTCCg                                                                                 >  2:221957/1‑149 (MQ=255)
                                                                 ggggTGGACTCGCTGGCCGGG‑GTCTGGCCCCAGAACACGCTGACAGGAATGATCTGCGCATCTTCTTCGGCATGCTGGCTGACGGCAGCCACCAGGCGCTCGAGCGTGGGCGGGGCGCCGCTTTTGTCCTGGCGGCCGAGCCAGTCCgg                                                                                <  4:75900/149‑1 (MQ=255)
                                                                      ggACTCGCTGGCCGGG‑GTCTGGCCCCAGAACACGCTGACAGGAATGATCTGCGCATCTTCTTCGGCATGCTGGCTGACGGCAGCCACCAGGCGCTCGAGCGTGGGCGGGGCGCCGCTTTTGTCCTGGCGGCCGAGCCAGTCCGGGTCgg                                                                           <  1:146067/149‑1 (MQ=255)
                                                                       gACTCGCTGGCCGGG‑GTCTGGCCCCAGAACACGCTGACAGGAATGATCTGCGCATCTTCTTCGGCATGCTGGCTGACGGCAGCCACCAGGCGCTCGAGCGTGGGCGGGGCGCCGCTTTTGTCCTGGCGGCCGAGCCAGTCCGGGTCggg                                                                          >  3:87218/1‑149 (MQ=255)
                                                                          tCGCTGGCCGGG‑GTCTGGCCCCAGAACACGCTGACAGGAATGATCTGCGCATCTTCTTCGGCATGCTGGCTGACGGCAGCCACCAGGCGCTCTAGCGTGGGCGGGGCGCCGCTTTTGTCCTGGCGGCCGAGCCAGTCCGGGTCGGGGGt                                                                       <  3:239021/149‑1 (MQ=255)
                                                                          tCGCTGGCCGGG‑GTCTGGCCCCAGAACACGCTGACAGGAATGATCTGCGCATCTTCTTCGGCATGCTGGCTGACGGCAGCCACCAGGCGCTCGAGCGTGGGCGGGGCGCCGCttt                                                                                                         >  5:168391/1‑115 (MQ=255)
                                                                          tCGCTGGCCGGG‑GTCTGGCCCCAGAACACGCTGACAGGAATGATCTGCGCATCTTCTTCGGCATGCTGGCTGACGGCAGCCACCAGGCGCTCGAGCGTGGGCGGGGCGCCGCttt                                                                                                         <  6:168391/115‑1 (MQ=255)
                                                                            gCTGGCCGGG‑GTCTGGCCCCAGAACACGCTGACAGGAATGATCTGCGCATCTTCTTCGGCATGCTGGCTGACGGCAGCCACCAGGCGCTCGAGCGTGGGCGGGGCGCCGCTTTTGTCCTGGCGGCCGAGCCAGTCCGGGTCGGGGGTCa                                                                     <  2:129080/149‑1 (MQ=255)
                                                                                     g‑gTCTGGCCCCAGAACACGCTGACAGGAATGATCTGCGCATCTTCTTCGGCATGCTGGCTGACGGCAGCCACCAGGCGCTCGAGCGTGGGCGGGGCGCCGCTTTTGTCCTGGCGGCCGAGCCAGTCCGGGTCGGGGGTCAGGTagaaga                                                            <  2:24868/149‑1 (MQ=255)
                                                                                            gCCCCAGAACACGCTGACAGGAATGATCTGCGCATCTTCTTCGGCATGCTGGCTGACGGCAGCCACCAGGCGCTCGAGCGTGGGCGGGGCGCCGCTTTTGTCCTGGCGGCCGAGCCAGTCCGGGTCGGGGGTCAGGTag                                                                <  6:52954/139‑1 (MQ=255)
                                                                                            gCCCCAGAACACGCTGACAGGAATGATCTGCGCATCTTCTTCGGCATGCTGGCTGACGGCAGCCACCAGGCGCTCGAGCGTGGGCGGGGCGCCGCTTTTGTCCTGGCGGCCGAGCCAGTCCGGGTCGGGGGTCAGGTag                                                                >  5:52954/1‑139 (MQ=255)
                                                                                              cccAGAACACGCTGACAGGAATGATCTGCGCATCTTCTTCGGCATGCTGGCTGACGGCAGCCACCAGGCGCTCGAGCGTGGGCGGGGCGCCGCTTTTGTCCTGGCGGCCGAGCCAGTCCGGGTCGGGGGTCAGGTAGAAGAACgcggcg                                                    <  2:86699/149‑1 (MQ=255)
                                                                                                     cacGCTGACAGGAATGATCTGCGCATCTTCTTCGGCATGCTGGCTGACGGCAGCCACCAGGCGCTCGAGCGTGGGCGGGGCGCCGCTTTTGTACTGGCGGCCGAGCCAGTCCGGGTCGGGGGTCAGGTAGAAGAACGCGGCGGGTTCCt                                             >  1:231998/1‑149 (MQ=255)
                                                                                                      acGCTGACAGGAATGATCTGCGCATCTTCTTCGGCATGCTGGCTGACGGCAGCCACCAGGCGCTCGAGCGTGGGCGGGGCGCCGCTTTTGTCCTGGCGGCCGAGCCAGTCCGGGTCGGGGGTCAGGTagaa                                                              >  2:82147/1‑131 (MQ=255)
                                                                                                      acGCTGACAGGAATGATCTGCGCATCTTCTTCGGCATGCTGGCTGACGGCAGCCACCAGGCGCTCGAGCGTGGGCGGGGCGCCGCTTTTGTCCTGGCGGCCGAGCCAGTCCGGGTCGGGGGTCAGGTagaa                                                              <  1:82147/131‑1 (MQ=255)
                                                                                                      acGCTGACAGGAATGATCTGCGCATCTTCTTCGGCATGCTGGCTGACGGCAGCCACCAGGCGCTCGAGCGTGGGCGGGGCGCCGCTTTTGTCCTGGCGGCCGAGCCAGTCCGGGTCAGGGGTCAGGTAGAAGAACGCGAAGCGTTCCt                                             >  6:188990/1‑148 (MQ=255)
                                                                                                        gCTGACAGGAATGATCTGCGCATCTTCTTCGGCATGCTGGCTGACGGCAGCCACCAGGCGCTCGAGCGTGGGCGGGGCGCCGCTTTTGTCCTggcggc                                                                                             <  1:87765/98‑1 (MQ=255)
                                                                                                        gCTGACAGGAATGATCTGCGCATCTTCTTCGGCATGCTGGCTGACGGCAGCCACCAGGCGCTCGAGCGTGGGCGGGGCGCCGCTTTTGTCCTggcggc                                                                                             >  2:87765/1‑98 (MQ=255)
                                                                                                        gCTGACAGGAATGATCTGCGCATCTTCTTCGGCATGCTGGCTGACGGCAGCCACCAGGCGCTCGAGCGTGGGCGGGGCGCCGCTTTTGTCCTGGCGGCCGAGCCAGTCCGGGTCGGGGGTCAGGTAGAAGAACGCGGCGGGTTCCTg                                            >  2:258263/1‑147 (MQ=255)
                                                                                                              aGGAATGATCTGCGCATCTTCTTCGGCATGCTGGCTGACGGCAGCCACCAGGCGCTCGAGCGTGGGCGGGGCGCCGCTTTTGTCCTGGCGGCCGAGCCAGTCCGGGTCGGGGGTCAGGTAGAAGAACGCGGCGGGTTCCTGCAGTGGGc                                    <  8:245452/149‑1 (MQ=255)
                                                                                                              aGGAATGATCTGCGCATCTTCTTCGGCATGCTGGCTGACGGCAGCCACCAGGCGCTCGAGCGTGGGCGGGGCGCCGCTTTTGTCCTGGCGGCCGAGCCAGTCCGGGTCGGGGGTCAGGTAGAAGAACGCGGCGGGTTCCTGCAGTGGGc                                    <  3:71286/149‑1 (MQ=255)
                                                                                                                        tGCGCATCTTCTTCGGCATGCTGGCTGACGGCAGCCACCAGGCGCTCGAGCGTGGGCGGGGCGCCGCTTTTGTCCTGGCGGCCGAGCCAGTCCGGGTCGGGGGTCAGGTAGAAGAACGCGGAGGGTTCCTGCAGTGGGCCGACGGCTAc                          >  3:246983/1‑149 (MQ=255)
                                                                                                                             aTCTTCTTCGGCATGCTGGCTGACGGCAGCCACCAGGCGCTCGAGCGTGGGCGGGGCGCCGCTTTTGTCCTGGCGGCCGAGCCAGTCCGGGTCGGGGGTCAGGTAGAAGAACGCGGCGGGTTCCTGCAGTGGGCCGACGGCTACCGgca                     <  2:231998/149‑1 (MQ=255)
                                                                                                                             aTCTTCTTCGGCATGCTGGCTGACGGCAGCCACCAGGCGCTCGAGCGTGGGCGGGGCGCCGCTTTTGTCCTGGCGGCCGAGCCAGTCCGGGTCGGGGGTCAGGTAGAAGAACGCGGCGGGTTCCTGCAGTGGGCAGACGGCTACCgg                       >  3:89370/1‑147 (MQ=255)
                                                                                                                                     cGGCATGCTGGCTGACGGCAGCCACCAGGCGCTCGAGCGTGGGCGGGGCGCCGCTTTTGTCCTGGCGGCCGAGCCAGTCCGGGTCGGGGGTCAGGTAGAAGAACGCGGCGGGTTCCTGCAGTGGGCCGACGGCTACCGGCAGCACCggg             >  8:42138/1‑149 (MQ=255)
                                                                                                                                          tGCTGGCTGACGGCAGCCACCAGGCGCTCGAGCGTGGGCGGGGCGCCGCTTTTGTCCTGGCGGCCGAGCCAGTCCGGGTCGGGGGTCAGGTAGAAGAACGCGGCGGGTTCCTGCAGTGGGCCGACGGCTACCGGCAGCACCGGGCGCgg        <  8:85401/149‑1 (MQ=255)
                                                                                                                                              ggctgACGGCAGCCACCAGGCGCTCGAGCGTGGGCGGGGCGCCGCTTTTGTCCTGGCGGCCGAGCCAGTCCGGGTCGGGGGTCAGGTAGAAGAACgcggc                                                     <  8:29468/100‑1 (MQ=255)
                                                                                                                                              ggctgACGGCAGCCACCAGGCGCTCGAGCGTGGGCGGGGCGCCGCTTTTGTCCTGGCGGCCGAGCCAGTCCGGGTCGGGGGTCAGGTAGAAGAACgcggc                                                     >  7:29468/1‑100 (MQ=255)
                                                                                                                                              ggctgACGGCAGCCACCAGGCGCTCGAGCGTGGGCGGGGCGCCGCTTTTGTCCTGGCGGCCGAGCCAGTCCGGGTCGGGGGTCAGGTAGAAGAACGCGGCGGGTTCCTGCAGTGGGCCGACGGCTACCGGCAGCACCGGGCGCGGCAGc    >  5:166779/1‑149 (MQ=255)
                                                                                                                                                ctgACGGCAGCCACCAGGCGCTCGAGCGTGGGCGGGGCGCCGCTTTTGTCCTGGCGGCCGAGCCAGTCCGGGTCGGGGGTCAGGTAGAAGAACGCGGCGGGTTCCTGCAGTGGGCCGACGGCTACCGGCAGCACCGGGCGCGGCAGccc  >  6:271822/1‑149 (MQ=255)
                                                                                                                                                 tgACGGCAGCCACCAGGCGCTCGAGCGTGGGCGGGGCGCCGCTTTTGTCCTGGCGGCCGAGCCAGTCCGGGTCGGGGGTCAGGTAGAAGAACGCGGCGGGTTCCTg                                            >  7:46554/1‑106 (MQ=255)
                                                                                                                                                 tgACGGCAGCCACCAGGCGCTCGAGCGTGGGCGGGGCGCCGCTTTTGTCCTGGCGGCCGAGCCAGTCCGGGTCGGGGGTCAGGTAGAAGAACGCGGCGGGTTCCTg                                            <  8:46554/106‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
GTACGGTCAGCAGCCGGCGCAGGCGCCCGGTTACCGCCCAGCTGTCGGCGAACAGCAGCTTCCAGGGGCTGGACTCGCTGGCCGGG‑GTCTGGCCCCAGAACACGCTGACAGGAATGATCTGCGCATCTTCTTCGGCATGCTGGCTGACGGCAGCCACCAGGCGCTCGAGCGTGGGCGGGGCGCCGCTTTTGTCCTGGCGGCCGAGCCAGTCCGGGTCGGGGGTCAGGTAGAAGAACGCGGCGGGTTCCTGCAGTGGGCCGACGGCTACCGGCAGCACCGGGCGCGGCAGCCC  >  NC_002947/1727680‑1727971

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