Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,726,485 A→T 7.4% L560Q (CTG→CAG)  plsB ← glycerol‑3‑phosphate acyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,726,4850AT 7.4% 186.6 / 3.5 64L560Q (CTG→CAG) plsBglycerol‑3‑phosphate acyltransferase
Reads supporting (aligned to +/- strand):  ref base A (30/29);  new base T (3/2);  total (33/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.72e-01

AGCACCGCGTTGCCTTCATCCAGGTAGAGGATGCGGCCCAGGGCGTCCTTCTGCTCGGCCACCAGGTTCATGCTGCGCACATGTTCGATCAGTGCCTGGCCGTCGCCCTCGGGCAGGGTGGTGTGCGGCGAGTAGGGCACTTGGCGCAGCAGCGCCAGGTACAGGTCGAGTACGCGGGTCAGGGCGCGCTCGTCCAGGGCAAGGCGGCTGGTGGACAGCAGTGCCAGGGCCACCAGGTTGACCGGGTTGATGGCGGCCGCCTCGTTGAGGTGGCGGGCCACGGTTTCTCCCA  >  NC_002947/1726338‑1726629
                                                                                                                                                   |                                                                                                                                                
aGCACCGCGTTGCCTTCATCCAGGTAGAGGATGCGGCCCAGGGCGTCCTTCTGCTCGGCCACCAGGTTCATGCTGCGCACATGTTCGATCAGTGCCTGGCCGTCGCCCTCGGGCAGGGTGGTGTGCGGCGAGTAGGGCACTTGGCgcag                                                                                                                                                 >  7:75132/1‑149 (MQ=255)
 gCACCGCGTTGCCTTCATCCAGGTAGAGGATGCGGCCCAGGGCGTCCTTCTGCTCGGCCACCAGGTTCATGCTGCGCACATGTTCGATCAGTGCCTGGCCGTCGCCCTCGGGCAGGGTGGTGTGCGGCGAGTAGGGCACTTGGCgcagc                                                                                                                                                <  2:263699/149‑1 (MQ=255)
    ccGCGTTGCCTTCATCCAGGTAGAGGATGCGGCCCAGGGCGTCCTTCTGCTCGGCCACCAGGTTCATGCTGCGCACATGTTCGATCAGTGCCTGGCCGTCGCCCTCGGGCAGGGTGGTGTGCGGCGAGTAGGGCACTTGGCGCagcagc                                                                                                                                             >  8:252547/1‑149 (MQ=255)
         ttGCCTTCATCCAGGTAGAGGATGCGGCCCAGGGCGTCCTTCTGCTCGGCCACCAGGTTCATGCTGCGCACATGTTCGATCAGTGCCTGGCCGTCGCCCTCGGGCAGGGTGGTGTGCGGCGAGTAGGGCACTTGGCGCTGCTGCGCCaa                                                                                                                                        <  3:291955/149‑2 (MQ=255)
         ttGCCTTCATCCAGGTAGAGGATGCGGCCCAGGGCGTCCTTCTGCTCGGCCACCAGGTTCATGCTGCGCACATGTTCGATCAGTGCCTGGCCGTCGCCCTCGGGCAGGGTGGTGTGCGGCGAGTAGGGCACTTGGCGCAGCAGCGCCAg                                                                                                                                        >  3:13602/1‑149 (MQ=255)
            ccTTCATCCAGGTAGAGGATGCGGCCCTGGGCGTCCTTCTGCTCGGCCACCAGGTTCATGCTGCGCACATGTTCGATCAGTGCCTGGCCGTCGCCCTCGGGCAGGGTGGTGTGCGGAGAGTAGGGCACTTGGCGCAGCAGCGCGAGGTa                                                                                                                                     >  6:7837/1‑149 (MQ=255)
                 aTCCAGGTAGAGGATGCGGCCCAGGGCGTCCTTCTGCTCGGCCACCAGGTTCATGCTGCGCACATGTTCGATCAGTGCCTGGCCGTCGCCCTCGGGCAGGGTGGTGTGCGGCGAGTAGGGCACTTGGCGCAGCAGCGCCAGGTACAGGt                                                                                                                                >  8:272983/1‑149 (MQ=255)
                  tCCAGGTAGAGGATGCGGCCCAGGGCGTCCTTCTGCTCGGCCACCAGGTTCATGCTGCGCACATGTTCGATCAGTGCCTGGCCGTCGCCCTCGGGGAGGGTGGTGTGCGGCGAGTAGGGCACTTGGCGCAGCAGCGCCAGGTACAGGTc                                                                                                                               >  1:268162/1‑149 (MQ=255)
                    cAGGTAGAGGATGCGGCCCAGGGCGTCCTTCTGCTCGGCCACCAGGTTCATGCTGCGCACATGTTCGATCAGTGCCTGGCCGTCGCCCTCGGGCAGGGTGGTGTGCGGCGAGTAGGGCACTTGGCGCAGCAGCGCCAGGTACAGGTc                                                                                                                               >  7:232802/1‑147 (MQ=255)
                                  ggCCCAGGGCGTCCTTCTGCTCGGCCACCAGGTTCATGCTGCGCACATGTTCGATCAGTGCCTGGCCGTCGCCCTCGGGCAGGGTGGTGTGCGGCGAGTAGGGCACTTGGCgcag                                                                                                                                                 >  6:276137/1‑115 (MQ=255)
                                  ggCCCAGGGCGTCCTTCTGCTCGGCCACCAGGTTCATGCTGCGCACATGTTCGATCAGTGCCTGGCCGTCGCCCTCGGGCAGGGTGGTGTGCGGCGAGTAGGGCACTTGGCgcag                                                                                                                                                 <  5:276137/115‑1 (MQ=255)
                                  ggCCCAGGGCGTCCTTCTGCTCGGCCACCAGGTTCATGCTGCGCACATGTTCGATCAGTGCCTGGCCGTCGCCCTCGGGCAGGGTGGTGTGCGGCGAGTAGGGCACTTGGCGCAGCAGCGCCAGGTACAGGTCGAGTACGCGGGTCAgg                                                                                                               <  2:34406/149‑1 (MQ=255)
                                    cccAGGGCGTCCTTCTGCTCGGCCACCAGGTTCATGCTGCGCACATGTTCGATCAGTGCCTGGCCGTCGCCCTCGGGCAGGGTGGTGTGCGGCGAGTAGGGCACTTGGCGCAGCAGCGCCAGGTACAGGTCGAGTACGCGGGTCAGGgc                                                                                                             <  5:173869/149‑1 (MQ=255)
                                    cccAGGGCGTCCTTCTGCTCGGCCACCAGGTTCATGCTGCGCACATGTTCGATCAGTGCCTGGCCGTCGCCCTCGGGCAGGGTGGTGTGCGGCGAGTAGGGCACTTGGCGCAGCAGCGCCAGGTACAGGTCGAGTACGCGGGTCAGGgc                                                                                                             >  8:97160/1‑149 (MQ=255)
                                      cAGGGCGTCCTTCTGCTCGGCCACCAGGTTCATGCTGCGCACATGTTCGATCAGTGCCTGGCCGTCGCCCTCGGGCAGGGTGGTGTGCGGCGAGTAGGGCACTTGGCGCAGCAGCGCCAGGTACAGGTCGAGTACGCGGGTCAGGgcgc                                                                                                           <  8:75132/149‑1 (MQ=255)
                                       aGGGCGTCCTTCTGCTCGGCCACCAGGTTCATGCTGCGCACATGTTCGATCAGTGCCTGGCCGTCGCCCTCGGGCAGGGTGGTGTGCGGCGAGTAGGGCACTTGGCGCAGCAGCGCCAGGTACAGGTCGAGTACGCGGGTCAg                                                                                                                >  1:325822/1‑143 (MQ=255)
                                       aGGGCGTCCTTCTGCTCGGCCACCAGGTTCATGCTGCGCACATGTTCGATCAGTGCCTGGCCGTCGCCCTCGGGCAGGGTGGTGTGCGGCGAGTAGGGCACTTGGCGCAGCAGCGCCAGGTACAGGTCGAGTACGCGGGTCAg                                                                                                                <  2:325822/143‑1 (MQ=255)
                                               cTTCTGCTCGGCCACCAGGTTCATGCTGAGCACATGTTCGATCAGTGCCTGGCCGTCGCCCTCGGGCAGGGTGGTGTGCGGCGAGTAGGGCACTTGGCGCAGCAGCGCCAGGTACAGGTCGAGTACGCGGGTCAGGGCGCGCTCGTc                                                                                                    >  1:19158/1‑147 (MQ=255)
                                                  cTGCTCGGCCACCAGGTTCATGCTGCGCACATGTTCGATCAGTGCCTGGCCGTCGCCCTCGGGCAGGGTGGTGTGCGGCGAGTAGGGCACTTGGCGCAGCAGCGCCAGGTACAGGTCGAGTACGCGGGTCAGGGCGCGCTCGTCCAgg                                                                                                <  6:135533/148‑1 (MQ=255)
                                                  cTGCTCGGCCACCAGGTTCATGCTGCGCACATGTTCGATCAGTGCCTGGCCGTCGCCCTCGGGCAGGGTGGTGTGCGGCGAGTAGGGCACTTGGCGCAGCAGCGCCAGGTACAGGTCGAGTACGCGGGTCAGGGCGCGCTCGTCCAg                                                                                                 <  7:72272/147‑1 (MQ=255)
                                                   tGCTCGGCCACCAGGTTCATGCTGCGCACATGTTCGATCAGTGCCTGGCCGTCGCCCTCGGGCAGGGTGGTGTGCGGCGAGTAGGGCACTTGGCGCAGCAGCGCCAGGTACAGGTCGAGTACGCGGGTCAGGGCGCGCTCGTCCAg                                                                                                 >  6:304413/1‑146 (MQ=255)
                                                   tGCTCGGCCACCAGGTTCATGCTGCGCACATGTTCGATCAGTGCCTGGCCGTCGCCCTCGGGCAGGGTGGTGTGCGGCGAGTAGGGCACTTGGCGCAGCAGCGCCAGGTACAGGTCGAGTACGCGGGTCAGGGCGCGCTCGTCCAg                                                                                                 <  5:304413/146‑1 (MQ=255)
                                                     cTCGGCCACCAGGTTCATGCTGCGCACATGTTCGATCAGTGCCTGGCCGTCGCCCTCGGGCAGGGTGGTGTGCGGCGAGTAGGGCACTTGGCGCAGCAGCGCCAGGTACAGGTCGAGTACGCGGGTCAGGGCGCGCTCGTCCAGGGc                                                                                              <  2:19158/147‑1 (MQ=255)
                                                        ggCCACCAGGTTCATGCTGCGCACATGTTCGATCAGTGCCTGGCGGTCGCCCTCGGGCAGGGTGGTGTGCGGCGAGTAGGGCACTTGGCGCAGCAGCGCCAGGTACAGGTCGAGTACGCGGGTCAGGGCGCGCTCGTCCAGGGCAAgg                                                                                          <  5:191797/148‑1 (MQ=255)
                                                           caccaGGTTCATGCTGCGCACATGTTCGATCAGTGCCTGGCCGTCGCCCTCGGGCAGGGTGGTGTGCGGCGAGTAGGGCACTTGGCGCAGCAGCGCCAGGTACAGGTCGAGTACGCGGGTCAGGGCGCGCTCGTCCAGGGCAAggcggc                                                                                      >  4:304867/1‑149 (MQ=255)
                                                            accaGGTTCATGCTGCGCACATGTTCGATCAGTGCCTGGCCGTCGCCCTCGGGCAGGGTGGTGTGCGGCGAGTAGGGCACTTGGCGCAGCAGCGCCAGGTACAGGTCGAGTACGCGGGTCAGGGCGCGCTCGTCCAGGGCAAGGCGGCt                                                                                     >  1:35532/1‑149 (MQ=255)
                                                             ccaGGTTCATGCTGCGCACATGTTCGATCAGTGCCTGGCCGTCGCCCTCGGGCAGGGTGGTGTGCGGCGAGTAGGGCACTTGGCGCAGCAGCGCCAGGTACAGGTCGAGTACGCGGGTCAGGGCGCGCTCGTCCAGGGCAAGGCGGCtg                                                                                    >  7:126811/1‑149 (MQ=255)
                                                              caGGTTCATGCTGCGCACATGTTCGATCAGTGCCTGGCCGTCGCCCTCGGGCAGGGTGGTGTGCGGCGAGTAGGGCACTTGGCGCAGCAGCGCCAGGTACAGGTCGAGTACGCGGGTCAGGgcgcgc                                                                                                         <  4:225595/127‑1 (MQ=255)
                                                              caGGTTCATGCTGCGCACATGTTCGATCAGTGCCTGGCCGTCGCCCTCGGGCAGGGTGGTGTGCGGCGAGTAGGGCACTTGGCGCAGCAGCGCCAGGTACAGGTCGAGTACGCGGGTCAGGGCGCGCTCGTCCAGGGCAAGGCGGCtgg                                                                                   <  7:252547/149‑1 (MQ=255)
                                                               aGGTTCATGCTGCGCACATGTTCGATCAGTGCCTGGCCGTCGCCCTCGGGCAGGGTGGTGTGCGGCGAGTAGGGCACTTGGCGCAGCAGCGCCAGGTACAGGTCGAGTACGCGGGTCAGGGCGCGCTCGTCCAGGGCAAGGCGGCtggt                                                                                  <  7:272983/149‑1 (MQ=255)
                                                                ggTTCATGCTGCGCACATGTTCGATCAGTGCCTGGCCGTCGCCCTCGGGCAGGGTGGTGTGCGGCGAGTAGGGCACTTGGCGCAGCAGCGCCAGGTACAGGTCGAGTACGCGGGTCAGGGCGCGCTCGTCCAGGGCAAGGCGGCtggtg                                                                                 <  4:131648/149‑1 (MQ=255)
                                                                           cgcACATCTTCGATCAGTGCCTGGCCGTCGCCCGCGGGCAGGGGGGGGTGCGGCGAGTCGGGCACTTGGCGCAGCAGCGCCAGGTACAGGTCGAGTACGCGGGTCAGGGCGCGCTCGTCCAGGGCAAGGCGGCTGGTGGACAGCAGt                                                                        <  5:145/147‑1 (MQ=255)
                                                                            gcACATGTTCGATCAGTGCCTGGCCGTCGCCCTCGGGCAGGGTGGTGTGCGGCGAGTAGGGCACTTGGCGCAGCAGCGCCAGGTACAGGTCGAGTACGCGGGTCAGGGCGCGCTCGTCCAGGGCAAGGCGGCTGGTGGACAGCAGTGc                                                                      >  3:84742/1‑148 (MQ=255)
                                                                              acaTGTTCGATCAGTGCCTGGCCGTCGCCCTCGGGCAGGGTGGTGTGCGGCGAGTAGGGCACATGGCGCAGCAGCGCCAGGTACAGGTCGAGTACGCGGGTCAGGGCGCGCTCGTCCAGGGCAAGGCGGCTGGTGGACAGCAGTGCCAg                                                                   >  1:95649/1‑149 (MQ=255)
                                                                               caTGTTCGATCAGTGCCTGGCCGTCGCCCTCGGGCAGGGTGGTGTGCGGCGAGTAGGGCACTTGGCGCAGCAGCGCCAGGTACAGGTCGAGTACGCGGGTCAGGGCGCGCTCGTCCAGGGCAAGGCGGCTGGTGGACAGCAGTGCCAgg                                                                  <  7:97160/149‑1 (MQ=255)
                                                                                  gTTCGATCAGTGCCTGGCCGTCGCCCTCGGGCAGGGTGGTGTGCGGCGAGTAGGGCACTTGGCGCAGCAGCGCCAGGTACAGGTCGAGTACGCGGGTCAGGGCGCGCTCGTCCAGGGCAAGGCGGCTGGTGGACAGCAGTGCCAGGGc                                                                >  2:87229/1‑148 (MQ=255)
                                                                                   ttCGATCAGTGCCTGGCCGTCGCCCTCGGGCAGGGTGGTGTGCGGCGAGTAGGGCACTTGGCGCAGCAGCGCCAGGTACAGGTCGAGTACGCGGGTCAGGGCGCGCTCGTCCAGGGCAAGGCGGCTGGTGGACAGCAGTGCCAGGGcca                                                              >  1:194901/1‑149 (MQ=255)
                                                                                      gATCAGTGCCTGGCCGTCGCCCTCGGGCAGGGTGGTGTGCGGCGAGTAGGGCACTTGGCGCAGCAGCGCCAGGTACAGGTCGAGTACGCGGGTCAAGGCGCGCTCGTCCAGGGCAAGGCGGCTGGTGGACAGCAGTGCCAGGGccacc                                                            <  1:87229/148‑1 (MQ=255)
                                                                                          aGTGCCTGGCCGTCGCCCTCGGGCAGGGTGGTGTGCGGCGAGTAGGGCACTTGGCGCAGCAGCGCCAGGTACAGGTCGAGTACGCGGGTCAGGGCGCGCTCGTCCAGGGCAAGGCGGCTGGTGGAcagca                                                                          >  4:267090/1‑130 (MQ=255)
                                                                                          aGTGCCTGGCCGTCGCCCTCGGGCAGGGTGGTGTGCGGCGAGTAGGGCACTTGGCGCAGCAGCGCCAGGTACAGGTCGAGTACGCGGGTCAGGGCGCGCTCGTCCAGGGCAAGGCGGCTGGTGGAcagca                                                                          <  3:267090/130‑1 (MQ=255)
                                                                                          aGTGCCTGGCCGTCGCCCTCGGGCAGGGTGGTGTGCGGCGAGTAGGGCACTTGGCGCAGCAGCGCCAGGTACAGGTCGAGTACGCGGGTCAGGGCGCGCTCGTCCAGGGCAAGGCGGCTGGTGGACAGCAGTGCCAGGGCCACCAGGtt                                                       <  4:13602/149‑1 (MQ=255)
                                                                                           gTGCCTGGCCGTCGCCCTCGGGCAGGGTGGTGTGCGGCGAGTAGGGCACTTGGCGCAGCAGCGCCAGGTACAGGTCGAGTACGCGGGTCAGGGCGCGCTCGTCCAGGGCAAGGCGGCTGGTGGACAGCAGTGCCAGGGCCACCAGGtt                                                       >  8:18777/1‑148 (MQ=255)
                                                                                            tGCCTGGCCGTCGCCCTCGGGCAGGGTGGTGTGCGGCGAGTAGGGCACTTGGCGCAGCAGCGCCAGGTACAGGTCGAGTACGCGGGTCAGGGCGCGCTCGTCCAGGGCAAGGCGGCTGGTGGACAGCAGTGCCAGGGCCACCAGGTTGa                                                     >  7:84328/1‑149 (MQ=255)
                                                                                            tGCCTGGCCGTCGCCCTCGGGCAGGGTGGTGTGCGGCGAGTAGGGCACTTGGCGCAGCAGCGCCAGGTACAGGTCGAGTACGCGGGTCAGGGCGCGCTCGTCCAGGGCAAGGCGGCTGGTGGACAGCAGTGCCAGGGCCACCAGGTTGa                                                     >  7:97520/1‑149 (MQ=255)
                                                                                                  gCCGTCGCCCTCGGGCAGGGTGGTGTGCGGCGAGTAGGGCACTTGGCGCAGCAGCGCCAGGTACAGGTCGAGTACGCGGGTCAGGGCGCGCTCGTCCAGGGCAAGGCGGCTGGTGGACAGCAGTGCCAGGGCCACCAGGTTGACCggg                                                <  8:232802/148‑1 (MQ=255)
                                                                                                   ccGTCGCCCTCGGGCAGGGTGGTGTGCGGCGAGTAGGGCACTTGGCGCAGCAGCGCCAGGTACAGGTCGAGTACGCGGGTCAGGGCGCGCTCGTCCAGGGCAAGGCGGCTGGTGGACAGCAGTGCCAGGGCCACCAGGTTGACCGGGtt                                              <  5:106899/149‑1 (MQ=255)
                                                                                                    cGTCGCCCTCGGGCAGGGTGGTGTGCGGCGAGTAGGGCACTTGGCGCAGCAGCGCCAGGTACAGGTCGAGTACGCGGGTCAGGGCGCGCTCGTCCAGGGCAAGGCGGCTGGTGGACAGCAGTGCCAGGGCCACCAGGTTGACCGGGtt                                              >  5:94291/1‑148 (MQ=255)
                                                                                                         cccTCGGGCAGGGTGGTGTGCGGCGAGTAGGGCACTTGGCGCAGCAGCGCCAGGTACAGGTCGAGTACGCGGGTCAGGGCGCGCTCGTCCAGGGCAAGGCGGCTGGTGGACAGCAGTGCCAGGGCCACCAGGTTGACCGGGTTGATggc                                        >  4:307559/1‑149 (MQ=255)
                                                                                                         cccTCGGGCAGGGTGGTGTGCGGCGAGTAGGGCACTTGGCGCAGCAGCGCCAGGTACAGGTCGAGTACGCGGGTCAGGGCGCGCTCGTCCAGGGCAAGGCGGCTGGTGGACAGCAGTGCCAGGGCCACCAGGTTGACCGGGTTGATggc                                        >  4:302786/1‑149 (MQ=255)
                                                                                                         cccTCGGGCAGGGTGGTGTGCGGCGAGTAGGGCACTTGGCGCAGCAGCGCCAGGTACAGGTCGAGTACGCGGGTCAGGGCGCGCTCGTCCAGGGCAAGGCGGCTGGTGGACAGCAGTGCCAGGGCCACCAGGTTGACCGGGTTGATggc                                        <  4:137045/149‑1 (MQ=255)
                                                                                                          ccTCGGGCAGGGTGGTGTGCGGCGAGTAGGGCACATGGCGCAGCAGCGCCAGGTACAGGTCGAGTACGCGGGTCAGGGCGCGCTCGTCCAGGGCAAGGCGGCTGGTGGACAGCAGTGCCAGGGCCACCAGGTTGACCGGGTTGATggcg                                       >  1:210129/1‑149 (MQ=255)
                                                                                                          acTCGGGCAGGGTGGTGTGCGGCGAGTAGGGCACATGGCGCAGCAGCGCCAGGTACAGGTCGAGTACGCGGGTCAGGGCGCGCTCGTACAGGGCAAGGCGGCTGGTGGACAGCAGTGCCAGGGCCACCAGGTTGACAGGGTTGATggcg                                       >  1:214725/2‑149 (MQ=255)
                                                                                                               ggCAGGGTGGTGTGCGGCGAGTAGGGCACTTGGCGCAGCAGCGCCAGGTACAGGTCGAGTACGCGGGTCAGGGCGCGCTCGTCCAGGGCAAGGCGGCTGGTGGACAGCAGTGCCAGGGCCACCAGGTTGACCGGGTTGATGGCGgccgc                                  <  3:304867/149‑1 (MQ=255)
                                                                                                                    ggtggtGTGCGGCGAGTAGGGCACTTGGCGCAGCAGCGCCAGGTACAGGTCGAGTACGCGGGTCAGGGCGCGCTCGTCCAGGGCAAGGCGGCTGGTGGACAGCAGTGCCAGGGCCACCAGGTTGACCGGGTTGATGGCGGCCGCCGCGt                             <  6:94291/149‑1 (MQ=255)
                                                                                                                     gtggtgTGCGGCGAGTAGGGCACTTGGCGCAGCAGCGCCAGGTACAGGTCGAGTACGCGGGTCAGGGCGCGCTCGTCCAGGGCAAGGCGGCTGGTGGACAGCAGTGCCAGGGCCACCAGGTTGACCGGGTTGATGGCGGCCGCCTCGtt                            <  8:126811/149‑1 (MQ=255)
                                                                                                                          gtgCGGCGAGTAGGGCACTTGGCGCAGCAGCGCCAGGTACAGGTCGAGTACGCGGGTCAGGGCGCGCTCGTCCAGGGCAAGGCGGCTGGTGGACAGCAGTGCCAGGGCCACCAGGTTGACCGGGTTGATGGCGGCCGCCTCGTTGAGGt                       >  6:246801/1‑149 (MQ=255)
                                                                                                                             cggcgAGTAGGGCACTTGGCGCAGCAGCGCCAGGTACAGGTCGAGTACGCGGGTCAGGGCGCGCTCGTCCAGGGCAAGGCGGCTGGTGGACAGCAGTGCCAGGGCCACCAGGTTGACCGGGTTGATGGCGGCCGCCTCGTTGAGGTGGc                    <  3:218370/149‑1 (MQ=255)
                                                                                                                               gctaGTAGGGCACTTGGCGCAGCAGCGCCAGGTACATGTCGAGTACGCGGGTCAGGGCGCGCTCGTCCAGGTCAAGTCGGCTGGTGGACAGCAGTGCCAGGGCCACCAGGTTGACCGGGTTGATGGCGGCCGACTCGTTGAGGTGGCgg                  <  2:214725/146‑1 (MQ=255)
                                                                                                                               gcgAGTAGGGCACTTGGCGCAGCAGCGCCAGGTACAGGTCGAGTACGCGGGTCAGGGCGCGCTCGTCCAGGGCAAGGCGGCTGGTGGACAGCAGTGCCAGGGCCACCAGGTTGACCGGGTTGATGGCGGCCGCCTCGTTGAGGTGGCgg                  <  2:210129/149‑1 (MQ=255)
                                                                                                                                      gggCACTTGGCGCAGCAGCGCCAGGTACAGGTCGAGTACGCGGGTCAGGGCGCGCTCGTCCAGGGCAAGGCGGCTGGTGGACAGCAGTGCCAGGGCCACCAGGTTGACCGGGTTGATGGCGGCCGCCTCGTTGAGGTGGCGGGCCACg            <  2:35532/148‑1 (MQ=255)
                                                                                                                                       ggCACTTGGCGCAGCAGCGCCAGGTACAGGTCGAGTACGCGGGTCAGGGCGCGCTCGTCCAGGGCAAGGCGGCTGGTGGACAGCAGTGCCAGGGCCACCAGGTTGACGGGGTTGATGGCGGCCGCCTCGTTGAGGTGGCGGGCCACGGt          >  5:283687/1‑149 (MQ=255)
                                                                                                                                         cACTTGGCGCAGCAGCGCCAGGTACAGGTCGAGTACGCGGGTCAGGGCGCGCTCGTCCAGGGCACGGCGGCTGGTGGACAGCAGTGCCATGGACACCAGGTTGACCGGGTTGATGGCGGCCGCCGCGTTGAGGTGCCGGGCCACGGttt        >  3:184998/1‑149 (MQ=255)
                                                                                                                                         cACTTGGCGCAGCAGCGCCAGGTACAGGTCGAGTACGCGGGTCAGGGCGCGCTCGTCCAGGGCAAGGCGGCTGGTGGACAGCAGTGCCAGGGCCACCAGGTTGACCGGGTTGATGGCGGCCGCCTCGTTGAGGTGGCGGGCCACGGttt        >  6:218006/1‑149 (MQ=255)
                                                                                                                                          aCTTGGCGCAGCAGCGCCAGGTACAGGTCGAGTACGCGGGTCAGGGCGCGCTCGTCCAGGGCAAGGCGGCTGGTGGACAGCAGTGCCAGGGCCACCAGGTTGACCGGGTTGATGGCGGCCGCCTCGTTGAGGTGGCGGGCCACGGTTtc       >  1:136518/1‑149 (MQ=255)
                                                                                                                                            ctGGCGCTGCTGCGCCAGGTACAGGTCGAGTACGCGGGTCAGGGCGCGCTCGTCCAGGGCAAGGCGGCTGGTGGACAGCAGTGCCAGGGCCACCAGGTTGACCGGGTTGATggcg                                       <  5:181885/114‑1 (MQ=255)
                                                                                                                                            ctGGCGCTGCTGCGCCAGGTACAGGTCGAGTACGCGGGTCAGGGCGCGCTCGTCCAGGGCAAGGCGGCTGGTGGACAGCAGTGCCAGGGCCACCAGGTTGACCGGGTTGATggcg                                       >  6:181885/2‑115 (MQ=255)
                                                                                                                                            ctGGCGCTGCTGCGCCAGGTACAGGTCGAGTACGCGGGTCAGGGCGCGCTCGTCCAGGGCAAGGCGGCTGGTGGACAGCAGTGCCAGGGCCACCAGGTTGACCGGGTTGATGGCGGCCGCCTCGTTGAGGTGGCGGGCCACGGTTtctc     >  4:75613/2‑149 (MQ=255)
                                                                                                                                            ctGGCGCTGCTGCGCCAGGTACAGGTCGAGTACGCGGGTCAGGGCGCGCTCGTCCAGGGCAAGGCGGCTGGTGGACAGCAGTGCCAGGGCCACCAGGTTGACCGGGTTGATGGCGGCCGCCTCGTTGAGGTGGCGGGCCACGGTTtctc     >  7:111829/2‑149 (MQ=255)
                                                                                                                                               gcgcAGCAGCGCCAGGTACAGGTCGAGTACGCGGGTCAGGGCGCGCTCGTCCAGGGCAAGGCGGCTGGTGGACAGCAGTGCCAGGGCCACCAGGTTGACCGGGTTGATGGCGGCCGCCTCGTTGAGGTGGCGGGCCACGGTTTCTCCCa  <  2:95649/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                
AGCACCGCGTTGCCTTCATCCAGGTAGAGGATGCGGCCCAGGGCGTCCTTCTGCTCGGCCACCAGGTTCATGCTGCGCACATGTTCGATCAGTGCCTGGCCGTCGCCCTCGGGCAGGGTGGTGTGCGGCGAGTAGGGCACTTGGCGCAGCAGCGCCAGGTACAGGTCGAGTACGCGGGTCAGGGCGCGCTCGTCCAGGGCAAGGCGGCTGGTGGACAGCAGTGCCAGGGCCACCAGGTTGACCGGGTTGATGGCGGCCGCCTCGTTGAGGTGGCGGGCCACGGTTTCTCCCA  >  NC_002947/1726338‑1726629

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: ‑