Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 447,157 A→T 8.6% intergenic (+34/‑268) 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,1570AT 8.6% 205.9 / 9.4 70intergenic (+34/‑268)PP_0368/PP_0369acyl‑CoA dehydrogenase/diguanylate cyclase
Reads supporting (aligned to +/- strand):  ref base A (39/25);  new base T (2/4);  total (41/29)
Fisher's exact test for biased strand distribution p-value = 2.24e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.33e-01

GCTTCTACTTCCAGCGCATCCTGCCACGTACCCGTACCCATGTGGCGACCATGCTGTCGGGCGCCAACAACCTGATGGCGATGGACGAAGAGCACTTCGGCCTGTCGTACTGATACTGGCGCTGTACCCAAAACCCGCCTGCCCTGACCGGGAGGCGGGTTTTTTACGTCTTCAGGTTCAGCGTGGTCCCTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCTGCTGGCCAATCTCCATCAGCCGTTCACCGTTGCCAGCGCTACCGTATTCGCGGGTAAACCCGCTC  >  NC_002947/447011‑447302
                                                                                                                                                  |                                                                                                                                                  
gcTTCTACTTCCAGCGCATCCTGCCACGTACCCGTACCCATGTGGCGACCATGCTGTCGGGCGCCAACAACCTGATGGCGATGGACGAAGAGCACTTCGGCCTGTCGTACTGATACTGGCGCTGTACCCAAAACCCGCCTGCCCTGACCg                                                                                                                                                 >  2:583058/1‑150 (MQ=255)
  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                                                                                                  >  1:914527/1‑126 (MQ=255)
                                                                      ccTGATGGCGATGGACGAAGAGCACTTCGGCCTGTCGTACTGATACTGGCGCTGTACCCAAAACCCGCCTGCCCTGACCGGGAGGCGGGTTTTTTACGTCTTCAGGTTCAGCGTGGTCCCTCGTAg                                                                                                  <  2:914527/126‑1 (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                                                                                                      <  2:963092/91‑1 (MQ=255)
                                                                                                     cTGTCGTACTGATACTGGCGCTGTACCCAAAACCCGCCTGCCCTGACCGGGAGGCGGGTTTTTTACGTCTTCAGGTTCAGCGTGGTCCCTc                                                                                                      >  1:963092/1‑91 (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                                    >  2:577149/1‑134 (MQ=255)
                                                                                                                            tACCCAAAACCCGCCTGCCCTGACCGGGAGGCGGGTTTTTTACGTCTTCAGGTTCAGCGTGGTCCCTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCTGCTGGCCAATCTCCATCAGCCGTTCACCGtt                                    <  1:577149/134‑1 (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                                               >  2:398649/1‑114 (MQ=255)
                                                                                                                                     cccGCCTGCCCTGACCGGGAGGCGGGTTTTTTACGTCTTCAGGTTCAGCGTGGTCCCTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCTGCTGGCCAATCTCCATCAGc                                               <  1:398649/114‑1 (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)
                                                                                                                                                  |                                                                                                                                                  
GCTTCTACTTCCAGCGCATCCTGCCACGTACCCGTACCCATGTGGCGACCATGCTGTCGGGCGCCAACAACCTGATGGCGATGGACGAAGAGCACTTCGGCCTGTCGTACTGATACTGGCGCTGTACCCAAAACCCGCCTGCCCTGACCGGGAGGCGGGTTTTTTACGTCTTCAGGTTCAGCGTGGTCCCTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCTGCTGGCCAATCTCCATCAGCCGTTCACCGTTGCCAGCGCTACCGTATTCGCGGGTAAACCCGCTC  >  NC_002947/447011‑447302

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: