Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,731,476 A→T 5.6% 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,731,4760AT 5.6% 1051.8 / 34.1 303L560Q (CTG→CAG) plsBglycerol‑3‑phosphate acyltransferase
Reads supporting (aligned to +/- strand):  ref base A (164/122);  new base T (8/9);  total (172/131)
Fisher's exact test for biased strand distribution p-value = 4.56e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.90e-01

Only 100 of 316 total aligned reads displayed.

TCAGCACCGCGTTGCCTTCATCCAGGTAGAGGATGCGGCCCAGGGCGTCCTTCTGCTCGGCCACCAGGTTCATGCTGCGCACATGTTCGATCAGTGCCTGGCCGTCGCCCTCGGGCAGGGTGGTGTGCGGCGAGTAGGGCACTTGGCGCAGCAGCGCCAGGTACAGGTCGAGTACGCGGGTCAGGGCGCGCTCGTCCAGGGCAAGGCGGCTGGTGGACAGCAGTGCCAGGGCCACCAGGTTGACCGGGTTGATGGCGGCCGCCTCGTTGAGGTGGCGGGCCACGGTTTCTCCCAGGC  >  NC_002947/1731327‑1731623
                                                                                                                                                     |                                                                                                                                                   
tcaGCACCGCGTTGCCTTCATCCAGGTAGAGGATGCGGCCCAGGGCGTCCTTCTGCTCGGCCACCAGGTTCATGCTGCGCACATGTTCGATCAGTGCCTGGCCGTCGCCCTCGGGCAGGGTGGTGTGCGGCGAGTAGGGCACTTGGCgca                                                                                                                                                     >  1:5824835/1‑150 (MQ=255)
  aGCACCGCGTTGCCTTCATCCAGGTAGAGGATGCGGTCCAGGGCGTCCTTCTGCTCGGCCACCAGGTTCATGCTGCGCACATGTTCGATCAGTGCCTGGCCGTCGCCCTCGGGCAGGGTGGTGTGCGGCGAGTAGGGCACTTGGCgcagc                                                                                                                                                   >  1:4748387/1‑150 (MQ=255)
   gCACCGCGTTGCCTTCATCCAGGTAGAGGATGCGGCCCAGGGCGTCCTTCTGCTCGGCCACCAGGTTCATGCTGCGCACATGTTCGATCAGTGCCTGGCCGTCGCCCTCGGGCAGGGTGGTGTGCGGCGAGTAGGGCACTTGGCgcagca                                                                                                                                                  <  1:4081521/150‑1 (MQ=255)
     aCCGCGTTGCCTTCATCCAGGTAGAGGATGTGGCCCAGGGCGTCCTTCTGCTCGGCCACCAGGTTCATGCTGCGCACATGTTCGATCAGTGCCTGGCCGTCGCCCTCGGGCAGGGTGGTGTGCGGCGAGTAGGGCACTTGGCGCagcagc                                                                                                                                                >  2:6227165/1‑150 (MQ=255)
     aCCGCGTTGCCTTCATCCAGGTAGAGGATGCGGCCCAGGGCGTCCTTCTGCTCGGCCACCAGGTTCATGCTGCGCACATGTTCGATCAGTGCCTGGCCGTCGCCCTCGGGCAGGGTGGTGTGCGGCGAGTAGGGCACTTGGCGCagcagc                                                                                                                                                >  1:965381/1‑150 (MQ=255)
     aCCGCGTTGCCTTCATCCAGGTAGAGGATGCGGCCCAGGGCGTCCTTCTGCTCGGCCACCAGGTTCATGCTGCGCACATGTTCGATCAGTGCCTGGCCGTCGCCCTCGGGCAGGGTGGTGTGCGGCGAGTAGGGCACTTGGCGCagcagc                                                                                                                                                <  2:2880731/150‑1 (MQ=255)
       cgcgTTGCCTTCATCCAGGTAGAGGATGCGGCCCAGGGCGTCCTTCTTCTCGGCCACCAGGTTCATGCTGCGCACATGTTCGATCAGTGCCTGGCCGTCGCCCTCGGGCAGGGTGGTGTGCGGCGAGTAGGGCACTTGGCGCAGCAgcgc                                                                                                                                              <  1:2661829/150‑1 (MQ=255)
        gcgTTGCCTTCATCCAGGTAGAGGATGCGGCCCAGGGCGTCCTTCTGCTCGGCCACCAGGTTCATGCTGCGCACATGTTCGATCAGTGCCTGGCCGTCGCCCTCGGGCAGGGTGGTGTGCGGCGAGTAGGGCACTTGGCGCAGCAGCGcc                                                                                                                                             >  2:4718570/1‑150 (MQ=255)
          gTTGCCTTCATCCAGGTAGAGGATGCGGCCCAGGGCGTCCTTCTGCTCGGCCACCAGGTTCATGCTGCGCACATGTTCGATCAGTGCCTGGCCGTCGCCCTCGGGCAGGGTGGTGTGCGGCGAGTAGGGCACTTGGCGCTGCTGCGCCaa                                                                                                                                           <  1:809667/150‑2 (MQ=255)
          gTTGCCTTCATCCAGGTAGAGGATGCGGCCCAGGGCGTCCTTCTGCTCGGCCACCAGGTTCATGCTGCGCACATGTTCGATCAGTGCCTGGCCGTCGCCCTCGGGCAGGGTGGTGTGCGGCGAGTAGGGCACTTGGCGCTGCTGCGCCaa                                                                                                                                           <  1:7196150/150‑2 (MQ=255)
          gTTGCCTTCATCCAGGTAGAGGATGCGGCCCAGGGCGTCCTTCTGCTCGGCCACCAGGTTCATGCTGCGCACATGTTCGATCAGTGCCTGGCCGTCGCCCTCGGGCAGGGTGGTGTGCGGCGAGTAGGGCACTTGGCGCAGCAGCGCCAg                                                                                                                                           >  1:5418908/1‑150 (MQ=255)
             gCCTTCATCCAGGTAGAGGATGCGGCCCAGGGCGTCCTTCTGCTCGGCCACCAGGTTCATGCTGCGCACATGTTCGATCAGTGCCTGGCCGTCGCCCTCGGGCAGGGTGGTGTGCGGCGAGTAGGGCACCTGGCGCTGCTGCGCCAAGtg                                                                                                                                        <  1:1038050/150‑2 (MQ=255)
               cTTCATCCAGGTAGAGGATGCGGCCCAGGGCGTCCTTCTGCTCGGCCACCAGGTTCATGCTGCGCACATGTTCGATCAGTGCCTGGCCGTCGCCCTCGGGCAGGGTGGTGTGCGGCGAGTAGGGCACTTGGCGCAGCAGCGCCAGGTACa                                                                                                                                      >  2:5384192/1‑150 (MQ=255)
               cTTCATCCAGGTAGAGGATGCGGCCCAGGGCGTCCTTCTGCTCGGCCACCAGGTTCATGCTGCGCACATGTTCGATCAGTGCCTGGCCGTCGCCCTCGGGCAGGGTGGTGTGCGGCGAGTAGGGCACTTGGCGCAGCAGCGCCAGGTACa                                                                                                                                      >  1:1474184/1‑150 (MQ=255)
                ttCATCCAGGTAGAGGATGCGGCCCAGGGCGTCCTTCTGCTCGGCCACCAGGTTCATGCTGCGCACATGTTCGATCAGTGCCTGGCCGTCGCCCTCGGGCAGGGTGGTGTGCGGCGAGTAGGGCACTTGGCGCAGCAGCGCCAGGTACAg                                                                                                                                     >  2:240670/1‑150 (MQ=255)
                ttCATCCAGGTAGAGGATGCGGCCCAGGGCGTCCTTCTGCTCGGCCACCAGGTTCATGCTGCGCACATGTTCGATCAGTGCCTGGCCGTCGCCCTCGGGCAGGGTGGTGTGCGGCGAGTAGGGCACTTGGCGCAGCAGCGCCAGGTACAg                                                                                                                                     <  1:6237762/150‑1 (MQ=255)
                ttCATCCAGGTAGAGGATGCGGCCCAGGGCGTCCTTCTGCTCGGCCACCAGGTTCATGCTGCGCACATGTTCGATCAGTGCCTGGCCGTCGCCCTCGGGCAGGGTGGTGTGCGGCGAGTAGGGCACTTGGCGCAGCAGCGCCAGGTACAg                                                                                                                                     >  2:4440064/1‑150 (MQ=255)
                ttCATCCAGGTAGAGGATGCGGCCCAGGGCGTCCTTCTGCTCGGCCACCAGGTTCATGCTGCGCACATGTTCGATCAGTGCCTGGCCGTCGCCCTCGGGCAGGGTGGTGTGCGGCGAGTAGGGCACTTGGCGCAGCAGCGCCAGGTACAg                                                                                                                                     >  2:2991009/1‑150 (MQ=255)
                   aTCCAGGTAGAGGATGCGGCCCAGGGCGTCCTTCTGCTCGGCCACCAGGTTCATGCTGCGCACATGTTCGATCAGTGCCTGGCCGTCGCCCTCGGGCAGGGTGGTGTGCGGCGAGTAGGGCACTTGGCGCAGCAGCGCCAGGTACAGGTc                                                                                                                                  <  1:2060526/150‑1 (MQ=255)
                     ccAGGTAGAGGATGCGGCCCAGGGCGTCCTTCTGCTCGGCCACCAGGTTCATGCTGCGCACATGTTCGATCAGTGCCTGGCCGTCGCCCTCGGGCAGGGTGGTGTGCGGCGAGTAGGGCACTTGGCGCAGCAGCGCCAGGTACAGGTCGa                                                                                                                                <  2:3560361/150‑1 (MQ=255)
                           agagGATGCGGCCCAGGGCGTCCTTCTGCTCGGCCACCAGGTTCATGCTGCGCACATGTTCGATCAGTGCCTGGCCGTCGCCCTCGGGCAGGGTGGTGTGCGGCGAGTAGGGCACTTGGCGCAGCAGCGCCAGGTACAGGTCGAGTAcgc                                                                                                                          <  1:3589300/150‑1 (MQ=255)
                             agGATGCGGCCCAGGGCGTCCTTCTGCTCGGCCACCAGGTTCATGCTGCGCACATGTTCGATCAGTGCCTGGCCGTCGCCCTCGGGCAGGGTGGTGTGCGGCGAGTAGGGCACTTGGCGCAGCAGCGCCAGGTACAGGTCGAGTACGCgg                                                                                                                        <  1:6137788/150‑1 (MQ=255)
                             agGATGCGGCCCAGGGCGTCCTTCTGCTCGGCCACCAGGTTCATGCTGCGCACATGTTCGATCAGTGCCTGGCCGTCGCCCTCGGGCAGGGTGGTGTGCGGCGAGTAGGGCACTTGGCGCAGCAGCGCCAGGTACAGGTCGAGTACGCgg                                                                                                                        >  1:6919399/1‑150 (MQ=255)
                                    ggCCCAGGGCGTCCTTCTGCTCGGCCACCAGGTTCATGCTGCGCACATGTTCGATCAGTGCCTGGCCGTCGCCCTCGGGCAGGGTGGTGTGCGGCGAGTAGGGCACTTGGCGCAGCAGCGCCAGGTACAGGTCGAGTACGCGGGTCAggg                                                                                                                 >  2:2646833/1‑150 (MQ=255)
                                    ggCCCAGGGCGTCCTTCTGCTCGGCCACCAGGTTCATGCTGCGCACATGTTCGATCAGTGCCTGGCCGTCGCCCTCGGGCAGGGTGGTGTGCGGCGAGTAGGGCACTTGGCGCAGCAGCGCCAGGTACAGGTCGAGTACGCGGGTCAggg                                                                                                                 >  2:1840611/1‑150 (MQ=255)
                                     gCCCAGGGCGTCCTTCTGCTCGGCCACCAGGTTCATGCTGCGCACATGTTCGATCAGTGCCTGGCCGTCGCCCTCGGGCAGGGTGGTGTGCGGCGAGTAGGGCACTTGGCGCAGCAGCGCCAGGTACAGGTCGAGTACGCGGGTCAGGgc                                                                                                                <  2:1165217/150‑1 (MQ=255)
                                     gCCCAGGGCGTCCTTCTGCTCGGCCACCAGGGTCATGCTGCGCACATGTTCGATCAGTGCCTGGCCGTCGCCCTCGGGCAGGGTGGTGTGCGGCGAGTAGGGCACTTGGCGCAGCAGCGCCAGGTACAGGTCGAGTACGCGGGTCAGGgc                                                                                                                <  1:3909517/150‑1 (MQ=255)
                                       ccAGGGCGTCCTTCTGCTCGGCCACCAGGTTCATGCTGCGCACATGTTCGATCAGTGCCTGGCCGTCGCCCTCGGGCAGGGTGGTGTGCGGCGAGTAGGGCACTTGGCGCAGCAGCGCCAGGTACAGGTCGAGTACGCGGGTCAGGgcgc                                                                                                              >  1:4120039/1‑150 (MQ=255)
                                         aGGGCGTCCTTCTGCTCGGCCACCAGGTTCATGCTGCGCACATGTTCGATCAGTGCCTGGCCGTCGCCCTCGGGCAGGGTGGTGTGCGGCGAGTAGGGCACTTGGCGCAGCAGCGCCAGGTACAGGTCGAGTACGCGGGTCAGGgcgcgc                                                                                                            <  2:4266529/150‑1 (MQ=255)
                                         aGGGCGTCCTTCTGCTCGGCCACCAGGTTCATGCTGCGCACATGTTCGATCAGTGCCTGGCCGTCGCCCTCGGGCAGGGTGGTGTGCGGCGAGTAGGGCACTTGGCGCAGCAGCGCCAGGTACAGGTCGAGTACGCGGGTCAGGgcgcgc                                                                                                            <  2:522222/150‑1 (MQ=255)
                                             cGTCCTTCTGCTCGGCCACCAGGTTCATGCTGCGCACATGTTCGATCAGTGCCTGGCCGTCGCCCTCGGGCAGGGTGGTGTGCGGCGAGTAGGGCACTTGGCGCAGCAGCGCCAGGTACAGGTCGAGTACGCGGGTCAGGGCGCGCTCGt                                                                                                        >  2:7186051/1‑150 (MQ=255)
                                                 cTTCTGCTCGGCCACCAGGTTCATGCTGCGCACATGTTCGATCAGTGCCTGGCCGTCGCCCTCGGGCAGGGTGGTGTGCGGCGAGTAGGGCACTTGGCGCAGCAGCGCCAGGTACAGGTCGAGTACGCGGGTCAGGGCGCGCTCGTCCAg                                                                                                    >  1:7371367/1‑150 (MQ=255)
                                                          ggCCACCAGGTTCATGCTGCGCACATGTTCGATCAGTGCCTGGCCGTCGCCCTCGGGCAGGGTGGTGTGCGGCGAGTAGGGCACTTGGCGCAGCAGCGCCAGGTACAGGTCGAGTACGCGGGTCAGGGCGCGCTCGTCCAGGGCAAggcg                                                                                           >  2:3331823/1‑150 (MQ=255)
                                                              accaGGTTCATGCTGCGCACATGTTCGATCAGTGCCTGGCCGTCGCCCTCGGGCAGGGTGGTGTGCGGCGAGTAGGGCACTTGGCGCAGCAGCGCCAGGTACAGGTCGAGTACGCGGGTCAGGGCGCGCTCGTCCAGGGCAAGGCGGCtg                                                                                       >  1:5098646/1‑150 (MQ=255)
                                                              accaGGTTCATGCTGCGCACATGTTCGATCAGTGCCTGGCCGTCGCCCTCGGGCAGGGTGGTGTGCGGCGAGTAGGGCACTTGGCGCAGCAGCGCCAGGTACAGGTCGAGTACGCGGGTCAGGGCGCGCTCGTCCAGGGCAAGGCGGCtg                                                                                       >  2:6483768/1‑150 (MQ=255)
                                                              accaGGTTCATGCTGCGCACATGTTCGATCAGTGCCTGGCCGTCGCCCTCGGGCAGGGTGGTGTGCGGCGAGTAGGGCACTTGGCGCAGCAGCGCCAGGTACAGGTCGAGTACGCGGGTCAGGGCGCGCTCGTCCAGGGCAAGGCGGCtg                                                                                       >  2:4314498/1‑150 (MQ=255)
                                                              accaGGTTCATGCTGCGCACATGTTCGATCAGTGCCTGGCCGTCGCCCTCGGGCAGGGTGGTGTGCGGCGAGTAGGGCACTTGGCGCAGCAGCGCCAGGTACAGGTCGAGTACGCGGGTCAGGGCGCGCTCGTCCAGGGCAAGGCGGCtg                                                                                       >  1:7296742/1‑150 (MQ=255)
                                                               ccaGGTTCATGCTGCGCACATGTTCGATCAGTGCCTGGCCGTCGCCCTCGGGCAGGGTGGTGTGCGGCGAGTAGGGCACTTGGCGCAGCAGCGCCAGGTACAGGTCGAGTACGCGGGTCAGGGCGCGCTCGTCCAGGGCAAGGCGGCtgg                                                                                      <  2:6297958/150‑1 (MQ=255)
                                                                   gTTCATGCTGCGCACATGTTCGATCAGTGCCTGGCCGTCGCCCTCGGGCAGGGTGGTGTGCGGCGAGTAGGGCACTTGGCGCAGCAGCGCCAGGTACAGGTCGAGTACGCGGGTCAGGGCGCGCTCGTCCAGGGCAAGGCGGCTGGTGGa                                                                                  <  1:5990057/150‑1 (MQ=255)
                                                                   gTTCATGCTGCGCACATGTTCGATCAGTGCCTGGCCGTCGCCCTCGGGCAGGGTGGTGTGCGGCGAGTAGGGCACTTGGCGCAGCAGCGCCAGGTACAGGTCGAGTACGCGGGTCAGGGCGCGCTCGTCCAGGGCAAGGCGGCTGGTGGa                                                                                  >  2:1686930/1‑150 (MQ=255)
                                                                   gTTCATGCTGCGCACATGTTCGATCAGTGCCTGGCCGTCGCCCTCGGGCAGGGTGGTGTGCGGCGAGTAGGGCACTTGGCGCAGCAGCGCCAGGTACAGGTCGAGTACGCGGGTCAGGGCGCGCTCGTCCAGGGCAAGGCGGCTGGTGGa                                                                                  <  2:5591393/150‑1 (MQ=255)
                                                                   gTTCATGCTGCGCACATGTTCGATCAGTGCCTGGCCGTCGCCCTCGGGCAGGGTGGTGTGCGGCGAGTAGGGCACTTGGCGCAGCAGCGCCAGGTACAGGTCGAGTACGCGGGTCAGGGCGCGCTCGTCCAGGGCAAGGCGGCTGGTGGa                                                                                  >  1:1132113/1‑150 (MQ=255)
                                                                     tCATGCTGCGCACATGTTCGATCAGTGCCTGGCCGTCGCCCTCGGGCAGGGTGGTGTGCGGCGAGTAGGGCACTTGGCGCAGCAGCGCCAGGTACAGGTCGAGTACGCGGGTCAGGGCGCGCTCGTCCAGGGCAAGGCGGCTGGTGGAca                                                                                >  2:3925649/1‑150 (MQ=255)
                                                                        tgctgcGCACATGTTCGATCAGTGCCTGGCCGTCGCCCTCGGGCAGGGTGGTGTGCGGCGAGTAGGGCACTTGGCGCAGCAGCGCCAGGTACAGGTCGAGTACGCGGGTCAGGGCGCGCTCGTCCAGGGCAAGGCGGCTGGTGGAcagca                                                                             <  2:6696811/150‑1 (MQ=255)
                                                                        tgctgcGCACATGTTCGATCAGTGCCTGGCCGTCGCCCTCGGGCAGGGTGGTGTGCGGCGAGTAGGGCACTTGGCGCAGCAGCGCCAGGTACAGGTCGAGTACGCGGGTCAGGGCGCGCTCGTCCAGGGCAAGGCGGCTGGTGGAcagca                                                                             <  1:5574746/150‑1 (MQ=255)
                                                                          ctgcGCACATGTTCGATCAGTGCCTGGCCGTCGCCCTCGGGCAGGGTGGTGTGCGGCGAGTAGGGCACTTGGCGCAGCAGCGCCAGGTACAGGTCGAGTACGCGGGTCAGGGCGCGCTCGTCCAGGGCAAGGCGGCTGGTGGACAGCAGt                                                                           <  1:2121209/150‑1 (MQ=255)
                                                                           tgcGCACATGTTCGATCAGTGCCTGGCCGTCGCCCTCGGGCAGGGTGGTGTGCGGCGAGTAGGGCACTTGGCGCAGCAGCGCCAGGTACAGGTCGAGTACGCGGGTCAGGGCGCGCTCGTCCAGGGCAAGGCGGCTGGTGGACAGCAGTg                                                                          >  1:904408/1‑150 (MQ=255)
                                                                            gcgcACATGTTCGATCAGTGCCTGGCCGTCGCCCTCGGGCAGGGTGGTGTGCGGCGAGTAGGGCACTTGGCGCAGCAGCGCCAGGTACAGGTCGAGTACGCGGGTCAGGGCGCGCTCGTCCAGGGCAAGGCGGCTGGTGGACAGCAGTGc                                                                         >  2:5825725/1‑150 (MQ=255)
                                                                                acctGTTCGATCAGTGCCTGGCCGTCGCCCTCGGGCAGGGTGGTGTGCGGCGAGTAGGGCACTTGGCGCAGCAGCGCCAGGTACAGGTCGAGTACGCGGGTCAGGGCGCGCTCGTCCAGGGCAAGGCGGCTGGTGGACAGCAGTGCCAgg                                                                     <  2:6409777/147‑1 (MQ=255)
                                                                                acaTGTTCGATCAGTGCCTGGCCGTCGCCCTCGGGCAGGGTGGTGTGCGGCGAGTAGGGCACTTGGCGCAGCAGCGCCAGGTACAGGTCGAGTACGCGGGTCAGGGCGCGCTCGTCCAGGGCAAGGCGGCTGGTGTACAGCAGTGCCAgg                                                                     >  2:1571634/1‑150 (MQ=255)
                                                                                  aTGTTCGATCAGTGCCTGGCCGTCGCCCTCGGGCAGGGTGGTGTGCGGCGAGTAGGGCACTTGGCGCAGCAGCGCCAGGTACAGGTCGAGTACGCGGGTCAGGGCGCGCTCGTCCAGGGCAAGGCGGCTGGTGGACAGCAGTGCCAGGGc                                                                   >  2:2426572/1‑150 (MQ=255)
                                                                                  aTGTTCGATCAGTGCCTGGCCGTCGCCCTCGGGCAGGGTGGTGTGCGGCGAGTAGGGCACTTGGCGCAGCAGCGCCAGGTACAGGTCGAGTACGCGGGTCAGGGCGCGCTCGTCCAGGGCAAGGCGGCTGGTGGACAGCAGTGCCAGGGc                                                                   >  2:372058/1‑150 (MQ=255)
                                                                                   tGTTCGATCAGTGCCTGGCCGTCGCCCTCGGGCAGGGTGGTGTGCGGCGAGTAGGGCACTTGGCGCAGCAGCGCCAGGTACAGGTCGAGTACGCGGGTCAGGGCGCGCTCGTCCAGGGCAAGGCGGCTGGTGGACAGCAGTGCCAGGGcc                                                                  >  2:804596/1‑150 (MQ=255)
                                                                                     ttCGATCAGTGCCTGGCCGTCGCCCTCGGGCAGGGTGGTGTGCGGCGAGTAGGGCACTTGGCGCAGCTGCGCCAGGTACAGGTCGAGTACGCGGGTCAGGGCGCGCTCGTCCAGGGCAAGGCGGCTGGTGGACAGCAGTGCCAGGGccac                                                                >  1:2191871/1‑150 (MQ=255)
                                                                                     ttCGATCAGTGCCTGGCCGTCGCCCTCGGGCAGGGTGGTGTGCGGCGAGTAGGGCACTTGGCGCAGCAGCGCCAGGTACAGGTCGAGTACGCGGGTCAGGgcgc                                                                                                              <  2:1114329/104‑1 (MQ=255)
                                                                                        gATCAGTGCCTGGCCGTCGCCCTCGGGCAGGGTGGTGTGCGGCGAGTAGGGCACTTGGCGCAGCAGCGCCAGGTACAGGTCGAGTACGCGGGTCAGGGCGCGCTCGTCCAGGGCAAGGCGGCTGGTGGACAGCAGTGCCAGGGCCACCAg                                                             >  2:5479816/1‑150 (MQ=255)
                                                                                         aTCAGTGCCTGGCCGTCGCCCTCGGGCAGGGTGGTGTGCGGCGAGTAGGGCACTTGGCGCAGCAGCGCCAGGTACAGGTCGAGTACGCGGGTCAGGGCGCGCTCGTCCAGGGCAAGGCGGCTGGTGGACAGCAGTGCCAGGGCCACCAgg                                                            <  1:4040212/150‑1 (MQ=255)
                                                                                            aGTGCCTGGCCGTCGCCCTCGGGCAGGGTGGTGTGCGGCGAGTAGGGCACTTGGCGCAGCAGCGCCAGGTACAGGTCGAGTACGCGGGTCAGGGCGCGCTCGTCCAGGGCAAGGCGGCTGGTGGACAGCAGTGCCAGGGCCACCAGGTTg                                                         <  1:1686930/150‑1 (MQ=255)
                                                                                              tGCCTGGCCGTCGCCCTCGGGCAGGGTGGTGTGCGGCGAGTAGGGCACTTGGCGCAGCAGCGCCAGGTACAGGTCGAGTACGCGGGTCAGGGCGCGCTCGTCCAGGGCAAGGCGGCTGGTGGACAGCAGTGCCAGGGCCACCAGGTTGAc                                                       <  1:4514170/150‑1 (MQ=255)
                                                                                                ccTGGCCGTCGCCCTCGGGCAGGGTGGTGTGCGGCGAGTAGGGCACTTGGCGCAGCAGCGCCAGGTACAGGTCGAGTACGCGGGTCAGGGCGCGCTCGTCCAGGGCAAGGCGGCTGGTGGAGAGCAGTGCCAGGGCCACCAGGTTGACCg                                                     >  2:1045303/1‑150 (MQ=255)
                                                                                                ccTGGCCGTCGCCCTCGGGCAGGGTGGTGTGCGGCGAGTAGGGCACTTGGCGCAGCAGCGCCAGGTACAGGTCGAGTACGCGGGTCAGGGCGCGCTCGTCCAGGGCAAGGCGGCTGGTGGACAGCAGTGCCAGGGCCACCAGGTTGACCg                                                     >  1:5288984/1‑150 (MQ=255)
                                                                                                  tGGCCGTCGCCCTCGGGCAGGGTGGTGTGCGGCGAGTAGGGCACTTGGCGCAGCAGCGCCAGGTACAGGTCGAGTACGCGGGTCAgg                                                                                                                  >  2:5781856/1‑87 (MQ=255)
                                                                                                  tGGCCGTCGCCCTCGGGCAGGGTGGTGTGCGGCGAGTAGGGCACTTGGCGCAGCAGCGCCAGGTACAGGTCGAGTACGCGGGTCAGGGCGCGCTCGTCCAGGGCAAGGCGGCTGGTGGACAGCAGTGCCAGGGCCACCAGGTTGACCggg                                                   >  1:2269043/1‑150 (MQ=255)
                                                                                                  tGGCCGTCGCCCTCGGGCAGGGTGGTGTGCGGCGAGTAGGGCACTTGGCGCAGCAGCGCCAGGTACAGGTCGAGTACGCGGGTCAGGGCGCGCTCGTCCAGGGCAAGGCGGCTGGTGGACAGCAGTGCCAGGGCCACCAGGTTGACCggg                                                   >  1:2974973/1‑150 (MQ=255)
                                                                                                  gggCCGTCGCCCTCGGGCAGGGTGGTGTGCGGCGAGTAGGGCACTTGGCGCAGCAGCGCCAGGTACAGGTCGAGTACGCGGGTCAGGGCGCGCTCGTCCAGGGCAAGGCGGCTGGTGGACAGCAGTGCCAGGGCCACCAGGTTGACCggg                                                   <  2:7244293/149‑1 (MQ=255)
                                                                                                    gCCGTCGCCCTCGGGCAGGGTGGTGTGCGGCGAGTAGGGCACTTGGCGCAGCAGCGCCAGGTACAGGTCGAGTACGCGGGTCAGGGCGCGCTCGTCCAGGGCAAGGCGGCTGGTGGACAGCAGTGCCAGGGCCACCAGGTTGACCGGGtt                                                 >  2:2752292/1‑150 (MQ=255)
                                                                                                           cccTCGGGCCGGGTGGTGTGCGGCGAGTAGGGCACTTGGCGCAGCAGCGCCAGGTACAGGTCGAGTACGCGGGTCAGGGCGCGCTCGTCCAGGGCAAGGCGGCTGGTGGACAGCAGTGCCAGGGCCACCAGGTTGACCGGGTTGATggcg                                          >  2:2113995/1‑150 (MQ=255)
                                                                                                           cccTCGGGCAGGGTGGTGTGCGGCGAGTAGGGCACTTGGCGCAGCAGCGCCAGGTACAGGTCGAGTACGCGGGTCAGGGCGCGCTCGTCCAGGGCAAGGCGGCTGGTGGACAGCAGTGCCAGGGCCACCAGGTTGACCGGGTTGATggcg                                          <  2:4090542/150‑1 (MQ=255)
                                                                                                            ccTCGGGCAGGGTGGTGTGCGGCGAGTAGGGCACTTGGCGCAGCAGCGCCAGGTACAGGTCGAGTACGCGGGTCAGGGCGCGCTCGTCCAGGGCAAGGCGGCTGGTGGACAGCAGTGCCAGGGCCACCAGGTTGACCGGGTTGATggcgg                                         >  1:2764287/1‑150 (MQ=255)
                                                                                                            ccTCGGGCAGGGTGGTGTGCGGCGAGTAGGGCACTTGGCGCAGCAGCGCCAGGTACAGGTCGAGTACGCGGGTCAGGGCGCGCTCGTCCAGGGCAAGGCGGCTGGTGGACAGCAGTGCCAGGGCCACCAGGTTGACCGGGTTGATggcgg                                         >  1:574638/1‑150 (MQ=255)
                                                                                                            ccTCGGGCAGGGTGGTGTGCGGCGAGTAGGGCACTTGGCGCAGCAGCGCCAGGTACAGGTCGAGTACGCGGGTCAGGGCGCGCTCGTCCAGGGCAAGGCGGCTGGTGGACAGCAGTGCCAGGGCCACCAGGTTGACCGGGTTGATggcgg                                         >  2:4864566/1‑150 (MQ=255)
                                                                                                             cTCGGGCAGGGTGGTGTGCGGCGAGTAGGGCACTTGGCGCAGCAGCGCCAGGTACAGGTCGAGTACGCGGGTCAGGGCGCGCTCGTCCAGGGCAAGGCGGCTGGTGGACAGCAGTGCCAGGGCCACCAGGTTGACCGGGTTGATggcggc                                        >  1:602670/1‑150 (MQ=255)
                                                                                                              tCGGGCAGGGTGGTGTGCGGCGAGTAGGGCACTTGGCGCAGCAGCGCCAGGTACAGGTCGAGTACGCGGGTCAGGGCGCGCTCGTCCAGGGCAAGGCGGCTGGTGGACAGCAGTGCCAGGGCCACCAGGTTGACCGGGTTGATGGCGgcc                                       <  1:3231315/150‑1 (MQ=255)
                                                                                                                gggCAGGGTGGTGTGCGGCGAGTAGGGCACTTGGCGCAGCAGCGCCAGGTACAGGTCGAGTACGCGGGTCAGGGCGCGCTCGTCCAGGGCAAGGCGGCTGGTGGACAGCAGTGCCAGGGCCACCAGGTTGACCGGGTTGATGGCGgccgc                                     >  1:4934650/1‑150 (MQ=255)
                                                                                                                      ggtggtGTGCGGCGAGTAGGGCACTTGGCGCAGCAGCGCCAGGTACAGGTCGAGTACGCGGGTCAGGGCGCGCTCGTCCAGGGCAAGGCGGCTGGTGGACAGCAGTGCCAGGGCCACCAGGTTGACCggg                                                   >  2:3833796/1‑130 (MQ=255)
                                                                                                                      ggtggtGTGCGGCGAGTAGGGCACTTGGCGCAGCAGCGCCAGGTACAGGTCGAGTACGCGGGTCAGGGCGCGCTCGTCCAGGGCAAGGCGGCTGGTGGACAGCAGTGCCAGGGCCACCAGGTTGACCggg                                                   <  1:3833985/130‑1 (MQ=255)
                                                                                                                      ggtggtGTGCGGCGAGTAGGGCACTTGGCGCAGCAGCGCCAGGTACAGGTCGAGTACGCGGGTCAGGGCGCGCTCGTCCAGGGCAAGGCGGCTGGTGGACAGCAGTGCCAGGGCCACCAGGTTGACCGGGTTGATGGCGGCCGCCTCGtt                               <  1:368112/150‑1 (MQ=255)
                                                                                                                       gtggtgTGCGGCGAGTAGGGCACGTGGCGCAGCAGCGCCAGGTACAGGTCGAGTACGCGGGTCAGGGCGCGCTCGTCCAGGGCAAGGCGGCTGGTGGACAGCAGTGCCAGGGCCACCAGGTTGACCGGGTTGATGGCGGCCGCCTCGTTg                              <  2:5969039/150‑1 (MQ=255)
                                                                                                                        tggtgTGCGGCGAGTAGGGCACTTGGCGCAGCAGCGCCAGGTACAGGTCGAGTACGCGGGTCAGGGCGCGCTCGTCCAGGGCAAGGCGGCTGGTGGACAGCAGTGCCAGGGCCACCAGGTTGACCGGGTTGATGGCGGCCGCCTCGTTGa                             <  2:6150268/150‑1 (MQ=255)
                                                                                                                         ggtgTGCGGCGAGTAGGGCACTTGGCGCAGCAGCGCCAGGTACAGGTCGAGTACGCGGGTCAGGGCGCGCTCGTCCAGGGCAAGGCGGCTGGTGGACAGCAGTGCCAGGGCCACCAGGTTGACCGGGTTGATGGCGGCCGCCTCGTTGAg                            <  1:1267698/150‑1 (MQ=255)
                                                                                                                         ggtgTGCGGCGAGTAGGGCACTTGGCGCAGCAGCGCCAGGTACAGGTCGAGTACGCGGGTCAGGGCGCGCTCGTCCAGGGCAAGGCGGCTGGTGGACAGCAGTGCCAGGGCCACCAGGTTGACCGGGTTGATGGCGGCCGCCTCGTTGAg                            <  1:6394503/150‑1 (MQ=255)
                                                                                                                            gtgCGGCGAGTAGGGCACTTGGCGCAGCAGCGCCAGGTACAGGTCGAGTACGCGGGTCAGGGCGCGCTCGTCCAGGGCAAGGCGGCTGGTGGACAGCAGTGCCAGGGCCACCAGGTTGACCGGGTTGATGGCGGCCGCCTCGTTGAGGTg                         <  2:741256/150‑1 (MQ=255)
                                                                                                                            gtgCGGCGAGTAGGGCACTTGGCGCAGCAGCGCCAGGTACAGGTCGAGTACGCGGGTCAGGGCGCGCTCGTCCAGGGCAAGGCGGCTGGTGGACAGCAGTGCCAGGGCCACCAGGTTGACCGGGTTGATGGCGGCCGCCTCGTTGAGGTg                         >  2:6533950/1‑150 (MQ=255)
                                                                                                                                  cgAGTAGGGCACTTGGCGCAGCAGCGCCAGGTACAGGTCGAGTACGCGGGTCAGGGCGCGCTCGTCCAGGGCAAGGCGGCTGGTGGACAGCAGTGCCAGGGCCACCAGGTTGACCGGGTTGATGGCGGCCGCCTCGTTGAGGTGGCGGGc                   >  1:3363048/1‑150 (MQ=255)
                                                                                                                                   gAGTAGGGCACTTGGCGCAGCAGCGCCAGGTACAGGTCGAGTACGCGGGTCa                                                                                                                    >  2:3091415/1‑52 (MQ=255)
                                                                                                                                        gggCACTTGGCGCAGCAGCGCCAGGTACAGGTCGAGTACGCGGGTCAGGGCGCGCTCGTCCAGGGCAAGGCGGCTGGTGGACAGCAGTGCCAGGGCCACCAGGTTGACCGGGTTGATGGCGGCCGCCTCGTTGAGGTGGCGGGCCACGGt             >  1:220213/1‑150 (MQ=255)
                                                                                                                                        gggCACTTGGCGCAGCAGCGCCAGGTACAGGTCGAGTACGCGGGTCAGGGCGCGCTCGTCCAGGGCAAGGCGGCTGGTGGACAGCAGTGCCAGGGCCACCAGGTTGACCGGGTTGATGGCGGCCGCCTCGTTGAGGTGGCGGGCCACGGt             <  1:7105065/150‑1 (MQ=255)
                                                                                                                                        ctgtacctGGCGCTGCTGCGCCAAGTACAGGTCGAGTACGCGGGTCAGGGCGCGCTCGTCCAGGGCAAGGCGGCTGGTGGACAGCAGTGCCAGGGCCACCAGGTTGACCGGGTTGATGGCGGCCGCCTCGTTGAGGTGGCGGGCCACGGt             >  2:708384/8‑150 (MQ=255)
                                                                                                                                        ctgCACTTGGCGCAGCAGCGCCAGGTACAGGTCGAGTACGCGGGTCAGGGCGCGCTCGTCCAGGGCAAGGCGGCTGGTGGACAGCAGTGCCAGGGCCACCAGGTTGACCGGGTTGATGGCGGCCGCCTCGTTGAGGTGGCGGGCCACGGt             >  1:3373343/3‑150 (MQ=255)
                                                                                                                                         ggCACTTGGCGCAGCAGCGCCAGGTACAGGTCGAGTACGCGGGTCAGGGCGCGCTCGTCCAGGGCAAGGCGGCTGGTGGACAGCAGTGCCAGGGCCACCAGGTTGACCGGGTTGATGGCGGCCGCCTCGTTGAGGt                          <  1:1396499/136‑1 (MQ=255)
                                                                                                                                         ggCACTTGGCGCAGCAGCGCCAGGTACAGGTCGAGTACGCGGGTCAGGGCGCGCTCGTCCAGGGCAAGGCGGCTGGTGGACAGCAGTGCCAGGGCCACCAGGTTGACCGGGTTGATGGCGGCCGCCTCGTGGAGGt                          >  2:1396499/1‑136 (MQ=255)
                                                                                                                                          gCACTTGGCGCAGCAGCGCCAGGTACAGGTCGAGTACGCGGGTCAGGGCGCGCTCGTCCAGGGCAAGGCGGCTGGTGGACAGCAGTGCCAGGGCCACCAGGTTGACCGGGTTGATGGCGGCCGCCTCGTTGAGGTGGCGGGCCACGGttt           <  1:186118/150‑1 (MQ=255)
                                                                                                                                            aCTTGGCGCAGCAGCGCCAGGTACAGGTCGAGTACGCGGGTCAGGGCGCGCTCGTCCAGGGCAAGGCGGCTGGTGGACAGCAGTGCCAGGGCCACCAGGTTGACCGGGTTGATGGCGGCCGCCTCGTTGAGGTGGCGGGCCACGGTTtct         <  1:6608975/150‑1 (MQ=255)
                                                                                                                                             cctGGCGCTGCTGCGCCAGGTACAGGTCGAGTACGCGGGTCAGGGCGCGCTCGTCCAGGGCAAGGCGGCTGGTGGACAGCAGTGCCAGGGCCACCAGGTTGACCGGGTTGATGGCGGCCGCCTCGTTGAGGTGGCGGGCCACGGTTtctc        >  2:3518746/3‑150 (MQ=255)
                                                                                                                                             cTTGGCGCAGCAGCGCCAGGTACAGGTCGAGTACGCGGGTCAGGGCGCGCTCGTCCAGGGCAAGGCGGCTGGTGGACAGCAGTGCCAGGGCCACCAGGTTGACCGGGTTGATGGCGGCCGCCTCGTTGAGGTGGCGGGCCACGGTTtctc        >  1:4286192/1‑150 (MQ=255)
                                                                                                                                              ttGGCGCAGCAGCGCCAGGTACAGGTCGAGTACGCGGGTCAGGGCGCGCTCGTCCAGGGCAAGGCGGCTGGTGGACAGCAGTGCCAGGGCCACCAGGTTGACCGGGTTGATGGCGGCCGCCTCGTTGAGGTGGCGGGCCACGGTTTCTcc       >  1:4421935/1‑150 (MQ=255)
                                                                                                                                              ttGGCGCAGCAGCGCCAGGTACAGGTCGAGTACGCGGGTCAGGGCGCGCTCGTCCAGGGCAAGGCGGCTGGTGGACAGCAGTGCCAGGGCCACCAGGTTGACCGGGTTGATGGCGGCCGCCTCGTTGAGGTGGCGGGCCACGGTTTCTcc       <  2:4966664/150‑1 (MQ=255)
                                                                                                                                                ggCGCAGCAGCGCCAGGTACAGGTCGAGTACGCGGGTCAGGGCGCGCTCGTCCAGGGCAAGGCGGCTGGTGGACAGCAGTGCCAGGGCCACCAGGTTGACCGGGTTGATGGCGGCCGCCTCGTTGAGGTGGCGGGCCACGGTTTCTCCCa     <  1:5918991/150‑1 (MQ=255)
                                                                                                                                                 gcgcAGCAGCGCCAGGTACAGGTCGAGTACGCGGGTCAGGGCGCGCTCGTCCAGGGCAAGGCGGCTGGTGGACAGCAGTGCCAGGGCCACCAGGTTGACCGGGTTGATGGCGGCCGCCTCGTTGAGGTGGCGGGCCACGGTTTCTCCCAg    >  1:74016/1‑150 (MQ=255)
                                                                                                                                                 gcgcAGCAGCGCCAGGTACAGGTCGAGTACGCGGGTCAGGGCGCGCTCGTCCAGGGCAAGGCGGCTGGTGGACAGCAGTGCCAGGGCCACCAGGTTGACCGGGTTGATGGCGGCCGCCTCGTTGAGGTGGCGGGCCACGGTTTCTCCCAg    <  1:1096658/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
TCAGCACCGCGTTGCCTTCATCCAGGTAGAGGATGCGGCCCAGGGCGTCCTTCTGCTCGGCCACCAGGTTCATGCTGCGCACATGTTCGATCAGTGCCTGGCCGTCGCCCTCGGGCAGGGTGGTGTGCGGCGAGTAGGGCACTTGGCGCAGCAGCGCCAGGTACAGGTCGAGTACGCGGGTCAGGGCGCGCTCGTCCAGGGCAAGGCGGCTGGTGGACAGCAGTGCCAGGGCCACCAGGTTGACCGGGTTGATGGCGGCCGCCTCGTTGAGGTGGCGGGCCACGGTTTCTCCCAGGC  >  NC_002947/1731327‑1731623

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