Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 447,162:1 +C 8.8% intergenic (+39/‑263) PP_0368 → / → PP_0369 acyl‑CoA dehydrogenase/diguanylate cyclase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947447,1621.C 8.8% 273.6 / 19.2 68intergenic (+39/‑263)PP_0368/PP_0369acyl‑CoA dehydrogenase/diguanylate cyclase
Reads supporting (aligned to +/- strand):  ref base . (37/25);  new base C (2/4);  total (39/29)
Fisher's exact test for biased strand distribution p-value = 3.90e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.03e-01

TTCTACTTCCAGCGCATCCTGCCACGTACCCGTACCCATGTGGCGACCATGCTGTCGGGCGCCAACAACCTGATGGCGATGGACGAAGAGCACTTCGGCCTGTCGTACTGATACTGGCGCTGTACCCAAAACCCGCCTGCCCTGACCGGGAGGCGGGTTTTTTACGTCTTCAGGTTCAGCGTGGTCCCTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCTGCTGGCCAATCTCCATCAGCCGTTCACCGTTGCCAGCGCTACCGTATTCGCGGGTAAACCCGCTCCTACAA  >  NC_002947/447013‑447308
                                                                                                                                                      |                                                                                                                                                  
ttCTACTTCCAGCGCATCCTGCCACGTACCCGTACCCATGTGGCGACCATGCTGTCGGGCGCCAACAACCTGATGGCGATGGACGAAGAGCACTTCGGCCTGTCGTACTGATACTGGCGCTGTACCCAAAACCCGCCTGCCCTGACCggg                                                                                                                                                     >  2:502562/1‑150 (MQ=255)
   tACTTCCAGCGCATCCTGCCACGTACCCGTACCCATGTGGCGACCATGCTGTCGGGCGCCAACAACCTGATGGCGATGGACGAAGAGCACTTCGGCCTGTCGTACTGATACTGGCGCTGTACCCAAAACCCGCCTGCCCTGACCGGGAgg                                                                                                                                                 >  2:176199/1‑150 (MQ=255)
    aCTTCCAGCGCATCCTGCCACGTACCCGTACCCATGTGGCGACCATGCTGTCGGGCGCCAACAACCTGATGGCGATGGACGAAGAGCACTTCGGCCTGTCGTACTGATACTGGCGCTGTACCCAAAACCCGCCTGCCCTGACCGGGAGGc                                                                                                                                                <  2:542755/150‑1 (MQ=255)
      ttCCGGCGCATCCTGCAACGTACCCGTACCCACGTGGCGACCATGCTGTCGGGCCCCAACGACCTGATGGGGATGGACGAAGAGCACTTCGGCGTGTCGTACTGAGACTGGCGCTGTACCCAAAACCCGCGTGCCCTGACCGGGAGGcag                                                                                                                                              >  2:718803/1‑148 (MQ=255)
      ttCCAGCGCATCCTGCCACGTACCCGTACCCATGTGGCGACCATGCTGTCGGGCGCCAACAACCTGATGGCGATGGACGAAGAGCACTTCGGCCTGTCGTACTGATACTGGCGCTGTACCCAAAACCCGCCTGCCCTGACCGGGAGGCgg                                                                                                                                              <  1:176199/150‑1 (MQ=255)
         cAGCGCATCCTGCCACGTACCCGTACCCATGTGGCGACCATGCTGTCGGGCGCCAACAACCTGATGGCGATGGACGAAGAGCACTTCGGCCTGTCGTACTGATACTGGCGCTGTACCCAAAACCCGCCTGCCCTGACCGGGAGGCGGGtt                                                                                                                                           <  1:959925/150‑1 (MQ=255)
            cgcATCCTGCCACGTACCCGTACCCATGTGGCGACCATGCTGTCGGGCGCCAACAACCTGATGGCGATGGACGAAGAGCACTTCGGCCTGTCGTACTGATACTGGCGCTGTACCCAAAACCCGCCTGCCCTGACCGGGAGGCGGGttttt                                                                                                                                        <  1:741490/150‑1 (MQ=255)
            cgcAGCCGGCCACGTACCCGGACCCATGTGGCGACCAGGGTGTCGGGCGCCAACAACCTGAGGGCGATGGACGAAGAGCACTTCGGCCTGTCGTACTGATACTGGCGCTGTACCCAAAACCCGCCTCCCGGTCAGGGCAGGCGGGttttg                                                                                                                                        <  1:1193778/150‑2 (MQ=255)
             gcATCCTGCCACGTACCCGTACCCATGTGGCGACCATGCTGTCGGGCGCCAACAACCTGATGGCGATGGACGAAGAGCACTTCGGCCTGTCGTACTGATACTGGCGCTGTACCCAAAACCCGCCTCCCGGTCAGGGCAGGCGGGttttgg                                                                                                                                       <  1:763967/150‑3 (MQ=38)
                     ccACGTACCCGTACCCATGTGGCGACCATGCTGTCGGGCGCCAACAACCTGATGGCGATGGACGAAGAGCACTTCGGCCTGTCGTACTGATACTGGCGCTGTACCCAAAACCCGCCTGCCCTGACCGGGAGGCGGGTTTTTTACGTCTTc                                                                                                                               <  1:1195522/150‑1 (MQ=255)
                      cACGTACCCGTACCCATGTGGCGACCATGCTGTCGGGCGCCAACAACCTGATGGCGATGGACGAAGAGCACTTCGGCCTGTCGTACTGATACTGGCGCTGTACCCAAAACCCGCCTGCCCTGACCGGGAGGCgg                                                                                                                                              >  1:1013515/1‑134 (MQ=255)
                      cACGTACCCGTACCCATGTGGCGACCATGCTGTCGGGCGCCAACAACCTGATGGCGATGGACGAAGAGCACTTCGGCCTGTCGTACTGATACTGGCGCTGTACCCAAAACCCGCCTGCCCTGACCGGGAGGCgg                                                                                                                                              <  2:1013515/134‑1 (MQ=255)
                       aCGTACCCGTACCCATGTGGCGACCATGCTGTCGGGCGCCAACAACCTGATGGCGATGGACGAAGAGCACTTCGGCCTGTCGTACTGATACTGGCGCTGTACCCAAAACCCGCCTGCCCTGACCGGGAGGCGGGTTTTTTACGTCTTCAg                                                                                                                             >  2:618456/1‑150 (MQ=255)
                         gtacccgtacccATGTGGCGACCATGCTGTCGGGCGCCAACAACCTGATGGCGATGGACGAAGAGCACTTCGGCCTGTCGTACTGATACTGGCGCTGTACCCAAAACCCGCCTGCCCTGACCGGGAGGCGGGTTTTTTACGTCTTCAGGt                                                                                                                           >  1:83288/1‑150 (MQ=255)
                                 acccATGTGGCGACCATGCTGTCGGGCGCCAACAACCTGATGGCGATGGACGAAGAGCACTTCGGCCTGTCGTACTGATACTGGCGCTGTACCCAAAACCCGCCTGCCCTGACCGGGAGGCGGGTTTTTTACGTCTTCAGGTTCAGCGTg                                                                                                                   <  1:583058/150‑1 (MQ=255)
                                            gACCATGCTGTCGGGCGCCAACAACCTGATGGCGATGGACGAAGAGCACTTCGGCCTGTCGTACTGATACTGGCGCTGTACCCAAAACCCGCCTGCCCTGACCGGGAGGCGGGTTTTTTACGTCTTCAGGTTCAGCGTGGTCCCTCGTAg                                                                                                        >  2:851093/1‑150 (MQ=255)
                                            gACCATGCTGTCGGGCGCCAACAACCTGATGGCGATGGACGAAGAGCACTTCGGCATGTCGTACTGATACTGGCGCTGTACCCAAAACCCGCCTGCCCTGACCGGGAGGCGGGGTTTTTACGTTTTCAGGTTCAGCGTGGTCCCTCGTAg                                                                                                        >  2:1132228/1‑150 (MQ=255)
                                                    tGTCGGGCGCCAACAACCTGATGGCGATGGACGAAGAGCACTTCGGCCTGTTGTACTGATACTGGCGCTGTACCCAAAACCCGCCTGCCCTGACCGGGAGGCGGGTTTTTTACGTCTTCAGGTTCAGCGTGGTCCCTCGTAGGAGCGGGt                                                                                                >  2:849453/1‑150 (MQ=255)
                                                     gTCGGGCGCCAACAACCTGATGGCGATGGACGAAGAGCACTTCGGCCTGTCGTACTGATACTGGCGCTGTACCCAAAACCCGCCTGCCCTGACCGGGAGGCGGGTTTTTTACGTCTTCAGGTTCAGCGTGGTCCCTCGTAGGAGCGGGtt                                                                                               >  1:588665/1‑150 (MQ=255)
                                                         ggCGCCAACAACCTGATGGCGATGGACGAAGAGCACTTCGGCCTGTCGTACTGATACTGGCGCTGTACCCAAAACCCGCCTGCCCTGACCGGGAGGCGGGTTTTTTACGTCTTCAGGTTCAGCGTGGTCCCTCGTAGGAGCGGGTTTAcc                                                                                           <  1:851093/150‑1 (MQ=255)
                                                            gcCAACAACCTGATGGCGATGGACGAAGAGCACTTCGGCCTGTCGTACTGATACTGGCGCTGTACCCAAAACCCGCCTGCCCTGACCGGGAGGCGGGTTTTTTACGTCTTCAGGTTCAGCGTGGTCCCTCGTAGGAGCGGGTTTACCcgc                                                                                        <  2:425780/150‑1 (MQ=255)
                                                              caccaacCTGATGGCGCTGGACGAAGAGCACTTCGGCCTGCCGTACTGATACTGGCGCTGTACCCAAAACCCGCCTGCCCTGACCGGGAGGCGGGTTTTTTACGTCTTCAGGTTCATCGTGGTCCCTCGTAGGAGCGGGTTTACCcgcgg                                                                                      >  2:1156278/4‑149 (MQ=255)
                                                              caacaaCCTGATGGCGATGGACGAAGAGCACTTCGGCCTGTCGTACTGATACTGGCGCTGTACCCAAAACCCGCCTGCCCTGACCGGGAGGCGGGATTTTTACGTCTTCAGGTTCAGCGTGGTCCCTCGTAGGAGCGGGTTTACCCGCGa                                                                                      >  2:436579/1‑150 (MQ=255)
                                                              caacaaCCTGATGGCGATGGACGAAGAGCACTTCGGCCTGTCGTACTGATACTGGCGCTGTACCCAAAACCCGCCTCCCGGTCAGGGCAGGCGGGtttt                                                                                                                                         <  1:755948/99‑1 (MQ=25)
                                                              caacaaCCTGATGGCGATGGACGAAGAGCACTTCGGCCTGTCGTACTGATACTGGCGCTGTACCCAAAACCCGCCTCCCGGTCAGGGCAGGCGGGtttt                                                                                                                                         >  2:755948/1‑99 (MQ=25)
                                                                    ccTGATGGCGATGGACGAAGAGCACTTCGGCCTGTCGTACTGATACTGGCGCTGTACCCAAAACCCGCCTGCCCTGACCGGGAGGCGGGTTTTTTACGTCTTCAGGTTCAGCGTGGTCCCTCGTAg                                                                                                        <  2:914527/126‑1 (MQ=255)
                                                                    ccTGATGGCGATGGACGAAGAGCACTTCGGCCTGTCGTACTGATACTGGCGCTGTACCCAAAACCCGCCTGCCCTGACCGGGAGGCGGGTTTTTTACGTCTTCAGGTTCAGCGTGGTCCCTCGTAg                                                                                                        >  1:914527/1‑126 (MQ=255)
                                                                      tGATGGCGATGGACGAAGAGCACTTCGGCCTGTCGTACTGATACTGGCGCTGTACCCAAAACCCGCCTGCCCTGACCGGGAGGCGGGTTTTTTACGTCTTCAGGTTCAGCGTGGTCCCTCGTAGGAGCGGGTTTACCCGCGAAGAggccg                                                                              >  1:112174/1‑150 (MQ=255)
                                                                                 gACGAAGAGCACTTCGGCCTGTCGTACTGATACTGGCGCTGTACCCAAAACCCGCCTGCCCTGACCGGGAGGCGGGTTTTTTACGTCTTCAGGTTCAGCGTGGTCCCTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCTGCTGGc                                                                   >  2:1035330/1‑150 (MQ=255)
                                                                                   cGAAGAGCACTTCGGCCTGTCGTACTGATACTGGCGCTGTACCCAAAACCCGCCTGCCCTGACCGGGAGGCGGGTTTTTTACGTCTTCAGGTTCAGCGTGGTCCCTCGTAGGAGCGGGTTTACCCGCGAAGAggcc                                                                               <  1:251003/136‑1 (MQ=255)
                                                                                   cGAAGAGCACTTCGGCCTGTCGTACTGATACTGGCGCTGTACCCAAAACCCGCCTGCCCTGACCGGGAGGCGGGTTTTTTACGTCTTCAGGTTCAGCGTGGTCCCTCGTAGGAGCGGGTTTACCCGCGAAGAggcc                                                                               >  2:251003/1‑136 (MQ=255)
                                                                                           aCTTCGGCCTGTCGTACTGATACTGGCGCTGTACCCAAAACCCGCCTGCCCTGACCGGGAGGCGGGTTTTTTACGTCTTCAGGTTCAGCGTGGTCCCTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCTGCTGGCCAATCTCCAt                                                         >  2:817448/1‑150 (MQ=255)
                                                                                             ttCGGCCTGTCGTACTGATACTGGCGCTGTACCCAAAACCCGCCTGCCCTGACCGGGAGGCGGGTTTTTTAc                                                                                                                                     <  2:111853/72‑1 (MQ=255)
                                                                                             ttCGGCCTGTCGTACTGATACTGGCGCTGTACCCAAAACCCGCCTGCCCTGACCGGGAGGCGGGTTTTTTAc                                                                                                                                     >  1:111853/1‑72 (MQ=255)
                                                                                             ttCGGCCTGTCGTACTGATACTGGCGCTGTACCCAAAACCCGCCTGCCCTGACCGGGAGGCGGGTTTTTTACGTCTTCAGGTTCAGCGTGGTCCCTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCTGCTGGCCAATCTCCATCa                                                       >  1:59992/1‑150 (MQ=255)
                                                                                              tCGGCCTGTTGTACTGATACTGGCGCTGTACCCAAAACCCGCCTGCCCTGACCGGGAGGCGGGTTTTTTACGTCTTCAGGTTCAGCGTGGTCCCTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCTGCTGGCCAATCTCCATCAg                                                      <  1:849453/150‑1 (MQ=255)
                                                                                                ggCCTGTCGTACTGATACTGGCGCTGTACCCAAAACCCGCCTGCCCTGACCGGGAGGCGGGTTTTTTACGTCTTCAGGTTCAGCGTGGTCCCTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCTGCTGGCCAATCTCCATCAGcc                                                    <  2:83288/150‑1 (MQ=255)
                                                                                                  ccTGTCGTACTGATACTGGCGCTGTACCCAAAACCCGCCTGCCCTGACCGGGAGGCGGGTTTTTTACGTCTTCAGGTTCAGCGTGGTCCCTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCTGCTGGCCAATCTCCATCAGCCGt                                                  <  1:1035330/150‑1 (MQ=255)
                                                                                                  ccTGTCGTACTGATACTGGCGCTGTACCCAAAACCCGCCTCCCGGTCAGGGCAGGCGGGtttt                                                                                                                                         <  2:910216/63‑1 (MQ=9)
                                                                                                  ccTGTCGTACTGATACTGGCGCTGTACCCAAAACCCGCCTCCCGGTCAGGGCAGGCGGGtttt                                                                                                                                         >  1:910216/1‑63 (MQ=9)
                                                                                                   cTGTCGTACTGATACTGGCGCTGTACCCAAAACCCGCCTGCCCTGACCGGGAGGCGGGTTTTTTACGTCTTCAGGTTCAGCGTGGTCCCTc                                                                                                            >  1:963092/1‑91 (MQ=255)
                                                                                                   cTGTCGTACTGATACTGGCGCTGTACCCAAAACCCGCCTGCCCTGACCGGGAGGCGGGTTTTTTACGTCTTCAGGTTCAGCGTGGTCCCTc                                                                                                            <  2:963092/91‑1 (MQ=255)
                                                                                                    tGTCGTACTGATACTGGCGCTGTACCCAAAACCCGCCTGCCCTGACCGGGAGGCGGGTTTTTTACGTCTTCAGGTTCAGCGTGGTCCCTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCTGCTGGCCAATCTCCATCAGCCGTTc                                                >  2:993521/1‑150 (MQ=255)
                                                                                                    tGTCGTACTGATACTGGCGCTGTACCCAAAACCCGCCTGCCCTGACCGGGAGGCGGGTTTTTTACGTCTTCAGGTTCAGCGTGGTCCCTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCTGCTGGCCAATCTCCATCAGCCGTTc                                                >  2:1494/1‑150 (MQ=255)
                                                                                                     gTCGTACTGATACTGGCGCTGTACCCAAAACCCGCCTGCCCTGACCGGGAGGCGGGGTTTTTACGTCTTCAGGTTCAGCGTGGTCTCTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCTGCTGGCCAATCTCCATCAGCCGTTCa                                               >  2:617103/1‑150 (MQ=255)
                                                                                                              aTACTGGCGCTGTACCCAAAACCCGCCTGCCCTGACCGGGAGGCGGGTTTTTTACGTCTTCAGGTTCAGCGTGGTCCCTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCTGCTGGCCAATCTCCATCAGCCGTTCACCGTTGCCa                                      >  1:377717/1‑150 (MQ=255)
                                                                                                                          tACCCAAAACCCGCCTGCCCTGACCGGGAGGCGGGTTTTTTACGTCTTCAGGTTCAGCGTGGTCCCTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCTGCTGGCCAATCTCCATCAGCCGTTCACCGtt                                          <  1:577149/134‑1 (MQ=255)
                                                                                                                          tACCCAAAACCCGCCTGCCCTGACCGGGAGGCGGGTTTTTTACGTCTTCAGGTTCAGCGTGGTCCCTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCTGCTGGCCAATCTCCATCAGCCGTTCACCGtt                                          >  2:577149/1‑134 (MQ=255)
                                                                                                                            cccAAAACCCGCCTGCCCTGACCGGGAGGCGGGATTTTTACGTCTTCAGGTTCAGCGTGGTCCCTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGTCCCTGCTGTCCAATTTCCATCAGCCGTTCACCGTTGCCAGCGCTACCGTATTc                        >  2:135969/1‑150 (MQ=255)
                                                                                                                             ccAAAACCCGCCTGCCCTGACCGGGAGGCGGGTTTTTTACGTCTTCAGGTTCAGCGTGGTCCCTCGTAGGAGCGGGTTTACCCGCTAAGAGGCCGGCCCTGCTGGCCAATCTCCATCAGCCGTTCAc                                              >  2:1189325/1‑127 (MQ=255)
                                                                                                                             ccAAAACCCGCCTGCCCTGACCGGGAGGCGGGTTTTTTACGTCTTCAGGTTCAGCGTGGTCCCTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCTGCTGGCCAATCTCCATCAGCCGtt                                                 <  1:893644/124‑1 (MQ=255)
                                                                                                                             ccAAAACCCGCCTGCCCTGACCGGGAGGCGGGTTTTTTACGTCTTCAGGTTCAGCGTGGTCCCTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCTGCTGGCCAATCTCCATCAGCCGtt                                                 >  2:893644/1‑124 (MQ=255)
                                                                                                                             ccAAAACCCGCCTGCCCTGACCGGGAGGCGGGTTTTTTACGTCTTCAGGTTCAGCGTGGTCCCTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCTGCTGGCCAATCTCCATCAGCCGTTCACCGTTGCCAGCGCTACCGTATTcg                       >  2:348420/1‑150 (MQ=255)
                                                                                                                             ccAAAACCCGCCTGCCCTGACCGGGAGGCGGGTTTTTTACGTCTTCAGGTTCAGCGTGGTCCCTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCTGCTGGCCAAGCTCCATCAGCCGTTCAc                                              <  1:1189325/127‑1 (MQ=255)
                                                                                                                              cAAAACCCGCCTGCCCTGACCGGGAGGCGGGTTTTTTACGTCTTCAGGTTCAGCGTGGTCCCTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCTGCTGGCCAATCTCCATCAGCCGTTCACCGTTGCCAGCGCTACCGTATTcgc                      >  2:493353/1‑150 (MQ=255)
                                                                                                                              cAAAACCCGCCTGCCCTGACCGGGAGGCGGGTTTTTTACGTCTTCAGGTTCAGCGTGGTCCCTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCTGCTGGCCAATCTCCATCAGCCGTTCACCGTTGCCAGCGCTACCGTATTcgc                      >  1:93732/1‑150 (MQ=255)
                                                                                                                               aagaCCCGCCTGCCCTGACCGGGAGGCGGGTTTTTTACGTCTTCAGGTTCAGCGTGGTCCCTCGTAGGAGCGGGTTTACCCGCGAAGAggc                                                                                >  1:1205813/4‑91 (MQ=255)
                                                                                                                               aaaaCCCGCCTGCCCTGACCGGGAGGCGGGTTTTTTACGTCTTCAGGTTCAGCGTGGTCCCTCGTAGGAGCGGGTTTACCCGCGAAGAggc                                                                                <  2:1205813/91‑1 (MQ=255)
                                                                                                                               aaaaCCCGCCTGCCCTGACCGGGAGGCGGGTTTTTTACGTCTTCAGGTTCAGCGTGGTCCCTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCTGCTGGCCAATCTCCATCAGCCGTTCACCGTTGCCAGCGCTACCGTATTcgcg                     >  1:862060/1‑150 (MQ=255)
                                                                                                                                aaaCCCGCCTGCCCTGACCGGGAGGCGGGTTTTTTACGTCTTCAGGTTCAGCGTGGTCCCTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCTGCTGGCCAATCTCCATCAGCCGTTCACCGTTGCCAGCGCTACCGTATTCGCgg                    <  2:59992/150‑1 (MQ=255)
                                                                                                                                aaaCCCGCCTGCCCTGACCGGGAGGCGGGTTTTTTACGTCTTCAGGTTCAGCGTGGTCCCTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCTGCTGGCCAATCTCCATCAGCCGTTCACCGTTGCCAGCGCTACCGTATTCGCgg                    <  1:817448/150‑1 (MQ=255)
                                                                                                                                 aaCCCGCCTGCCCTGACCGGGAGGCGGGTTTTTTACGTCTTCAGGTTCAGCGTGGTCCCTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCTGCTGGCCAATCTCCATCAGCCGTTCACCGTTGCCAGCGCTACCGTATTCGCggg                   <  1:993521/150‑1 (MQ=255)
                                                                                                                                   cccGCCTGCCCTGACCGGGAGGCGGGTTTTTTACGTCTTCAGGTTCAGCGTGGTCCCTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCTGCTGGCCAATCTCCATCAGc                                                     <  1:398649/114‑1 (MQ=255)
                                                                                                                                   cccGCCTGCCCTGACCGGGAGGCGGGTTTTTTACGTCTTCAGGTTCAGCGTGGTCCCTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCTGCTGGCCAATCTCCATCAGc                                                     >  2:398649/1‑114 (MQ=255)
                                                                                                                                    ccGCCTGCCCTGACCGGGAGGCGGGTTTTTTACGTCTTCAGGTTCAGCGTGGTCCCTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCTGCTGGCCAATCTCCATCAGCCGTTCGCCGTTGCCAGCGCTACCGTATTCGCGGGTaa                >  1:809986/1‑150 (MQ=255)
                                                                                                                                      gCCTGCCCTGACCGGGAGGCGGGTTTTTTACGTCTTCAGGTTCAGCGTGGTCCCTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCTGCTGGCCAATCTCCATCAGCCGTTCACCGTTGCCAGCGCTACCGTATTCGCGGGTAAAc              <  1:618456/150‑1 (MQ=255)
                                                                                                                                         tGCCCTGACCGGGAGGCGGGTTTTTTACGTCTTCAGGTTCAGCGTGGTCCCTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCTGCTGGCCAATCTCCATCAGCCGTTCACCGTTGCCAGCGCTACCGTATTCGCGGGTAAACCCg           >  1:26654/1‑150 (MQ=255)
                                                                                                                                           cccTGACCGGGAGGCGGGTTTTTTACGTCTTCAGGTTCAGCGTGGTCCCTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCTGCTGGCCAATCTCCATCAGCCGTTCACCGTTGCCAGCGCTACCGTATTCGCGGGTAAACCCGCt         >  2:439676/1‑150 (MQ=255)
                                                                                                                                            ccTGACCGGGAGGCGGGTTTTTTACGTCTTCAGGTTCAGCGTGGTCCCTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCTGCTGGCCAATCTCCATCAGCCGTTCACCGTTGCCAGCGCTACCGTATTCGCGGGTAAACCAGCTc        >  2:5402/1‑150 (MQ=255)
                                                                                                                                                  ggggAGGCGGGTTTTTTACGTCTACAGGTACACCGTGGTCCCTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCTGCTGGCCAATCTACATCAGCCGTTCACCGTTGCCAGCGCTACCGTATTCGCGGGTAAACCCGCTCCTACaa  <  2:862060/149‑1 (MQ=255)
                                                                                                                                                     gAGGCGGGTTTTTTACGTCTTCAGGTTCAGCGTGGTCCCTCGTAGGAGCGGGTTTACCCGCGaa                                                                                     >  1:499126/1‑64 (MQ=255)
                                                                                                                                                     gAGGCGGGTTTTTTACGTCTTCAGGTTCAGCGTGGTCCCTCGTAGGAGCGGGTTTACCCGCGaa                                                                                     <  2:499126/64‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                  
TTCTACTTCCAGCGCATCCTGCCACGTACCCGTACCCATGTGGCGACCATGCTGTCGGGCGCCAACAACCTGATGGCGATGGACGAAGAGCACTTCGGCCTGTCGTACTGATACTGGCGCTGTACCCAAAACCCGCCTGCCCTGACCGGGAGGCGGGTTTTTTACGTCTTCAGGTTCAGCGTGGTCCCTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCTGCTGGCCAATCTCCATCAGCCGTTCACCGTTGCCAGCGCTACCGTATTCGCGGGTAAACCCGCTCCTACAA  >  NC_002947/447013‑447308

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