Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,522,792 T→A 10.0% intergenic (+26/+14) PP_4856 → / ← yhjG Dps family ferritin/protein YhjG

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,522,7920TA10.0% 221.1 / 12.7 80intergenic (+26/+14)PP_4856/yhjGDps family ferritin/protein YhjG
Reads supporting (aligned to +/- strand):  ref base T (41/31);  new base A (2/6);  total (43/37)
Fisher's exact test for biased strand distribution p-value = 1.35e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.53e-01

CCGCGAGATCATTCAATACATCGGGGATAAAGACCCGACTACCCGTCGTATCTTCGAAGACATCCTGGCGCAGGAAGAGGAGCACGCGGACGATATGTCCGACCTGTTGCAAGGCCTTTGATTTCTGACACCCTCAAGGGGCCGCTTTTGCGGCCCCTTGCTATTTCCGTGCTTTCACCGCCGTCGGCGCCTTGCCTGCCTTCATCTGCTGCAACAACGGTGTGCACTGGTTCGGGTCATCCCCGCTGCTGGGCGCAATCAGCGCCAGCAACCCCGCCGCCGGCCCTGCCAC  >  NC_002947/5522646‑5522937
                                                                                                                                                  |                                                                                                                                                 
ccGCGAGATCATTCAATACATCGGGGATAAAGACCCGACTACCCGTCGTATCTTCGAAGACATCCTGGCGCAGGAAGAGGAGCACGCGGACGATATGTCCGACCTGTTGCAAGGCCTTTGATTTCTGACACCCTCAAGGGGCCGCtt                                                                                                                                                   <  7:239324/147‑1 (MQ=255)
   cgAGATCATTCAATACATCGGGGATAAAGACCCGACTACCCGTCGTATCTTCGAAGACATCCTGGCGCAGGAAGAGGAGCACGCGGACGATATGTCCGACCTGTTGCAAGGCCTTTGATTTCTGACACCCTCAAGGGGCCGCTTTTGCg                                                                                                                                              >  8:22720/1‑149 (MQ=255)
      gaTCATTCAATACATCGGGGATAAAGACCCGACTACCCGTCGTATCTTCGAAGACATCCTGGCGCAGGAAGAGGAGCACGCGGACGATATGTCCGACCTGTTGCAAGGCCTTTGATTTCTGACACCCTCAAGGGGCCGCTTTTGCGGcc                                                                                                                                           >  3:85111/1‑149 (MQ=255)
         cATTCAATACATCGGGGATAAAGACCCGACTACCCGTCGTATCTTCGAAGACATCCTGGCGCAGGAAGAGGAGCACGCGGACGATATGTCCGACCTGTTGCAAGGCCTTTGATTTCTGACACCCTCAAGGGGCCGCTTTTGCGGCCCCt                                                                                                                                        >  3:363540/1‑149 (MQ=255)
          aTTCAATACATCGGGGATAAAGACCCGACTACCCGTCGTATCTTCGAAGACATCCTGGCGCAGGAAGAGGAGCACGCGGACGATATGTCCGACCTGTTGCAAGGCCTTTGATTTCTGACACCCTCAAGGGGCCGCtttt                                                                                                                                                 >  4:332505/1‑139 (MQ=255)
           ttCAATACATCGGGGATAAAGACCCGACTACCCGTCGTATCTTCGAAGACATCCTGGCGCAGGAAGAGGAGCACGCGGACGATATGTCCGACCTGTTGCAAGGCCTTTGATTTCTGACACCCTCAAGGGGCCGCtttt                                                                                                                                                 >  5:325774/1‑138 (MQ=255)
           ttCAATACATCGGGGATAAAGACCCGACTACCCGTCGTATCTTCGAAGACATCCTGGCGCAGGAAGAGCAGCACGCGGACGATATGTCCGACCTGTTGCAAGGCCTTTGATTTCTGACACCCTCAAGGGGCCGCtttt                                                                                                                                                 <  6:325774/138‑1 (MQ=255)
             cAATACATCGGGGATAAAGACCCGACTACCCGTCGTATCTTCGAAGACATCCTGGCGCAGGAAGAGGAGCACGCGGACGATATGTCCGACCTGTTGGAAGGCCTTTGATTTCTGACACCCTCAAGGGGCCGCAAAAGCGGCCCCTTga                                                                                                                                     <  1:359859/148‑2 (MQ=255)
             cAATACATCGGGGATAAAGACCCGACTACCCGTCGTATCTTCGAAGACATCCTGGCGCAGGAAGAGGAGCACGCGGACGATATGTCCGACCTGTTGCAAGGCCTTTGATTTCTGACACCCTCAAGGGGCCGCAAAAGCGGCCCCTTga                                                                                                                                     <  8:339843/148‑2 (MQ=255)
             cAATACATCGGGGATAAAGACCCGACTACCCGTCGTATCTTCGAAGACATCCTGGCGCAGGAAGAGGAGCACGCGGACGATATGTCCGACCTGTTGCAAGGCCTTTGATTTCTGACACCCTCAAGGGGCCGCAAAAGCGGCCCCTTga                                                                                                                                     <  4:189065/148‑2 (MQ=255)
              aaTACATCGGGGATAAAGACCCGACTACCCGTCGTATCTTCGAAGACATCCTGGCGCAGGAAGAGGAGCACGCGGACGATATGTCCGACCTGTTGCAAGGCCTTTGATTTCTGACACCCTCAAGGGGCCGCTTTTGCGGCCCCTTGCTa                                                                                                                                   >  6:232211/1‑149 (MQ=255)
              aaTACATCGGGGATAAAGACCCGACTACCCGTCGTATCTTCGAAGACATCCTGGCGCAGGAAGAGGAGCACGCGGACGATATGTCCGACCTGTTGCAAGGCCTTTGATTTCTGACACCCTCAAGGGGCCGCAAAAGCGGCCCCTTgagg                                                                                                                                   <  5:70125/149‑4 (MQ=255)
              aaTACATCGGGGATAAAGACCCGACTACCCGTCGTATCTTCGAAGACATCCTGGCGCAGGAAGAGGAGCACGCGGACGATATGTCCGACCTGTTGCAAGGCCTTTGATTTCTGACACCCTCAAGGGGCCGCAAAAGCGGCCCCTTgagg                                                                                                                                   <  8:25030/149‑4 (MQ=255)
              aaTACATCGGGGATAAAGACCCGACTACCCGTCGTATCTTCGAAGACATCCTGGCGCAGGAAGAGGAGCACGCGGACGATATGTCCGACCTGTTGCAAGGCCTTTGATTTCGGACACCCTCAAGGTGCCGCTTTTGCGGCCCGTTGCTa                                                                                                                                   >  8:130133/1‑149 (MQ=255)
                   aTCGGGGATAAAGACCCGACTACCCGTCGTATCTTCGAAGACATCCTGGCGCAGGAAGAGGAGCACGCGGACGATATGTCCGACCTGTTGCAAGGCCTTTGATTTCTGACACCCTCAAGGGGCCGCTTTTGCGGCCCCTTGCTATTTcc                                                                                                                              <  3:339518/149‑1 (MQ=255)
                    tCGGGGATAAAGACCCGACTACCCGTCGTATCTTCGAAGCCATCCTGGCGCAGGAAGAGGAGCACGCGGACGATATGTCCGACCTGTTGCAAGGCCTTTGATTTCTGACACCCTCAAGGGGCCGCTTTTGCGGCCCCTTGCTATTTCCg                                                                                                                             >  3:89043/1‑149 (MQ=255)
                    tCGGGGATAAAGACCCGACTACCCGTCGTATCTTCGAAGACATCCTGGCGCAGGAAGAGGAGCACGCGGACGATATGTCCGACCTGTTGCAAGGCCTTTGATTTCTGACACCCTCAAGGGGCCGCTTTTGCGGCCCCTTGCTATTTCCg                                                                                                                             >  1:133898/1‑149 (MQ=255)
                    tCGGGGATAAAGACCCGACTACCCGTCGTATCTTCGAAGACATCCTGGCGCAGGAAGAGGAGCACGCGGACGATATGTCCGACCTGTTGCAAGGCCTTTGATTTCTGACACCCTCAAGGGGCCGCAAAAGCGGCCCCTTgag                                                                                                                                    >  6:275885/1‑140 (MQ=255)
                    tCGGGGATAAAGACCCGACTACCCGTCGTATCTTCGAAGACATCCTGGCGCAGGAAGAGGAGCACGCGGACGATATGTCCGACCTGTTGCAAGGCCTTTGATTTCTGACACCCTCAAGGGGCCGCAAAAGCGGCCCCTTgag                                                                                                                                    <  5:275885/142‑3 (MQ=255)
                        ggATAAAGACCCGACTACCCGTCGTATCTTCGAAGACATCCTGGCGCAGGAAGAGGAGCACGCGGACGATATGTCCGACCTGTTGCAAGGCCTTTGATTTCTGACACCCTCAAGGGGCCGCTTTTGCGGCCCCTTGCTATTTCCGTGCt                                                                                                                         >  6:322684/1‑149 (MQ=255)
                        ggATAAAGACCCGACTACCCGTCGTATCTTCGAAGACATCCTGGCGCAGGAAGAGGAGCACGCGGACGATATGTCCGACCTGTTGCAAGGCCTTTGATTTCTGACACCCTCAAGGGGCCGCTTTTGCGGCCCCTTGCTATTTCCGTGCt                                                                                                                         >  6:126971/1‑149 (MQ=255)
                          aTAAAGACCCGACTACCCGTCGTATCTTCGAAGACATCCTGGCGCAGGAAGAGGAGCACGCGGACGATATGTCCGACCTGTTGCAAGGCCTTTGATTTCTGACACCCTCAAGGGGCCGCTTTTGCGGCCCCTTGCTATTTCCGTGCttt                                                                                                                       <  2:133898/149‑1 (MQ=255)
                            aaaGACCCGACTACCCGTCGTATCTTCGAAGACATCCTGGCGCAGGAAGAGGAGCACGCGGACGATATGTCCGACCTGTTGCAAGGCCTTTGATTTCTGACACCCTCAAGGGGCCGCTTTTGCGGCCCCTTGCTATTTCCGTGCTTTCa                                                                                                                     >  2:331017/1‑149 (MQ=255)
                            aaaGACCCGACTACCCGTCGTATCTTCGAAGACATCCTGGCGCAGGAAGAGGAGCACGCGGACGATATGTCCGACCTGTTGCAAGGCCTTTGATTTCTGACACCCTCAAGGGGCCGCTTTTGCGGCCCCTTGCTATTTCCGTGCTTTCa                                                                                                                     >  8:241303/1‑149 (MQ=255)
                             aaGACCCGACTACCCGTCGTATCTTCGAAGACATCCTGGCGCAGGAAGAGGAGCACGCGGACGATATGTCCGACCTGTTGCAAGGCCTTTGATTTCTGACACCCTCAAGGGGCCGCTTTTGCGGCCCCTTGCTATTTCCGTGCTTTCAc                                                                                                                    <  7:22720/149‑1 (MQ=255)
                                   cGACTACCCGTCGTATCTTCGAAGACATCCTGGCGCAGGAAGAGGAGCACGCGGACGATATGTCCGACCTGTTGCAAGGCCTTTGATTTCTGACACCCTCAAGGGGCCGCTTTTGCGGCCCCTTGCTATTTCCGTGCTTTCACCGCCGt                                                                                                              >  7:342974/1‑149 (MQ=255)
                                         cccGTCGTATCTTCGAAGACATCCTGGCGCAGGAAGAGGAGCACGCGGACGATATGTCCGACCTGTTGCAAGGCCTTTGATTTCTGACACCCTCAAGGGGCCGCTTTTGCGGCCCCTTGCTATTTCCGTGCTTTCACCGCCGTCGgcgc                                                                                                        >  6:41779/1‑149 (MQ=255)
                                         cccGTCGTATCTTCGAAGACATCCTGGCGCAGGAAGAGGAGCACGCGGACGATATGTCCGACCTGTTGCAAGGCCTTTGATTTCTGACACCCTCAAGGGGCCGCTTTTGCGGCCCCTTGCTATTTCCGTGCTTTCACCGCCGTCGgcgc                                                                                                        >  1:132107/1‑149 (MQ=255)
                                          ccGTCGTATCTTCGAAGACATCCTGGCGCAGGAAGAGGAGCACGCGGACGATATGTCCGACCTGTTGCAAGGCCTTTGATTTCTGACACCCTCAAGGGGCCGCTTTTGCGGCCCCTTGCTATTTCCGTGCTTTCACCGCCGTCGGCGcc                                                                                                       >  7:226823/1‑149 (MQ=255)
                                           cgtcgtATCTTCGAAGACATCCTGGCGCAGGAAGAGGAGCACGCGGACGATATGTCCGACCTGTTGCAAGGCCTTTGATTTCTGACACCCTCAAGGGGCCGCTTTTGCGGCCCCTTGCTATTTCCGTGCt                                                                                                                         <  8:330410/130‑1 (MQ=255)
                                           cgtcgtATCTTCGAAGACATCCTGGCGCAGGAAGAGGAGCACGCGGACGATATGTCCGACCTGTTGCAAGGCCTTTGATTTCTGACACCCTCAAGGGGCCGCTTTTGCGGCCCCTTGCTATTTCCGTGCt                                                                                                                         >  7:330410/1‑130 (MQ=255)
                                                tATCTTCGAAGACATCCTGGCGCAGGAAGAGGAGCACGCGGACGATATGTCCGACCTGTTGCAAGGCCTTTGATTTCTGACACCCTCAAGGGGCCGCTTTTGCGGCCCCTTGCTATTTCCGTGCTTTCACCGCCGTCGGCGCCTtgcct                                                                                                 >  4:157080/1‑149 (MQ=255)
                                                         aGACATCCTGGCGCAGGAAGAGGAGCACGCGGACGATATGTCCGACCTGTTGCAAGGCCTTTGATTTCTGACACCCTCAAGGGGCCGCTTTTGCGGCCCCTTGCTATTTCCGTGCTTTCACCGCCGTCGGCGCCTTGCCTGCCTTCATc                                                                                        <  1:331017/149‑1 (MQ=255)
                                                         aGACATCCTGGCGCAGGAAGAGGAGCACGCGGACGATATGTCCGACCTGTTGCAAGGCCTTTGATTTCTGACACCCTCAAGGGGCCGCTTTTGCGGCCCCTTGCTATTTCCGTGCTTTCACCGCCGTCGGCGCCTTGCCTGCCTTCATc                                                                                        >  5:181733/1‑149 (MQ=255)
                                                            cATCCTGGCGCAGGAAGAGGAGCACGCGGACGATATGTCCGACCTGTTGCAAGGCCTTTGATTTCTGACACCCTCAAGGGGCCGCTTTTGCGGCCCCTTGCTATTTCCGTGCTTTCACCGCCGTCGGCGCCTTGCCTGCCTTCATctgc                                                                                     <  4:265664/149‑1 (MQ=255)
                                                              tCCTGGCGCAGGAAGAGGAGCACGCGGACGATATGTCCGACCTGTTGCAAGGCCTTTGATTTCTGACACCCTCAAGGGGCCGCTTTTGCGGCCCCTTGCTATTTCCGTGCTTTCACCGCCGTCGGCGCCTTGCCTGCCTTCATctgctg                                                                                   <  3:157080/149‑1 (MQ=255)
                                                              tCCTGGCGCAGGAAGAGGAGCACGCGGACGATATGTCCGACCTGTTGCAAGGCCTTTGATTTCTGACACCCTCAAGGGGCCGCTTTTGCGGCCCCTTGCTATTTCCGTGCTTTCACCGCCGTCGGCGCCTTGCCTGCCTTCATctgctg                                                                                   <  6:181733/149‑1 (MQ=255)
                                                               ccTGGCGCAGGAAGAGGAGCACGCGGACGATATGTCCGACCTGTTGCAAGGCCTTTGATTTCTGACACCCTCAAGGGGCCGCTTTTGCGGCCCCt                                                                                                                                        <  6:44266/95‑1 (MQ=255)
                                                               ccTGGCGCAGGAAGAGGAGCACGCGGACGATATGTCCGACCTGTTGCAAGGCCTTTGATTTCTGACACCCTCAAGGGGCCGCTTTTGCGGCCCCt                                                                                                                                        >  5:44266/1‑95 (MQ=255)
                                                                cTGGCGCAGGAAGAGGAGCACGCGGACGATATGTCCGACCTGTTGCAAGGCCTTTGATTTCTGACACCCTCAAGGGGCCGCTTTTGCGGCCCCTTGCTATTTCCGTGCTTTCACCGCCGTCGGCGCCTTGCCTGCCTTCATCtgctgc                                                                                  <  4:85111/148‑1 (MQ=255)
                                                                cTGGCGCAGGAAGAGGAGCACGCGGACGATATGTCCGACCTGTTGCAAGGCCTTTGATTTCTGACACCCTCAAGGGGCCGCTTTTGCGGCCCCTTGCTATTTCCGTGCTTTCACCGCCGTCGGCGCCTGGCCTGCATTCATCtgctgc                                                                                  >  8:188449/1‑148 (MQ=255)
                                                                   gcgcAGGAAGAGGAGCACGCGGACGATATGTCCGACCTGTTGCAAGGCCTTTGATTTCTGACACCCTCAAGGGGCCGCTTTTGCGGCCCCTTGCTATTTCCGTGCTTTCACCGCCGTCGGCGCCTTGCCTGCCTTCATCTGCTGcaac                                                                               <  8:39702/148‑1 (MQ=255)
                                                                       aGGAAGAGGAGCACGCGGACGATATGTCCGACCTGTTGCAAGGCCTTTGATTTCTGACACCCTCAAGGGGCCGCTTTTGCGGCCCCTTGCTATTTCCGTGCTTTCACCGCCGTCGGCGCCTTGCCTGCCTTCATCTGCTGCAACAACgg                                                                          >  2:226894/1‑149 (MQ=255)
                                                                       aGGAAGAGGAGCACGCGGACGATATGTCCGACCTGTTGCAAGGCCTTTGATTTCTGACACCCTCAAGGGGCCGCTTTTGCGGCCCCTTGCTATTTCCGTGCTTTCACCGCCGTCGGCGCCTTGCCTGCCTTCATCTGCTGCAACAACgg                                                                          >  7:191829/1‑149 (MQ=255)
                                                                        ggAAGAGGAGCACGCGGACGATATGTCCGACCTGTTGCAAGGCCTTTGATTTCTGACACCCTCAAGGGGCCGCTTTTGCGGCCCCTTGCTATTTCCGTGCTTTCACCGCCGTCGGCGCCTTGCCTGCCTTCATCTGCTGCAACAACGgt                                                                         <  8:226823/149‑1 (MQ=255)
                                                                         gAAGAGGAGCACGCGGACGATATGTCCGACCTGTTGCAAGGCCTTTGATTTCTGACACCCTCAAGGGGCCGCTTTTGCGGCCCCTTGCTATTTCCGTGCTTTCACCGCCGTCGGCGCCTTGCCTGCCTTCATCTGCTGCAACAACGgtg                                                                        >  5:312055/1‑149 (MQ=255)
                                                                                agCACGCGGACGATATGTCCGACCTGTTGCAAGGCCTTTGATTTCTGACACCCTCAAGGGGCCGCTTTTGCGGCCCCTTGCTATTTCCGTGCTTTCACCGCCGTCGGCGCCTTGCCTGCCTTCATCTGCTGCAACAACGGTGTGCACTg                                                                 >  3:58001/1‑149 (MQ=255)
                                                                                 gCACGCGGACGATATGTCCGACCTGTTGCAAGGCCTTTGATTTCTGACACCCTCAAGGGGCCGCTTTTGCGGCCCCTTGCTATTTCCGTGCTTTCACCGCCGTCGGCGCCTTGCCTGCCTTCATCTGCTGCAACAACGGTGTGCACTgg                                                                >  5:163134/1‑149 (MQ=255)
                                                                                    cgcgGACGATATGTCCGACCTGTTGCAAGGCCTTTGATTTCTGACACCCTCAAGGGGCCGCTTTTGCGGCCCCTTGCTATTTCCGTGCTTTCACCGCCGTCGGCGCCTTGCCTGCCTTCATCTGCTGCAACAACGGTGTGCACTGGTTc                                                             <  5:232211/149‑1 (MQ=255)
                                                                                     gcgGACGATATGTCCGACCTGTTGCAAGGCCTTTGATTTCTGACACCCTCAAGGGGCCGCTTTTGCGGCCCCTTGCTATTTCCGTGCTTTCACCGCCGTCGGCGCCTTGCCTGCCTTCATCTGCTGCAACAACGGTGTGCACTGGTTCg                                                            >  5:403440/1‑149 (MQ=255)
                                                                                       ggACGATATGTCCGACCTGTTGCAAGGACTTTGATTTCTGACACCCTCAAGGGGCCGCTTTTGCGGCCCCTTGCTATTTCCGTGCTTTCACCGCCGTCGGCGCCTTGCCTGCCTTCATCTGCTGCAACAACGGTGTGCACTGGTTCggg                                                          <  7:241303/149‑1 (MQ=255)
                                                                                              tcttacGATCTGTTGCAAGGCCTTTGATTTCTGACACCCTCAAGGGGCCGCTTTTGCGGCCCCTTGCTATTTCCGTGCTTTCACCTCCGTCGGCGCCTTGCCTGCCTTCATCTGCTGCAACAACGGTGTGCACTGGTTCGGGTCATccc                                                   <  7:339532/144‑1 (MQ=255)
                                                                                               tGTCCGACCTGTTGCAAGGCCTTTGATTTCTGACACCCTCAAGGGGCCGCTTTTGCGGCCCCTTGCTATTTCCGTGCTTTCACCGCCGTCGGCGCCTTGCCTGCCTTCATCTGCTGCAACAACGGTGTGCACTGGTTCGGGTCATcccc                                                  <  8:191829/149‑1 (MQ=255)
                                                                                                              aaGGCCTTTGATTTCTGACACCCTCAAGGGGCCGCTTTTGCGGCCCCTTGCTATTTCCGTGCTTTCAcc                                                                                                                   <  2:130328/69‑1 (MQ=255)
                                                                                                               aGGCCTTTGATTTCTGACACCCTCAAGGGGCCGCTTTTGCGGCCCCTTGCTATTTCCGTGCTTTCACCGCCGTCGGCGCCTTGCCTGCCTTCATCTGCTGCAACAACGgt                                                                         <  2:311783/110‑1 (MQ=255)
                                                                                                               aGGCCTTTGATTTCTGACACCCTCAAGGGGCCGCTTTTGCGGCCCCTTGCTATTTCCGTGCTTTCACCGCCGTCGGCGCCTTGCCTGCCTTCATCTGCTGCAACAACGGTGTGCACTGGTACGGGTCATCCCCGCTGCTGGGCGCAATc                                  <  1:226894/149‑1 (MQ=255)
                                                                                                                 gCCTTTGATTTCTGACACCCTCAAGGGGCCGCTTTTGCGGCCCCTTGCTATTTCCGTGCTTTCACCGCCGTCGGCGCCTTGCCTGCCTTCATCTGCTGCAACAACGGTGTGCACTGGTTCGGGTCATCCCCGCTGCTGGGCGCAATCAg                                >  6:141866/1‑149 (MQ=255)
                                                                                                                  ccTTTGATTTCTGACACCCTCAAGGGGCCGCTTTTGCGGCCCCTTGCTATTTCCGTGCTTTAACCGCCGTCGGCGCCTTGCCTTCCTTCATCTGCTGCAACAACGGTGTGCACTGATTCTGGACATCACCGCTGCTGGGAGCAATCAgc                               >  8:103959/1‑149 (MQ=255)
                                                                                                                   cTTTGATTTATGACACCCTCAAGGGGCCGCTTTTGCGGCCCCTTGCTATTTCCGTGCATTAACCGCCGTCGGAGCCTTGCCTGCCATCCTCTGCTGCAACAAAGGTGTGCACTGGTTCGGGTAATCCCACCTGCTGGTCGCAATCAgcg                              >  1:367144/1‑149 (MQ=255)
                                                                                                                      tGATTTCTGACACCCTCAAGGGGCCGCTTTTGCGGCCCCTTGCTATTTCCGTGCTTTCACCGCCGTCGGCGCCTTGCCTGCCTTCATCTGCTGCAACAACGGTGTGCACTGGTTCGGGTCATCCCCGCTGCTGGGCGCAATCAgcgc                             <  6:312055/147‑1 (MQ=255)
                                                                                                                      tGATTTCTGACACCCTCAAGGGGCCGCTTTTGCGGCCCCTTGCTATTTCCGTGCTTTCACCGCCGTCGGCGCCTTGCCTGCCTTCATCTGCTGCAACAACGGTGTGCACTGGTTCGGGTCATCCCCGCTGCTGGGCGCAATCAgcgc                             <  5:126971/147‑1 (MQ=255)
                                                                                                                      tGATTTCTGACACCCTCAAGGGGCCGCTTTTGCGGCCCCTTGCTATTTCCGTGCTTTCACCGCCGTCGGCGCCTTGCCTGCCTTCATCTGCTGCAACAACGGTGTGCACTGGTTCGGGTCATCCCCGCTGCTGGGCGCAATCAGCGCCa                           >  8:307048/1‑149 (MQ=255)
                                                                                                                       gATTTCTGACACCCTCAAGGGTCCGCTTTTGCGGCCCCTTGCTATTTCCGTGCTTTCACCGCCGTCGGCGCCTTGCCTGCCTTCATCTGCTGCAACAACGGTGTGCACTGGTTCGGGTCATCCCCGCTGCTGGGCGCAATCAGCGCCa                           <  2:132107/148‑1 (MQ=255)
                                                                                                                       cagttCTGACACCCTCACGGGGCCGCTTTTGCCGCACCTTGCTATTTCCGTGCCTTCACCGCCGTCGGCGCCTTGCCTGCCTTCATCTGCTGCAACAACCGTGTGCACTGTTTCGGCTCATCCCCGCTGCTGGGCGCAATCAGCGCCa                           <  8:69363/145‑1 (MQ=255)
                                                                                                                        aTTTCTGACACCCTCAAGGGGCCGCTTTTGCGGCCCCTTGCTATTTCCGTGCTTTCACCGCCGTCgg                                                                                                           >  6:84304/1‑67 (MQ=255)
                                                                                                                        aTTTCTGACACCCTCAAGGGGCCGCTTTTGCGGCCCCTTGCTATTTCCGTGCTTTCACCGCCGTCgg                                                                                                           <  5:84304/67‑1 (MQ=255)
                                                                                                                         tttCTGACACCCTCAAGGGGCCGCTTTTGCGGCCCCTTGCTATTTCCGTGCTTTCACCGCCGTCGGCGCCTTGCCTGCCTTCATCTGCTGCAACAACGGTGTGCACTGGTTCGGGTCATCCCCGCTGCTGGGCGCAATCAgcgc                             >  5:211698/1‑144 (MQ=255)
                                                                                                                         tttCTGACACCCTCAAGGGGCCGCTTTTGCGGCCCCTTGCTATTTCCGTGCTTTCACCGCCGTCGGCGCCTTGCCTGCCTTCATCTGCTGCAACAACGGTGTGCACTGGTTCGGGTCATCCCCGCTGCTGGGCGCAATCAgcgc                             <  6:211698/144‑1 (MQ=255)
                                                                                                                           tCTGACACCCTCAAGGGGCCGCTTTTGCGGCCCCTTGCTATTTCCGTGCTTTCACCGCCGTCGGCGCCTTGCCTGCCTTCATCTGCTGCAACAACGGTGTGCACTGGTTCGGGTCATCCCCGCTGCTGGGCGCAATCAGCGCCAGCaa                       <  1:53014/148‑1 (MQ=255)
                                                                                                                           tCTGACACCCTCAAGGGGCCGCTTTTGAGGCCCCTTGCTATTTCCGTGCTTTCACCGCCGTCGGCGCCTTGCCTGCCTTCATCTGCTGCAACAACGGTGTGCACTGGTTCGGGTCATCCCCGCTGCTGGGCGCAATCAGCGCCAGCAAc                      >  1:300530/1‑149 (MQ=255)
                                                                                                                            cTGACACCCTCAAGGGGCCGCTTTTGCGGCCCCTTGCTATTTCCGTGCTTTCACCGCCGTCGGCGCCTTGCCTGCCTTCATCTGCTGCAACCACGGTGTGCACTGGTTCGGGTCATCCCCGCTGCTGGGCGCAATCAGCGCCAGCAAcc                     >  1:255790/1‑149 (MQ=255)
                                                                                                                            cTGACACCCTCAAGGGGCCGCTTTTGCGGCCCCTTGCTATTTCCGTGCTTTCACCGCCGTCGGCGCCTTGCCTGCCTTCATCTGCTGCAACAACGGTGTGCACTGGTTCGGGTCATCCCCGCTGCTGGGCGCAATCAGCGCCAGCAAc                      <  5:41779/148‑1 (MQ=255)
                                                                                                                             tGACACCCTCAAGGGGCCGCTTTTGCGGCCCCTTGCTATTTCCGTGCTTTCACCGCCGTCGGCGCCTTGCCTGCCTTCATCTGCTGCAACAACGGTGTGCACTGGTTCGGATCATCCCCGCTACTGGGCGCAATCAGCGCCAGCAAccc                    >  6:309251/1‑149 (MQ=255)
                                                                                                                               acacCCTCAAGGGGCCGCTTTTGCGGCCCCTTGCTATTTCCGTGCTTTCACCGCCGTCGGCGCCTTGCCTGCCTTCATCTGCTGCAACAACGGTGTGCACTGGTTCGGGTCATCCCCGCTGCTGGGCGCAATc                                  <  6:312287/133‑1 (MQ=255)
                                                                                                                               acacCCTCAAGGGGCCGCTTTTGCGGCCCCTTGCTATTTCCGTGCTTTCACCGCCGTCGGCGCCTTGCCTGCCTTCATCTGCTGCAACAACGGTGTGCACTGGTTCGGGTCATCCCCGCTGCTGGGCGCAATc                                  >  5:312287/1‑133 (MQ=255)
                                                                                                                                   ccTCAAGGGGCCGCTTTTGCGGCCCCTTGCTATTTCCGTGCTTTCACCGCCGTCGGCGCCTTGCCTGCCTTCATCTGCTGCAACAACGGTGTGCACTGGTTCGGGTCATCCCCGCTGCTGGGCGCAATCAGCGCCAGCAACCccgccgc              >  1:253615/1‑149 (MQ=255)
                                                                                                                                   ccTCAAGGGGCCGCTTTTGCGGCCCCTTGCTATTTCCGTGCTTTCACCGCCGTCGGCGCCTTGCCTGCCTTCATCTGCTGCAACAACGGTGTGCACTGGTTCGGGTCATCCCCGCTGCTGGGCGCAATCAGCGCCAGCAACCccgccgc              >  4:99267/1‑149 (MQ=255)
                                                                                                                                      cAAGGGGCCGCAAAAGCGGCCCCTTGCTATTTCCGTGCTTTCACCGCCGTCGGCGCCTTGCCTGCCTTCATCTGCTGCAACAACGGTGTGCACTGGTTCGGGTCATCCCCGCTGCTGGGCGCAATCAGCGCCAGCAACCCCGCCGCCgg           >  7:255688/1‑149 (MQ=255)
                                                                                                                                       aaGGGGCCGCTTTTGCGGCCCCTTGCTATTTCCGTGCTTTCACCGCCGTCGGCGCCTTGCCTGCCTTCATCTGCTGCAACAACGGTGTGCACTGGTTCGGGTCATCCCCGCTGCTGGGCGCAATCAGCGCCAGCAACCCCGCCGCCgg           <  5:141866/148‑1 (MQ=255)
                                                                                                                                         ggggCCGCTTTTGCGGCCCCTTGCTATTTCCGTGCTTTCACCGCCGTCGGCGCCTTGCCTGCCTTCATCTGCTGCAACAACGGTGTGCACTGGTTCGGGTCATCCCCGCTGCTGGGCGCAATCAGCGCCAGCAACCCCGCCGCCGGccc        >  1:189143/1‑149 (MQ=255)
                                                                                                                                          gggCCGCTTTTGCGGCCCCTTGCTATTTCCGTGCTTTCACCGCCGTCGGCGCCTTGCCTGCCTTCATCTGCTGCAACAACGGTGTGCACTGGTTCGGGTCATCCCCGCTGCTGGGCGCAATCAGCGCCAGCAACCCCGCCGCCGGCCCt       <  4:363540/149‑1 (MQ=255)
                                                                                                                                           ggCCGCTTTTGCGGCCCCTTGCTATTTCCGTGCTTTCACCGCCGTCGGCGCCTTGCCTGCCTTCATCTGCTGCAACAACGGTGTGCACTGGTTCGGGTCATCCCCGCTGCTGGGCGCAATCAGCGCCAGCAACCCCGCCGCCGGCCCTg      <  5:322684/149‑1 (MQ=255)
                                                                                                                                            gCCGCTTTTGCGGCCCCTTGCTATTTCCGTGCTTTCACCGCCGTCGGCGCCTTGCCTGCCTTCATCTGCTGCAACAACGGTGTGCACTGGTTCGGGTCATCCCCGCTGCTGGGCGCAATCAGCGCCAGCAACCCCGCCGCCGGCCCTGc     >  2:388219/1‑149 (MQ=255)
                                                                                                                                                cTTTTGCGGCCCCTTGCTATTTCCGTGCTTTCACCGCCGTCGGCGCCTTGCCTGCCTTCATCTGCTGCAACAACGGTGTGCACTGGTTCGGGTCATCCCCGCTGCTGGGCGCAATCAGCGCCAGCAACCCCGCCGCCGGCCCTGccac  >  1:406320/1‑148 (MQ=255)
                                                                                                                                                  |                                                                                                                                                 
CCGCGAGATCATTCAATACATCGGGGATAAAGACCCGACTACCCGTCGTATCTTCGAAGACATCCTGGCGCAGGAAGAGGAGCACGCGGACGATATGTCCGACCTGTTGCAAGGCCTTTGATTTCTGACACCCTCAAGGGGCCGCTTTTGCGGCCCCTTGCTATTTCCGTGCTTTCACCGCCGTCGGCGCCTTGCCTGCCTTCATCTGCTGCAACAACGGTGTGCACTGGTTCGGGTCATCCCCGCTGCTGGGCGCAATCAGCGCCAGCAACCCCGCCGCCGGCCCTGCCAC  >  NC_002947/5522646‑5522937

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: