Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,962,562 T→C 6.5% N68S (AAC→AGC) ‡ ppdC ← putative prepilin peptidase‑dependent protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,962,5620TC 6.5% 348.4 / 11.2 107N68S (AAC→AGC) ‡ppdCputative prepilin peptidase‑dependent protein
Reads supporting (aligned to +/- strand):  ref base T (51/49);  new base C (2/5);  total (53/54)
Fisher's exact test for biased strand distribution p-value = 4.37e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.85e-01

Only 100 of 109 total aligned reads displayed.

TAGACCCTTAACATAGCGGCTCCTGACTACTGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCGCCAGAGCTGTTGGTACTGGTTTCTGCTGGCAAGACTGTTCATTAATGTTCGCTGGAAACCCGATAACGCAGTGACAATCATCACCATCAACACCATCG  >  NC_000913/2962415‑2962709
                                                                                                                                                   |                                                                                                                                                   
tAGACCCTTAACATAGCGGCTCCTGACTACTGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGtt                                                                                                                                                    <  2:128463/149‑1 (MQ=255)
tAGACCCTTAACATAGCGGCTCCTGACTACTGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGtt                                                                                                                                                    <  4:489183/149‑1 (MQ=255)
   aCCCTTAACATAGCGGCTCCTGACTACTGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAACTGGc                                                                                                                                                 <  1:366766/149‑1 (MQ=255)
     ccTTAACATAGCGGCTCCTGACTACTGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACATGCCAGTTGGCAg                                                                                                                                               >  4:286678/1‑149 (MQ=255)
       ttAACATAGCGGCTCCTGACTACTGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGt                                                                                                                                             >  4:364645/1‑149 (MQ=255)
       ttAACATAGCGGCTCCTGACTACTGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAACTGGCAGGt                                                                                                                                             <  3:222701/149‑1 (MQ=255)
        tAACATAGCGGCTCCTGACTACTGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTg                                                                                                                                            <  1:189277/149‑1 (MQ=255)
        tAACATAGCGGCTCCTGACTACTGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAACTGGCAGGtc                                                                                                                                            <  2:599776/149‑2 (MQ=255)
         aaCATAGCGGCTCCTGACTACTGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTgg                                                                                                                                           <  3:122373/149‑1 (MQ=255)
          aCATAGCGGCTCCTGACTACTGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGc                                                                                                                                          <  1:5545/149‑1 (MQ=255)
            aTAGCGGCTCCTGACTACTGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGa                                                                                                                                        >  2:195360/1‑149 (MQ=255)
             tAGCGGCTCCTGACTACTGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGaa                                                                                                                                       >  4:147852/1‑149 (MQ=255)
             tAGCGGCTCCTGACTACTGACGATTCGGGCAATGCAGGCGGGACATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGaa                                                                                                                                       >  3:8280/1‑149 (MQ=255)
                  gCTCCTGACTACTGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATcgc                                                                                                                                  <  4:54476/149‑1 (MQ=255)
                  gCTCCTGACTACTGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCGAAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATcgc                                                                                                                                  <  3:250201/149‑1 (MQ=255)
                   cTCCTGACTACTGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATcgcg                                                                                                                                 <  3:318393/149‑1 (MQ=255)
                        gACTACTGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAg                                                                                                                            >  1:565516/1‑149 (MQ=255)
                          ctactGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGtt                                                                                                                          <  1:65439/149‑1 (MQ=255)
                           tactGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTg                                                                                                                         >  4:162231/1‑149 (MQ=255)
                           tactGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTg                                                                                                                         >  2:289436/1‑149 (MQ=255)
                              tGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGc                                                                                                                        >  1:241976/1‑147 (MQ=255)
                              tGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGc                                                                                                                        <  2:241976/147‑1 (MQ=255)
                              tGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGt                                                                                                                      >  2:386844/1‑149 (MQ=255)
                              agACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGt                                                                                                                      >  3:112829/2‑149 (MQ=255)
                                  gATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTtgc                                                                                                                  >  4:519590/1‑149 (MQ=255)
                                   aTTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTtgct                                                                                                                 <  1:559465/149‑1 (MQ=255)
                                    ttCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTtgctg                                                                                                                >  2:541282/1‑149 (MQ=255)
                                        ggCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTgccag                                                                                                            <  3:162231/149‑1 (MQ=255)
                                           aaTGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTgccagcca                                                                                                         >  1:561827/1‑149 (MQ=255)
                                            aTGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTtgct                                                                                                                 >  1:310666/1‑140 (MQ=255)
                                            aTGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTtgct                                                                                                                 <  2:310666/140‑1 (MQ=255)
                                               cAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATgcc                                                                                                     <  1:136181/149‑1 (MQ=255)
                                                      gTCATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCGCCagag                                                                                              <  2:197278/149‑1 (MQ=255)
                                                         aTCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCGCCAGAGCTg                                                                                           >  2:369734/1‑149 (MQ=255)
                                                         aTCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCGCCAGAGCTg                                                                                           >  4:177112/1‑149 (MQ=255)
                                                            tcGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCGCCAGAGCTGTTg                                                                                        >  2:669148/1‑149 (MQ=255)
                                                            tcGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCGCCAGAGCTGTTg                                                                                        >  2:644934/1‑149 (MQ=255)
                                                             tgCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCGCCAGAGCTGTTgg                                                                                       <  2:157182/148‑1 (MQ=255)
                                                               ccTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCGCCAGAGCTGTTGGTa                                                                                     <  1:203297/149‑1 (MQ=255)
                                                                 ttCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTtgctgc                                                                                                               >  4:353278/1‑121 (MQ=255)
                                                                 ttCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCGCCAGAGCTGTTGGTACt                                                                                   >  3:371662/1‑149 (MQ=255)
                                                                 ttCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCGCCAGAGCTGTTGGTACt                                                                                   >  3:527156/1‑149 (MQ=255)
                                                                 ttCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCGCCAGAGCTGTTGGTACt                                                                                   >  2:615021/1‑149 (MQ=255)
                                                                 ttCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCGCCAGAGCTGTTGGTACt                                                                                   >  1:550525/1‑149 (MQ=255)
                                                                      tGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCGCCAGAGCTGTTGGTACTGGttt                                                                              <  1:289436/149‑1 (MQ=255)
                                                                      tGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCGCCAGAGCTGTTGGTACTGGttt                                                                              <  1:369734/149‑1 (MQ=255)
                                                                       gCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCGCCAGAGCTGTTGGTACTGGTTTc                                                                             <  2:13321/149‑1 (MQ=255)
                                                                             gggTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCGCCAGAGCTGTTGGTACTGGTTTCTGCTgg                                                                       <  4:445391/149‑1 (MQ=255)
                                                                                   aaCTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCGCCAGAGCTGTTGGTACTGGTTTCTGCTGGCAAGAc                                                                 <  4:205980/149‑1 (MQ=255)
                                                                                    aCTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATc                                                                                                                                    >  2:352763/1‑81 (MQ=255)
                                                                                        gCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCGCCAGAGCTGTTGGTACTgg                                                                                 >  2:472095/1‑128 (MQ=255)
                                                                                        gCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCGCCAGAGCTGTTGGTACTgg                                                                                 <  1:472095/128‑1 (MQ=255)
                                                                                         cGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCGCCAGAGCTGTTGGTACTGGTTTCTGCTGGCAAGACTGTTCa                                                           >  1:266947/1‑149 (MQ=255)
                                                                                           tAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATcgc                                                                                                                                  >  4:251824/1‑76 (MQ=255)
                                                                                           tAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATcgc                                                                                                                                  <  3:251824/76‑1 (MQ=255)
                                                                                           tAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCGCCAGAGCTGTTGGTACTGGTTTCTGCTGGCAAGACTGTTCAtt                                                         <  3:529672/149‑1 (MQ=255)
                                                                                           tAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCGCCAGAGCTGTTGGTACTGGTTTCTGCTGGCAAGACTGTTCAtt                                                         >  3:465978/1‑149 (MQ=255)
                                                                                             aCGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCGCCAGAGCTGTTGGTACTGGTTTCTGCTGGCAAGACTGTTCATTaa                                                       <  3:177219/149‑1 (MQ=255)
                                                                                              cGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCGCCAGAGCTGTTGGTACTGGTTTCTGCTGGCAAGACTGTTCATTAAt                                                      <  2:690107/149‑1 (MQ=255)
                                                                                                 tGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCGCCAGAGCTGTTGGTACTGGTTTCTGCTGGCAAGACTGTTCATTAATGtt                                                   >  1:648553/1‑149 (MQ=255)
                                                                                                  gATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCGCCAGAGCTGTTg                                                                                        <  4:514155/111‑1 (MQ=255)
                                                                                                  gATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCGCCAGAGCTGTTg                                                                                        >  3:514155/1‑111 (MQ=255)
                                                                                                    tGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCGCCAGAGCTGTTGGTACTGGTTTCTGCTGGCAAGACTGTTCATTAATGTTCGc                                                <  2:298946/149‑1 (MQ=255)
                                                                                                       tGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCGCCAGAGCTGTTGGTACTGGTTTCTGCTGGCAAGACTGt                                                              >  1:292476/1‑132 (MQ=255)
                                                                                                       tGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCGCCAGAGCTGTTGGTACTGGTTTCTGCTGGCAAGACTGt                                                              <  2:292476/132‑1 (MQ=255)
                                                                                                       tGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCGCCAGAGCTGTTGGTACTGGTTTCTGCTGGCAAGACTGTTCATTAATGTTCGCTgg                                             >  1:650345/1‑149 (MQ=255)
                                                                                                        gACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCGCCAGAGCTGTTGGTACTGGTTTCTGCTGGCAAGACTGTTCATTAATGTTCGCTGGa                                            >  4:513286/1‑149 (MQ=255)
                                                                                                         acacATCCCGCCTGCGATGTCTGCATTTGGTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCGCCAGAGCTGTTGGTACTGGTTTCTGCTGGCAAGACTGTTCATTAATGTTCGCTGGaa                                           >  4:635560/1‑149 (MQ=255)
                                                                                                         acacATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCGCCAGAGCTGTTGGTACTGGTTTCTGCTGGCAAGACTGTTCATTAATGTTCGCTGGaa                                           >  2:79359/1‑149 (MQ=255)
                                                                                                           acaTCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCGCCAGAGCTGTTGGTACTGGTTTCTGCTGGCAAGACTGTTCATTAATGTTCGCTGGAAAc                                         >  4:170564/1‑149 (MQ=255)
                                                                                                             aTCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCGCCAGAGCTGTTGGTACTGGTTTCTGCTGGCAAGACTGTTCATTAATGTTCGCTGGAAAccc                                       <  3:286678/149‑1 (MQ=255)
                                                                                                              tCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCGCCAGAGCTGTTGGTACTGGTTTCTGCTGGCAAGACTGTTCATTAATGTTCGCTGGAAACCCg                                      >  1:591626/1‑149 (MQ=255)
                                                                                                               cccGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCGCCAGAGCTGTTGGTACTGGTTTCTGCTGGCAAGACTGTTCATTAATGTTCGCTGGAAACCCGa                                     <  3:437913/149‑1 (MQ=255)
                                                                                                               cccGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCGCCAGAGCTGTTGGTACTGGTTTCTGCTGGCAAGACTGTTCATTAATGTTCGCTGGAAACCCGa                                     <  1:79359/149‑1 (MQ=255)
                                                                                                                ccGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCGCCAGAGCTGTTGGTACTGGTTTCTGCTGGCAAGACTGTTCATTAATGTTCGCTGGAAACCCGAt                                    <  2:650345/149‑1 (MQ=255)
                                                                                                                      gCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCGCCAGAGCTGTTGGTACTGGTTTCTGCTGGCAAGACTGTTCATTAATGTTCGCTGGAAACCCGATAACGCa                              <  2:612180/149‑1 (MQ=255)
                                                                                                                      gCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCGCCAGAGCTGTTGGTACTGGTTTCTGCTGGCAAGACTGTTCATTAATGTTCGCTGGAAACCCGATAACGCa                              <  1:386844/149‑1 (MQ=255)
                                                                                                                        gATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCGCCAGAGCTGTTGGTACTGGTTTCTGCTGGCAAGACTGTTCATTAATGTTCGCTGGAAACCCGATAACGCAGt                            <  3:177112/149‑1 (MQ=255)
                                                                                                                        gATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCGCCAGAGCTGTTGGTACTGGTTTCTGCTGGCAAGACTGTTCATTAATGTTCGCTGGAAACCCGATAACGCAGt                            <  1:195360/149‑1 (MQ=255)
                                                                                                                            tCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCGCCAGAGCTGTTGGTACTGGTTTCTGCTGGCAAGACTGTTc                                                            >  4:611537/1‑113 (MQ=255)
                                                                                                                            tCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCGCCAGAGCTGTTGGTACTGGTTTCTGCTGGCAAGACTGTTc                                                            <  3:611537/113‑1 (MQ=255)
                                                                                                                              tGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCGCCAGAGCTGTTGGTACTGGTTTCTGCTGGCAAGACTGtt                                                             >  4:534072/1‑110 (MQ=255)
                                                                                                                              tGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCGCCAGAGCTGTTGGTACTGGTTTCTGCTGGCAAGACTGTTCATTAATGTTCGCTGGAAACCCGATa                                   >  1:292042/1‑136 (MQ=255)
                                                                                                                              tGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCGCCAGAGCTGTTGGTACTGGTTTCTGCTGGCAAGACTGTTCATTAATGTTCGCTGGAAACCCGATAACGCAGTGACAat                      >  3:156050/1‑149 (MQ=255)
                                                                                                                              tGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCGCCAGAGCTGTTGGTACTGGTTACTGCTGGCAAGACTGTTCATTAATGTTCGCTGGAAACCCGATa                                   <  2:292042/136‑1 (MQ=255)
                                                                                                                              tGCATTCGCTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCGCCAGAGCTGTTGGTACTGGTTTCTGCTGGCAAGACTGtt                                                             <  3:534072/110‑1 (MQ=255)
                                                                                                                               gCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCGCCAGAGCTGTTGGTACTGGTTTCTGCTGGCAAGACTGTTCATTAATGTTCGCTGGAAACCCGATAACGCAGTGACAatc                     <  4:238639/149‑1 (MQ=255)
                                                                                                                                cATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCGCCAGAGCTGTTGGTACTGGTTTCTGCTGGCAAGACTGTTCATTAATGTTCGCTGGAAACCCGATAACGCAg                             <  2:129623/140‑1 (MQ=255)
                                                                                                                                cATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCGCCAGAGCTGTTGGTACTGGTTTCTGCTGGCAAGACTGTTCATTAATGTTCGCTGGAAACCCGATAACGCAg                             >  1:129623/1‑140 (MQ=255)
                                                                                                                                cATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCGCCAGAGCTGTTGGTACTGGTTTCTGCTGGCAAGACTGTTCATTAATGTTCGCTGGAAACCCGATAACGCAGTGACAatca                    >  3:37774/1‑149 (MQ=255)
                                                                                                                                 aTTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCGCCAGAGCTGTTGGTACTGGTTTCTGCTGGCAAGACTGTTCATTAATGTTCGCTGGAAACCCGATAACGCAGTGACAatcat                   >  3:687737/1‑149 (MQ=255)
                                                                                                                                      gTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCGCCAGAGCTGTTGGTACTGGTTTCTGCTGGCAAGACTGTTCATTAATGTTCGCTGGAAACCCGATAACGCAGTGACAat                      >  4:206691/1‑141 (MQ=255)
                                                                                                                                      gTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCGCCAGAGCTGTTGGTACTGGTTTCTGCTGGCAAGACTGTTCATTAATGTTCGCTGGAAACCCGATAACGCAGTGACAat                      <  3:206691/141‑1 (MQ=255)
                                                                                                                                         caCCTGCCAACTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCGCCAGAGCTGTTGGTACTGGTTTCTGCTGGCAAGACTGTTCATTAATGTTCGCt                                               >  1:257835/2‑113 (MQ=255)
                                                                                                                                         caCCTGCCAACTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCGCCAGAGCTGTTGGTACTGGTTTCTGCTGGCAAGACTGTTCATTAATGTTCGCTGGAAACCCGATAACGCAGt                            <  2:473451/131‑1 (MQ=255)
                                                                                                                                         caCCTGCCAACTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCGCCAGAGCTGTTGGTACTGGTTTCTGCTGGCAAGACTGTTCATTAATGTTCGCTGGAAACCCGATAACGCAGt                            >  1:473451/2‑132 (MQ=255)
                                                                                                                                           ccTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCGCCAGAGCTGTTGGTACTGGTTTCTGCTGGCAAGACTGTTCATTAATGTTCGCTGGAAACCCGATAACGCAGTGACAATCATCACCATCAac         >  4:461743/1‑149 (MQ=255)
                                                                                                                                           ccTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCGCCAGAGCTGTTGGTACTGGTTTCTGCTGGCAAGACTGTTCATTAATGTTCGCTGGAAACCCGATAACGCAGTGACAATCATCACCATCAac         <  3:170564/149‑1 (MQ=255)
                                                                                                                                                 aGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCGCCAGAGCTGTTGGTACTGGTTTCTGCTGGCAAGACTGTTCATTAATGTTCGCTGGAAACCCGATAACGCAGTGACAATCATCACCATCAACACCATc   <  4:527156/149‑1 (MQ=255)
                                                                                                                                                  gTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCGCCAGAGCTGTTGGTACTGGTTTCTGCTGGCAAGACTGTTCATTAATGTTCGCTGGAAACCCGATAACGCAGTGACAATCATCACCATCAACACCATCg  >  2:613376/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
TAGACCCTTAACATAGCGGCTCCTGACTACTGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCGCCAGAGCTGTTGGTACTGGTTTCTGCTGGCAAGACTGTTCATTAATGTTCGCTGGAAACCCGATAACGCAGTGACAATCATCACCATCAACACCATCG  >  NC_000913/2962415‑2962709

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: