Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,395,523 G→C 7.2% A84G (GCC→GGC)  PP_3871 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,395,5230GC 7.2% 205.3 / 4.2 66A84G (GCC→GGC) PP_3871hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (28/33);  new base C (3/2);  total (31/35)
Fisher's exact test for biased strand distribution p-value = 6.59e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.81e-01

AGCCGAGAAGCTGAAGCGGTAGACCACCCGCGCACGGTTGATCAGTACCAACCCGCCACCGTCGTAGACGATGGGGTCGTACTCATCACCAGGCCGCCATCCGACCAAGGCCCGCCACAACTCCTTGCGAAAGTCGTGCAGCAGGTCGGCGGCCTGCTGGCCGCGCTCGTCCTTGGTGTCCATGGCGAACACCACCTCAAAGCGATCCTCGATGTCCTGGGTAACTGCGTTCTGCGCCTTGTTCTCGCCGGCTTTGTCATCAGCGGTGATCACATAGGCCGAGGGGTTCTTCAGT  >  NC_002947/4395375‑4395669
                                                                                                                                                    |                                                                                                                                                  
agccGAGAAGCTGAAGCGGTAGACCACCCGCGCACGGTTGATCAGTACCAACCCGCCACCGTCGTAGACGATGGGGTCGTACTCATCACCAGGCCGCCATCCGACCAAGGCCCGCCACAACTCCTTGCGAAAGTCGTGCAGCAGGTcgg                                                                                                                                                    >  5:319933/1‑149 (MQ=255)
         gCTGAAGCGGTAGACCACCCGCGCACGGTTGATCAGTACCAACCCGCCACCGTCGTAGACGATGGGGTCGTACTCATCACCAGGCCGCCATCCGACCAAGGCCCGCCACAACTCCTTGCGAAAGTCGTGCAGCAGGTCGGCGGCctgct                                                                                                                                           >  1:326626/1‑149 (MQ=255)
             aaGCGGTAGACCACCCGCGCACGGTTGATCAGTACCAACCCGCCACCGTCGTAGACGATGGGGTCGTACTCATCACCAGGCCGCCATCCGACCAAGGCCCGCCACAACTCCTTGCGAAAGTCGTGCAGCAGGTCGGCGGCCTGCTgg                                                                                                                                         <  8:306822/147‑1 (MQ=255)
                  gTAGACCACCCGCGCACGGTTGATCAGTACCAACCCGCCACCGTCGTAGACGATGGGGTCGTACTCATCACCATGCCGCCATCCGACCAAGGCCCGCCACAACTCCTTGCGAAAGTCGTGCAGCAGGTCGGCGGCCTGCTGGCCgcgc                                                                                                                                   >  6:275319/1‑148 (MQ=255)
                  gTAGACCACCCGCGCACGGTTGATCAGTACCAACCCGCCACCGTCGTAGACGATGGGGTCGTACTCATCACCAGGCCGCCATCCGACCAAGGCCCGCCACAACTCCTTGCGAAAGTCGTGCAGCAGGTCGGCGGCCTGCTGGCCgcgc                                                                                                                                   <  5:275319/148‑1 (MQ=255)
                  gTAGACCACCCGCGCACGGTTGATCAGTACCAACCCGCCACCGTCGTAGACGATGGGGTCGTACTCATCACCAGGCCGCCATCCGACCAAGGCCCGCCACAACTCCTTGCGAAAGTCGTGCAGCAGGCCGCCGACctgctgca                                                                                                                                        >  6:130069/1‑141 (MQ=255)
                  gTAGACCACCCGCGCACGGTTGATCAGTACCAACCCGCCACCGTCGTAGACGATGGGGTCGTACTCATCACCAGGCCGCCATCCGACCAAGGCCCGCCACAACTCCTTGCGAAAGTCGTGCAGCAGGCCGCCGACctgctgca                                                                                                                                        <  5:124984/143‑3 (MQ=255)
                  gTAGACCACCCGCGCACGGTTGATCAGTACCAACCCGCCACCGTCGTAGACGATGGGGTCGTACTCATCACCAGGCCGCCATCCGACCAAGGCCCGCCACAACTCCTTGCGAAAGTCGTGCAGCAGGCCGCCGACctgctgca                                                                                                                                        <  5:130069/143‑3 (MQ=255)
                  gTAGACCACCCGCGCACGGTTGATCAGTACCAACCCGCCACCGTCGTAGACGATGGGGTCGTACTCATCACCAGGCCGCCATCCGACCAAGGCCCGCCACAACTCCTTGCGAAAGTCGTGCAGCAGGCCGCCGACctgctgca                                                                                                                                        >  6:124984/1‑141 (MQ=255)
                   tAGACCACCCGCGCACGGTTGATCAGTACCAACCCGCCACCGTCGTAGACGATGGGGTCGTACTCATCACCAGGCCGCCATCCGACCAAGGCCCGCCACAACTCCTTGCGAAAGTCGTGCAGCAGGTCGGCGGCctgctgc                                                                                                                                         >  4:229224/1‑140 (MQ=255)
                   tAGACCACCCGCGCACGGTTGATCAGTACCAACCCGCCACCGTCGTAGACGATGGGGTCGTACTCATCACCAGGCCGCCATCCGACCAAGGCCCGCCACAACTCCTTGCGAAAGTCGTGCAGCAGGTCGGCGGCctgctgc                                                                                                                                         <  3:229224/141‑2 (MQ=255)
                   tAGACCACCCGCGCACGGTTGATCAGTACCAACCCGCCACCGTCGTAGACGATGGGGTCGTACTCATCACCAGGCCGCCATCCGACCAAGGCCCGCCACAACTCCTTGCGAAAGTCGTGCAGCAGGTCGGCGGCCTGCTGGCCGCGCTc                                                                                                                                 <  6:288808/149‑1 (MQ=255)
                        caccCGCGCACGGTTGATCAGTACCAACCCGCCACCGTCGTAGACGATGGGGTCGTACTCATCACCAGGCCGCCATCCGACCAAGGCCCGCCACAACTCCTTGCGAAAGTCGTGCAGCAGGTCGGCGGCCTGCTGGCCGCGCTCGTCCt                                                                                                                            >  6:145440/1‑149 (MQ=255)
                           ccGCGCACGGTTGATCAGTACCAACCCGCCACCGTCGTAGACGATGGGGTCGTACTCATCACCAGGCCGCCATCCGACCAAGGCCCGCCACAACTCCTTGCGAAAGTCGTGCAGCAGGTCGGCGGCCTGCTGGCCGCGCTCGTCCTTgg                                                                                                                         >  2:235382/1‑149 (MQ=255)
                            cgcgCACGGTTGATCAGTACCAACCCGCCACCGTCGTAGACGATGGGGTCGTACTCATCACCAGGCCGCCATCCGACCAAGGCCCGCCACAACTCCTTGCGAAAGTCGTGCAGCAGGTCGGCGGCCTGCTGGCCGCGCTCGTCCTTGgt                                                                                                                        <  7:161432/149‑1 (MQ=255)
                               gcACGGTTGATCAGTACCAACCCGCCACCGTCGTAGACGATGGGGTCGTACTCATCCCCAGGCCGCCATCCGACCAAGGCCAGCCACAACTCCATGCGAAAGTCGTGCAGCAGGTCGGCGGCCTGCTGGCCGCGCTCGTCCTTGGTGTc                                                                                                                     >  2:227687/1‑149 (MQ=255)
                                   ggTTGATCAGTACCAACCCGCCACCGTCGTAGACGATGGGGTCGTACTCATCACCAGGCCGCCATCCGACCAAGGCCCGCCACAACTCCTTGCGAAAGTCGTGCAGCAGGTCGGCGGCCTGCTGGCCGCGCTc                                                                                                                                 <  6:47890/133‑1 (MQ=255)
                                   ggTTGATCAGTACCAACCCGCCACCGTCGTAGACGATGGGGTCGTACTCATCACCAGGCCGCCATCCGACCAAGGCCCGCCACAACTCCTTGCGAAAGTCGTGCAGCAGGTCGGCGGCCTGCTGGCCGCGCTc                                                                                                                                 >  5:47890/1‑133 (MQ=255)
                                     ttGATCAGTACCAACCCGCCACCGTCGTAGACGATGGGGTCGTACTCATCACCAGGCCGCCATCCGACCAAGGCCCGCCACAACTCCTTGCGAAAGTCGTGCAGCAGGTCGGCGGCCTGCTGGCCGCGCTCGTCCTTGGTGTCCATGGc                                                                                                               <  8:15160/149‑1 (MQ=255)
                                       gATCAGTACCAACCCGCCACCGTCGTAGACGATGGGGTCGTACTCATCACCAGGCCGCCATCCGACCAAGGCCCGCCACAACTCCTTGCGAAAGTCGTGCAGCAGGTCGGCGGCCTGCTGGCCGCGCTCGTCCTTGGTGTCCATGGCg                                                                                                              >  6:331752/1‑148 (MQ=255)
                                        aTCAGTACCAACCCGCCACCGTCGTAGACGATGGGGTCGTACTCATCACCATGCCCCAATCCGACCAAGGCCCGCCACAACTCCTTGCGAAAGTCTTGAAGAAGGACGGCGGCCTGCTGGACGCGCTCGTCCTTGGTGTCCAAGGCGaa                                                                                                            >  4:247125/1‑149 (MQ=255)
                                         tCAGTACCAACCCGCCACCGTCGTAGACGATGGGGTCGTACTCATCACCAGGCCGCCATCCGACCAAGGCCCGCCACAACTCCTTGCGAAAGTCGTGCAGCAGGTCGGCGGCCTGCTGGCCGCGCTCGTCCTTGGTGTCCATGGCGAac                                                                                                           <  5:145440/149‑1 (MQ=255)
                                         tCAGTACCAACCCGCCACCGTCGTAGACGATGGGGTCGTACTCATCACCAGGCCCCCATCCGCCCAACGACCGCCACAACTCCTTGCGAAAGTCGTGCACCAGGTCGCCGGCCAGCTCGCCCCGCACGTCCTTGGTGGCCAACGCGaa                                                                                                            >  5:225165/1‑148 (MQ=255)
                                                cAACCCGCCACCGTCGTAGACGATGGGGTCGTACTCATCACCAGGCCGCCATCCGACCAAGGCCCGCCACAACTCCTTGCGAAAGTCGTGCAGCAGGTCGGCGGCCTGCTGGCCGCGCTCGTCCTTGGTGTCCATGGCGAACACCACCt                                                                                                    >  7:297896/1‑149 (MQ=255)
                                                cAACCCGCCACCGTCGTAGACGATGGGGTCGTACTCATCACCAGGCCGCCATCCGACCAAGGCCCGCCACAACTCCTTGCGAAAGTCGTGCAGCAGGTCGGCGGCCTGCTGGCCGCGCTCGTCCTTGGTGTCCATGGCGAACACCACCt                                                                                                    >  2:179953/1‑149 (MQ=255)
                                                     cGCCACCGTCGTAGACGATGGGGTCGTACTCATCACCAGGCCGCCATCCGACCAAGGCCCGCCACAACTCCTTGCGAAAGTCGTGCAGCAGGTCGGCGGCCTGCTGGCCGCGCTCGTCCTTTGTGTCCATGGCGAACACCACCTCaaa                                                                                                <  8:8729/148‑1 (MQ=255)
                                                     cGCCACCGTCGTAGACGATGGGGTCGTACTCATCACCAGGCCGCCATCCGACCAAGGCCCGCCACAACTCCTTGCGAAAGTCGTGCAGCAGGTCGGCGGCCTGCTGCCCGCGCTCGTCCTTGGTGTCCATGGCGAACACCACCTCa                                                                                                  >  7:8729/1‑146 (MQ=255)
                                                          ccGTCGTAGACGATGGGGTCGTACTCATCACCAGGCCGCCATCCGACCAAGGCCCGCCACAACTCCTTGCGAAAGTCGTGCAGCAGGTCGGCGGCCTGCTGGCCGCGCTCGTCCTTGGTGTCCATGGCGAACACCACCTCAAAGCGATc                                                                                          <  5:331752/149‑1 (MQ=255)
                                                               gtAGACGATGGGGTCGTACTCATCACCAGGCCGCCATCCGACCAAGGCCCGCCACAACTCCTTGCGAAAGTCGTGCAGCAGGTCGGCGGCCTGCTGGCCGCGCTCGTCCTTGGTGTCCATGGCGAACACCACCTCAAAGCGATCCTCGa                                                                                     <  8:297896/149‑1 (MQ=255)
                                                                tAGACGATGGGGTCGTACTCATCACCAGGCCGCCATCCGACCAAGGCCCGCCACAACTCCTTGCGAAAGTCGTGCAGCAGGTCGGCGGCCTGCTGGCCGCGCTCGTCCTTGGTGTCCATgg                                                                                                                >  7:175553/1‑121 (MQ=255)
                                                                tAGACGATGGGGTCGTACTCATCACCAGGCCGCCATCCGACCAAGGCCCGCCACAACTCCTTGCGAAAGTCGTGCAGCAGGTCGGCGGCCTGCTGGCCGCGCTCGTCCTTGGTGTCCATgg                                                                                                                <  8:175553/121‑1 (MQ=255)
                                                                tAGACGATGGGGTCGTACTCATCACCAGGCCGCCATCCGACCAAGGCCCGCCACAACTCCTTGCGAAAGTCGTGCAGCAGGTCGGCGGCCTGCTGGCCGCGCTCGTCCTTGGTGTCCATGGCGAACACCACCTCAAAGCGATCCTCGa                                                                                     <  6:279529/148‑1 (MQ=255)
                                                                tAGACGATGGGGTCGTACTCATCACCAGGCCGCCATCCGACCAAGGCCCGCCACAACTCCTTGCGAAAGTCGTGCAGCAGGTCGGCGGCCTGCTGGCCGCGCTCGTCCTTGGTGTCCATGGCGAACACCACCTCAAAGCGATCCTCGa                                                                                     >  5:279529/1‑148 (MQ=255)
                                                                       tGGGGTCGTACTCATCACCAGGCCGCCATCCGACCAAGGCCCGCCACAACTCCTTGCGAAAGTCGTGCAGCAGGTCGGCGGCCTGCTGGCCGCGCTCGTCCTTGGTGTCCATGGCGAACACCACCTCAAAGCGATCCTCGATGTCCTgg                                                                             <  6:205215/149‑1 (MQ=255)
                                                                       tGGGGTCGTACTCATCACCAGGCCGCCATCCGACCAAGGCCCGCCACAACTCCTTGCGAAAGTCGTGCAGCAGGTCGGCGGCCTGCTGGCCGCGCTCGTCCTTGGTGTCCATGGCGAACACCACCTCAAAGCGATCCTCGATGTCCTgg                                                                             <  2:326626/149‑1 (MQ=255)
                                                                           gTCGTACTCATCACCAGGCCGCCATCCGACCAAGGCCCGCCACAACTCCTTGCGAAAGTCGTGCAGCAGGTCGGCGGCCTGCTGGCCGCGCTCGTCCTTGGTGGCCATGGCGAACACCACCTCAAAGCGATCCTCGATGTCCTGGGTaa                                                                         <  1:235382/149‑1 (MQ=255)
                                                                            tCGTACTCATCACCAGGCCGCCATCCGACCAAGGCCCGCCACAACTCCTTGCGAAAGTCGTGCAGCAGGTCGGCGGCCTGCTGGCCGCGCTCGTCCTTGGTGTCCATg                                                                                                                 <  2:285959/108‑1 (MQ=255)
                                                                            tCGTACTCATCACCAGGCCGCCATCCGACCAAGGCCCGCCACAACTCCTTGCGAAAGTCGTGCAGCAGGTCGGCGGCCTGCTGGCCGCGCTCGTCCTTGGTGTCCATg                                                                                                                 >  1:285959/1‑108 (MQ=255)
                                                                             cGTACTCATCACCAGGCCGCCATCCGACCAAGGCCCGCCACAACTCCTTGCGAAAGTCGTGCAGCAGGTCGGCGGCCTGCTGGCCGCGCTCGTCCTTGGTGTCCATGGCGAACACCACCTCAAAGCGATCCTCGATGTCCTGGGTAACt                                                                       <  3:88910/149‑1 (MQ=255)
                                                                                  tcatcaCCAGGCCGCCATCCGACCAAGGCCTGCCACAACTCCTTGCGAAAGTCGTGCAGCAGGTCGGCGGCCTGCTGGCCGCGCTCGTCCTTGGTGTCCATGGCGAACACCACCTCAAAGCGATCCTCGATGTCCTGGGTAACTGCGtt                                                                  >  1:197096/1‑149 (MQ=255)
                                                                                  tcatcaCCAGGCCGCCATCCGACCAAGGCCCGCCACAACTCCTTGCGAAAGTCGTGCAGCAGGTCGGCGGCCTGCTGGCCGCGCTCGTCCTTg                                                                                                                          <  5:143467/93‑1 (MQ=255)
                                                                                  tcatcaCCAGGCCGCCATCCGACCAAGGCCCGCCACAACTCCTTGCGAAAGTCGTGCAGCAGGTCGGCGGCCTGCTGGCCGCGCTCGTCCTTg                                                                                                                          >  6:143467/1‑93 (MQ=255)
                                                                                  tcatcaCCAGGCCGCCATCCGACCAAGGCCCGCCACAACTCCTTGCGAAAGTCGTGCAGCAGGTCGGCGGCCTGCTGGCCGCGCTCGTCCTTGGTGTCCATGGCGAACACCACCTCAAAGCGATCCTCGATGTCCTGGGTAACTGCGtt                                                                  >  5:173720/1‑149 (MQ=255)
                                                                                  tcatcaCCAGGCCGCCATCCGACCAAGGCCCGCCACAACTCCTTGCGAAAGTCGTGCAGCAGGTCGGCGGCCTGCTGGCCGCGCTCGTCCTTGGTGTCCATGGCGAACACCACCTCAAAGCGATCCTCGATGTCCTGGGTAACTGCGtt                                                                  >  1:310922/1‑149 (MQ=255)
                                                                                    atcaCCAGGCCGCCATCCGACCAAGGCCCGCCACAACTCCTTGCGAAAGTCGTGCAGCAGGTCGGCGGCCTGCTGGCCGCGCTCGTCCTTGGTGTCCATGGCGAACACCACCTCAAAGCGATCCTCGATGTCCTGGGTAACTGCGTTCt                                                                <  3:247125/149‑1 (MQ=255)
                                                                                           ggCCGCCATCCGACCAAGGCCCGCCACAACTCCTTGCGAAAGTCGTGCAGCAGGTCGGCGGCCTGCTGGCCGCGCTCGTCCTTGGTGTCCATGGCGAACACCACCTCAAAGCGATGCTCGATGTCCTGGGTAACTGCGTTCTGCGCCtt                                                         <  2:310922/149‑1 (MQ=255)
                                                                                           ggCCGCCATCCGACCAAGGCCCGCCACAACTCCTTGCGAAAGTCGTGCAGCAGGTCGGCGGCCTGCTGGCCGCGCTCGTCCTTGGTGTCCATGGCGAACACCACCTCAAAGCGATCCTCGATGTCCTGGGTAACTGCGTTCTGCGCCtt                                                         <  2:68479/149‑1 (MQ=255)
                                                                                               gccATCCGACCAAGGCCCGCCACAACTCCTTGCGAAAGTCGTGCAGCAGGTCGGCGGCCTGCTGGCCGCGCTCGTCCTTGGTGTCCATGGCGAACACCACCTCAAAGCGATCCTCGATGTCCTGGGt                                                                           >  4:224546/1‑127 (MQ=255)
                                                                                               gccATCCGACCAAGGCCCGCCACAACTCCTTGCGAAAGTCGTGCAGCAGGTCGGCGGCCTGCTGGCCGCGCTCGTCCTTGGTGTCCATGGCGAACACCACCTCAAAGCGATCCTCGATGTCCTGGGt                                                                           <  3:224546/127‑1 (MQ=255)
                                                                                                 cATCCGACCAAGGCCCGCCACAACTCCTTGCGAAAGTCGTGCAGCAGGTCGGCGGCCTGCTGGCCGCGCTCGTCCTTGGTGTCCATGGCGAACACCACCTCAAAGCGATCCTCGATGTCCTgg                                                                             <  6:266592/123‑1 (MQ=255)
                                                                                                 cATCCGACCAAGGCCCGCCACAACTCCTTGCGAAAGTCGTGCAGCAGGTCGGCGGCCTGCTGGCCGCGCTCGTCCTTGGTGTCCATGGCGAACACCACCTCAAAGCGATCCTCGATGTCCTgg                                                                             >  5:266592/1‑123 (MQ=255)
                                                                                                  aTCCGACCAATGCCCGACACAACTCCTTGCGAAAGTCGTGCAGCAGGTCTGCGGCCTGCTGGCCGCGCTCGTCCTTGTTGTCCATGGCGAACACCACCTCAAAGAGATCATCGATGTCCTGGTTAACTGCGTTATGCGCCTTGTtcac                                                   >  4:42901/1‑146 (MQ=255)
                                                                                                  aTCCGACCAAGGCCCGCCACAACTCCTTGCGAAAGTCGTGCAGCAGGTCGGCGGCCTGCTGGCCGCGCTCGTCCTTGGTGTCCATGGCGAACACCACCTCAAAGCGATCCTCGATGTCCTGGGTAACTGCGTTCTGCGCCTTGTTCTCg                                                  >  4:173155/1‑149 (MQ=255)
                                                                                                  aTCCGACCAAGGCCCGCCACAACTCCTTGCGAAAGTCGTGCAGCAGGTCGGCGGCCTGCTGGCCGCGCTCGTCCTTGGTGTCCATGGCGAACACCACCTCAAAGCGATCCTCGATGTCCTGGGTAACTGCGTTCTGCGCCTTGTTCTCg                                                  >  6:298045/1‑149 (MQ=255)
                                                                                                   tCCGACGAAGGCGCGCCACAACTGCTTGCGTAAGTCGTGCAGCAGGTCGGCGGCCTGCTGGCCGCGGTCGTCCTTGGTGTCCATGGCGAACACCACCTCAAAGCGATCCTCGATGTCCTGGGTAAGTGCGTTCTGCGCCTTGTTCTCGc                                                 <  6:225165/149‑1 (MQ=255)
                                                                                                   tCCGACCAAGGCCCGCCACAACTCCTTGCGAAAGTCGTGCAGCAGGTCGGCGGCCTGCTGGCCGCGCTCGTCCTTGGTGTCCATGGCGAACACCACCTCAAAGCGATCCTCGATGTCCTGGGTAACTGCGTTCTGCGCCTTGTTCTCGc                                                 >  6:93966/1‑149 (MQ=255)
                                                                                                      gACCAAGGCCCGCCACAACTCCTTGCGAAAGTCGTGCAGCAGGTCGGCGGCCTGCTGGCCGCGCTCGTCCTTGGTGTCCATGGCGAACACCACCTCAAAGCGATCCTCGATGTCCTGGGTAACTGCGTTCTGCGCCTTGTTCTCGCCgg                                              >  8:263679/1‑149 (MQ=255)
                                                                                                        ccAAGGCCCGCCACAACTCCTTGCGAAAGTCGTGCAGCAGGTCGGCGGCCTGCTGGCCGCGCTCGTCCTTGGTGTCCATGGCGAACACCACCTCAAAGCGATCCTCGATGTCCTGGGTAACTGCGTTCTGCGCCTTGTTCTCGCCGGCt                                            <  3:173155/149‑1 (MQ=255)
                                                                                                         cAAGGCCCGCCACAACTCCTTGCGAAAGTCGTGCAGCAGGTCGGCGGCCTGCTGGCCGCGCTCGTCCTTGGTGTCCATGGCGAACACCACCTCAAAGCGATCCTCGATGTCCTGGGTAACTGCGTTCTGCGCCTTGTTCTCGCCGGCtt                                           >  1:79870/1‑149 (MQ=255)
                                                                                                            ggCCCGCCACAACTCCTTGCGAAAGTCGTGCAGCAGGTCGGCGGCCTGCTGGCCGCGCTCGTCCTTGGTGTCCATGGCGAACACCACCTCAAAGCGATCCTCGATGTCCTGGGTAACTGCGTTCTGCGCCTTGTTCTCGCCGGCTTTGt                                        <  3:75673/149‑1 (MQ=255)
                                                                                                             gCCCGCCACAACTCCTTGCGAAAGTCGTGCAGCAGGTCGGCGGCCTGCTGGCCGCGCTCGTCCTTGGTGTCCATGGCGAACACCACCTCAAAGCGATCCTCGATGTCCTGGGTAACTGCGTTCTGCGCCTTGTTCTCGCCGGCTTTGtc                                       <  7:263679/149‑1 (MQ=255)
                                                                                                                     caACTCCTTGCGAAAGTCGTGCAGCAGGTCGGCGGCCTGCTGGCCGCGCTCGTCCTTGGTGTCCATGGCGAACACCACCTCAAAGCGATCCTCGATGTCCTGGGTAACTGCGTTCTGCGCCTTGTTCTCGCCGGCTTTGTCATCAGCgg                               >  1:43168/1‑149 (MQ=255)
                                                                                                                                gAAAGTCGTGCAGCAGGTCGGCGGCCTGCTGGCCGCGCTCGTCCTTGGTGTCCATGGCGAACACCACCTCAAAGCGATCCTCGATGTCCTGGGTAACTGCGTTCTGCGCCTTGTTCTCGCCGGCTTTGTCATCAGCGGTGATCACATAg                    <  2:79870/149‑1 (MQ=255)
                                                                                                                                       gTGCAGCAGGTCGGCGGCCTGCTGGCCGCGCTCGTCCTTGGTGTCCATGGCGAACACCACCTCAAAGCGATCCTCGATGTCCTGGGTAACTGCGTTCTGCGCCTTGTTCTCGCCGGCTTTGTCATCAGCGGTGATCACATAGGCCGAg              <  7:115737/148‑1 (MQ=255)
                                                                                                                                          cagcagGTCGGCGGCCTGCTGGCCGCGCTCGTCCTTGGTGTCCATGGCGAACACCACCTCAAAGCGATCCTCGATGTCCTGGGTAACTGCGTTCTGCGCCTTGTTCTCGCCGGCTTTGTCATCAGCGGTGATCACATAGGCCGAGGGGt          <  2:197096/149‑1 (MQ=255)
                                                                                                                                              agGTCGGCGGCCTGCTGGCCGCGCTCGTCCTTGGTGTCCATGGCGAACACCACCTCAACGCGATCCTCGATGTCCTGGGTAACTGCGTTCTGCGCCTTGTTCTCGCCGGCTTTGTCATCAGCGGTGATCACATAGGCCGAGGGGttctt      <  6:173720/149‑1 (MQ=255)
                                                                                                                                              agGTCGGCGGCCTGCTGGCCGCGCTCGTCCTTGGTGTCCATGGCGAACACCACCTCAAAGCGATCCTCGATGTCCTGGGTAACTGCGTTCTGCGCCTTGTTCTCGCCGGCTTTGTCATCAGCGGTGATCACATAGGCCGAGGGGttctt      <  2:43168/149‑1 (MQ=255)
                                                                                                                                                   ggcggcCTGCTGGCCGCGCTCGTCCTTGGTGTCCATGGCGAACACCACCTCAAAGCGATCCTCGATGTCCTGGGTAACTGCGTTCTGCGCCTTGTTCTCGCCGGCTTTGtcat                                     <  7:121817/113‑1 (MQ=255)
                                                                                                                                                   ggcggcCTGCTGGCCGCGCTCGTCCTTGGTGTCCATGGCGAACACCACCTCAAAGCGATCCTCGATGTCCTGGGTAACTGCGTTCTGCGCCTTGTTCTCGCCGGCTTTGtcat                                     >  8:121817/1‑113 (MQ=255)
                                                                                                                                                   ggcggcCTGCTGGCCGCGCTCGTCCTTGGTGTCCATGGCGAACACCACCTCAAAGCGATCCTCGATGTCCTGGGTAACTGCGTTCTGCGCCTTGTTCTCGCCGGCTTTGTCATCAGCGGTGATCACATAGGCCGAGGGGTTCTTCAGt  <  1:179953/148‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
AGCCGAGAAGCTGAAGCGGTAGACCACCCGCGCACGGTTGATCAGTACCAACCCGCCACCGTCGTAGACGATGGGGTCGTACTCATCACCAGGCCGCCATCCGACCAAGGCCCGCCACAACTCCTTGCGAAAGTCGTGCAGCAGGTCGGCGGCCTGCTGGCCGCGCTCGTCCTTGGTGTCCATGGCGAACACCACCTCAAAGCGATCCTCGATGTCCTGGGTAACTGCGTTCTGCGCCTTGTTCTCGCCGGCTTTGTCATCAGCGGTGATCACATAGGCCGAGGGGTTCTTCAGT  >  NC_002947/4395375‑4395669

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: ‑