Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,013,762 C→G 6.0% R207R (CGG→CGC)  potF‑I ← putrescine‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,013,7620CG 6.0% 262.0 / 5.9 83R207R (CGG→CGC) potF‑Iputrescine‑binding protein
Reads supporting (aligned to +/- strand):  ref base C (37/41);  new base G (2/3);  total (39/44)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.06e-01

TCGAATAGCCCACCTTCACGCCATTGTTCGCCTCGTCGGCGCGGTTCTTCGCCTGCAGCACATCCC‑‑‑CT‑GAGAACCC‑‑‑‑‑‑CACTACCACGCAGATATCGCCGTTGG‑CCAGGTCGCCCGTGTACT‑TGGACGAGTGGAAATAGCTGATGTAAGGCCGCACTTTCATCAGCAGGTCCTTGGCCTTCTCGTAGTCCTTGGGGTCTTGGCTGTGATGGGGCAAGCCCAGGTA‑GTGCAGGGCGATCGGCAGCAGTTCGGGGCCATTGTCCAGCACCGCGACGCCACAGCTTTTCAGCTTGC  >  NC_002947/1013614‑1013904
                                                                                                                                                                |                                                                                                                                               
tCGAATAGCCCACCTTCACGCCATTGTTCGCCTCGTCGGCGCGGTTCTTCGCCTGCAGCACATCCC‑‑‑CT‑GAGAACCC‑‑‑‑‑‑CACTACCACGCAGATATCGCCGTTGG‑ACAGGTCGCCCGTGTACT‑TGGACGAGTGGAAATAGCTGATGTAAGGc                                                                                                                                                 >  7:27565/1‑149 (MQ=255)
  gAATAGCCCACCTTCACGCCATTGTTCGCCTCGTCGGCGCGGTTCTTCGCCTGCAGCACATCCC‑‑‑CT‑GAGAACCC‑‑‑‑‑‑CACTACCACGCAGATATCGCCGTTGG‑CCAGGTCGCCCGTGTACT‑TGGACGAGTGGAAATAGCTGATGTATGGCCg                                                                                                                                               >  4:331629/1‑149 (MQ=255)
         ccACCTTCACGCCATTGTTCGCCTCGTCGGCGCGGTTCTTCGCCTGCAGCACATCCC‑‑‑CT‑GAGAACCC‑‑‑‑‑‑CACTACCACGCAGATATCGCCGTTGG‑CCAGGTCGCCCGTGTACT‑TGGACGAGTGGAAATAGCTGATGTAAGGCCGCACTTtc                                                                                                                                        <  1:9143/149‑1 (MQ=255)
         ccACCTTCACGCCATTGTTCGCCTCGTCGGCGCGGTTCTTCGCCTGCAGCACATCCC‑‑‑CT‑GAGAACCC‑‑‑‑‑‑CACTACCACGCAGATATCGCCGTTGG‑CCAGGTCGCCCGTGTACT‑TGGACGAGTGGAAATAGCTGATGTAAGGCCGCACTTtc                                                                                                                                        <  6:259063/149‑1 (MQ=255)
          cACCTTCACGCCATTGTTCGCCTCGTCGGCGCGGTTCTTCGCCTGCAGCACATCCC‑‑‑CT‑GAGAACCC‑‑‑‑‑‑CACTACCACGCAGATATCGCCGTTGG‑CCAGGTCGCCCGTGTACT‑TGGACGAGTGGAAATAGCTGATGTAAGGCCGCACTTtca                                                                                                                                       >  7:318673/1‑149 (MQ=255)
           aCCTTCACGCCATTGTTCGCCTCGTCGGCGCGGTTCTTCCCCTGCAGCACATCCC‑‑‑CT‑GAGAACCC‑‑‑‑‑‑CACTACCACGCAGATATCGCCGTTGG‑CCAGGTCGCCCGTGTACT‑TGGACGAGTGGAAATAGCTGATGTAAGGCCGCACTTtcat                                                                                                                                      >  6:168941/1‑149 (MQ=255)
             cTTCACGCCATTGTTCGCCTCGTCGGCGCGGTTCTTCGCCTGCAGCACATCCC‑‑‑CT‑GAGAACCC‑‑‑‑‑‑CACTACCACGCAGATATCGCCGTTGG‑CCAGGTCGCCCGTGTACT‑TGGACGAGTGGAAATAGCTGATGAAAGTGCGGCCTTacatca                                                                                                                                    <  1:96982/149‑1 (MQ=255)
                cACGCCTTTGTTCGCCTCGTCGGGGGGGTTCTTCGCCTGCAGCACATCCCTTTCTCTTCATCCCTTTTAACTCTACCACGCAGATATCGCCGTTGG‑CCAGGTCGCCCGTGTACT‑TGGACGAGTGGAAATAGCTGATGTAAGGCCGCa                                                                                                                                             <  1:334597/147‑1 (MQ=16)
                     cATTGTTCGCCTCGTCGGCGCGGTTCTTCGCCTGCAGCACATCCC‑‑‑CT‑GAGAACCC‑‑‑‑‑‑CACTACCACGCAGATATCGCCGTTGG‑CCAGGTCGCCCGTGTACTCTGGACGAGTGGAAATAGCTGATGTAAGGCCGCACTTTCATCAGCAGGTc                                                                                                                             <  1:299876/149‑1 (MQ=255)
                     cATTGTTCGCCTCGTCGGCGCGGTTCTTCGCCTGCAGCACATCCC‑‑‑CT‑GAGAACCC‑‑‑‑‑‑CACTACCACGCAGATATCGCCGTTGG‑CCAGGTCGCCCGTGTACT‑TGGACGAGTGGAAATAGCTGATGTAAGGCCGCACTTTCATCAGCACGTcc                                                                                                                            >  2:394520/1‑149 (MQ=255)
                      aTTGTTCGCCTCGTCGGCGCGGTTCTTCGCCTGCAGCACATCCC‑‑‑CT‑GAGAACCC‑‑‑‑‑‑CACTACCACGCAGATATCGCCGTTGG‑CCAGGTCGCCCGTGTACT‑TGGACGAGTGGAAATAGCTGATGAAAGTGCGGCCTTACATcagctatttcc                                                                                                                           <  4:354127/149‑8 (MQ=38)
                         gTTCGCCTCGTCGGCGCGGTTCTTCGCCTGCAGCACATCCC‑‑‑CT‑GAGAACCC‑‑‑‑‑‑CACTACCACGCAGATATCGCCGTTGG‑CCAGGTCGCCCGTGTACT‑TGGACGAGTGGAAATAGCTGATGTAAGGCCGCACTTTCATCAGCAGGTCCTTg                                                                                                                         <  1:326040/148‑1 (MQ=255)
                         gTTCGCCTCGTCGGCGCGGTTCTTCGCCTGCAGCACATCCC‑‑‑CT‑GAGAACCC‑‑‑‑‑‑CACTACCACGCAGATATCGCCGTTGG‑CCAGGTCGCCCGTGTACT‑TGGACGAGTGGAAATAGCTGATCTAAGGCCGCACTTTCATCAGCAGGTCCTTg                                                                                                                         <  1:34315/148‑1 (MQ=255)
                             gCCTCGTCGGCGCGGTTCTTCGCCTGCAGCACATCCC‑‑‑CT‑GAGAACCC‑‑‑‑‑‑CACTACCACGCAGATATCGCCGTTGG‑CCAGGTCGCCCGTGTACT‑TGGACGAGTGGAAATAGCTGATGTAAGGCCGCACTTTCATCAGCAGGTCCTTGGCCtt                                                                                                                    >  3:10924/1‑149 (MQ=255)
                              ccTCGTCGGCGCGGTTCTTCGCCTGCAGCACATCCC‑‑‑CT‑GAGAACCC‑‑‑‑‑‑CACTACCACGCAGATATCGCCGTTGG‑CCAGGTCGCCCGTGTACT‑TGGACGAGTGGAAATAGCTGATGTAAGGCCGCACTTTCATCAGCAGGTCCTTGGCCTtc                                                                                                                   >  1:195625/1‑149 (MQ=255)
                              ccTCGTCGGCGCGGTTCTTCGCCTGCAGCACAAACC‑‑‑CT‑GAGAACCC‑‑‑‑‑‑CACTACCACGCAGATATCGCCGTTGG‑CCAGGTCGCCCGTGTACT‑TGGAGGAGTGGAAATAGCTGATGTAAGGCCGCACTTTCATCAGCCGGTCCTTGGCCTtc                                                                                                                   >  6:136518/1‑149 (MQ=255)
                                  gtcgGCGCGGTTCTTCGCCTGCAGCACATCCC‑‑‑CT‑GAGAACCC‑‑‑‑‑‑CACTACCACGCAGATATCGCCGTTGG‑CCAGGTCGCCCGTGTACT‑TGGACGAGTGGAAATAGCTGATGTAAGGCCGCACTTTCATCAGCAGGTCCTTGGCCTtctc                                                                                                                 <  1:81382/147‑1 (MQ=255)
                                  gtcgGCGCGGTTCTTCGCCTGCAGCACATCCC‑‑‑CT‑GAGAACCC‑‑‑‑‑‑CACTACCACGCAGATATCGCCGTTGG‑CCAGGTCGCCCGTGTACT‑TGGACGAGTGGAAATAGCTGATGTAAGGCCGCACTTTCATCAGCAGGTCCTTGGCCTtctc                                                                                                                 <  4:400854/147‑1 (MQ=255)
                                  gtcgGCGCGGTTCTTCGCCTGCAGCACATCCC‑‑‑CT‑GAGAACCC‑‑‑‑‑‑CACTAACACGCAGATATCGCCGTTGG‑CCAGGTCGCCCGTGTACT‑TGGACGAGTGGAAATAGCTGATGTAAGGCCGCACTTTCATCAGCAGGTCCTTGGCCTTCTCg                                                                                                                <  1:310840/148‑1 (MQ=255)
                                   tcgGCGCGGTTCTTCGCCTGCAGCACATCCC‑‑‑CT‑GAGAACCC‑‑‑‑‑‑CACTACCACGCAGATATCGCCGTTGG‑CCAGGTCGCCCGTGTACT‑TGGACGAGTGGAAATAGCTGATGTAAGGCCGCACTTTCATCAGCAGGTCCTTGGCCTTCTCGTa                                                                                                              <  3:1871/149‑1 (MQ=255)
                                     ggCGCGGTTCTTCGCCTGCAGCACATCCC‑‑‑CT‑GAGAACCC‑‑‑‑‑‑CACTACCACGCAGATATCGCCGTTGG‑CCAGGTCGCCCGTGTACT‑TGGACGAGTGGAAATAGCTGATGTAAGGCCGCACTTTCATCAGCAGGTCCTTg                                                                                                                         >  1:387528/1‑136 (MQ=255)
                                     ggCGCGGTTCTTCGCCTGCAGCACATCCC‑‑‑CT‑GAGAACCC‑‑‑‑‑‑CACTACCACGCAGATATCGCCGTTGG‑CCAGGTCGCCCGTGTACT‑TGGACGAGTGGAAATAGCTGATGTAAGGCCGCACTTTCATCAGCAGGTCCTTg                                                                                                                         <  2:387528/136‑1 (MQ=255)
                                           gTTCTTCGCCTGCAGCACATCCC‑‑‑CT‑GAGAACCC‑‑‑‑‑‑CACTACCACGCAGATATCGCCGTTGG‑CCAGGTCGCCCGTGTACT‑TGGACGAGTGGAAATGGCTGATGTAAGGCCGCACTTTCATTAGCAGGTCCTTGGCCTTCTCGTAGTCCTTgg                                                                                                      >  7:222488/1‑149 (MQ=255)
                                           gTTCTTCGCCTGCAGCACATCCC‑‑‑CT‑GAGAACCC‑‑‑‑‑‑CACTACCACGCAGATATCGCCGTTGG‑CCAGGTCGCCCGTGTACT‑TGGACGAGTGGAAATAGCTGATGTAAGGCCGCACTTTCATCAGCAGGTCCTTg                                                                                                                         <  7:133244/130‑1 (MQ=255)
                                           gTTCTTCGCCTGCAGCACATCCC‑‑‑CT‑GAGAACCC‑‑‑‑‑‑CACTACCACGCAGATATCGCCGTTGG‑CCAGGTCGCCCGTGTACT‑TGGACGAGTGGAAATAGCTGATGTAAGGCCGCACTTTCATCAGCAGGTCCTTg                                                                                                                         >  8:133244/1‑130 (MQ=255)
                                           gTTCTTCGCCTGCAGCACATCCC‑‑‑CT‑GAGAACCC‑‑‑‑‑‑CACTACCACGCAGATATCGCCGTTGG‑CCAGGTCGCCCATGTACT‑TGGACGAGTGGAAATAGCTGATGTAAGGCCGCACTTTCATCAGCAGGTCCTTGGCCTTCTCGTAGTCCTTgg                                                                                                      >  1:191904/1‑149 (MQ=255)
                                            ttcttcGCCTGCAGCACATCCC‑‑‑CT‑GAGAACCC‑‑‑‑‑‑CACTACCACGCAGATATCGCCGTTGG‑CCAGGTCGCCCGTGTACT‑TGGACGAGTGGAAATAGCTGATGTAAGGCCGCACTTTCATCAGCAGGTCCTTGGCCTTCTCGTAGTCCTTggg                                                                                                     <  8:318673/149‑1 (MQ=255)
                                                   ccTGCAGCACATCCC‑‑‑CT‑GAGAACCC‑‑‑‑‑‑CACTACCACGCAGATATCGCCGTTGGACCAGGTCGCCCGTGTACT‑TGGACGAGTGGAAATAGCTGATGTAAGGCCGCACTTTCATCAGCAGGTCCTTGGCCTTCTCGTAGTCCTTGGGGTCTTg                                                                                               >  1:305118/1‑149 (MQ=255)
                                                         gcaCATCCC‑‑‑CT‑GAGAACCC‑‑‑‑‑‑CACTACCACGCAGATATCGCCGTTGG‑CCAGGTCGCCCGTGTACT‑TGGACGAGTGGAAATAGCTGATGTAAGGCCGCACTTTCATCAGCAGGTCCTTGGCCTTCTCGTAGTCCTTGGGGTCTTGGCTGTGa                                                                                        <  5:244972/149‑1 (MQ=255)
                                                           acaTCCC‑‑‑CT‑GAGAACCC‑‑‑‑‑‑CACTACCACGCAGATATCGCCGTTGG‑CCAGGTCGCCCGTGTACT‑TGGACGAGTGGAAATAGCTGATGTAAGGCCGCACTTTCATCAGCAGGTCCTTGGCCTTCTCGTAGTCCTTGGGGTCTTGGCTGTGATg                                                                                      >  4:306675/1‑149 (MQ=255)
                                                            caTCCC‑‑‑CT‑GAGAACCC‑‑‑‑‑‑CACTACCACGCAGATATCGCCGTTGG‑CCAGGTCGCCCGTGTACT‑TGGACGAGTGGAAATAGCTGATGTAAGGCCGCACTTGCATCAGCAGGTCCTTGGCCTTCTCGTAGTCCTTGGGGTCTTGGCTGTGATgg                                                                                     >  8:155826/1‑149 (MQ=255)
                                                             aTCCC‑‑‑CT‑GAGAACCC‑‑‑‑‑‑CACTACCACGCAGATATCGCCGTTGG‑CCAGGTCGCCCGTGTACT‑TGGACGAGTGGAAATAGCTGATGTAAGGCCGCACTTTCATCAGCAGGTCCTTGGCCTTCTCGTAGTCCTTGGGGTCTTGGCTGTGATggg                                                                                    <  2:195625/149‑1 (MQ=255)
                                                             aTCCC‑‑‑CT‑GAGAACCC‑‑‑‑‑‑CACTACAACGCAGATATCGCCGTTGG‑CCAGGTCGCCCGTGTACT‑TGGACGAGTGGAAATAGCTGATGTAAGGCCGCACTTTCATCAGCAGGTCCTTGGCCTTCTCGTAGTCCTTGGGGTCTTGGCTGTGATggg                                                                                    >  8:138047/1‑149 (MQ=255)
                                                                cc‑‑‑cT‑GAGAACCC‑‑‑‑‑‑CACTACCACGCAGATATCGCCGTTGG‑CCAGGTCGCCCGTGTACT‑TGGACGAGTGGAAATAGCTGATGTAAGGCCGCACTTTCATCAGCAGGTCCTTGGCCTTCTCGTAGTCCTTGGGGTCTTGGCTGTGATGGGGc                                                                                  >  4:157662/1‑148 (MQ=255)
                                                                        gagaACCC‑‑‑‑‑‑CACTACCACGCAGATATCGCCGTTGG‑CCAGGTCGCCCGTGTACT‑TGGACGAGTGGAAATAGCTGATGTAAGGCCGCACTTTCATCAGCAGGTCCTTGGCCTTCTCGTAGTCCTTGGGGTCTTGGCTGTGATGGGGCAAGcc                                                                             >  6:355361/1‑149 (MQ=255)
                                                                            aCCC‑‑‑‑‑‑CACTACCACGCAGATATCGCCGTTGG‑CCAGGTCGCCCGTGTACT‑TGGACGAGTGGAAATAGCTGATGTAAGGCCGCACTTTCATCAGCAGGTCCTTGGCCTTCTCGTAGTCCTTGGGGTCTTGGCTGTGATGGGGCAAGCCCAgg                                                                         <  4:68296/149‑1 (MQ=255)
                                                                            aCCC‑‑‑‑‑‑CACTACCACGCAGATATCGCCGTTGG‑CCAGGTCGCCCATGTACT‑TGGACGAGTGGAAATAGCTGATGTAAGGCCGCACTTTCATCAGCAGGTCCTTGGCCTTCTCGTAGTCCTTGGGGTCTTGGCTGTGATGGGGCAAGCCCAgg                                                                         <  2:191904/149‑1 (MQ=255)
                                                                             ccc‑‑‑‑‑‑cACTACCACGCAGATATCGCCGTTGG‑CCAGGTCGCCCGTGTACT‑TGGACTAGTGGAAATAGCTGATGTAAGGCCGCACTTTCATCAGCAGTTCCTTGGCCTTCTCGTAGTCCTTGGGGTCTTGGCTGTGATGGGGCAAGCCCAGGt                                                                        <  1:210433/149‑1 (MQ=255)
                                                                              cc‑‑‑‑‑‑cACTACCACGCAGATATCGCCGTTGG‑CCAGGTCGCCCGTGTACT‑TGGACGAGTGGAAATAGCTGATGTAAGGCCGCACTTTCATCAGCAGGTCCTTGGCCTTCTCGTAGTCCTTGGGGTCTTGGCTGTGATGGGGCAAGCCCAGGTa                                                                       >  2:246668/1‑149 (MQ=255)
                                                                               c‑‑‑‑‑‑cACTACCACGCAGATATCGCCGTTGG‑CCAGGTCGCCCGTGTACT‑TGGACGAGTGGAAATAGCTGATGTAAGGCCGCACTTTCATCAGCAGGTCCTTGGCCTTCTCGTAGTCCTTGGGGTCTTGGCTGTGATGGGGCAAGCCCAGGTA‑g                                                                     >  3:317647/1‑149 (MQ=255)
                                                                                        cTACCACGCAGATATCGCCGTTGG‑CCAGGTCGCCCGTGTACT‑TGGACGAGTGGAAATAGCTGATGTAAGGCCGCACTTTCATCAGCAGGTCCTTGGCCTTCTCGTAGTCCTTGGGGTCTTGGCTGTGATGGGGCAAGCCCAGGTA‑GTGc                                                                  <  6:390894/149‑1 (MQ=255)
                                                                                           ccACGCAGATATCGCCGTTGG‑CCAGGTCGCCCGTGTACT‑TGGACGAGTGGAAATGGCTGATGTAAGGCCGCACTTTCATTAGCAGGTCCTTGGCCTTCTCGTAGTCCTTGGGGTCTTGGCTGTGATGGGGCAAGCCCAGGTA‑GTGCAg                                                                <  8:222488/148‑1 (MQ=255)
                                                                                            cACGCAGATATCGCCGTTGG‑CCAGGTCGCCCGTGTACT‑TGGACGAGTGGAAATAGCTGATGTAAGGCCGCACTTTCATCAGCAGGTCCTTGGCCTTCTCGTAGTCCTTGGGGTCTTGGCTGTGATGGGGCAAGCCCAGGTA‑GTGCAggg                                                              <  4:10924/149‑1 (MQ=255)
                                                                                            cACGCAGATATCGCCGTTGG‑CCAGGTCGCCCGTGTACT‑TGGACGAGTGGAAATAGCTGATGTAAGGCCGCACTTTCATCAGCAGGTCCTTGGCCTTCTCGTAGTCCTTGGGGTCTTGGCTGTGATGGGGCAAGCCCAGGTA‑GTGCAggg                                                              <  5:168941/149‑1 (MQ=255)
                                                                                             aCGCAGATATCGCCGTTGG‑CCAGGTCGCCCGTGTACT‑TGGACGAGTGGAAATAGCTGATGTAAGGCCGCACTTTCATCAGCAGGTCCTTGGCCTTCTCGTAGTCCTTGGGGTCTTGGCTGTGATGGGGCAAGCCCAg                                                                          <  1:377681/137‑1 (MQ=255)
                                                                                             aCGCAGATATCGCCGTTGG‑CCAGGTCGCCCGTGTACT‑TGGACGAGTGGAAATAGCTGATGTAAGGCCGCACTTTCATCAGCAGGTCCTTGGCCTTCTCGTAGTCCTTGGGGTCTTGGCTGTGATGGGGCAAGCCCAg                                                                          >  2:377681/1‑137 (MQ=255)
                                                                                             aCGCAGATATCGCCGTTCG‑CCAGGTCGCCCGTGTACT‑TGGACGAGTGGAAATAGCTGATGTAAGGCCGCACTTTCATCAGCAGGTCCTTGGCCTTCTCGTAGTCCTAGGGGTCTTGGCTGTGATGGGGCAAGCCCAGGTA‑CTGCAGGGc                                                             >  7:106501/1‑149 (MQ=255)
                                                                                                 aGATATCGCCGTTGG‑CCAGGTCGCCCGTGTACT‑TGGACGAGTGGAAATAGCTGATGTAAGGCCGCACTTTCATCAGCAGGTCCTTGGCCTTCTCGTAGTCCTTGGGGTCTTGGCTGTGATTGGGCAAGCCCAGGTA‑GTGCAgagatcc                                                          >  1:143120/1‑142 (MQ=255)
                                                                                                 aGATATCGCCGTTGG‑CCAGGTCGCCCGTGTACT‑TGGACGAGTGGAAATAGCTGATGTAAGGCCGCACTTTCATCAGCAGGTCCTTGGCCTTCTCGTAGTCCTTGGGGTCTTGGCTGTGATGGGGCAAGCCCAGGTA‑GTGCa                                                                 <  2:143120/141‑1 (MQ=255)
                                                                                                 aGATATCGCCGTTGG‑CCAGGTCGCCCGTGTACT‑TGGACGAGTGGAAATAGCTGATGTAAGGCCGCACTTTCATCAGCAGGTCCTTGGCCTTCTCGTAGTCCTTGGGGTCTTGGCTGTGATGGGGCAAGCCCAGGTA‑GTGCAGGGCGATc                                                         <  7:246205/149‑1 (MQ=255)
                                                                                                    tatCGCCGTTGG‑CCAGGTCGCCCGTGTACT‑TGGACGAGTGGAAATAGCTGATGTAAGGCCGCACTTTCATCAGCAGGTCCTTGGCCTTCTCGTAGTCCTTGGGGTCTTGGCTGTGATGGGGCAAGCCCAGGTA‑GTGCAGGGCGATCgg                                                       >  6:379801/1‑148 (MQ=255)
                                                                                                     atCGCCGTTGG‑CCAGGTCGCCCGTGTACT‑TGGACGAGTGGAAATAGCTGATGTAAGGCCGCACTTTCATCAGCAGGTCCTTGGCCTTCTCGTAGTCCTTGGGGTCTTGGCTGTGATGGGGCAAGCCCAGGTA‑GTGCAGGGCGATCgg                                                       >  6:50163/1‑147 (MQ=255)
                                                                                                      tCGCCGTTGG‑CCAGGTCGCCCGTGTACT‑TGGACGAGTGGAAATAGCTGATGTAAGGCCGCACTTTCATCAGCAGGTCCTTGGCCTTCTCGTAGTCCTTGGGGTCTTGGCTGTGATGGGGCAAGCCCAGGTA‑GTGCAGGGCGATCGgcag                                                    <  4:317647/149‑1 (MQ=255)
                                                                                                        gCCGTTGG‑CCAGGTCGCCCGTGTACT‑TGGACGAGTGGAAATAGCTGATGTAAGGCCGCACTTTCATCAGCAGGTCCTTGGCCTTCTCGTAGTCCTTGGGGTCTTGGCTGTGATGGGGCAAGCCCAGGTA‑GTGCAGGGCGATCGgcag                                                    >  6:29087/1‑147 (MQ=255)
                                                                                                          cGTTGG‑CCAGGTCGCCCGTGTACT‑TGGACGAGTGGAAATAGCTGATGTAAGGCCGCACTTTCATCAGCAGGTCCTTGGCCTTCTCGTAGTCCTTGGGGTCTTg                                                                                               >  5:107621/1‑103 (MQ=255)
                                                                                                          cGTTGG‑CCAGGTCGCCCGTGTACT‑TGGACGAGTGGAAATAGCTGATGTAAGGCCGCACTTTCATCAGCAGGTCCTTGGCCTTCTCGTAGTCCTTGGGGTCTTg                                                                                               <  6:107621/103‑1 (MQ=255)
                                                                                                           gTTGG‑CCAGGTCGCCCGTGTACT‑TGGACGAGTGGAAATAGCTGATGTAAGGCCGCACTTTCATCAGCAGGTCCTTGGCCTTCTCGTAGTCCTTGGGGTCTTGGCTGTGATGGGGCAAGCCCAGGTA‑GTGCAGGGCGATCGGCAGCAGtt                                               <  5:29087/149‑1 (MQ=255)
                                                                                                            ttGG‑CCAGGTCGCCCGTGTACT‑TGGACGAGTGGAAATAGCTGATGTAAGGCCGCACTTTCATCAGCAGGTCCTTGGCCTTCTCGTAGTCCTTGGGGTCTTGGCTGTGATGGGGCAAGCCCAGGTA‑GTGCAGGGCGATCGGCAGCAGTTc                                              <  2:280211/149‑1 (MQ=255)
                                                                                                            ttGG‑ACAGGTCGCCCGTGTACT‑TGGACGAGTGGAAATAGCTGATGTAAGGCCGCACTTTCATCAGCAGGTCCTTGGCCTTCTCGTAGTCCTTGGGGTCTTGGCTGTGATGGGGCAAGCCCAGGTA‑GTGCAGGGCGATCGGCAGCAGTTc                                              <  8:27565/149‑1 (MQ=255)
                                                                                                                 ccAGGTCGCCCGTGTACT‑TGGACGAGTGGAAATAGCTGATGTAAGGCCGCACTTTCATCAGCAGGTCCTTGGCCTTCTCGTAGTCCTTGGGGTCTTGGCTGTGATGGGGCAAGCCCAGGTA‑GTGCAGGGCGATCGGCAGCAGTTCgggg                                          >  2:2724/1‑149 (MQ=255)
                                                                                                                        gCCCGTGTACT‑TGGACGAGTGGAAATAGCTGATGTAAGGCCGCACTTTCATCAGCAGGTCCTTGGCCTTCTCGTAGTCCTTGGGGTCTTGGTTGTGATGGGGCAAGCCCAGGTA‑GTGCAGGGCGATCGGCAGCAGTTCGGGGCCATTGt                                   >  7:8920/1‑149 (MQ=255)
                                                                                                                           cGTGTACT‑TGGACGAGTGGAAATAGCTGATGTAAGGCCGCACTTTCATCAGCAGGTCCTTGGCCTTCTCGTAGTCCTTGGGGTCTTGGCTGTGATGGGGCAAGCCCAGGTA‑GTGCAGGGCGATCGGCAGCAGTTCGGGGCCATTGTCCa                                <  5:355361/149‑1 (MQ=255)
                                                                                                                           cGTGTACT‑TGGACGAGTGGAAATAGCTGATGTAAGGCCGCACTTTCATCAGCAGGTCCTTGGCCTTCTCGTAGTCCTTGGGGTCTTGGCTGTGATGGGGCAAGCCCAGGTA‑GTGCAGGGCGATCGGCAGCAGTTCGGGGCCATTGTCCa                                <  2:305118/149‑1 (MQ=255)
                                                                                                                            gtgtACT‑TGGACGAGTGGAAATAGCTGATGTAAGGCCGCACTTTCATCAGCAGGTCCTTGGCCTTCTCGTAGTCCTTGGGGTCTTGGCTGGGATGTGGCCAGGCCAGGTA‑GTGCAGGGGGATCGGCAGCAGTTCGGGGCCATTGtgaag                               >  6:337366/1‑145 (MQ=255)
                                                                                                                                  t‑tGGACGAGTGGAAATAGCTGATGTAAGGCCGCACTTTCATCAGCAGGTCCTTGGCCTTCTCGTAGTCCTTGGGGTCTTGGCTGTGATGGGGCAAGCCCAGGTACGTGCAGGGCGATCGGCAGCAGTTCGGGGCCAtt                                     <  1:259556/138‑1 (MQ=255)
                                                                                                                                  t‑tGGACGAGTGGAAATAGCTGATGTAAGGCCGCACTTTCATCAGCAGGTCCTTGGCCTTCTCGTAGTCCTTGGGGTCTTGGCTGTGATGGGGCAAGCCCAGGTA‑GTGCAGGGCGATCGGCAGCAGTTCGGGGCCAtt                                     >  2:259556/1‑137 (MQ=255)
                                                                                                                                          aGTGGAAATAGCTGATGTAAGGCCGCACTTTCATCAGCAGGTCCTTGGCCTTCTCGTAGTCCTTGGGGTCTTGGCTGTGATGGGGCAAGCCCAGGTA‑GTGCAGGGCGATCGGCAGCAGTTCGGGGCCATTGTCCAGCACCGCGACGCca                  >  3:8365/1‑149 (MQ=255)
                                                                                                                                          aGTGGAAATAGCTGATGTAAGGCCGCACTTTCATCAGCAGGTCCTTGGCCTTCTCGTAGTCCTTGGGGTCTTGGCTGTGATGGGGCAAGCCCAGGTA‑GTGCAGGGCGATCGGCAGCAGTTCGGGGCCATTGTCCAGCACCGCGACGCca                  <  1:246668/149‑1 (MQ=255)
                                                                                                                                           gTGGAAATAGCTGATGTAAGGCCGCACTTTCATCAGCAGGTCCTTGGCCTTCTCGTAGTCCTTGGGGTCTTGGCTGTGATGGGGCAAGCCCAGGTA‑GTGCAGGGCGATCGGCAGCAGTTCGGGGCCATTGTCCAGCACCGCGACGCca                  >  4:225200/1‑148 (MQ=255)
                                                                                                                                              gAAATAGCTGATGTAAGGCCGCACTTTCATCAGCAGGTCCTTGGCCTTCTCGTAGTCCTTGGGGTCTTGGCTGTGATGGGGCAAGCCCa                                                                           >  6:200888/1‑89 (MQ=255)
                                                                                                                                              gAAATAGCTGATGTAAGGCCGCACTTTCATCAGCAGGTCCTTGGCCTTCTCGTAGTCCTTGGGGTCTTGGCTGTGATGGGGCAAGCCCa                                                                           <  5:200888/89‑1 (MQ=255)
                                                                                                                                              gAAATAGCTGATGTAAGGCCGCACTTTCATCAGCAGGTCCTTGGCCTTCTCGTAGTCCTTGGGGTCTTGGCTGTGATGGGGCAAGCCCAGGTA‑GTGCAGGGCGATCGGCAGGAGTTCGGGGCCATTGTCCAGCACCGCGACGCCACAGc              >  2:276205/1‑149 (MQ=255)
                                                                                                                                              gAAATAGCTGATGTAAGGCCGCACTTTCATCAGCAGGTCCTTGGCCTTCTCGTAGTCCTTGGGGTCTTGGCTGTGATGGGGCAAGCCCAGGTA‑GTGCAGGGCGATCGGCAGCAGTTCGGGGCCATTGTCCAGCACCGCGACGCCACAGc              >  2:167334/1‑149 (MQ=255)
                                                                                                                                                 aTAGCTGATGTAAGGCCGCACTTTCATCAGCAGGTCCTTGGCCTTCTCGTAGTCCTTGGGGTCTTGGCTGTGATGGGGCAAGCCCAGGTA‑GTGCAGGGCGATCGGCAGCAGTTCGGGGc                                         <  3:65738/119‑1 (MQ=255)
                                                                                                                                                 aTAGCTGATGTAAGGCCGCACTTTCATCAGCAGGTCCTTGGCCTTCTCGTAGTCCTTGGGGTCTTGGCTGTGATGGGGCAAGCCCAGGTA‑GTGCAGGGCGATCGGCAGCAGTTCGGGGc                                         >  4:65738/1‑119 (MQ=255)
                                                                                                                                                 aTAGCTGATGTAAGGCCGCACTTTCATCAGCAGGTCCTTGGCCTTCTCGTAGTCCTTGGGGTCTTGGCTGTGATGGGGCAAGCCCAGGTA‑GTGCAGGGCGATCGGCAGCAGTTCGGGGCCATTGTCCAGCACCGCGACGCCACAGCttt           <  5:379801/149‑1 (MQ=255)
                                                                                                                                                  ctgCTGATGAAAGTGCGGCCTTACATCAGCAGGTCCTTGGCCTTCTCGTAGTCCTTGGGGTCTTGGCTGTGATGGGGCAAGCCCAGGTA‑GTGCAGGGCGATCGGCAGCAGTTCGGGGCCATTGTCCAGCACCGCGACGCCACAGCttt           >  8:104510/3‑148 (MQ=255)
                                                                                                                                                  ctgCTGATGAAAGTGCGGCCTTACATCAGCAGGTCCTTGGCCTTCTCGTAGTCCTTGGGGTCTTGGCTGTGATGGGGCAAGCCCAGGTA‑GTGCAGGGCGATCGGCAGCAGTTCGGGGCCATTGTCCAGCACCGCGACGCCACAGCtt            <  8:372513/145‑1 (MQ=255)
                                                                                                                                                  ctgCTGATGAAAGTGCGGCCTTACATCAGCAGGTCCTTGGCCTTCTCGTAGTCCTTGGGGTCTTGGCTGTGATGGGGCAAGCCCAGGTA‑GTGCAGGGCGATCGGCAGCAGTTCGGGGCCATTGTCCAGCACCGCGACGCCACAGCtt            >  7:372513/3‑147 (MQ=255)
                                                                                                                                                    gCTGATGTAAGGCCGCACTTTCATCAGCAGGTCCTTGGCCTTCTCGTAGTCCTTGGGGTCTTGGCTGTGATGGGGCAAGCCCAGGTA‑GTGCAGGGCGATCGGCAGCAGTTCGGGGCCATTGTCCAGCACCGCGACGCCACAGCTTTTCa        <  5:136518/149‑1 (MQ=255)
                                                                                                                                                    gCTGATGTAAGGCCGCACTTTCATCAGCAGGTCCTTGGCCTTCTCGTAGTCCTTGGGGTCTTGGCTGTGATGGGGCAAGCCCAGGTA‑GTGCAGGGCGATCGGCAGCAGTTCGGGGCCATTGTCCAGCACCGCGACGCCACAGCTTTTCa        <  3:306675/149‑1 (MQ=255)
                                                                                                                                                    gCTGATGTAAGGCCGCACTTTCATCAGCAGGTCCTTGGCCTTCTCGTAGTCCTTGGGGTCTTGGCTGTGATGGGGCAAGCCCAGGTA‑GTGCAGGGCGATCGGCAGCAGTTCGGGGCCATTGTCCAGCACCGCGACGCCACAGCTTTTCa        <  3:157662/149‑1 (MQ=255)
                                                                                                                                                     cTGATGTAAGGCCGCACTTTCATCAGCAGGTCCTTGGCCTTCTCGTAGTCCTTGGGGTCTTGGCTGTGATGGGGCAAGCCCAGGTA‑GTGCAGGGCGATCGGCAGCAGTTCGGGGCCATTGTCCAGCACCGCGACGCCACAGCTTTTCAg       >  6:261174/1‑149 (MQ=255)
                                                                                                                                                           tAAGGCCGCACTTTCATCAGCAGGTCCTTGGCCTTCTCGTAGTCCTTGGGGTCTTGGCTGTGATGGGGCAAGCCCAGGTA‑GTGCAGGGCGATCGGCAGCAGTTCGGGGCCATTGTCCAGCACCGCGACGCCACAGCTTTTCAGCTTGc  <  5:41620/148‑1 (MQ=255)
                                                                                                                                                                |                                                                                                                                               
TCGAATAGCCCACCTTCACGCCATTGTTCGCCTCGTCGGCGCGGTTCTTCGCCTGCAGCACATCCC‑‑‑CT‑GAGAACCC‑‑‑‑‑‑CACTACCACGCAGATATCGCCGTTGG‑CCAGGTCGCCCGTGTACT‑TGGACGAGTGGAAATAGCTGATGTAAGGCCGCACTTTCATCAGCAGGTCCTTGGCCTTCTCGTAGTCCTTGGGGTCTTGGCTGTGATGGGGCAAGCCCAGGTA‑GTGCAGGGCGATCGGCAGCAGTTCGGGGCCATTGTCCAGCACCGCGACGCCACAGCTTTTCAGCTTGC  >  NC_002947/1013614‑1013904

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