Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,522,794 T→A 10.0% intergenic (+28/+12) 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,7940TA10.0% 222.9 / 12.8 80intergenic (+28/+12)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 = 9.05e-01

CGAGATCATTCAATACATCGGGGATAAAGACCCGACTACCCGTCGTATCTTCGAAGACATCCTGGCGCAGGAAGAGGAGCACGCGGACGATATGTCCGACCTGTTGCAAGGCCTTTGATTTCTGACACCCTCAAGGGGCCGCTTTTGCGGCCCCTTGCTATTTCCGTGCTTTCACCGCCGTCGGCGCCTTGCCTGCCTTCATCTGCTGCAACAACGGTGTGCACTGGTTCGGGTCATCCCCGCTGCTGGGCGCAATCAGCGCCAGCAACCCCGCCGCCGGCCCTGCCAC  >  NC_002947/5522649‑5522937
                                                                                                                                                 |                                                                                                                                               
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                                                                                                                                   <  8:25030/149‑4 (MQ=255)
           aaTACATCGGGGATAAAGACCCGACTACCCGTCGTATCTTCGAAGACATCCTGGCGCAGGAAGAGGAGCACGCGGACGATATGTCCGACCTGTTGCAAGGCCTTTGATTTCTGACACCCTCAAGGGGCCGCAAAAGCGGCCCCTTgagg                                                                                                                                   <  5:70125/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                                                                                                                         >  7:330410/1‑130 (MQ=255)
                                        cgtcgtATCTTCGAAGACATCCTGGCGCAGGAAGAGGAGCACGCGGACGATATGTCCGACCTGTTGCAAGGCCTTTGATTTCTGACACCCTCAAGGGGCCGCTTTTGCGGCCCCTTGCTATTTCCGTGCt                                                                                                                         <  8:330410/130‑1 (MQ=255)
                                             tATCTTCGAAGACATCCTGGCGCAGGAAGAGGAGCACGCGGACGATATGTCCGACCTGTTGCAAGGCCTTTGATTTCTGACACCCTCAAGGGGCCGCTTTTGCGGCCCCTTGCTATTTCCGTGCTTTCACCGCCGTCGGCGCCTtgcct                                                                                                 >  4:157080/1‑149 (MQ=255)
                                                      aGACATCCTGGCGCAGGAAGAGGAGCACGCGGACGATATGTCCGACCTGTTGCAAGGCCTTTGATTTCTGACACCCTCAAGGGGCCGCTTTTGCGGCCCCTTGCTATTTCCGTGCTTTCACCGCCGTCGGCGCCTTGCCTGCCTTCATc                                                                                        >  5:181733/1‑149 (MQ=255)
                                                      aGACATCCTGGCGCAGGAAGAGGAGCACGCGGACGATATGTCCGACCTGTTGCAAGGCCTTTGATTTCTGACACCCTCAAGGGGCCGCTTTTGCGGCCCCTTGCTATTTCCGTGCTTTCACCGCCGTCGGCGCCTTGCCTGCCTTCATc                                                                                        <  1:331017/149‑1 (MQ=255)
                                                         cATCCTGGCGCAGGAAGAGGAGCACGCGGACGATATGTCCGACCTGTTGCAAGGCCTTTGATTTCTGACACCCTCAAGGGGCCGCTTTTGCGGCCCCTTGCTATTTCCGTGCTTTCACCGCCGTCGGCGCCTTGCCTGCCTTCATctgc                                                                                     <  4:265664/149‑1 (MQ=255)
                                                           tCCTGGCGCAGGAAGAGGAGCACGCGGACGATATGTCCGACCTGTTGCAAGGCCTTTGATTTCTGACACCCTCAAGGGGCCGCTTTTGCGGCCCCTTGCTATTTCCGTGCTTTCACCGCCGTCGGCGCCTTGCCTGCCTTCATctgctg                                                                                   <  6:181733/149‑1 (MQ=255)
                                                           tCCTGGCGCAGGAAGAGGAGCACGCGGACGATATGTCCGACCTGTTGCAAGGCCTTTGATTTCTGACACCCTCAAGGGGCCGCTTTTGCGGCCCCTTGCTATTTCCGTGCTTTCACCGCCGTCGGCGCCTTGCCTGCCTTCATctgctg                                                                                   <  3:157080/149‑1 (MQ=255)
                                                            ccTGGCGCAGGAAGAGGAGCACGCGGACGATATGTCCGACCTGTTGCAAGGCCTTTGATTTCTGACACCCTCAAGGGGCCGCTTTTGCGGCCCCt                                                                                                                                        >  5:44266/1‑95 (MQ=255)
                                                            ccTGGCGCAGGAAGAGGAGCACGCGGACGATATGTCCGACCTGTTGCAAGGCCTTTGATTTCTGACACCCTCAAGGGGCCGCTTTTGCGGCCCCt                                                                                                                                        <  6:44266/95‑1 (MQ=255)
                                                             cTGGCGCAGGAAGAGGAGCACGCGGACGATATGTCCGACCTGTTGCAAGGCCTTTGATTTCTGACACCCTCAAGGGGCCGCTTTTGCGGCCCCTTGCTATTTCCGTGCTTTCACCGCCGTCGGCGCCTTGCCTGCCTTCATCtgctgc                                                                                  <  4:85111/148‑1 (MQ=255)
                                                             cTGGCGCAGGAAGAGGAGCACGCGGACGATATGTCCGACCTGTTGCAAGGCCTTTGATTTCTGACACCCTCAAGGGGCCGCTTTTGCGGCCCCTTGCTATTTCCGTGCTTTCACCGCCGTCGGCGCCTGGCCTGCATTCATCtgctgc                                                                                  >  8:188449/1‑148 (MQ=255)
                                                                gcgcAGGAAGAGGAGCACGCGGACGATATGTCCGACCTGTTGCAAGGCCTTTGATTTCTGACACCCTCAAGGGGCCGCTTTTGCGGCCCCTTGCTATTTCCGTGCTTTCACCGCCGTCGGCGCCTTGCCTGCCTTCATCTGCTGcaac                                                                               <  8:39702/148‑1 (MQ=255)
                                                                    aGGAAGAGGAGCACGCGGACGATATGTCCGACCTGTTGCAAGGCCTTTGATTTCTGACACCCTCAAGGGGCCGCTTTTGCGGCCCCTTGCTATTTCCGTGCTTTCACCGCCGTCGGCGCCTTGCCTGCCTTCATCTGCTGCAACAACgg                                                                          >  7:191829/1‑149 (MQ=255)
                                                                    aGGAAGAGGAGCACGCGGACGATATGTCCGACCTGTTGCAAGGCCTTTGATTTCTGACACCCTCAAGGGGCCGCTTTTGCGGCCCCTTGCTATTTCCGTGCTTTCACCGCCGTCGGCGCCTTGCCTGCCTTCATCTGCTGCAACAACgg                                                                          >  2:226894/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                             <  5:126971/147‑1 (MQ=255)
                                                                                                                   tGATTTCTGACACCCTCAAGGGGCCGCTTTTGCGGCCCCTTGCTATTTCCGTGCTTTCACCGCCGTCGGCGCCTTGCCTGCCTTCATCTGCTGCAACAACGGTGTGCACTGGTTCGGGTCATCCCCGCTGCTGGGCGCAATCAgcgc                             <  6:312055/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                                                                                                           <  5:84304/67‑1 (MQ=255)
                                                                                                                     aTTTCTGACACCCTCAAGGGGCCGCTTTTGCGGCCCCTTGCTATTTCCGTGCTTTCACCGCCGTCgg                                                                                                           >  6:84304/1‑67 (MQ=255)
                                                                                                                      tttCTGACACCCTCAAGGGGCCGCTTTTGCGGCCCCTTGCTATTTCCGTGCTTTCACCGCCGTCGGCGCCTTGCCTGCCTTCATCTGCTGCAACAACGGTGTGCACTGGTTCGGGTCATCCCCGCTGCTGGGCGCAATCAgcgc                             <  6:211698/144‑1 (MQ=255)
                                                                                                                      tttCTGACACCCTCAAGGGGCCGCTTTTGCGGCCCCTTGCTATTTCCGTGCTTTCACCGCCGTCGGCGCCTTGCCTGCCTTCATCTGCTGCAACAACGGTGTGCACTGGTTCGGGTCATCCCCGCTGCTGGGCGCAATCAgcgc                             >  5:211698/1‑144 (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)
                                                                                                                                                 |                                                                                                                                               
CGAGATCATTCAATACATCGGGGATAAAGACCCGACTACCCGTCGTATCTTCGAAGACATCCTGGCGCAGGAAGAGGAGCACGCGGACGATATGTCCGACCTGTTGCAAGGCCTTTGATTTCTGACACCCTCAAGGGGCCGCTTTTGCGGCCCCTTGCTATTTCCGTGCTTTCACCGCCGTCGGCGCCTTGCCTGCCTTCATCTGCTGCAACAACGGTGTGCACTGGTTCGGGTCATCCCCGCTGCTGGGCGCAATCAGCGCCAGCAACCCCGCCGCCGGCCCTGCCAC  >  NC_002947/5522649‑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: