Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,632,777 C→T 5.8% A175T (GCC→ACC)  pgsA ← CDP‑diacylglycerol‑glycerol‑3‑phosphate 3‑phosphatidyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,632,7770CT 5.8% 202.2 / 2.0 69A175T (GCC→ACC) pgsACDP‑diacylglycerol‑glycerol‑3‑phosphate 3‑phosphatidyltransferase
Reads supporting (aligned to +/- strand):  ref base C (33/32);  new base T (2/2);  total (35/34)
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.37e-01

CCCCGACCTTGGCAAGGTCGTGCTCTACCAACTGAGCTATTCCCGCGTCGTGATGGGTGCCATTCTAGCGATTTCTAAAACAGCGTCAACCCCTTGATTCAAAAAACTTTTATTTTTTTTCCGAGCCTTCGCGCAGGTGCGGCCAGGCGGCCAGCAGGTAGTGCACCATCGACCACAGCGTCAGCCCCGCTGCCACCAGCAGCAGCCCGTAGCCCAGGATCACCCAGAAGGTCACCACCGGCGGGTTGCCCAGCAGGATCACCAGCGCCAGCATCTGCGCCGCAGTCTTCCACTTA  >  NC_002947/4632630‑4632925
                                                                                                                                                   |                                                                                                                                                    
