Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,725,784 C→G 9.3% E794Q (GAG→CAG)  plsB ← glycerol‑3‑phosphate acyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,725,7840CG 9.3% 221.2 / 11.5 75E794Q (GAG→CAG) plsBglycerol‑3‑phosphate acyltransferase
Reads supporting (aligned to +/- strand):  ref base C (36/31);  new base G (3/4);  total (40/35)
Fisher's exact test for biased strand distribution p-value = 7.01e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.60e-01

CACGGACTTTAGCGGATTTTGCCGACCGATGAACCGTTAGAGGGTCGAAGTTTCCAGCCCGTCGTCACGGTGCAGGGCCACCTGGCGGATCGACAGGCGCAGCTCGGCCGACAGTACTCGCTTGGCCACGCCCTCGGCCAGCTCGCCGAGCTTGTCGTGATAACTCAGCTTGCCTTGCGCGTCGGCGTGCAGCACGCCTTGCTGAAGCAAGGTCTGAATGAAGTGGCGGAACAGCGTCTTGTCGAAGAACTCCGGTGCGTTCAGGCCATGCAGGATCGACAGGCGCTGGG  >  NC_002947/1725641‑1725930
                                                                                                                                               |                                                                                                                                                   
cACGGACTTTAGCGGATTTTGCCGACCGATGAACCGTTAGAGGGTCGAAGTTTCCAGCCCGTCGTCACGGTGCAGGGCCACCTGGCGGATCGACAGGCGCAGCTCGGCCGACAGTACTCGCTTGGCCACGCCCTCGGCCAGCTCGcc                                                                                                                                                  >  2:31325/1‑147 (MQ=255)
 aCGGACTTTAGCGGATTTTGCCGACCGATGAACCGTTAGAGGGTCGAAGTTTCCAGCCCGTCGTCACGGTGCAGGGCCACCTGGCGGATCGACAGGCGCAGCTCGGCCGACAGTACTCGCTTGGCCACGCCCTCGGCGAGCTGGCCGAg                                                                                                                                               >  3:220362/1‑149 (MQ=255)
  cGGACTTTAGCGGATTTTGCCGACCGATGAACCGTTAGAGGGTCGAAGTTTCCAGCCCGTCTTCACGGTGCAGGGCCACCTGGCGGATCGACAGGCGCAGCTCGGCCGACAGTACTCGCTTGGCCACGCCCTCGGCGAGCTGGCCGAgg                                                                                                                                              <  7:247182/149‑2 (MQ=255)
       tttAGCGGATTTTGCCGACCGATGAACCGTTAGAGGGTCGAAGTTTCCAGCCCGTCGTCACGGTGCAGGGCCACCTGGCGGATCGACAGGCGCAGCTCGGCCGACAGTACTCGCTTGGCCACGCCCTCGGCCAGCTCGCCGAGCTTGTc                                                                                                                                         >  8:103061/1‑149 (MQ=255)
       tttAGCGGATTTTGCCGACCGATGAACCGTTAGAGGGTCGAAGTTTCCAGCCCGTCGTCACGGTGCAGGGCCACCTGGCGGATCGACAGGCGCAGCTCGGCCGACAGTACTCGCTTGGCCACGCCCTCGGCCAGCTCGCCGAGCTTGTc                                                                                                                                         >  8:108105/1‑149 (MQ=255)
         tAGCGGATTTTGCCGACCGATGAACCGTTAGAGGGTCGAAGTTTCCAGCCCGTCGTCACGGTGCAGGGCCACCTGGCGGATCGACAGGCGCAGCTCGGCCGACAGTACTCGCTTGGCCACGCCCTCGGCGAGCTGGCCGAgggcgtg                                                                                                                                         <  4:220362/147‑7 (MQ=255)
            cGGATTTTGCCGACCGATGAACCGTTAGAGGGTCGAAGTTTCCAGCCCGTCGTCACGGTGCAGGGCCACCTGGCGGATCGACAGGCGCAGCTCGGCCGACAGTACTCGCTTGGCCACGCCCTCGGCCAGCTCGCCGAGCTTGTCGTGa                                                                                                                                     >  5:195752/1‑148 (MQ=255)
              gATTTTGCCGACCGATGAACCGTTAGAGGGTCGAAGTTTCCAGCCCGTCGTCACGGTGCAGGGCCACCTGGCGGATCGACAGGCGCAGCTCGGCCGACAGTACTCGCTTGGCCACGCCCTCGGCCAGCTCGCCGAGCTTGTCGTGATaa                                                                                                                                  <  5:305788/149‑1 (MQ=255)
                ttttGCCGACCGATGAACCGTTAGAGGGTCGAAGTTTCCAGCCCGTCGTCACGGTGCAGGGCCACCTGGCGGATCGACAGGCGCAGCTCGGCCGACAGTACTCGCTTGGCCACGCCCTAGGCCAGCTCGCCGAGCTTGTCGTGATAACt                                                                                                                                <  1:31325/149‑1 (MQ=255)
                        accgaTGAACCGTTAGAGGGTCGAAGTTTCCAGCCCGTCGTCACGGTGCAGGGCCACCTGGCGGATCGACAGGCGCAGCTCGGCCGACAGTACTCGCTTGGCCACGCCCTCGGCCAGCTCGCCGAGCTTGTCGTGATAACTCAGcttgc                                                                                                                        >  2:262342/1‑149 (MQ=255)
                        accgaTGAACCGTTAGAGGGTCGAAGTTTCCAGCCCGTCGTCACGGTGCAGGGCCACCTGGCGGATCGACAGGCGCAGCTCGGCCGACAGTACTCGCTTGGCCACGCCCTCGGCCAGCTCGCCGAGCTTGGCGTGATAACTCAGcttgc                                                                                                                        >  4:106984/1‑149 (MQ=255)
                            aTGAACCGTTAGAGGGTCGAAGTTTCCAGCCCGTCGTCACGGTGCAGGGCCACCTGGCGGATCGACAGGCGCAGCTCGGCCGACAGTACTCGCTTGGCCACGCCCTCGGCCAGCTCGCCGAGCTTGTCGTGATAACTCAGcttgccttg                                                                                                                    >  6:46365/1‑149 (MQ=255)
                                aCCGTTAGAGGGTCGAAGTTTCCAGCCCGTCGTCACGGTGCAGGGCCACCTGGCGGATCGACAGGCGCAGCTCGGCCGACAGTACTCGCTTGGCCACGCCCTCGGCCAGCTCGCCGAGCTTGTCg                                                                                                                                        <  1:90604/125‑1 (MQ=255)
                                aCCGTTAGAGGGTCGAAGTTTCCAGCCCGTCGTCACGGTGCAGGGCCACCTGGCGGATCGACAGGCGCAGCTCGGCCGACAGTACTCGCTTGGCCACGCCCTCGGCCAGCTCGCCGAGCTTGTCg                                                                                                                                        >  2:90604/1‑125 (MQ=255)
                                aCCGTTAGAGGGTCGAAGTTTCCAGCCCGTCGTCACGGTGCAGGGCCACCTGGCGGATCGACAGGCGCAGCTCGGCCGACAGTACTCGCTTGGCCACGCCCTCGGCCAGCTCGCCGAGCTTGTCGTGATAACTCAGCTTGCCTTGcgcg                                                                                                                >  6:29501/1‑149 (MQ=255)
                                  cGTTAGAGGGTCGAAGTTTCCAGCCCGTCGTCACGGTGCAGGGCCACCTGGCGGATCGACAGGCGCAGCTCGGCCGACAGTACTCGCTTGGCCACGCCCTCGGCCAGCTCGCCGAGCTTGTCGTGATAACTCAGCTTGCCTTGCGCGTc                                                                                                              >  6:143040/1‑149 (MQ=255)
                                       gagGGTCGAAGTTTCCAGCCCGTCGTCACGGTGCAGGGCCACCTGGCGGATCGACAGGCGCAGCTCGGCCGACAGTACTCGCTTGGCCACGCCCTCGGCCAGCTCGCCGAGCTTGTCGTGATAACTCAGCTTGCCTTGCGCGTCGGCGt                                                                                                         >  1:91545/1‑149 (MQ=255)
                                        agGGTCGAAGTTTCCAGCCCGTCGTCACGGTGCAGGGCCACCTGGCGGATCGACAGGCGCAGCTCGGCCGACAGTACTCGCTTGGCCACGCCCTCGGCCAGCTCGCCGAGCTTGTCGTGATAACTCAGCTTGCCTTGCGCGTCGGCGTg                                                                                                        >  1:84821/1‑149 (MQ=255)
                                            tCGAAGTTTCCAGCCCGTCGTCACGGTGCAGGGCCACCTGGCGGATCGACAGGCGCAGCTCGGCCGACAGTACTCGCTTGGCCACGCCCTCGGCCAGCTCGCCGAGCTTGTCGTGATAACTCAGcttgccttgc                                                                                                                   >  5:246278/1‑134 (MQ=255)
                                            tCGAAGTTTCCAGCCCGTCGTCACGGTGCAGGGCCACCTGGCGGATCGACAGGCGCAGCTCGGCCGACAGTACTCGCTTGGCCACGCCCTCGGCCAGCTCGCCGAGCTTGTCGTGATAACTCAGcttgccttgc                                                                                                                   >  8:268215/1‑134 (MQ=255)
                                            tCGAAGTTTCCAGCCCGTCGTCACGGTGCAGGGCCACCTGGCGGATCGACAGGCGCAGCTCGGCCGACAGTACTCGCTTGGCCACGCCCTCGGCCAGCTCGCCGAGCTTGTCGTGATAACTCAGcttgccttgc                                                                                                                   <  7:268215/134‑1 (MQ=255)
                                            tCGAAGTTTCCAGCCCGTCGTCACGGTGCAGGGCCACCTGGCGGATCGACAGGCGCAGCTCGGCCGACAGTACTCGCTTGGCCACGCCCTCGGCCAGCTCGCCGAGCTTGTCGTGATAACTCAGcttgccttgc                                                                                                                   <  6:246278/134‑1 (MQ=255)
                                            tCGAAGTTTCCAGCCCGTCGTCACGGTGCAGGGCCACCTGGCGGATCGACAGGCGCAGCTCGGCCGACAGTACTCGCTTGGCCACGCCCTCGGCCAGCTCGCCGAGCTTGTCGTGATAACTCAGCTTGCCTTGCGCGTCGGCGTgcag                                                                                                     >  1:113552/1‑148 (MQ=255)
                                             cGAAGTTTCCAGCCCGTCGTCACGGTGCAGGGCCACCTGGCGGATCGACAGGCGCAGCTCCGCCGACAGTACTCGCTTGGCCACGCCCTCGGCCAGCTCGCCGAGCTTGTCGTGATAACTCAGCTTGCCTTGCGCGTCGGCGTgcagca                                                                                                   >  6:146045/1‑149 (MQ=255)
                                               aaGTTTCCAGCCCGTCGTCACGGTGCAGGGCCACCTGGCGGATCGACAGGCGCAGCTCGGCCGACAGTACTCGCTTGGCCACGCCCTCGGCCAGCTCGCCGAGCTTGTCGTGATAACTCAGCTTGCCTTGCGCGTCGGCGTGCAGCACg                                                                                                 >  5:149729/1‑149 (MQ=255)
                                                      cAGCCCGTCGTCACGGTGCAGGGCCACCTGGCGGATCGACAGGCGCAGCTCGGCCGACAGTACTCGCTTGGCCACGCCCTCGGCCAGCTCGCCGAGCTTGCCGTGATAACTCAGCTTGCCTTGCGCGTCGGCGTGCAGCACGCCTTGc                                                                                           <  3:106984/148‑1 (MQ=255)
                                                      cAGCCCGTCGTCACGGTGCAGGGCCACCTGGCGGATCGACAGGCGCAGCTCCGCCGACAGTAATCGCTTGGCCACGCCCTCGGCCAGCTAGGCGCTCTTGTGGTGATAACTCAGCTTGCCTTGCGCGTAGGCGTGAAGCACGCCTTGCt                                                                                          >  6:341162/1‑149 (MQ=255)
                                                         cccGTCGTCACGGTGCAGGGCCACCTGGCGGATCGACAGGCGCAGCTCGGCCGACAGTACTCGCTTGGCCACGCCCTCGGCCAGCTCGCCGAGCTTGTCGTGATAACTCAGCTTGCCTTGCGCGTCGGCGTGCAGCACGCCTTGCTGaa                                                                                       <  5:29501/149‑1 (MQ=255)
                                                         cccGTCGTCACGGTGCAGGGCCACCTGGCGGATCGACAGGCGCAGCTCGGCCGACAGTACTCGCTTGGCCACGCCCTCGGCCAGCTCGCCGAGCTTGTCGTGATAACTCAGCTTGCCTTGCGCGTCGGCGTGCAGCACGCCTTGCTGaa                                                                                       >  2:63532/1‑149 (MQ=255)
                                                           cgtcgtCACGGTGCAGGGCCACCTGGCGGATCGACAGGCGCAGCTCGGCCGACAGTACTCGCTTGGCCACGCCCTCGGCCAGCTCGCCGAGCTTGTCGTGATAACTCAGCTTGCCTTGCGCGTCGGCGTGCAGCACGCCTTGCTGAAg                                                                                      >  8:340139/1‑148 (MQ=255)
                                                               gtcACGGTGCAGGGCCACCTGGCGGATCGACAGGCGCAGCTCGGCCGACAGTACTCGCTTGGCCACGCCCTCGGCCAGCTCGCCGAGCTTGTCGTGATAACTCAGCTTGCCTTGCGCGTCGGCGTGCAGCACGCCTTGCTGAAGCAAgg                                                                                 <  7:108105/149‑1 (MQ=255)
                                                               gtcACGGTGCAGGGCCACCTGGCGGATCGACAGGCGCAGCTCGGCCGACAGTACTCGCTTGGCCACGCCCTCGGCCAGCTCGCCGAGCTTGTCGTGATAACTCAGCTTGCCTTGCGCGTCGGCGTGCAGCACGCCTTGCTGAAGCAAgg                                                                                 <  7:103061/149‑1 (MQ=255)
                                                               gtcACGGTGCAGGGCCACCTGGCGGATCGACAGGCGCAGCTCGGCCGACAGTACTCGCTTGGCCACGCCCTCGGCCAGCTCGCCGAGCTTGTCGTGATAACTCAGCTTGCCTTGCGCGTCGGCGTGCAGCACGCCTTGCTGAAGCAAg                                                                                  <  1:63532/148‑1 (MQ=255)
                                                                     gTGCAGGGCCACCTGGCGGATCGACAGGCGCAGCTCGGCCGACAGTACTCGCTTGGCCACGCCCTCGGCCAGCTCGCCGAGCTTGTCGTGATAACTCAGCTTGCCTTGCGCGTCGGCGTGCAGCACGCCTTGCTGAAGCAAGGTCtgaa                                                                           <  7:170565/149‑1 (MQ=255)
                                                                        cAGGGCCACCTGGCGGATCGACAGGCGCAGCTCGGCCGACAGTACTCGCTTGGCCACGCCCTCGGCCAGCTCGCCGAGCTTGTCGTGATAACTCAGCTTGCCTTGCGCGTCGGCGTGCAGCACGCCTTGCTGAAGCAAGGTCtgaatga                                                                        <  5:46365/149‑1 (MQ=255)
                                                                         aGGGCCACCTGGCGGATCGACAGGCGCAGCTCGGCCGACAGTACTCGCTTGGCCACGCCCTCGGCCAGCTCGCCGAGCTTGTCGTGATAACTCAGCTTGCCTTGCGCGTCGGCGTgca                                                                                                      >  7:297082/1‑118 (MQ=255)
                                                                         aGGGCCACCTGGCGGATCGACAGGCGCAGCTCGGCCGACAGTACTCGCTTGGCCACGCCCTCGGCCAGCTCGCCGAGCTTGTCGTGATAACTCAGCTTGCCTTGCGCGTCGGCGTgca                                                                                                      <  8:297082/118‑1 (MQ=255)
                                                                              cACCTGGCGGATCGACAGGCGCAGCTCGGCCGACAGTACTCGCTTGGCCACGCCCTCGGCCAGCTCGCCGAGCTTGTCGTGATAACTCAGCTTGCCTTGCGCGTCGGCGTGCAGCACGCCTTGCTGAAGCAAGGTCTGAATGAAGTGGc                                                                  <  2:84821/149‑1 (MQ=255)
                                                                                  tGGCGGATCGACAGGCGCAGCTCGGCCGACAGTACTCGCTTGGCCACGCCCTCGGCCAGCTCGCCGAGCTTGTCGTGATAACTCAGCTTGCCTTGCGCGTCGGCGTg                                                                                                        >  5:340043/1‑107 (MQ=255)
                                                                                  tGGCGGATCGACAGGCGCAGCTCGGCCGACAGTACTCGCTTGGCCACGCCCTCGGCCAGCTCGCCGAGCTTGTCGTGATAACTCAGCTTGCCTTGCGCGTCGGCGTg                                                                                                        <  6:340043/107‑1 (MQ=255)
                                                                                    gCGGATCGACAGGCGCAGCTCGGCCGACAGTACTCGCTTGGCCACGCCCTCGGCCAGCTCGCCGAGCTTGTCGTGATAACTCAGCTTGCCTTGCGCGTCGGCGTGCAGCACGCCTTGCTGAAGCAAGGTCTGAATGAAGTGGCGGAAc                                                             <  2:298060/148‑1 (MQ=255)
                                                                                    gCGGATCGACAGGCGCAGCTCGGCCGACAGTACTCGCTTGGCCACGCCCTCGGCCAACTCGCCGAGCTTGTCGTGATAACTCAGCTTGCCTTGCGCGTCGGCGTGCAGCACGCCTTGCTGAAGCAAGGTCTGAATGAAGTGGCGGaa                                                              >  8:270945/1‑147 (MQ=255)
                                                                                       gATCGACAGGCGCAGCTCGGCCGACAGTACTCGCTTGGCCACGCCCTCGGCCAGCTCGCCGAGCTTGTCGTGATAACTCAGCTTGCCTTGCGCGTCGGCGGGCCGCACGCCTTGCTGCAGCAAGGTCGGAATGAAGTGGCGGAACAGCg                                                         >  8:308750/1‑149 (MQ=255)
                                                                                          cGACAGGCGCAGCTCGGCCGACAGTACTCGCTTGGCCACGCCCTCGGCCAGCTCGCCGAGCTTGTCGTGATAACTCAGCTTGCCTTGCGCGTCGGCGTGCAGCACGCCTTGCTGAAGCAAGGTCTGAATGAAGTGGCGGAACAGCGTCt                                                      >  8:140788/1‑149 (MQ=255)
                                                                                             cAGGCGCAGCTCGGCCGACAGTACTCGCTTGGCCACGCCCTCGGCCAGCTCGCCGAGCTTGTCGTTATAACTCAGCTTGCCTTGCGCGTCGGCGTGCAGCACGCCTTGCTGAAGCAAGGTCTGAATGAAGTGGCGGAACAGCGTCTTGt                                                   <  6:195752/149‑1 (MQ=255)
                                                                                              aGGCGCAGCTCGGCCGACAGTACTCGCTTGGCCACGCCCTCGGCCAGCTCGCCGAGCTTGTCGTGATAACTCAGCTTGCCTTGCGCGTCGGCGTGCAGCACGCCTTGCTGAAGCAAGGTCt                                                                              <  5:205872/121‑1 (MQ=255)
                                                                                              aGGCGCAGCTCGGCCGACAGTACTCGCTTGGCCACGCCCTCGGCCAGCTCGCCGAGCTTGTCGTGATAACTCAGCTTGCCTTGCGCGTCGGCGTGCAGCACGCCTTGCTGAAGCAAGGTCt                                                                              >  6:205872/1‑121 (MQ=255)
                                                                                                gcgcAGCTCGGCCGACAGTACTCGCTTGGCCACGCCCTCGGCCAGCTCGCCGAGCTTGTCGTGATAACTCAGCTTGCCTTGCGCGTCGGTGTGCAGCACGCCTTGCTGAAGCAAGGTCTGAATGAAGTGGCGGAACAGCGTCTTGTCg                                                 >  6:61216/1‑148 (MQ=255)
                                                                                                gcgcAGCTCGGCCGACAGTACTCGCTTGGCCACGCCCTCGGCCAGCTCGCCGAGCTTGTCGTGATAACTCAGCTTGCCTTGCGCGTCGGCGTGCAGCACGCCTTGCTGAAGCAAGGTCTGAATGAAGTGGCGGAACAGCGTCTTGTCga                                                <  6:149729/149‑1 (MQ=255)
                                                                                                 cgcAGCTCGGCCGACAGTACTCGCTTGGCCACGCCCTCGGCCAGCTCGCCGAGCTTGTCGTGATAACTCAGCTTGCCTTGCGCGTCGGCGTGCAGCACGCCTTGCTGAAGCAAGGTCTGAATGAAGTGGCGGAACAGCGTCTTGTCgaa                                               <  2:91545/149‑1 (MQ=255)
                                                                                                  gcAGCTCGGCCGACAGTACTCGCTTGGCCACGCCCTCGGCCAGCTCGCCGAGCTTGTCGTGATAACTCAGCTTGCCTTGCGCGTCGGCGTGCAGCACGCCTTGCTGAAGCAAGGTCTGAATGAAGTGGCGGAACAGCGTCTTGTCgaag                                              <  7:140788/149‑1 (MQ=255)
                                                                                                         ggCCGACAGTACTCGCTTGGCCACGCCCTCGGCCAGCTCGCCGAGCTTGTCGTGATAACTCAGCTTGCCTTGCGCGTCGGCGTGCAGCACGCCTTGCTGAAGCAAGGTCTGAATGAAGTGGCGGAACAGCGTCTTGTCGAAGAACTcc                                        <  1:262342/148‑1 (MQ=255)
                                                                                                         ggCCGACAGTACTCGCTTGGCCACGCCCTCGGCCAGCTCGCCGAGCTTGTCGTGATAACTCAGCTTGCCTTGCGCGTCGGCGTGCAGCACGCCTTGCTGAAGCAAGGTCTGAATGAAGTGGCGGAACAGCGTCTTGTCGAAGAACTCCg                                       <  7:340139/149‑1 (MQ=255)
                                                                                                              aCAGTACTCGCTTGGCCACGCCCTCGGCCAACTCGCCGAGCTTGTCGTGATAACTCAGCTTGCCTTGCGCGTCGGCGTGCAGCACGCCTTGCTGAAGCAAGGTCTGAATGAAGTGGCGGAACAGCGTCTTGTCGAAGAACTCCGGTGCg                                 <  7:270945/149‑1 (MQ=255)
                                                                                                                    cTCGCTTGGCCACGCCCTCGGCCAGCTCGCCGAGCTTGTCGTGATAACTCAGCTTGCCTTGCGCGTCGGCGTGCAGCACGCCTTGCTGAAGCAAGGTCTGAATGAAGTGGCGGAACAGCGTCTTGTCGAAGAACTCCGGTGCGTTCAgg                           >  8:34686/1‑149 (MQ=255)
                                                                                                                    cTCGCTTGGCCACGCCCTCGGCCAGCTCGCCGAGCTTGTCGTGATAACTCAGCTTGCCTTGCGCGTCGGCGTGCAGCACGCCTTGCTGAAGCAAGGTCTGAATGAAGTGGCGGAACAGCGTCTTGTCGAAGAACTCCGGTGCGTTCAgg                           >  6:365190/1‑149 (MQ=255)
                                                                                                                    cTCGCTTGGCCACGCCCTCGGCCAGCTCGCCGAGCTTGTCGTGATAACTCAGCTTGCCTTGCGCGTCGGCGTGCAGCACGCCTTGCTGAAGCAAGGTCTGAATGAAGTGGCGGAACAGCGTCTTGTCGAAGAACTCCGGTGCGTTCAgg                           <  5:143040/149‑1 (MQ=255)
                                                                                                                    cTCGCTTGGCCACGCCCTCGGCCAGCTCGCCGAGCTTGTCGTGATAACTCAGCTTGCCTTGCGCGTCGGCGTGCAGCACGCCTTGCTGAAGCAAGGTCTGAATGAAGTGGCGGAACAGCGTCTTGTCGAAGAACTCCGGTGCGTTCAg                            <  5:146045/148‑1 (MQ=255)
                                                                                                                       gCTTGGCCACGCCCTCGGCCAGCTCGCCGAGCTTGTCGTGATAACTCAGCTTGCCTTGCGCGTCGGCGTGCAGCACGCCTTGCTGAAGCAAGGTCTGAATGAAGTGGCGGAACAGCGTCTTGTCgaa                                               <  1:66780/127‑1 (MQ=255)
                                                                                                                       gCTTGGCCACGCCCTCGGCCAGCTCGCCGAGCTTGTCGTGATAACTCAGCTTGCCTTGCGCGTCGGCGTGCAGCACGCCTTGCTGAAGCAAGGTCTGAATGAAGTGGCGGAACAGCGTCTTGTCgaa                                               >  2:66780/1‑127 (MQ=255)
                                                                                                                       gCTTGGCCACGCCCTCGGCCAGCTCGCCGAGCTTGTCGTGATAACTCAGCTTGCCTTGCGCGTCGGCGTGCAGCACGCCTTGCTGAAGCAAGGTCTGAATGAAGTGGCGGAACAGCGTCTTGTCGAAGAACTCCGGTGCGTTCAGGc                          <  7:34686/147‑1 (MQ=255)
                                                                                                                             ccACGCCCTCGGCCAGCTCGCCGAGCTTGTCGTGATAACTCAGCTTGCCTTGCGCGTCGGCGTGCAGCACGCCTTGCTGAAGCAAGGTCTGAATGAAGTGGCGGAACAGCGTCTTGTCgaagaa                                            >  4:122954/1‑124 (MQ=255)
                                                                                                                             ccACGCCCTCGGCCAGCTCGCCGAGCTTGTCGTGATAACTCAGCTTGCCTTGCGCGTCGGCGTGCAGCACGCCTTGCTGAAGCAAGGTCTGAATGAAGTGGCGGAACAGCGTCTTGTCgaagaa                                            <  3:122954/124‑1 (MQ=255)
                                                                                                                             ccACGCCCTCGGCCAGCTCGCCGAGCTTGTCGTGATAACTCAGCTTGCCTTGCGCGTCGGCGTGCAGCACGCCTTGCTGAAGCAAGGTCTGAATGAAGTGGCGGAACAGCGTCTTGTCGAAGAACTCCGGTGCGTTCAGGCCATGCAg                   >  3:242551/1‑148 (MQ=255)
                                                                                                                                caagcTCGGCGAGCTGGCCGAGCTTGTCGTGATAACTCAGCTTGCCTTGCGCGTCGGCGTGCAGCACGCCTTGCTGAAGCAAGGTCTGAATGAAGTGGCGGAACAGCGTCTTGTCGAAGAACTCCGGTGCGTTCAGGcc                         <  1:49086/135‑1 (MQ=255)
                                                                                                                                caagcTCGGCGAGCTGGCCGAGCTTGTCGTGATAACTCAGCTTGCCTTGCGCGTCGGCGTGCAGCACGCCTTGCTGAAGCAAGGTCTGAATGAAGTGGCGGAACAGCGTCTTGTCGAAGAACTCCGGTGCGTTCAGGcc                         >  2:49086/5‑139 (MQ=255)
                                                                                                                                caagcTCGGCGAGCTGGCCGAGCTTGTCGTGATAACTCAGCTTGCCTTGCGCGTCGGCGTGCAGCACGCCTTGCTGAAGCAAGGTCTGAATGAAGTGGCGGAACAGCGTCTTGTCGAAGAACTCCGGTGCGTTCAGGCCATGCa                    <  6:337651/140‑1 (MQ=255)
                                                                                                                                caagcTCGGCGAGCTGGCCGAGCTTGTCGTGATAACTCAGCTTGCCTTGCGCGTCGGCGTGCAGCACGCCTTGCTGAAGCAAGGTCTGAATGAAGTGGCGGAACAGCGTCTTGTCGAAGAACTCCGGTGCGTTCAGGCCATGCa                    >  5:337651/5‑144 (MQ=255)
                                                                                                                                cGCCCTCGGCCAGCTCGCCGAGCTTGTCGTGATAACTCAGCTTGCCTTGCGCGTCGGCGTGCAGCACGCCTTGCTGAAGCAAGGTCTGAATGAAGTGGCGGAACAGCGTCTTGTCGAAGAACTCCGGTGCGTTCAGGCCATGCAGGATc               >  4:39530/1‑149 (MQ=255)
                                                                                                                                    cTCGGCCAGCTCGCCGAGCTTGTCGTGATAACTCAGCTTGCCTTGCGCGTCGGCGTGCAGCACGCCTTGCTGAAGCAAGGTCTGAATGAAGt                                                                     <  2:51766/92‑1 (MQ=255)
                                                                                                                                    cTCGGCCAGCTCGCCGAGCTTGTCGTGATAACTCAGCTTGCCTTGCGCGTCGGCGTGCAGCACGCCTTGCTGAAGCAAGGTCTGAATGAAGt                                                                     >  1:51766/1‑92 (MQ=255)
                                                                                                                                    cTCGGCCAGCTCGCCGAGCTTGTCGTGATAACTCAGCTTGCCTTGCGCGTCGGCGTGCAGCACGCCTTGCTGAAGCAAGGTCTGAATGAAGTGGCGGAACAGCGTCTTGTCGAAGAACTCCGGTGCGTTCAGGCCATGCAGGATCGa             >  8:110335/1‑147 (MQ=255)
                                                                                                                                    cTCGGCCAGCTCGCCGAGCTTGTCGTGATAACTCAGCTTGCCTTGCGCGTCGGCGTGCAGCACGCCTTGCTGAAGCAAGGTCTGAATGAAGTGGCGGAACAGCGTCTTGTCGAAGAACTCCGGTGCGTTCAGGCCATGCAGGATCGa             >  3:41721/1‑147 (MQ=255)
                                                                                                                                           aGCTCGCCGAGCTTGTCGTGATAACTCAGCTTGCCTTGCGCGTCGGCGTGCAGCACGCCTTGCTGAAGCAAGGTCTGAATGAAGTGGCGGAACAGCGTCTTGTCGAAGAACTCCGGTGCGTTCAGGCCATGCAGGATCGACAGGCGCTg    >  2:234147/1‑149 (MQ=255)
                                                                                                                                              tCGCCGAGCTTGTCGTGATAACTCAGCTTGCCTTGCGCGTCGGCGTGCAGCACGCCTTGCTGAAGCAAGGTCTGAATGAAGTGGCGGAACAGCGTCTTGTCGAAGAACTCCGGCTGCGTTCAGGCCATGCAGGATCGACAGGCGCTggg  <  1:234147/149‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                                   
CACGGACTTTAGCGGATTTTGCCGACCGATGAACCGTTAGAGGGTCGAAGTTTCCAGCCCGTCGTCACGGTGCAGGGCCACCTGGCGGATCGACAGGCGCAGCTCGGCCGACAGTACTCGCTTGGCCACGCCCTCGGCCAGCTCGCCGAGCTTGTCGTGATAACTCAGCTTGCCTTGCGCGTCGGCGTGCAGCACGCCTTGCTGAAGCAAGGTCTGAATGAAGTGGCGGAACAGCGTCTTGTCGAAGAACTCCGGTGCGTTCAGGCCATGCAGGATCGACAGGCGCTGGG  >  NC_002947/1725641‑1725930

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: