Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 467,098 A→G 7.7% intergenic (‑31/+22) PP_5445 ← / ← PP_t71 hypothetical protein/tRNA‑Met

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947467,0980AG 7.7% 184.8 / 4.6 65intergenic (‑31/+22)PP_5445/PP_t71hypothetical protein/tRNA‑Met
Reads supporting (aligned to +/- strand):  ref base A (32/28);  new base G (3/2);  total (35/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.12e-01

GCTGGGTGACCGGCTTCTTGGTGGGGCTGATTTCGCTGAGATACCGTTTTAGCGCTACTTCCAAGGTCATGCGCTCGGCGCCGCTGCGCTGGATGTACGCGCCCCTGACCATTTATTTCTTGCGGGCATGAAAAAGCCGCGCAATCGCGCGGCTTTGTCGTATCTGGTGGGCCCACACGGA‑‑CTCGAACCGTGGACCAAAGGATTATGAGTCCTCTGCTCTAACCAACTGAGCTATAGGCCCTCAGTAAGTGGGCGGATTATAACGAGGGTTTCGCACCTGTGCCATCTGAA  >  NC_002947/466956‑467246
                                                                                                                                              |                                                                                                                                                      
gcTGGGTGACCGGCTTCTTGGTGGGGCTGATTTCGCTGAGATACCGTTTTAGCGCTACTTCCAAGGTCATGCGCTCGGCGCCGCTGCGCTGGATGTACGCGCCCCTGACCATTTATTTCTTGCGGGCATGAAAAAGCCGCGCAATcgcg                                                                                                                                                  <  6:197935/149‑1 (MQ=255)
 cTGGGTGACCGGCTTCTTGGTGGGGCTGATTTCGCTGAGATACCGTTTTAGCGCTACTTCCAAGGTCATGCGCTCGGCGCCGCTGCGCTGGATGTACGCGCCCCTGACCATTTATTTCTTGCGGGCATGAAAAAGCCGCTCAATcgc                                                                                                                                                   >  4:222774/1‑147 (MQ=255)
      tGACCGGCTTCTTGGTGGGGCTGATTTCGCTGAGATACCGTTTTAGCGCTACTTCCAAGGTCATGCGCTCGGCGCCGCTGCGCTGGATGTACGCGCCCCTGACCATTTATTTCTTGCGGGCATGAAAAAGCCGCGCAATCGCGCGGCtt                                                                                                                                            >  1:140872/1‑149 (MQ=255)
          cGGCTTCTTGGTGGGGCTGATTTCGCTGAGATACCGTTTTAGCGCTACTTCCAAGGTCATGCGCTCGGCGCCGCTGCGCTGGATGTACGCGCCCCTGACCATTTATTTCTTGCGGGCATGAAAAAGCCGCGCAATCGCGCGGCTTTGTc                                                                                                                                        >  5:79017/1‑149 (MQ=255)
          cGGCTTCTTGGTGGGGCTGATTTCGCTGAGATACCGTTTTAGCGCTACTTCCAAGGTCATGCGCTCGGCGCCGCTGCGCTGGATGTACGCGCCCCTGACCATTTATTTCTTGCGGGCATGAAAAAGCCGCGCAATCGCGCGGCTTTGTc                                                                                                                                        <  2:61123/149‑1 (MQ=255)
                 ttGGTGGGGCTGATTTCGCTGAGATACCGTTTTAGCGCTACTTCCAAGGTCATGCGCTCGGCGCCGCTGCGCTGGATGTACGCGCCCCTGACCATTTATTTCTTGCGGGCATGAAAAAGCCGCGCAATCGCGCGGCttt                                                                                                                                           >  4:235158/1‑139 (MQ=255)
                 ttGGTGGGGCTGATTTCGCTGAGATACCGTTTTAGCGCTACTTCCAAGGTCATGCGCTCGGCGCCGCTGCGCTGGATGTACGCGCCCCTGACCATTTATTTCTTGCGGGCATGAAAAAGCCGCGCAATCGCGCGGCttt                                                                                                                                           <  3:235158/139‑1 (MQ=255)
                    gtggGGCTGATTTCGCTGAGATACCGTTTTAGCGCTACTTCCAAGGTCATGCGCTCGGCGCCGCTGCGCTGGATGTACGCGCCCCTGACCATTTATTTCTTGCGGGCATGAAAAAGCCGCGCAATCGCGCGGCTTTGTCGTATCtggt                                                                                                                               <  7:193312/148‑1 (MQ=255)
                         gCTGATTTCGCTGAGATACCGTTTTAGCGCTACTTCCAAGGTCATGCGCTCGGCGCCGCTGCGCTGGATGTACGCGCCCCTGACCATTTATTTCTTGCGGGCATGAAAAAGCCGCGCAATCGCGCGGCTTTGTCGTATCTGGTGGGcc                                                                                                                          >  8:83034/1‑148 (MQ=255)
                             aTTTCGCTGAGATACCGTTTTAGCGCTACTTCCAAGGTCATGCGCTCGGCGCCGCTGCGCTGGATGTACGCGCCCCTGACCATTTATTTCTTGCGGGCATGAAAAAGCCGCGCAATCGCGCGGCTTTGTCGTATCTGGTGGGCCCacac                                                                                                                     <  4:168710/149‑1 (MQ=255)
                             aTTTCGCTGAGATACCGTTTTAGCGCTACTTCCAAGGTCATGCGCTCGGCGCCGCTGCGCTGGATGTACGCGCCCCTGACCATTTATTTCTTGCGGGCATGAAAAAGCCGCGCAATCGCGCGGCTTTGTCGTATCTGGTGGGCCCacac                                                                                                                     <  7:193062/149‑1 (MQ=255)
                               ttCGCTGAGTTAGCGTTTTAGCGCTACTTCCAAGGTCATGCGGTTGGCGCCGCTGCGGTGGATGTACGCGCCCCTGACCATTTATTTCTTGCGGGCATGAATAAGCCGCGCGATTGCGCGGCtttttg                                                                                                                                        <  4:59672/128‑4 (MQ=255)
                               ttCGCTGAGATACAGTTTTAGCGCTACTTCCAAGCTCAAGCGCTCGGCGCCGCTGCCCTGGATGTACGAGCCCCTGACCATTAAATTCTTGCGGGCATGAAAAACACGCCCGATTGCGCGGCtatatg                                                                                                                                        >  3:59672/1‑123 (MQ=255)
                                  gCTGAGATACCGTTTTAGCGCTACTTCCAAGGTCATGCGCTCGGCGCCGCTGCGCTGGATGTACGCGCCCCTGACCATTTATTTCTTGCGGGCATGAAAAAGCCGCGCAATCGCGCGGCTTTGTCGTATCTGGTGGGCCCACACGGA‑‑c                                                                                                               <  2:177507/148‑1 (MQ=255)
                                  gCTGAGATACCGTTTTAGCGCTACTTCCAAGGTCATGCGCTCGGCGCCGCTGCGCTGGATGTACGCGCCCCTGACCATTTATTTCTTGCGGGCATGAAAAAGCCGCGCAATCGCGCGGCTTTGTCGTATCTGGTGGGCCCACACGGA‑‑Ct                                                                                                              >  2:111783/1‑149 (MQ=255)
                                   cTGAGATACCGTTTTAGCGCTACTTCCAAGGTCATGCGCTCGGCGCCGCTGCGCTGGATGTACGCGCCCCTGACCATTTATTTCTTGCGGGCATGAAAAAGCCGCGCAATCGCGCGGCTTTGTCGTATCTGGTGGGCCCACACGGA‑‑CTc                                                                                                             <  2:140872/149‑1 (MQ=255)
                                    tGAGATACCGTTTTAGCGCTACTTCCAAGGTCATGCGCTCGGCGCCGCTGCGCTGGATGTACGCGCCCCTGACCATTTATTTCTTGCGGGCATGAAAAAGCCGCGCAATCGCGCGGCTTTGTCGTATCTGGTGGGCCCACACGGA‑‑CTc                                                                                                             >  7:217164/1‑148 (MQ=255)
                                    tGACATACCGTTTTAGCGCTACTTCCACGGGAATGCGCTCGGCGCCGCTGCGCTGGATGTACGCGCCCCTGACCATTTATTTCTTGCGGTCATGAAAAAGCCGCGCAATCGCGCGGCTTTGTCGTATCTGGTGGGCCCACACGGA‑‑CTCg                                                                                                            <  8:152796/149‑1 (MQ=255)
                                      agaTACCGTTCTAGCGCTACTTCCACGGTCATGCGCTCGGCGCCGCTGCGCTGGATGTACGCGCCCCTGACCATTTATTTCTTGCGGGCATGAAAAAGCCGCGCAATCGCGCGGCTTTGTCGTATCTGGTGGGCCCACACGGA‑‑CTCg                                                                                                            <  4:349097/147‑1 (MQ=255)
                                           ccGTTTTAGCGCTACTTCCAAGGTCATGCGCTCGGCGCCGCTGCGCTGGATGTACGCGCCCCTGACCATTTATTTCTTGCGGGCATGAAAAAGCCGCGCAATCGCGCGGCTTTGTCGTATCTGGTGGGCCCACACGGA‑‑CTCGAACCGt                                                                                                      <  8:211192/148‑1 (MQ=255)
                                                   gcgcTACTTCCCAGGTCATGCGCTCGGCGCCGCTGCGCTGGATGTACGCGCCCCTGACCATTTATTTCTTGCGGGCATGAAAAAGCCGCGCAATCGCGCGGCTTTGTCGTATCTGGTGGGCCCACACGGA‑‑CTCGAACCGTGGACCaaa                                                                                              <  5:8507/148‑1 (MQ=255)
                                                   gcgcTACTTCCAAGGTCATGCGCTCGGCGCCGCTGCGCTGGATGTACGCGCCCCTGACCATTTATTTCTTGCGGGCATGAAAAAGCCGCGCAATCGCGGGGCTTTGTCGTATCTGGTGGGCCCACACGGA‑‑CTCGAACCGTGGAc                                                                                                  >  7:198344/1‑144 (MQ=255)
                                                   gcgcTACTTCCAAGGTCATGCGCTCGGCGCCGCTGCGCTGGATGTACGCGCCCCTGACCATTTATTTCTTGCGGGCATGAAAAAGCCGCGCAATCGCGCGGCTTTGTCGTATCTGGTGGGCCCACACGGA‑‑CTCGAACCGTGGAc                                                                                                  <  8:198344/144‑1 (MQ=255)
                                                          ttCCAAGGTCATGCGCTCGGCGCCGCTGCGCTGGATGTACGCGCCCCTGACCATTTATTTCTTGCGGGCATGAAAAAGCCGCGCAATCGCGCGGCTTTGTCGTATCTGGTGGGCCCACACGGA‑‑CTCGAACCGTGGACCAAAGGATTAt                                                                                       >  5:154099/1‑148 (MQ=255)
                                                           tCCAAGGTCATGCGCTCGGCGCCGCTGCGCTGGATGTACGCGCCCCTGACCATTTATTTCTTGCGGGCATGAAAAAGCCGCGCAATCGCGCGGCTTTGTCGTATCTGGTGGGCCCACACGGA‑‑CTCGAACCGTGGACCAAAGGATTAt                                                                                       >  8:326870/1‑147 (MQ=255)
                                                             cAAGGTCATGCGCTCGGCGCCGCTGCGCTGGATGTACGCGCCCCTGACCATTTATTTCTTGCGGGCATGAAAAAGCCGCGCAATCGCGCGGCTTTGTCGTATCTGGTGGGCCCACACGGA‑‑CTCGAACCGTGGACCAAAGGATTATGAGt                                                                                   >  3:180685/1‑149 (MQ=255)
                                                             cAAGGTCATGCGCTCGGCGCCGCTGCGCTGGATGTACGCGCCCCTGACCATTTATTTCTTGCGGGCATGAAAAAGCCGCGCAATCGCACGGCTTTGTCGTATCTGGTGGGCCCACACGGA‑‑CTCGAACCGTGGACCAAAGGATTATGAGt                                                                                   >  5:190119/1‑149 (MQ=255)
                                                                  tCATGTGCTCGGCGCCGCTGCGCTGGATGTACGCGCCCCTGACCATTTATTTCTTGCGGGCATGAAAAAGCCGCGCAATCGCGCGGCTTTGTCGTATCTGGTGGGCCCACACGGA‑‑CTCGAACCGTGGACCAAAGGATTATGAGTCctct                                                                              >  5:254210/1‑149 (MQ=255)
                                                                    aTGCGCTCGGCGCCGCTGCGCTGGATGTACGCGCCCCTGACCATTTATTTCTTGCGGGCATGAAAAAGCCGCGCAATCGCGCGGCTTTGTCGTATCTGGTGGGCCCACACGGA‑‑CTCGAACCGTGGACCAAAGGATTATGAGTCCTCTg                                                                             <  3:222774/148‑1 (MQ=255)
                                                                      gcgcTCGGCGCCGCTGCGCTGGATGTACGAGCCCCTGACCATTTATTTCTTGCGGGCATGAAAAAGCCGCGCAATCGAGCGGCTTTGTCGAATCTGCTGGGCCCACACGGA‑‑CTAGAACCGTGGACCATAGGATTATGCGTACTCTGgtc                                                                          >  8:336195/1‑149 (MQ=255)
                                                                       cgcTCGGCGCCGCTGCGCTTGATGTACGCGCCCCTGACCATTAATTTCTTGCGGGCATGAAAAAGCCGCGCAATCGCGCGGCTTTGTCGTATCTGGTGGGCCCACACGGA‑‑CTCGAACCGTGGACCAAAGGATTATGAGTCCTCTGctct                                                                         <  8:310597/149‑1 (MQ=255)
                                                                        gcTCGGCGCCGCTGCGCTGGATGTACGCGCCCCTGACCATTTATTTCTTGCGGGCATGAAAAAGCCGCGCAATCGCGCGGCTTTGTCGTATCTGGTGGGCCCACACGGA‑‑CTCGAACCGTGGACCAAAGGATTATGAGTCCTCTGctct                                                                         >  7:89892/1‑148 (MQ=255)
                                                                         cTCGGCGCCGCTGCGCTGGATGTACGCGCCCCTGACCATTTATTTCTTGCGGGCATGAAAAAGCCGCGCAATCGCGCGGCTTTGTCGTATCTGGTGGGCCCACACGGA‑‑CTCGAACCGTGGACCAAAGGATTATGAGTCCTCTGCTCTaa                                                                       <  8:89892/149‑1 (MQ=255)
                                                                         cTCGGCGCCGCTGCGCTGGATGTACGCGCCCCTGACCATTTATTTCTTGCGGGCATGAAAAAGCCGCGCAATCGCGCGGCTTTGTCGTATCTGGTGGGCCCACACGGA‑‑CTCGAACCGTGGACCAAAGGATTATGAGTCCTCTGCTCTa                                                                        >  4:147281/1‑148 (MQ=255)
                                                                          tCGGCGCCGCTGCGCTGGATGTACGCGCCCCTGACCATTTATTTCTTGCGGGCATGAAAAAGCCGCGCGATTGCGCGGCttttt                                                                                                                                         <  8:320534/84‑3 (MQ=255)
                                                                          tCGGCGCCGCTGCGCTGGATGTACGCGCCCCTGACCATTTATTTCTTGCGGGCATGAAAAAGCCGCGCGATTGCGCGGCttttt                                                                                                                                         >  7:320534/1‑82 (MQ=255)
                                                                                 cgctgcgctgGATGTACGCGCCCCTGACCATTTATTTCTTGCGGGCATGAAAAAGCCGCGCAATCGCGCGGCTTTGTCGTATCTGGTGGGCCCACACGGA‑‑CTCGAACCGTGGACCAAAGGATTATGAGTCCTCTGCTCTAACCAACTGa                                                               <  7:326870/149‑1 (MQ=255)
                                                                                  gctgcgctgGATGTACGCGCCCCTGACCATTTATTTCTTGCGGGCATGAAAAAGCCGCGCAATCGCGCGGCTTTGTCGTATCTGGTGGGCCCACACGGA‑‑CTCGAACCGTGGACCAAAGGATTATGAGTCCTCTGCTCTAACCAACTGAg                                                              <  7:83034/149‑1 (MQ=255)
                                                                                       gctgGATGTACGCGCCCCTGACCATTTATTTCTTGCGGGCATGAAAAAGCCGCGCAATCGCGCGGCTTTGTCGTATCTGGTGGGCCCACACGGA‑‑CTCGAACCGTGGACCAAAGGATTATGAGTCCTCTGCTCTAACCAACTGAGCta                                                           >  5:190312/1‑147 (MQ=255)
                                                                                         tgGATGTACGCGCCCCTGACCATTTATTTCTTGCGGGCATGAAAAAGCCGCGCAATCGCGCGGCTTTGTCGTATCTGGTGGGCCCACACGGA‑‑CTCGAACCGTGGACCAGAGGATTATGAGTCCTCTGCTCTAAc                                                                      >  6:136834/1‑134 (MQ=255)
                                                                                         tgGATGTACGCGCCCCTGACCATTTATTTCTTGCGGGCATGAAAAAGCCGCGCAATCGCGCGGCTTTGTCGTATCTGGTGGGCCCACACGGA‑‑CTCGAACCGTGGACCAGAGGATTATGAGTCCTCTGCTCTAAc                                                                      <  5:136834/134‑1 (MQ=255)
                                                                                         tgGATGTACGCGCCCCTGACCATTTATTTCTTGCGGGCATGAAAAAGCCGCGCAATCGCGCGGCTTTGTCGTATCTGGTGGGCCCACACGGA‑‑CTCGAACCGTGGACCAAAGGATTATGAGTCCTCTGCTCTAACCAACTGAGCTATAgg                                                       <  1:158812/149‑1 (MQ=255)
                                                                                             tGTACGCGCCCCTGACCATTTATTTCTTGCGGGCATGAAAAAGCCGCGCAATCGCGCGGCTTTGTCGTATCTGGTGGGCCCACACGGA‑‑CTCGAACCGTGGACCAAAGGATTa                                                                                        <  8:289557/112‑1 (MQ=255)
                                                                                             tGTACGCGCCCCTGACCATTTATTTCTTGCGGGCATGAAAAAGCCGCGCAATCGCGCGGCTTTGTCGTATCTGGTGGGCCCACACGGA‑‑CTCGAACCGTGGACCAAAGGATTa                                                                                        >  7:289557/1‑112 (MQ=255)
                                                                                                 cgcgCCCCTGACCATTTATTTCTTGCGGGCATGAAAAAGCCGCGCAATCGCGCGGCTTTGTCGTATCTGGTGGGCCCACACGGA‑‑CTCGAACCGTGGACCAAAGGATTATGAGTCCTCTGCTCTAACCAACTGAGCTATAGGCCCTCAGt                                               <  6:190312/149‑1 (MQ=255)
                                                                                                  gcgcCCCTGACCATTTATTTCTTGCGGGCATGAAAAAGCCGCGCAATCGCGCGGCTTTGTCGTATCTCGTGGGCCCACACGGA‑‑CTCGAACCGTGGCCCATAGGCTTATGCGTCCTCTGCTCTTACCCACTGAGCTATAGGCCGTCAGt                                               >  8:343007/1‑148 (MQ=255)
                                                                                                   cgcCCCTGACCATTTATTTCTTGCGGGCATGAAAAAGCCGCGCAATCGCGCGGCTTTGTCGTATCTGGTGGGCCCACACGGACTCTCGAACCGTGGACCAAAGGATTATGAGTCCTCTGCTCTAACCAACTGAGCTATAGGCCCTCAGt                                               <  1:111783/149‑1 (MQ=255)
                                                                                                      cccTGACCATTTATTTCTTGCGGGCATGAAAAAGCCGCGCAATCGCGCGGCTTTGTCGTATCTGGTGGGCCCACACGGA‑‑CTCGAACCGTGGACCAAAGGATTATGAGTCCTCTGCTCTAACCAACTGAGCTATAGGCCCTCAGTAAGTg                                          >  5:353292/1‑149 (MQ=255)
                                                                                                       ccTGACCATTTATTTCTTGCGGGCATGAAAAAGCCGCGCAATCGCGCGGCTTTGTCGTATCTGGTGGGCCCACACGGA‑‑CTCGAACCGTGGACCAAAGGATTATGAGTCCTCTGCTCTAACCAACTGAGCTATAGGCCCTCAGTAAGTgg                                         >  4:371396/1‑149 (MQ=255)
                                                                                                            ccATTTATTTCTTGCGGGCATGAAAAAGCCGAGCAATCGCGCGGCTTTGTCGTTTCTGGTGGGCCCAAACGGA‑‑CTCGAACCGTGGACCAAAGGATTATTAGTCATATGATCTAAACAAATGAGCTATAGGCCATAAGTAAGTGGGCgca                                    >  5:89831/1‑147 (MQ=255)
                                                                                                             cATTTATTTCTTGCGGGCATGAAAAAGCCGCGCAATCGCGCGGCTTTGTCGTATCTGGTGGGCCCACACGGA‑‑CTCGAACCGTGGACCAAAGGATTATGAGTCCTCTGCTCTAACCAACTGAGCTATAGGCCCTCAGTAAGTGGGCGGAt                                   <  6:154099/149‑1 (MQ=255)
                                                                                                              atttatatCTTGCGGGCATGAAAAAGCCGCGCAATCGCGCGGCTTTGTCGTATCTGGTGGGACCAAACGGA‑‑ATCGAACCGTGGACCAAAGGATTATGAGTCCTCTCGACTACACAACTGAGCTATACGCCCACATGAAGTGGGCGGAt                                   >  8:304881/1‑148 (MQ=255)
                                                                                                                  atttCTTGCGGGCATGAAAAAGCCGCGCAATCGCGCGGCTTTGTCGTATCTGGTGGGCCCACACGGA‑‑CTCGAACCGTGGACCAAAGGATTATGAGTCCTCTGCTCTAACCAACTGAGCTATAGGCCCTCAGTAAGTGGGCGGATTATaa                              >  1:251614/1‑149 (MQ=255)
                                                                                                                   tttCTTGCGGGCATGACAAAGCCGCGCAATCGCGCGGGTTTGTCGTATCTGGGGGGCCCACACGGA‑‑CTCGAACCGTGGACCAAAGGATTATGAGTCCTCTGCTCTAACCAACTGAGCTATAGGCCCTCAGTAAGTGGGCGGATTATaa                              <  6:295230/148‑1 (MQ=255)
                                                                                                                    ttCTTGCGGGCATGAAAAAGCCGCGCAATCGCGCGGCTTTGTCGTATCTGGTGGGCCCACACGGA‑‑CTCGAACCGTGGACCAAAGGATTATGAGTCCTCTGCTCTAACCAACTGAGCTATAGGCCCTCAGTAAGTGGGCGGATTATAAc                             >  7:345376/1‑148 (MQ=255)
                                                                                                                    ttCTTGCGGGCATGAAAAAGCCGCGCAATCGCGCGGCTTTGTCGTATCTGGTGGGCCCACACGGA‑‑CTCGAACCGTGGACCAAAGGATTATGAGTCCTCTGCTCTAACAAACTGAGCTATAGGCCCTCAGTAAGTGGGCGGATTATAAc                             >  1:149085/1‑148 (MQ=255)
                                                                                                                     tCTTGCGGGCATGAAAAAGCCGCGCAATCGCGCGGCTTTGTCGTATCTGGTGGGCCCACACGGA‑‑CTCGAACCGTGGACCAAAGGATTATGAGTCCTCTGCTCTAACCAACTGAGCTATAGGCCCTCAGTAAGTGGGCGGATTATaa                              >  3:330529/1‑146 (MQ=255)
                                                                                                                     tCTTGCGGGCATGAAAAAGCCGCGCAATCGCGCGGCTTTGTCGTATCTGGTGGGCCCACACGGA‑‑CTCGAACCGTGGACCAAAGGATTATGAGTCCTCTGCTCTAACCAACTGAGCTATAGGCCCTCAGTAAGTGGGCGGATTATaa                              <  4:330529/146‑1 (MQ=255)
                                                                                                                       ttGCGGGCATGAAAAAGCCGCGCAATCGCGCGGCTTTGTCGTATCTGGTGGGCCCACACGGA‑‑CTCGAACCGTGGACCAAAGGATTATGAGTCCTCTGCTCTAACCAACTGAGCTATAGGCCCTCAGTAAGTGGGCGGATTATAACGAgg                         >  7:295753/1‑149 (MQ=255)
                                                                                                                           gggCATGAAAAAGCCGCGCAATCGCGCGGCTTTGTCGTATCTGGTGGGCCCACACGGA‑‑CTCGAACCGTGGACCAAAGGATTATGAGTCGTCTGCTCTAAGCAACTGAGCTAAAGGCCCTCAGTAAGTGGGCGGATTATAACGAGGGttt                     >  7:71796/1‑149 (MQ=255)
                                                                                                                                 gAAAAAGCCGCGCAATCGCGCGGCTTTGTCGTATCTGGTGGGCCCACACGGA‑‑CTCGAACCGTGGACCAAAGGATTATGAGTCCTCTGCTCTAACCAACTGAGCTATAGGCCCTCAGTAAGTGGGCGGATTATAACGAGGGTTTCGCAcc               >  8:222029/1‑149 (MQ=255)
                                                                                                                                   aaaaGCCGCGCAATCGCGCGGCTTTGTCGTATCTGGTGGGCCCACACGGA‑‑CTCGAACCCTGGACCAAAGGATTATGAGTCCGCGGCTCTAACCACCTGAGCTATAGGCCCTGAGTAAGTGGGCGGATTATAACGAGGGTTTCGCACCtg             >  4:269672/1‑149 (MQ=255)
                                                                                                                                    aaaGCCGCGCAATCGCGCGGCTTTGTCGTATCTGGTGGGCCCACACGGA‑‑CTCGAACCGTGGACCAAAGGATTATGAGTCCTCTGCTCTAACCAACTGAGCTATAGGCCCTCAGTAAGTGGGCGGATTATAACGAGGGTTTCGCACCtgt            <  6:79017/149‑1 (MQ=255)
                                                                                                                                       gCCGCGCGATTGCGCGGCTTTGTCGTATCTGGTGGGCCCACACGGA‑‑CTCGAACCGTGGACCAAAGGATTATGAGTCCTCTGCTCTAACCAACTGAGCTATAGGCCCTCAGTAAGTGGGCGGATTATAACGAGGGTTTCGCACCTGTGc          >  2:67222/1‑148 (MQ=255)
                                                                                                                                           cgcAATCGCGCGGCTTTGTCGTATCTGGTGGGCCCACACGGA‑‑CTCGAACCGTGGACCAAAGGATTATGAGTCCTCTGCTCTAACCAACTGAGCTATAGGCCCTCAGTAAGTGGGCGGATTATAACGAGGGTTTCGCACCTGTGCCATCt     <  3:147281/149‑1 (MQ=255)
                                                                                                                                              aaTCGCGCGGCTTTGTCGTATCTGGTGGGCCCACACGGA‑‑CTCGAACCGTGGACCAAAGGATTATGAGTCCTCTGCTCTAACCAACTGAGCTATAGGCCCTCAGTAAGTGGGCGGATTATAACGAGGGTTTCGCACCTGTGCCATCTGaa  <  8:217164/149‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                                      
GCTGGGTGACCGGCTTCTTGGTGGGGCTGATTTCGCTGAGATACCGTTTTAGCGCTACTTCCAAGGTCATGCGCTCGGCGCCGCTGCGCTGGATGTACGCGCCCCTGACCATTTATTTCTTGCGGGCATGAAAAAGCCGCGCAATCGCGCGGCTTTGTCGTATCTGGTGGGCCCACACGGA‑‑CTCGAACCGTGGACCAAAGGATTATGAGTCCTCTGCTCTAACCAACTGAGCTATAGGCCCTCAGTAAGTGGGCGGATTATAACGAGGGTTTCGCACCTGTGCCATCTGAA  >  NC_002947/466956‑467246

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