ccccGACCTTGGCAAGGTCGTGCTCTACCAACTGAGCTATTCCCGCGTCGTGATGGGTGCCATTCTAGCGATTTCTAAAACAGCGTCAACCCCTTGATTCAAAAAACTTTTATTTTTTTTCCGAGCCTTCGCGCAGGTgcggccaggcgg                                                                                                                                                    >  2:1599570/1‑150 (MQ=255)
 cccGACCTTGGCAAGGTCGTGCTCTACCAACTGAGCTCTAGCCTCGTCGTGATGGGTGCCATTCTAGCGATTTCTAAAACAGCGTCAACCCCTTGATTCAAAAAACTTTTATTTTTTTTCCGAGCCTTCGCGCAGGTgcggccaggcggc                                                                                                                                                   <  1:1452839/150‑1 (MQ=255)
 cccGACCTTGGCAAGGTCGTGCGCTACCAACTGAGCTATTCCCGCGTCGTGATGGGTGCCATTCTAGCGATTTCTAAAACAGCGTCAACCCCTTGATTCAAAAAACTTTTATTTTTTTTCCGAGCCTTCGCGCAGGTgcggccaggcggc                                                                                                                                                   >  2:1115380/1‑150 (MQ=255)
 cacGACCTGGCCAAGGCCTTGCTCTACCACGTGATCTATTCCCGCGTCGCGATGGGTGCCATACTAGCGATTTCTAAAACAGCGTCAACCGCTTGATTCAACAAGCCTTTCTTTTTTTTCCGGGCCTTCTCGCCGGTgtggccagccggc                                                                                                                                                   >  2:307142/3‑150 (MQ=255)
        ttGGCCAGGTCGGTCTCTACCAACTGAGCTATTCCCGCGTCGTGATGGGTGCCATTCTAGCGATTTCTAAAACAGCGTCAACCCCTTCATTCAAAAAACTTTTATTTTTTTTCCGAGCCTTCGCGCAGGTGCGGCCAGGCGGCCAGCAgg                                                                                                                                            <  2:567188/150‑1 (MQ=255)
          ggCAAGGTCGTGCTCTACCAACTGAGCTATTCCCGCGTCGTGATGGGTGCCATTCTAGCGATTTCTAAAACAGCGTCAACCCCTTGATTCAAAAAACTTTTATTTTTTTTCCGAGCCTTCGCGCAGGTGCGGCCAGGCGGCCAGCAGGTa                                                                                                                                          >  1:1706499/1‑150 (MQ=255)
            cAAGGTCGTGCTCTACCAACTGAGCTATTCCCGCGTCGTGATGGGTGCCATTCTAGCGATTTCTAAAACAGCGTCAACCCCTTGATTCAAAAAACTTTTATTTTTTTTCCGAGCCTTCGCGCAGGTGCGGCCAGGCGGCCAGCAGGTAGt                                                                                                                                        >  1:31450/1‑150 (MQ=255)
            cAAGGTCGTGCTCTACCAACTGAGCTATTCCCGCGTCGTGATGGGTGCCATTCTAGCGATTTCTAAAACAGCGTCAACCCCTTGATTCAAAAAACTTTTATTTTTTTTCCGAGCCTTCGCGCAGGTGCGGCCAGGCGGCCAGCAGGTAGt                                                                                                                                        >  1:313725/1‑150 (MQ=255)
                 tCGTGCTCTACCAACTGAGCTATTCCCGGGTCGTGATGGGTGCCATTCTAGCGATTTCTAAAACAGCGTCAACCCCTTGATTCAAAAAACTTTTATTTTTTTTCCGAGCCTTCGCGCAGGTGCGGCCAGGCGGCCAGCAGGTAGTGCAcc                                                                                                                                   >  1:991702/1‑150 (MQ=255)
                   gTGCTCTACCAACTGAGCTATTCCCGCGTCGTGATGGGTGCCATTCTAGCGATTTCTAAAACAGCGTCAACCCCTTGATTCAAAAAACTTTTATTTTTTTTCCGATCCTTCGCGCAGGTGCGGCCAGGCGGCCAGCAGGTAGTGCACCAt                                                                                                                                 >  2:545179/1‑150 (MQ=255)
                    tGCTCTACCAACTGAGCTATTCCCGCGTCGTGATGGGTGCCATTCTAGCGATTTCTAAAACAGCGTCAACCCCTTGATTCAAAAAACTTTTATTTTTTTTCCGAGCCTTCGCGCAGGTGCGGCCAGGCGGCCAGCAGGTAGTGCACCATc                                                                                                                                <  2:555900/150‑1 (MQ=255)
                       tctACCAACTGAGCTATTCCCGCGTCGTGATGGGTGCCATTCTAGCGATTTCTAAAACAGCGTCAACCCCTTGATTCAAAAAACTTTTATTTTTTTTCCGAGCCTTCGCGCAGGTGCGGCCAGGCGGCCAGCAgg                                                                                                                                            <  1:775163/135‑1 (MQ=255)
                       tctACCAACTGAGCTATTCCCGCGTCGTGATGGGTGCCATTCTAGCGATTTCTAAAACAGCGTCAACCCCTTGATTCAAAAAACTTTTATTTTTTTTCCGAGCCTTCGCGCAGGTGCGGCCAGGCGGCCAGCAgg                                                                                                                                            >  2:775163/1‑135 (MQ=255)
                          aCCAACTGAGCTATTCCCGCGTCGTGATGGGTGCCATTCTAGCGATTTCTAAAACAGCGTCAACCCCTTGATTCAAAAAACTTTTATTTTTTTTCCGAGCCTTCGCGCAGGTGCGGCCAGGCGGCCAGCAGGTAGTGCACCATCGAc                                                                                                                             <  1:71756/147‑1 (MQ=255)
                          aCCAACTGAGCTATTCCCGCGTCGTGATGGGTGCCATTCTAGCGATTTCTAAAACAGCGTCAACCCCTTGATTCAAAAAACTTTTATTTTTTTTCCGAGCCTTCGCGCAGGTGCGGCCAGGCGGCCAGCAGGTAGTGCACCATCGAc                                                                                                                             >  2:71756/1‑147 (MQ=255)
                                tGAGCTATTCCCGCGTCGTGATGGGTGCCATTCTAGCGATTTCTAAAACAGCGTCAACCCCTTGATTCAAAAAACTTTTATTTTTTTTCCGAGCCTTCGCGCAGGTGCGGCCAGGCGGCCAGCAGGTAGTGCACCATCGACCACAGCGTc                                                                                                                    >  2:929208/1‑150 (MQ=255)
                                  aGCTATTCCCGCGTCGTGATGGGTGCCATTCTAGCGATTTCTAAAACAGCGTCAACCCCTTGATTCAAAAAACTTTTATTTTTTTTCCGAGCCTTCGCGCAGGTGCGGCCAGGCGGCCAGCAGGTAGTGCACCATCGACCACAGCGTCAg                                                                                                                  <  2:1150643/150‑1 (MQ=255)
                                     tATTCCCGCGTCGTGATGGGTGCCATTCTAGCGATTTCTAAAACAGCGTCAACCCCTTGATTCAAAAAACTTTTATTTTTTTTCCGAGCCTTCGCGCAGGTGCGTCCAGGCGGCCAGCAGGTAGTGCACCATCGACCACAGCGTCAGccc                                                                                                               >  1:545886/1‑150 (MQ=255)
                                      aTTCCCGCGTCGTGATGGGTGCCATTCTAGCGATTTCTAAAACAGCGTCAACCCCTTGATTCAAAAAACTTTTATTTTTTTTCCGAGCCTTCGCGCAGGTGCGGCCAGGCGGCCAGCAGGTAGTGCACCATCGACCACAGCGTCAGcccc                                                                                                              <  2:94103/150‑1 (MQ=255)
                                      aTTCCCGCGTCGTGATGGGTGCCATTCTAGCGATTTCTAAAACAGCGTCAACCCCTTGATTCAAAAAACTTTTATTTTTTTTCCGAGCCTTCGCGCAGGTGCGGCCAGGCGGCCAGCAGGTAGTGCACCATCGACCACAGCGTCAGcccc                                                                                                              <  1:545025/150‑1 (MQ=255)
                                       ttCCCGCGTCGTGATGGGTGCCATTCTAGCGATTTCTAAAACAGCGTCAACCCCTTGATTCAAAAAACTTTTATTTTTTTTCCGAGCCTTCGCGCAGGTGCGGCCAGGCGGCCAGCAGGTAGTGCACCATCGACCACAGCGTCAGCCCCg                                                                                                             >  2:813195/1‑150 (MQ=255)
                                        tCCCGCGTCGTGATGGGTGCCATTCTAGCGATTTCTAAAACAGCGTCAACCCCTTGATTCAAAAAACTTTTATTTTTTTTCCGAGCCTTCGCGCAGGTGCGGCCAGGCGGCCAGCAGGTAGTGCACCATCGACCACAGCGTCAGCCCCGc                                                                                                            <  2:1874743/150‑1 (MQ=255)
                                                  tGATGGGTGCCATTCTAGCGATTTCTAAAACAGCGTCAACCCCTTGATTCAAAAAACTTTTATTTTTTTTCCGAGCCTTCGCGCAGGTgcggccaggcggc                                                                                                                                                   <  2:34753/101‑1 (MQ=255)
                                                  tGATGGGTGCCATTCTAGCGATTTCTAAAACAGCGTCAACCCCTTGATTCAAAAAACTTTTATTTTTTTTCCGAGCCTTCGCGCAGGTgcggccaggcggc                                                                                                                                                   >  1:34753/1‑101 (MQ=255)
                                                  tGATGGGTGCCATTCTAGCGATTTCTAAAACAGCGTCAACCCCTTGATTCAAAAAACTTTTATTTTTTTTCCGAGCCTTCGCGCAGGTGCGGCCAGGCGGCCAGCAGGTAGTGCACCATCGACCACAGCGTCAGCCCCGCTGCCACcagc                                                                                                  >  2:1030880/1‑150 (MQ=255)
                                                  tGATGGGTGCCATTCTAGCGATTTCTAAAACAGCGTCAACCCCTTGATTCAAAAAACTTTTATTTTTTTTCCGAGCCTTCGCGCAGGTGCGGCCAGGCGGCCAGCAGGTAGTGCACCATCGACCACAGCGTCAGCCCCGCTGCCACcagc                                                                                                  >  2:1033244/1‑150 (MQ=255)
                                                         tGCCATTCTAGCGATTTCTAAAACAGCGTCAACCCCTTGATTCAAAAAACTTTTATTTTTTTTCCGAGCCTTCGCGCAGGTGCGGCCAGGCGGCCAGCAGGTAGTGCACCATCGACCACAGCGTCAGCCCCGCTGCCACCAGCAGCAGcc                                                                                           >  1:1893714/1‑150 (MQ=255)
                                                             aTTCTAGCGATTTCTAAAACAGCGTCAACCCCTTGATTCAAAAAACTTTTATTTTTTTTCCGAGCCTTCGCGCAGGTGCGGCCAGGCGGCCAGCAGGTAGTGCACCATCGACCAGAGCGTCAGCCCCGCTGCCACCTGCAGCAGCCCGTa                                                                                       <  1:1094829/150‑1 (MQ=255)
                                                               tCTAGCGATTTCTAAAACAGCGTCAACCACTTGCACAAAGAAACTTTTATTTTTTTCCGAGGCTTCGCGCAGGTGCGGCCAGGCGGCCAGCAGGTAGTGCACCATCGACCACAGCGTCAGCCCCGCTGCCACCAGCAGCAGCCCGTAGcc                                                                                    <  2:1401071/150‑1 (MQ=255)
                                                                 tAGCGATTTCTAAAACAGCGTCAACCCCTTGATTCAAAAAACTTTTATTTTTTTTCCGAGCCTTCGCGCAGGTGCGGCCAGGCGGCCAGCAGGTAGTGCACCATCGACCACAGCGTCAGCCCCGCTGCCACCAGCAGCAGCCCGTAGccc                                                                                   >  2:560434/1‑150 (MQ=255)
                                                                  aGCGATTTCTAAAACAGCGTCAACCCCTTGATTCACAAAACTTTTATTTTTTTTCCGAGCCTTCGCGCAGGTGCGGCCAGGCGGCCAGCAGGTAGTGCACCATCGACCACAGCGTCAGCCCCGCTGCCACCAGCAGCAGCCCGTAGCCCa                                                                                  <  2:1893714/150‑1 (MQ=255)
                                                                  aGCGATTTCTAAAACAGCGTCAACCCCTTGATTCAAAAAACTTTTATTTTTTTTCCGAGCCTTCGCGCAGGTGCGGCCAGGCGGCCAGCAGGTAGTGCACCATCGACCACAGCGTCAGCCCCGCTGCCACCAGCAGCAGCCCGTAGCCCa                                                                                  <  1:1878479/150‑1 (MQ=255)
                                                                      aTTTCTAAAACAGCGTCAACCCCTTGATTCAAAAAACTTTTATTTTTTTTCCGAGCCTTCGCGCAGGTGCGGCCAGGCGGCCAGCAGGTAGTGCACCATCGACCACAGCGTCAGCCCCGCTGCCACCAGCAGCAGCCCGTAGCCCAGGAt                                                                              >  1:1898568/1‑150 (MQ=255)
                                                                          cTAAAACAGCGTCAACCCCTTGATTCAAAAAACTTTTACTTTTTCCTGTAACCATCGCGCAGGTGCGGCCAGGCGGCCAGCAGGTAGTGCACCATCGACCACAGCGTCAGCCCCGCTGCCACCAGCAGCAGCCCGTAGCCCAGGATCAcc                                                                          >  1:331105/1‑150 (MQ=255)
                                                                           tAAAACAGCGTCAACCCCTTGATTCAAAAAACTTTTATTTTTTTTCCGAGCCTTCGCGCAGTTGCGGCCAGGCGGCCAGCAGGTAGTGCACCATCGACCACAGCGTCAGCCCCGCTGCCACCAGCAGCAGCCCGTAGCCCAGGATCAccc                                                                         >  2:695072/1‑150 (MQ=255)
                                                                             aaaCAGCGTCAACCCCTTGATTCAAAAAACTTTTATTTTTTTTCCGAGCCTTCGCGCAGGTGCGGCCAGGCGGCCAGCAGGTAGTGCACCATCGACCACAGCGTCAGCCCCGCTGCCACCCGCAGCAGCCCGTAGCCCAGGATCAccc                                                                         >  1:1368748/1‑148 (MQ=255)
                                                                             aaaCAGCGTCAACCCCTTGATTCAAAAAACTTTTATTTTTTTTCCGAGCCTTCGCGCAGGTGCGGCCAGGCGGCCAGCAGGTAGTGCACCATCGACCACAGCGTCAGCCCCGCTGCCACCAGCAGCAGCCCGTAGCCCAGGATCAccc                                                                         <  2:1368748/148‑1 (MQ=255)
                                                                               aCAGCGTCAACCCCTTGATTCAAAAAACTTTTATTTTTTTTCCGAGCCTTCGCGCAGGTGCGGCCAGGCGGCCAGCAGGTAGTGCa                                                                                                                                     >  1:297734/1‑86 (MQ=255)
                                                                               aCAGCGTCAACCCCTTGATTCAAAAAACTTTTATTTTTTTTCCGAGCCTTCGCGCAGGTGCGGCCAGGCGGCCAGCAGGTAGTGCa                                                                                                                                     <  2:297734/86‑1 (MQ=255)
                                                                                cAGCGTCAACCCCTTGATTCAAAAAACTTTTATTTTTTTTCCGAGCCTTCGCGCAGGTGCGGCCAGGCGGCCAGCAGGTAGTGCACCATCGACCACAGCGTCAGCCCCGCTGCCACCAGCAGCAGCCCGTAGCCCAGGATCACCCAGAAg                                                                    <  1:1318968/150‑1 (MQ=255)
                                                                                cAGCGTCAACCCCTTGATTCAAAAAACTTTTATTTTTTTTCCGAGCCTTCGCGCAGGTGCGGCCAGGCGGCCAGCAGGTAGTGCACCATCGACCACAGCGTCAGCCCCGCTGCCACCAGCAGCAGCCCGTAGCCCAGGATCACCCAGAAg                                                                    <  1:545179/150‑1 (MQ=255)
                                                                                       aaCCCCTTGATTCAAAAAACTTTTATTTTTTTTCCGAGCCTTCGCGCAGGTGCGGCCAGGCGGCCAGCAGGTAGTGCACCATCGACCACAGCGTCAGCCCCGCTGCCACCAGCAGCAGCCCGTAGCCCAGGATCACCCAGAAGGTcacca                                                             <  2:313725/150‑1 (MQ=255)
                                                                                       aaCCCCTTGATTCAAAAAACTTTTATTTTTTTTCCGAGCCTTCGCGCAGGTGCGGCCAGGCGGCCAGCAGGTAGTGCACCATCGACCACAGCGTCAGCCCCGCTGCCACCAGCAGCAGCCCGTAGCCCAGGATCACCCAGAAGGTcacca                                                             >  2:906884/1‑150 (MQ=255)
                                                                                       aaCCCCTTGATTCAAAAAACTTTTATTTTTTTTCCGAGCCTTCGCGCAGGTGCGGCCAGGCGGCCAGCAGGTAGTGCACCATCGACCACAGCGTCAGCCCCGCTGCCACCAGCAGCAGCCCGTAGCCCAGGATCACCCAGAAGGTcacca                                                             <  2:31450/150‑1 (MQ=255)
                                                                                         ccccTTGATTCAAAAAACTTTTATTTTTTTTCCGAGCCTTCGCGCAGGTGCGGCCAGGCGGCCAGCAGGTAGTGCACCATCGACCACAGCGTCAGCCCCGCTGCCACCAGCAGCAGCCCGTAGCCCAGGATCACCCAGAAGGTCaccacc                                                           >  1:1616560/1‑150 (MQ=255)
                                                                                          cccTTGATTCAAAAAACTTTTATTTTTTTTCCGAGCCTTCGCGCAGGTGCGGCCAGGCGGCCAGCAGGTAGTGCACCATCGACCACAGCGTCAGCCCCGCTGCCACCAGCAGCAGCCCGTAGCCCAGGATCACCCAGAAGGTCACCACcg                                                          >  2:1185196/1‑150 (MQ=255)
                                                                                            cTTGATTCAAAAAACTTTTATTTTTTTTCCGAGCCTTCGCGCAGGTGCGGCCAGGCGGCCAGCAGGTAGTCCACCATCGACCACAGCGTCAGCCCCGCTGCCACCAGCAGCAGCCCGTAGCCCAGGATCACCCAGAAGGTCACCACcggc                                                        >  1:1610236/1‑150 (MQ=255)
                                                                                                 ttCAAAAAACTTTTATTTTTTTTCCGAGCCTTCGCGCAGGTGCGGCCAGGCGGCCAGCAGGTAGTGCACCATCGAc                                                                                                                             >  1:1544680/1‑76 (MQ=255)
                                                                                                 ttCAAAAAACTTTTATTTTTTTTCCGAGCCTTCGCGCAGGTGCGGCCAGGCGGCCAGCAGGTAGTGCACCATCGAc                                                                                                                             <  2:1544680/76‑1 (MQ=255)
                                                                                                   cAAAAAACTTTTATTTTTTTTCCGAGCCTTCGCGCAGGTGCGGCCAGGCGGCCAGCAGGTAGTGCACCATCGACCACAGCGTCAGCCCCGCTGCCACCAGCAGCAGCCCGTAGCCCAGGATCACCCAGAAGGTCACCACCGGCGGGTTGc                                                 <  2:895293/150‑1 (MQ=255)
                                                                                                    aaaaaaCTTTTATTTTTTTTCCGAGCCTTCGCGCAGGTGCGTCCAGGCGGCCAGCAGGTAGTGCACCATCGACCACAGCGTCAGCCCCGCTGCCACCAGCAGCAGCCCGTAGCCCAGGATCACCCAGAAGGTCACCACCGGCGGGTTGcc                                                <  2:545886/150‑1 (MQ=255)
                                                                                                    aaaaaaCTTTTATTTTTTTTCCGAGCCTTCGCGCAGGTGCGGCCAGGCGGCCAGCAGGTAGTGCACCATCGACCACAGCGTCAGCCCCGCTGCCACCAGCAGCAGCCCGTAGCCCAGGATCACCCATAAGGTCACCACCGGCGGGTTGcc                                                >  2:713096/1‑150 (MQ=255)
                                                                                                    aaaaaaCTTTTATTTTTTTTCCGAGCCTTCGCGCAGGTGCGGCCAGGCGGCCAGCAGGTAGTGCACCATCGACCACAGCGTCAGCCCCGCTGCCACCAGCAGCAGCCCGTAGCCCAGGATCACCCAGAAGGTCACCACCGGCGGGTTGcc                                                <  1:929208/150‑1 (MQ=255)
                                                                                                     aaaaaCTTTTATTTTTTTTCCGAGCCTTCGCGCAGGTGCGGCCAGGTGGCCAGCAGGTAGTGCACCATCGACCACAGCGTCAGCCCCGCTGCCACCAGCAGCAGCCCGTAGCCCAGGATCACCCAGAAGGTCACCACcggcgg                                                      >  1:1144338/1‑143 (MQ=255)
                                                                                                     aaaaaCTTTTATTTTTTTTCCGAGCCTTCGCGCAGGTGCGGCCAGGTGGCCAGCAGGTAGTGCACCATCGACCACAGCGTCAGCCCCGCTGCCACCAGCAGCAGCCCGTAGCCCAGGATCACCCAGAAGGTCACCACcggcgg                                                      >  1:1146019/1‑143 (MQ=255)
                                                                                                     aaaaaCTTTTATTTTTTTTCCGAGCCTTCGCGCAGGTGCGGCCAGGTGGCCAGCAGGTAGTGCACCATCGACCACAGCGTCAGCCCCGCTGCCACCAGCAGCAGCCCGTAGCCCAGGATCACCCAGAAGGTCACCACcggcgg                                                      <  2:1144338/143‑1 (MQ=255)
                                                                                                     aaaaaCTTTTATTTTTTTTCCGAGCCTTCGCGCAGGTGCGGCCAGGTGGCCAGCAGGTAGTGCACCATCGACCACAGCGTCAGCCCCGCTGCCACCAGCAGCAGCCCGTAGCCCAGGATCACCCAGAAGGTCACCACcggcgg                                                      <  2:1146019/143‑1 (MQ=255)
                                                                                                      aaaaCTTTTATTTTTTTTCCGAGCCTTCGCGCAGGTGCGGCCAGGCGGCCAGCAGGTAGTGCACCATCGACCACAGCGTCAGCCCCGCTGCCACCAGCAGCAGCCCGTAGCCCAGGATCAc                                                                           <  1:1631569/121‑1 (MQ=255)
                                                                                                      aaaaCTTTTATTTTTTTTCCGAGCCTTCGCGCAGGTGCGGCCAGGCGGCCAGCAGGTAGTGCACCATCGACCACAGCGTCAGCCCCGCTGCCACCAGCAGCAGCCCGTAGCCCAGGATCAc                                                                           >  2:1631569/1‑121 (MQ=255)
                                                                                                       aaaCTTTTATTTTTTTTCCGAGCCTTCGCGCAGGTGCGGCCAGGCGGCCAGCAGGTAGTGCACCATCGACCACCGAGTCAGCCCCGCTGCCACCAGTAAAGGCTTCCATGCTATAAGCACCCAGAAGGTCACCACCGGCGGGTTGCCcag                                             <  1:1773814/150‑1 (MQ=255)
                                                                                                         aCTTTTATTTTTTTTCCGAGCCTTCGCGCAGGTGCGGCCAGGCGGCCAGCAGGTAGTGCACCATCGACCACAGCGTCAGCCCCGCTGCCACCAGCAGCAGCCCGTAGCCCAGGATCACCCAGAAGGTCACCACCGGCGGGTTGCCcagca                                           <  2:1706499/150‑1 (MQ=255)
                                                                                                         aCTTTTATTTTTTTTCCGAGCCTTCGCGCAGGTGCGGCCAGGCGGCCAGCAGGTAGTGCACCATCGACCACAGCGTCAGCCCCGCTGCCACCAGCAGCAGCCCGTAGCCCAGGATCACCCAGAAGGTCACCACCGGCGGGTTGCCcagca                                           <  1:813195/150‑1 (MQ=255)
                                                                                                           ttttATTTTTTTTCCGAGCCTTCGCGCAGGTGCGGCCAGGCGGCCAGCAGGTAGTGCACCATCGACCACAGCGTCAGCCCCGCTGCCACCAGCAGCAGCCCGTAGCCCAGGATCACCCAGAAGGTCACCACCGGCGGGTTGCCCAGCAgg                                         <  1:1293432/150‑1 (MQ=255)
                                                                                                            tttATTTTTTTTCCGAGCCTTCGCGCAGGTGCGGCCAGGCGGCCAGCAGGTAGTGCACCATCGACCACAGCGTCAGCCCCGCTGCCACCAGCAGCAGCCCGTAGCCCAGGATCACCCAGAAGGTCACCACCGGCGGGTTGCCCAGCAGGa                                        <  2:991702/150‑1 (MQ=255)
                                                                                                             ttATTTTTTTTCCGAGCCTTCGCGCAGGTGCGGCCAGGCGGCCAGCAGGTGTTGCACCATCGACCACAGCGTCAGCCCCGCTGCCACCAGCAGCAGCCCGTAGCCCAGGATCACCCAGAAGGTCACCACCGGCGGGTTGCCCAGCAGGAt                                       >  1:1787740/1‑150 (MQ=255)
                                                                                                               aTTTTTTTGCCGAGCCTTCGCGCAGGTGCGGCCAGGCGGCCAGCAGGTAGTGCACCATCGACCACAGCGTCAGCCCCGCTGCCACCAGCAGCAGCCCGTAGCCCAGGATCACCCAGAAGGTCACCACCGGCGGGTTGCCCAGCAGGATCa                                     <  2:786728/150‑1 (MQ=255)
                                                                                                                ttttttttCCGAGCCTTCGCGCAGGTGCGGCCAGGCGGCCAGCAGGTAGTGCACCATCGACCACAGCGTCAGCCCCGCTGCCACCAGCAGCAGCCCGTAGCCCAGGATCACCCAGAAGGTCACCACCGGCGGGTTGCCCAGCAGGATCAc                                    <  1:1030880/150‑1 (MQ=255)
                                                                                                                          gAGCCTTCGCGCAGGTGCGGCCAGGCGGCCAGCAGGTAGTGCACCATCGACCACAGCGTCAGCCCCGCTGCCACCAGCAGCAGCCCGTAGCCCAGGATCACCCAGAAGGTCACCACCGGCGGGTTGCCCAGCAGGATCACCAGCGCCAGc                          >  1:854020/1‑150 (MQ=255)
                                                                                                                           aGCCTTCGCGCAGGTGCGGCCAGGCGGCCAGCAGGTAGTGCACCATCGACCTCATCGTCAGCCCCGCTGCCACCAGCAGCAGCCCGTAGCCCAGGATCACCCAGAAGGTCACCACCGGCGGGTTGCCCAGCAGGATCACCAGCGCCAgta                         >  1:1079656/1‑148 (MQ=255)
                                                                                                                           aGCCTTCGCGCAGGTGCGGCCAGGCGGCCAGCAGGTAGTGCACCATCGACCACAGCGTCAGCCCCGCTGCCACCAGCAGCAGCCCGTAGCCCAGGATCACCCAGAAGGTCACCACCGGCGGGTTGCCCAGCAGGATCACCAGCGCCAGCa                         >  2:1105167/1‑150 (MQ=255)
                                                                                                                             ccTTCGCGCAGGTGCGGCCAAGCGGCCAGCAGGTAGTGCACCATCGACCACAGCGTCAGCCCCGCTGCCACCAGCAGCAGCCCATAGCCCAGGATCACCCAGAAGGTCACCACCGGCGGCTTGCCCAGCGGGCTCACCAGCGCCAGCATc                       >  2:417297/1‑150 (MQ=255)
                                                                                                                              cTTCGCGCAGGTGCGGCCAGGCGGCAAGATTGTAGTGCACCATCGACCACAGCGTCAGCCCCGCTGCCACCAGCAGCAGCCCGTAGCCCAGGATCACCCAGAAGGTCACCACCGGCGGGTTGCCCAGCAGGATCACCAGCGCCAGCATCt                      >  1:1297179/1‑150 (MQ=255)
                                                                                                                              cTTCGCGCAGGGGCGGCCAGGCGGCCAGCAGGTAGTGCACCATCGACCACAGCGTCAGCCCCGCTGCCACCAGCAGCAGCCCGCAGCCCAGGATCACCCAGAAGGTCACCACCGGCGGGTTGCCCAGCAGGATCACCAGCGCCAGCATCt                      <  1:401921/150‑1 (MQ=255)
                                                                                                                                  gcgcAGGTGCGGCCAGGCGGCCAGCAGGTAGTGCACCATCGACCACAGCGTCAGCCCCGCTGCCACCAGCAGCAGCCCGTAGCCCAGGATCACCCAGAAGGTCACCACCGGCGGGTTGCCCAGCAGGATCACCAGCGCCAGCATCTgcgc                  >  2:153251/1‑150 (MQ=255)
                                                                                                                                         tGCGGCCAGGCGGCCAGCAGGTAGTGCACCATCGACCACAGCGTCAGCCCCGCTGCCACCAGCAGCAGCCCGTAGCCCAGGATCACCCAGAAGGTCACCACCGGCGGGTTGCCCAGCAGGATCACCAGCGCCAGCATCTGCGCCGCAGTc           <  1:1061274/150‑1 (MQ=255)
                                                                                                                                            ggccaggcggccagCAGGTAGTGCACCATCGACCACAGCGTCAGCCCCGCTGCCACCAGCAGCAGCCCGTAGCCCAGGATCACCCAGAAGGTCACCACCGGCGGGTTGCCCAGCAGGACCACCAGCGCCAGCATCTGCGCCGCAGTCTTc        >  1:1759741/1‑150 (MQ=255)
                                                                                                                                                aggcggccagCAGGTAGTGCACCATCGACCACAGCGTCAGCCCCGCTGCCACCAGCAGCAGCCCGTAGCCCAGGATCACCCAGAAGGTCACCACCGGCGGGTTGCCCAGCAGGATCACCAGCGCCAGCATCTGCGCCGCAGTCTTCCACt    >  1:827220/1‑150 (MQ=255)
                                                                                                                                                  gcggccagCAGGTAGTGCACCATCGACCACAGCGTCAGCCCCGCTGCCACCAGCAGCAGCCCGTAGCCCAGGATCACCCAGAAGGTCACCACCGGCGGGTTGCCCAGCGGGATCACCAGCGCCAGCATCTGCGCCGCAGTCTTCCACTTa  <  1:1599570/150‑1 (MQ=255)
                                                                                                                                                  gcggccagCAGGTAGTGCACCATCGACCACAGCGTCAGCCCCGCTGCCACAAGCAGCAGCCCGTAGCCCAGGATCACCCAGAAGGTCACCACCGGCGGGTTGCCCAGCAGGATCACCAGCGCCAGCATCTGCGCCGCAGTCTTCCACTTa  <  2:1297179/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
CCCCGACCTTGGCAAGGTCGTGCTCTACCAACTGAGCTATTCCCGCGTCGTGATGGGTGCCATTCTAGCGATTTCTAAAACAGCGTCAACCCCTTGATTCAAAAAACTTTTATTTTTTTTCCGAGCCTTCGCGCAGGTGCGGCCAGGCGGCCAGCAGGTAGTGCACCATCGACCACAGCGTCAGCCCCGCTGCCACCAGCAGCAGCCCGTAGCCCAGGATCACCCAGAAGGTCACCACCGGCGGGTTGCCCAGCAGGATCACCAGCGCCAGCATCTGCGCCGCAGTCTTCCACTTA  >  NC_002947/4632630‑4632925

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