Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,616,399 T→G 5.7% V73G (GTG→GGG)  nikC → nickel ABC transporter permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,616,3990TG 5.7% 368.3 / 13.3 106V73G (GTG→GGG) nikCnickel ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base T (49/51);  new base G (2/4);  total (51/55)
Fisher's exact test for biased strand distribution p-value = 6.80e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.36e-01

Only 100 of 110 total aligned reads displayed.

TGCGCTCACCAGCCAGTGGTGGCTGCCGTATGACCCACAGGCGATTGATTTGCCGTCGCGCCTGCTTTCGCCGGATGCGCAGCACTGGCTGGGCACCGATCACTTAGGTCGCGATATTTTCTCGCGGCTGATGGCAGCGACCCGCGTGTCGCTCGGTTCGGTAATGGCCTGCCTGCTGCTGGTGCTGACATTAGGGCTGGTTATTGGCGGCAGCGCCGGGTTGATTGGCGGGCGCGTTGATCAGGCCACCATGCGCGTCGCCGATATGTTTATGACCTTCCCGACCTCGATTCTG  >  NC_000913/3616253‑3616547
                                                                                                                                                  |                                                                                                                                                    
tGCGCTCACCAGCCAGTGGTGGCTGCCGTATGACCCACAGGCGATTGATTTGCCGTCGCGCCTGCTTTCGCCGGATGCGCAGCACTGGCTGGGCACCGATCACTTAGGTCGCGATATTTTCTCGCGGCTGATGGCAGCGACCCGCgtgt                                                                                                                                                    <  4:838583/149‑1 (MQ=255)
 gcgcTCACCAGCCAGTGGTGGCTGCCGTATGACCCACAGGCGATTGATTTGCCGTCGCGCCTGCTTTCGCCGGATGCGCAGCACTGGCTGGGCACCGATCACTTAGGTCGCGATATTTTCTCGCGGCTGATGGCAGCGACCCGCGTGTc                                                                                                                                                   <  3:168195/149‑1 (MQ=255)
    cTCACCAGCCAGTGGTGGCTGCCGTATGACCCACAGGCGATTGATTTGCCGTCGCGCCTGCTTTCGCCGGATGCGCAGCACTGGCTGGGCACCGATCACTTAGGTCGCGATATTTTCTCGCGGCTGATGGCAGCGACCCGCGTGTCGc                                                                                                                                                 >  1:850490/1‑148 (MQ=255)
      cACCAGCCAGTGGTGGCTGCCGTATGACCCACAGGCGATTGATTTGCCGTCGCGACTGCTTTCGCCGGATGCGCAGCACTGGCTGGGCACCGATCACTTAGGTCGCGATATTTTCTCGCGGCTGATGGCAGCGACCCGCGTGTCGCtc                                                                                                                                               <  1:693810/148‑1 (MQ=255)
       aCCAGCCAGTGGTGGCTGCCGTATGACCCACAGGCGATTGATTTGCCGTCGCGCCTGCTTTCGCCGGATGCGCAGCACTGGCTGGGCACCGATCACTTAGGTCGCGATATTTTCTCGCGGCTGATGGCAGCGACACGCGGGTCGCtgc                                                                                                                                              <  2:707164/148‑3 (MQ=255)
          agccagTGGTGGCTGCCGTATGACCCACAGGCGATTGATTTGCCGTCGCGCCTGCTTTCGCCGGATGCGCAGCACTGGCTGGGCACCGATCACTTAGGTCGCGATATTTTCTCGCGGCTGATGGCAGCGACACGCGGGTCGCtgc                                                                                                                                              <  4:354063/145‑3 (MQ=255)
           gccagTGGTGGCTGCCGTATGACCCACAGGCGATTGATTTGCCGTCGCGCCTGCTTTCGCCGGATGCGCAGCACTGGCTGGGCACCGATCACTTAGGTCGCGATATTTTCTCGCGGCTGATGGCAGCGACCCGCGTGTCGCtcggttc                                                                                                                                          >  1:825264/1‑148 (MQ=255)
                tggtggCTGCCGTATGACCCACAGGCGATTGATTTGCCGTCGCGCCTGCTTTCGCCGGATGCGCAGCACTGGCTGGGCACCGATCACTTAGGTCGCGATATTTTCTCGCGGCTGATGGCAGCGACCCGCGTGTCGCTCGGTTCGGTaa                                                                                                                                     <  1:530893/148‑1 (MQ=255)
                 ggtggCTGTCGTATGACCCACAGGCGATTGATTTGCCGTCGCGTCTGCTTTCGCCGGATGCGCAGCACTGGCTGGGCAGCGATCGCTTATGTCGCGATATTTTCACGCGGCTGATGGCAGCGACCCGCGTGTCGCTCGGTTCGGTAATg                                                                                                                                   <  4:565122/149‑1 (MQ=255)
                 ggtggCTGCCGTATGACCCACAGGCGATTGATTTGCCGTCGCGCCTGCTTTCGCCGGATGCGCAGCACTGGCTGGGCACCGATCACTTAGGTCGCGATATTTTCTCGCGGCTGATGGCAGCGACCCGCGTGTCGCtcggt                                                                                                                                            >  1:416566/1‑140 (MQ=255)
                 ggtggCTGCCGTATGACCCACAGGCGATTGATTTGCCGTCGCGCCTGCTTTCGCCGGATGCGCAGCACTGGCTGGGCACCGATCACTTAGGTCGCGATATTTTCTCGCGGCTGATGGCAGCGACCCGCGTGTCGCtcggt                                                                                                                                            <  2:416566/140‑1 (MQ=255)
                 ggtggCTGCCGTATGACCCACAGGCGATTGATTTGCCGTCGCGCCTGCTTTCGCCGGATGCGCAGCACTGGCTGGGCACCGATCACTTAGGTCGCGATATTTTCTCGCGGCTGATGGCAGCGACCCGCGTGTCGCTCGGTTCGGTAATg                                                                                                                                   <  4:533162/149‑1 (MQ=255)
                   tggCTGCCGTATGACCCACAGGCGATTGATTTGCCGTCGCGCCTGCTTTCGCCGGATGCGCAGCACTGGCTGGGCACCGATCACTTAGGTCGCGATATTTTCTCGCGGCTGATGGCAGCGACCCGCGTGTCGCTCGGTTCGGTAATgg                                                                                                                                  >  1:290784/1‑148 (MQ=255)
                    ggCTGCCGTATGACCCACAGGCGATTGATTTGCCGTCGCGCCTGCTTTCGCCGGATGCGCAGCACTGGCTGGGCACCGATCACTTAGGTCGCGATATTTTCTCGCGGCTGATGGCAGCGACCCGCGTGTCGCTCGGTTCGGTAATGgcc                                                                                                                                >  4:610937/1‑149 (MQ=255)
                        gCCGTATGACCCACAGGCGATTGATTTGCCGTCGCGCCTGCTTTCGCCGGATGCGCAGCACTGGCTGGGCACCGATCACTTAGGTCGCGATATTTTCTCGCGGCTGATGGCAGCGACCCGCGTGTCGCTCGGTTCGGTAATGgcctgc                                                                                                                             >  1:197234/1‑148 (MQ=255)
                        gCCGTATGACCCACAGGCGATTGATTTGCCGTCGCGCCTGCTTTCGCCGGATGCGCAGCACTGGCTGGGCACCGATCACTTAGGTCGCGATATTTTCTCGCGGCTGATGGCAGCGACCCGCGTGTCGCTCGGTTCGGTAATGgcctgc                                                                                                                             >  2:734451/1‑148 (MQ=255)
                        gCCGTATGACCCACAGGCGATTGATTTGCCGTCGCGCCTGCTTTCGCCGGATGCGCAGCACTGGCTGGGCACCGATCACTTAGGTCGCGATATTTTCTCGCGGCTGATGGCAGCGACCCGCGTGTCGCTCGGTTCGGTAATGgcctgc                                                                                                                             >  1:808583/1‑148 (MQ=255)
                         ccGTATGACCCACAGGCGATTGATTTGCCGTCGCGCCTGCTTTCGCCGGATGCGCAGCACTGGCTGGGCACCGATCACTTAGGTCGCGATATTTTCTCGCGGCTGATGGCAGCGACCCGCGTGTCGCTCGGTTCGGTAATGgcctgcct                                                                                                                           >  4:126694/1‑149 (MQ=255)
                           gTATGACCCACAGGCGATTGATTTGCCGTCGCGCCTGCTTTCGCCGGATGCGCAGCACTGGCTGGGCACCGATCACTTAGGTCGCGATATTTTCTCGCGGCTGATGGCAGCGACCCGCGTGTCGCTCGGTTCGGTAATGGcctgcctg                                                                                                                          <  1:544957/148‑1 (MQ=255)
                            tATGACCCACAGGCGATTGATTTGCCGTCGCGCCTGCTTTCGCCGGATGCGCAGCACTGGCTGGGCACCGATCACTTAGGTCGCGATATTTTCTCGCGGCTGATGGCAGCGACCCGCGTGTCGCTCGGTTCGGTAATGGCCTGCctgct                                                                                                                        <  3:927810/149‑1 (MQ=255)
                              tGACCCACAGGCGATTGATTTGCCGTCGCGCCTGCTTTCGCCGGATGCGCAGCACTGGCTGGGCACCGATCACTTAGGTCGCGATATTTTCTCGCGGCTGATGGCAGCGACCCGCGTGTCGCTCGGTTCGGTAATGGCCTGCctgctg                                                                                                                       <  1:598283/148‑1 (MQ=255)
                                 cccACAGGCGATTGATTTGCCGTCGCGCCTGCTTTCGCCGGATGCGCAGCACTGGCTGGGCACCGATCACTTAGGTCGCGATATTTTCTCGCGGCTGATGGCAGCGACCCGCGTGTCGCTCGGTTCGGTAATgg                                                                                                                                  >  1:353122/1‑134 (MQ=255)
                                 cccACAGGCGATTGATTTGCCGTCGCGCCTGCTTTCGCCGGATGCGCAGCACTGGCTGGGCACCGATCACTTAGGTCGCGATATTTTCTCGCGGCTGATGGCAGCGACCCGCGTGTCGCTCGGTTCGGTAATgg                                                                                                                                  <  2:353122/134‑1 (MQ=255)
                                 cccACAGGCGATTGATTTGCCGTCGCGCCTGCTTTCGCCGGATGCGCAGCACTGGCTGGGCACCGATCACTTAGGTCGCGATATTTTCTCGCGGCTGATGGCAGCGACCCGCGTGTCGCTCGGTTCGGTAATGGCCTGCctgctgctg                                                                                                                    <  1:809093/148‑1 (MQ=255)
                                 cccACAGGCGATTGATTTGCCGTCGCGCCTGCTTTCGCCGGATGCGCAGCACTGGCTGGGCACCGATCACTTAGGTCGCGATATTTTCTCGCGGCTGATGGCAGCGACCCGCGTGTCGCTCGGTTCGGTAATGGCCTGCctgctgctg                                                                                                                    >  2:2717/1‑148 (MQ=255)
                                 cccACAGGCGATTGATTTGCCGTCGCGCCTGCTTTCGCCGGATGCGCAGCACTGGCTGGGCACCGATCACTTAGGTCGCGATATTTTCTCGCGGCTGATGGCAGCGACCCGCGTGTCGCTCGGTTCGGTAATGGCCTGCCTGCTGCTgg                                                                                                                   >  4:386571/1‑149 (MQ=255)
                                   cacaGGCGATTGATTTGCCGTCGCGCCTGCTTTCGCCGGATGCGCAGCACTGGCTGGGCACCGATCACTTAGGTCGCGATATTTTCTCGCGGCTGATGGCAGCGACCCGCGTGTCGCTCGGTTCGGTAATGGCCTGCCTGCTGCTGGt                                                                                                                  >  1:812443/1‑148 (MQ=255)
                                   cacaGGCGATTGATTTGCCGTCGCGCCTGCTTTCGCCGGATGCGCAGCACTGGCTGGGCACCGATCACTTAGGTCGCGATATTTTCTCGCGGCTGATGGCAGCGACCCGCGTGTCGCTCGGTTCGGTAATGGCCTGCCTGCTGCTGGTg                                                                                                                 <  3:602307/149‑1 (MQ=255)
                                    acaGGCGATTGATTTGCCGTCGCGCCTGCTTTCGCCGGATGCGCAGCACTGGCTGGGCACCGATCACTTAGGTCGCGATATTTTCTCGCGGCTGATGGCAGCGACCCGCGTGTCGCTCGGTTCGGTAATGGCCTGCCTGCTGCTGGTGc                                                                                                                <  4:607471/149‑1 (MQ=255)
                                     caGGCGATTGATTTGCCGTCGCGCCTGCTTTCGCCGGATGCGCAGCACTGGCTGGGCACCGATCACTTAGGTCGCGATATTTTCTCGCGGCTGATGGCAGCGACCCGCGTGTCGCTCGGTTCGGTAATGGCCTGCCTGCTGCTGGTGc                                                                                                                <  2:57943/148‑1 (MQ=255)
                                       ggCGATTGATTTGCCGTCGCGCCTGCTTTCGCCGGATGCGCAGCACTGGCTGGGCACCGATCACTTAGGTCGCGATATTTTCTCGCGGCTGATGGCAGCGACCCGCGTGTCGCTCGGTTCGGTAATGGCCTGCCTGCTGCTGGTGCTg                                                                                                              <  1:112940/148‑1 (MQ=255)
                                         cGATTGATTTGCCGTCGCGCCTGCTTTCGCCGGATGCGCAGCACTGGCTGGGCACCGATCACTTAGGTCGCGATATTTTCTCGCGGCTGATGGCAGCGACCCGCGTGTCGCTCGGTTCGGTAATGGCCTGCCTGCTGCTGGTGCTGAc                                                                                                            >  1:584051/1‑148 (MQ=255)
                                           attgattTGCCGTCGCGCCTGCTTTCGCCGGATGCGCAGCACTGGCTGGGCACCGATCACTTAGGTCGCGATATTTTCTCGCGGCTGATGGCAGCGACCCGCGTGTCGCTCGGTTCGGTAATGGCCTGCCTGCTGCTGGTGCTGACAtt                                                                                                         >  3:110808/1‑149 (MQ=255)
                                             tgattTGCCGTCGCGCCTGCTTTCGCCGGATGCGCAGCACTGGCTGGGCACCGATCACTTAGGTCGCGATATTTTCTCGCGGCTGATGGCAGCGACCCGCGTGTCGCTCGGTTCGGTAATGGCCTGCCTGCTGCTGGTGCTGACATTAg                                                                                                       >  3:45494/1‑149 (MQ=255)
                                              gattTGCCGTCGCGCCTGCTTTCGCCGGATGCGCAGCACTGGCTGGGCACCGATCACTTAGGTCGCGATATTTTCTCGCGGCTGATGGCAGCGACCCGCGTGTCGCTCGGTTCGGTAATGGCCTGCCTGCTGCTGGTGCTGACATTAg                                                                                                       <  2:584051/148‑1 (MQ=255)
                                                tttGCCGTCGCGCCTGCTTTCGCCGGATGCGCAGCACTGGCTGGGCACCGATCACTTAGGTCGCGATATTTTCTCGCGGCTGATGGCAGCGACCCGCGTGTCGCTCGGTTCGGTAATGGCCTGCCTGCTGCTGGTGCTGACATTAggg                                                                                                     >  1:83291/1‑148 (MQ=255)
                                                 ttGCCGTCGCGCCTGCTTTCGCCGGATGCGCAGCACTGGCTGGGCACCGATCACTTAGGTCGCGATATTTTCTCGCGGCTGATGGCAGCGACCCGCGTGTCGCTCGGTTCGGTAATGGCCTGCCTGCTGCTGGTGCTGACATTAGGGc                                                                                                    >  1:347523/1‑148 (MQ=255)
                                                 ttGCCGTCGCGCCTGCTTTCGCCGGATGCGCAGCACTGGCTGGGCACCGATCACTTAGGTCGCGATATTTTCTCGCGGCTGATGGCAGCGACCCGCGTGTCGCTCGGTTCGGTAATGGCCTGCCTGCTGCTGGTGCTGACATTAGGGCt                                                                                                   <  3:49409/149‑1 (MQ=255)
                                                   gCCGTCGCGCCTGCTTTCGCCGGATGCGCAGCACTGGCTGGGCACCGATCACTTAGGTCGCGATATTTTCTCGCGGCTGATGGCAGCGACCCGCGTGTCGCTCGGTTCGGTAATGGCCTGCCTGCTGCTGGTGCTGACATTAGGGCTgg                                                                                                 <  4:45494/149‑1 (MQ=255)
                                                   gCCGTCGCGCCTGCTTTCGCCGGATGCGCAGCACTGGCTGGGCACCGATCACTTAGGTCGCGATATTTTCTCGCGGCTGATGGCAGCGACCCGCGTGTCGCTCGGTTCGGTAATGGCCTGCCTGCTGCTGGTGCTGACATTAGGGCTgg                                                                                                 <  4:527128/149‑1 (MQ=255)
                                                   accGTCGCGCCTGCTTTCGCCGGATGCGCAGCACTGGCTGGGCACCGATCACTTAGGTCGCGATATTTTCTCGCGGCTGATGGCAGCGACCCGCGTGTCGCTCGGTTCGGTAATGGCCTGCCTGCTGCTGGTGCTGACATTAGGGCTgg                                                                                                 <  4:615251/148‑1 (MQ=255)
                                                       tCGCGCCTGCTTTCGCCGGATGCGCAGCACTGGCTGGGCACCGATCACTTAGGTCGCGATATTTTCTCGCGGCTGATGGCAGCGACCCGCGTGTCGCTCGGTTCGGTAATGGCCTGCCTGCTGCTGGTGCTGACATTAGGGCTGGTTa                                                                                              >  2:18/1‑148 (MQ=255)
                                                       tCGCGCCTGCTTTCGCCGGATGCGCAGCACTGGCTGGGCACCGATCACTTAGGTCGCGATATTTTCTCGCGGCTGATGGCAGCGACCCGCGTGTCGCTCGGTTCGGTAATGGCCTGCCTGCTGCTGGTGCTGACATTAGGGCTGGTTa                                                                                              >  1:919633/1‑148 (MQ=255)
                                                         gcgcCTGCTTTCGCCGGATGCGCAGCACTGGCTGGGCACCGATCACTTAGGTCGCGATATTTTCTCGCGGCTGATGGCAGCGACCCGCGTGTCGCTCGGTTCGGTAATGGCCTGCCTGCTGCTGGTGCTGACATTAGGGCTGGTTATTg                                                                                           <  4:110808/149‑1 (MQ=255)
                                                           gcCTGCTTTCGCCGGATGCGCAGCACTGGCTGGGCACCGATCACTTAGGTCGCGATATTTTCTCGCGGCTGATGGCAGCGACCCGCGTGTCGCTCGGTTCGGTAATGGCCTGCCTGCTGCTGGTGCTGACATTAGGGCTGGTTATTggc                                                                                         <  3:610937/149‑1 (MQ=255)
                                                           gcCTGCTTTCGCCGGATGCGCAGCACTGGCTGGGCACCGATCACTTAGGTCGCGATATTTTCTCGCGGCTGATGGCAGCGACCCGCGTGTCGCTCGGTTCGGTAATGGCCTGCCTGCTGCTGGTGCTGACATTAGGGCTGGTTATTggc                                                                                         >  4:124086/1‑149 (MQ=255)
                                                           gcCTGCTTTCGCCGGATGCGCAGCACTGGCTGGGCACCGATCACTTAGGTCGCGATATTTTCTCGCGGCTGATGGCAGCGACCCGCGTGTCGCTCGGTTCGGTAATGGCCTGCCTGCTGCTGGTGCTGACATTAGGGCTGGTTATTgg                                                                                          <  1:490569/148‑1 (MQ=255)
                                                            ccTGCTTTCGCCGGATGCGCAGCACTGGCTGGGCACCGATCACTTAGGTCGCGATATTTTCTCGCGGCTGATGGCAGCGACCCGCGTGTCGCTCGGTTCGGTAAt                                                                                                                                    >  2:380883/1‑105 (MQ=255)
                                                            ccTGCTTTCGCCGGATGCGCAGCACTGGCTGGGCACCGATCACTTAGGTCGCGATATTTTCTCGCGGCTGATGGCAGCGACCCGCGTGTCGCTCGGTTCGGTAATGGCCTGCCTGCTGCTGGTGCTGACATTAGGGCTGGTTATTggc                                                                                         >  1:170954/1‑148 (MQ=255)
                                                             cTGCTTTCGCCGGATGCGCAGCACTGGCTGGGCACCGATCACTTAGGTCGCGATATTTTCTCGCGGCTGATGGCAGCGACCCGCGTGTCGCTCGGTTCGGTAATGGCCTGCCTGCTGCTGGTGCTGACATTAGGGCTGGTTATTagcga                                                                                       >  3:698715/1‑148 (MQ=255)
                                                              tGCTTTCGCCGGATGCGCAGCACTGGCTGGGCACCGATCACTTAGGTCGCGATATTTTCTCGCGGCTGATGGCAGCGACCCGCGTGTCGCTCGGTTCGGTAATGGCCTGCCTGCTGCTGGTGCTGACATTAGGGCTGGTTATTggcggc                                                                                      >  3:925665/1‑149 (MQ=255)
                                                               gCTTTCGCCGGATGCGCAGCACTGGCTGGGCACCGATCACTTAGGTCGCGATATTTTCTCGCGGCTGATGGCAGCGACCCGCGTGTCGCTCGGTTCGGTAATGGCCTGCCTGCTGCTGGTGCTGACATTAGGGCTGGTTATTggcggc                                                                                      >  1:879323/1‑148 (MQ=255)
                                                               gCTTTCGCCGGATGCGCAGCACTGGCTGGGCACCGATCACTTAGGTCGCGATATTTTCTCGCGGCTGATGGCAGCGACCCGCGTGTCGCTCGGTTCGGTAATGGCCTGCCTGCTGCTGGTGCTGACATTAGGGCTGGTTATTGGCGGCa                                                                                     <  4:522551/149‑1 (MQ=255)
                                                                 tttCGCCGGATGCGCAGCACTGGCTGGGCACCGATCACTTAGGTCGCGATATTTTCTCGCGGCTGATGGCAGCGACCCGCGTGTCGCTCGGTTCGGTAATGGCCTGCCTGCTGCTGGTGCTGACATTAGGGCTGGTTATTGGCGGCAg                                                                                    <  2:175649/148‑1 (MQ=255)
                                                                        ggATGCGCAGCACTGGCTGGGCACCGATCACTTAGGTCGCGATATTTTCTCGCGGCTGATGGCAGCGACCCGCGTGTCGCTCGGTTCGGTAATGGCCTGCCTGCTGCTGGTGCTGACATTAGGGCTGGTTATTGGCGGCAGCGCCgg                                                                              >  3:756821/1‑147 (MQ=255)
                                                                        ggATGCGCAGCACTGGCTGGGCACCGATCACTTAGGTCGCGATATTTTCTCGCGGCTGATGGCAGCGACCCGCGTGTCGCTCGGTTCGGTAATGGCCTGCCTGCTGCTGGTGCTGACATTAGGGCTGGTTATTGGCGGCAGCGCCgg                                                                              <  4:756821/147‑1 (MQ=255)
                                                                         gATGCGCAGCACTGGCTGGGCACCGATCACTTAGGTCGCGATATTTTCTCGCGGCTGATGGCAGCGACCCGCGTGTCGCTCGGTTCGGTAATGGCCTGCCTGCTGCTGGTGCTGACATTAGGGCTGGTTATTGGCGGCAGCGCCGGGtt                                                                           >  3:443039/1‑149 (MQ=255)
                                                                               cagcaCTGGCTGGGCACCGATCACTTAGGTCGCGATATTTTCTCGCGGCTGATGGCAGCGACCCGCGTGTCGCTCGGTTCGGTAATGGCCTGCCTGCTGCTGGTGCTGACATTAGGGCTGGTTATTGGCGGCAGCGCCGGGTTGATTgg                                                                     >  4:154003/1‑149 (MQ=255)
                                                                                agcaCTGGCTGGGCACCGATCACTTAGGTCGCGATATTTTCTCGCGGCTGATGGCAGCGACCCGCGTGTCGCTCGGTTCGGTAATGGCCTGCCTGCTGCTGGTGCTGACATTAGGGCTGGTTATTGGCGGCAGCGCCGGGTTGATTggc                                                                    <  4:105335/149‑1 (MQ=255)
                                                                                agcaCTGGCTGGGCACCGATCACTTAGGTCGCGATATTTTCTCGCGGCTGATGGCAGCGACCCGCGTGTCGCTCGGTTCGGTAATGGCCTGCCTGCTGCTGGTGCTGACATTAGGGCTGGTTATTGGCGGCAGCGCCGGGTTGATTgg                                                                     <  2:879323/148‑1 (MQ=255)
                                                                                agcaCTGGCTGGGCACCGATCACTTAGGTCGCGATATTTTCTCGCGGCTGATGGCAGCGACCCGCGTGTCGCTCGGTTCGGTAATGGCCTGCCTGCTGCTGGTGCTGACATTAGGGCTGGTTATTGGCGGCAGCGCCGGGTTGATTgg                                                                     <  2:290784/148‑1 (MQ=255)
                                                                                    ctggctggGCACCGATCACTTAGGTCGCGATATTTTCTCGCGGCTGATGGCAGCGACCCGCGTGTCGCTCGGTTCGGTAATGGCCTGCCTGCTGCTGGTGCTGACATTAGGGCTGGTTATTGGCGGCAGCGCCGGGTTGATTggcgggc                                                                <  3:126694/149‑1 (MQ=255)
                                                                                    ctggctggGCACCGATCACTTAGGTCGCGATATTTTCTCGCGGCTGATGGCAGCGACCCGCGTGTCGCTCGGTTCGGTAATGGCCTGCCTGCTGCTGGTGCTGACATTAGGGCTGGTTATTGGCGGCAGCGCCGGGTTGATTggcgggc                                                                >  3:211946/1‑149 (MQ=255)
                                                                                      ggctggGCACCGATCACTTAGGTCGCGATATTTTCTCGCGGCTGATGGCAGCGACACGCGGGTCGCtgcc                                                                                                                                             >  3:324630/1‑67 (MQ=38)
                                                                                      ggctggGCACCGATCACTTAGGTCGCGATATTTTCTCGCGGCTGATGGCAGCGACACGCGGGTCGCtgcc                                                                                                                                             <  4:324630/70‑4 (MQ=38)
                                                                                       gctggGCACCGATCACTTAGGTCGCGATATTTTCTCGCGGCTGATGGCAGCGACCCGCGTGTCGCTCGGTTCGGTAATGGCCTGCctgctgc                                                                                                                      <  3:616054/92‑1 (MQ=255)
                                                                                       gctggGCACCGATCACTTAGGTCGCGATATTTTCTCGCGGCTGATGGCAGCGACCCGCGTGTCGCTCGGTTCGGTAATGGCCTGCctgctgc                                                                                                                      >  4:616054/1‑92 (MQ=255)
                                                                                         tggGCACCGATCACTCAGGTCGCGATATTTTCTCGCGGCTGATGGCAGCGACCCGCGTGTCGCTCGGTTCGGTAATGGCCTGCCTGCTGCTGGCGCTGACATTAGGGCTGGTTATTGGCGGCAGCGCCGGGTTGATTGGCGGGCGCGt                                                            <  2:935722/148‑1 (MQ=255)
                                                                                          gggCACCGATCACTTAGGTCGCGATATTTTCTCGCGGCTGATGGCAGCGACCCGCGTGTCGCTCGGTTCGGTAATGGCCTGCCTGCTGCTGGTGCTGACATTAGGGCTGGTTATTGGCGGCAGCGCCGGGTTGATTGGCGGGCGCGtt                                                           >  2:921430/1‑148 (MQ=255)
                                                                                              aCCGATCACTTAGGTCGCGATATTTTCTCGCGGCTGATGGCAGCGACCCGCGTGTCGCTCGGTTCGGTAATGGCCTGCCTGCTGCTGGTGCTGACATTAGGGCTGGTTATTGGCGGCAGCGCCGGGTTGATTGGCGGGCGCGTTGATCa                                                      <  4:325496/149‑1 (MQ=255)
                                                                                               ccGATCACTTAGGTCGCGATATTTTCTCGCGGCTGATGGCAGCGACCCGCGTGTCGCTCGGTTCGGTAATGGCCTGCCTGCTGCTGGTGCTGACATTAGGGCTGGTTATTGGCGGCAGCGCCGGGTTGATTGGCGGGCGCGTTGATCa                                                      >  1:322764/1‑148 (MQ=255)
                                                                                               ccGATCACTTAGGTCGCGATATTTTCTCGCGGCTGATGGCAGCGACCCGCGTGTCGCTCGGTTCGGTAATGGCCTGCCTGCTGCTGGTGCTGACATTAGGGCTGGTTATTGGCGGCAGCGCCGGGTTGATTGGCGGGCGCGTTGATCa                                                      >  1:838853/1‑148 (MQ=255)
                                                                                               ccGATCACTTAGGTCGCGATATTTTCTCGCGGCTGATGGCAGCGACCCGCGTGTCGCTCGGTTCGGTAATGGCCTGCCTGCTGCTGGTGCTGACATTAGGGCTGGTTATTGGCGGCAGCGCCGGGTTGATTGGCGGGCGCGTTGATCAg                                                     <  3:386571/149‑1 (MQ=255)
                                                                                                    cACTTAGGTCGCGATATTTTCTCGCGGCTGATGGCAGCGACCCGCGTGTCGCTCGGTTCGGTAATGGCCTGCCTGCTGCTGGTGCTGACATTAGGGCTGGTTATTGGCGGCAGCGCCGGGTTGATTGGCGGGcgcg                                                             <  1:390107/136‑1 (MQ=255)
                                                                                                        tAGGTCGCGATATTTTCTCGCGGCTGATGGCAGCGACCCGCGTGTCGCTCGGTTCGGTAATGGCCTGCCTGCTGCTGGTGCTGACATTAGGGCTGGTTATTGGCGGCAGCGCCGGGTTGATTGGCGGGCGCGTTGATCAGGCCACCAt                                             >  1:245641/1‑148 (MQ=255)
                                                                                                          ggTCGCGATATTTTCTCGCGGCTGATGGCAGCGACCCGCGTGTCGCTCGGTTCGGTAATGGCCTGCCTGCTGCTGGTGCTGACATTGGGGCTGGTTATTGGCGGCAGCGCCGGGTTGCTTGGCGGGCGCGTTGATCAGGCCACCATgc                                           >  2:635940/1‑148 (MQ=255)
                                                                                                           gTCGCGATATTTTCTCGCGGCTGATGGCAGCGACCCGCGTGTCGCTCGGTTCGGTAATGGCCTGCCTGCTGCTGGTGCTGACATTAGGGCTGGTTATTGGCGGCAGCGCCGGGTTGATTGGCGGGCGCGTTGATCAGGCCACCATgcgc                                         >  4:758949/1‑149 (MQ=255)
                                                                                                           gTCGCGATATTTTCTCGCGGCTGATGGCAGCGACCCGCGTGTCGCTCGGTTCGGTAATGGCCTGCCTGCTGCTGGTGCTGACATTAGGGCTGGTTATTGGCGGCAGCGCCGGGTTGATTGGCGGGCGCGTTGATCAGGCCACCATgcgc                                         >  3:873038/1‑149 (MQ=255)
                                                                                                                gATATTTTCTCGCGGCTGATGGCAGCGACCCGCGTGTCGCTCGGTTCGGTAATGGCCTGCCTGCTGCTGGTGCTGACATTAGGGCTGGTTATTGGCGGCAGCGCCGGGTTGATTGGCGGGCGCGTTGATCAGGCCACCATGCGCGTCg                                     <  2:682635/148‑1 (MQ=255)
                                                                                                                gATATTTTCTCGCGGCTGATGGCAGCGACCCGCGTGTCGCTCGGTTCGGTAATGGCCTGCCTGCTGCTGGTGCTGACATTAGGGCTGGTTATTGGCGGCAGCGCCGGGTTGATTGGCGGGCGCGTTGATCAGGCCACCATGCGCGTCg                                     >  2:379177/1‑148 (MQ=255)
                                                                                                                   aatccctcGCGGCTGATGGCAGCGACCCGCGTGTCGCTCGGTTCGGTAATGGCCTGCCTGCTGCTGGTGCTGACATTAGGGCTGGTTATTGGCGGCAGCGCCGGGTTGATTGGCGGGCGCGTTGATCAGGCCACCATGCGCGTCGCCg                                  <  1:445630/143‑1 (MQ=255)
                                                                                                                    ttttCTCGCGGCTGATGGCAGCGACCCGCGTGTCGCTCGGTTCGGTAATGGCCTGCCTGCTGCTGGTGCTGACATTAGGGCTGGTTATTGGCGGCAGCGCCGGGTTGATTGGCGGGCGCGTTGATCAGGCCACCATGCGCGTCGCCGat                                >  3:742704/1‑149 (MQ=255)
                                                                                                                         tcGCGGCTGATGGCAGCGACCCGCGTGTCGCTCGGTTCGGTAATGGCCTGCctgctgctg                                                                                                                    >  4:234653/1‑60 (MQ=255)
                                                                                                                         tcGCGGCTGATGGCAGCGACCCGCGTGTCGCTCGGTTCGGTAATGGCCTGCctgctgctg                                                                                                                    <  3:234653/60‑1 (MQ=255)
                                                                                                                         tcGCGGCTGATGGCAGCGACCCGCGTGTCGCTCGGTTCGGTAATGGCCTGCCTGCTGCTGGTGCTGACATTAGGGCTGGTTATTGGCGGCAGCGCCGGGTTGATTGGCGGGCGCGTTGATCAGGCCACCATGCGCGTCGCCGATATGtt                           >  3:274057/1‑149 (MQ=255)
                                                                                                                            cgGCTGATGGCAGCGACCCGCGTGTCGCTCGGTTCGGTAATGGCCTGCCTGCTGCTGGTGCTGACATTAGGGCTGGTTATTGGCGGCAGCGCCGGGTTGATTGGCGGGCGCGTTGATCAGGCCACCATGCGCGTCGCCGATATGTTTa                         <  1:2717/148‑1 (MQ=255)
                                                                                                                            cgGCTGATGGCAGCGACCCGCGTGTCGCTCGGTTCGGTAATGGCCTGCCTGCTGCTGGTGCTGACATTAGGGCTGGTTATTGGCGGCAGCGCCGGGTTGATTGGCGGGCGCGTTGATCAGGCCACCATGCGCGTCGCCGATATGTTTa                         <  2:825264/148‑1 (MQ=255)
                                                                                                                               cTGATGGCAGCGACCCGCGTGTCGCTCGGTTCGGTAATGGCCTGCCTGCTGCTGGTGCTGACATTAGGGCTGGTTATTGGCGGCAGCGCCGGGTTGATTGGCGGGCGCGTTGATCAGGCCACCATGCGCGTCGCCGATATGTTTATGAc                     <  4:443039/149‑1 (MQ=255)
                                                                                                                                   tGGCAGCGACCCGCGTGTCGCTCGGTTCGGTAATGGCCTGCCTGCTGCTGGTGCTGACATTAGGGCTGGTTATTGGCGGCAGCGCCGGGTTGATTGGCGGGCGCGTTGATCAGGCCACCATGCGCGTCGCCGATATGTTTATGACCtt                  <  2:347523/148‑1 (MQ=255)
                                                                                                                                     gCAGCGACCCGCGTGTCGCTCGGTTCGGTAATGGCCTGCCTGCTGCTGGTGCTGACATTAGGGCTGGTTATTGGCGGCAGCGCCGGGTTGATTGGCGGGCGCGTTGATCAGGCCACCATGCGCGTCGCCGATATGTTTATGACCTTccc               >  3:446626/1‑149 (MQ=255)
                                                                                                                                     gCAGCGACCCGCGTGTCGCTCGGTTCGGTAATGGCCTGCCTGCTGCTGGTGCTGACATTAGGGCTGGTTATTGGCGGCAGCGCCGGGTTGATTGGCGGGCGCGTTGATCAGGCCACCATGCGCGTCGCCGATATGTTTATGACCTTcc                >  2:39307/1‑148 (MQ=255)
                                                                                                                                        gCGACCCGCGTGTCGCTCGGTTCGGTAATGGCCTGCCTGCTGCTGGTGCTGACATTAGGGCTGGTTATTGGCGGCAGCGCCGGGTTGATTGGCGGGCGCGTTGATCAGGCCACCATGCGCGTCGCCGATATGTTTATGACCTTCCCGa             >  1:238827/1‑148 (MQ=255)
                                                                                                                                         cGACCCGCGTGTCGCTCGGTTCGGTAATGGCCTGCCTGCTGCTGGTGCTGACATTAGGGCTGGTTATTGGCGGCAGCGCCGGGTTGATTGGCGGGCGCGTTGATCAGGCCACCATGCGCGTCGCCGATATGTTTATGACCTTCCCGAc            <  2:179856/148‑1 (MQ=255)
                                                                                                                                           aCCCGCGTGTCGCTCGGTTCGGTAATGGCCTGCCTGCTGCTGGTGCTGACATTAGGGCTGGTTATTGGCGGCAGCGCCGGGTTGATTGGCGGGCGCGTTGATCAGGCCACCATGCGCGTCGCCGATATGTTTATGACCTTCCCGACCt          >  2:111858/1‑148 (MQ=255)
                                                                                                                                              cgcgTGTCGCTCGGTTCGGTAATGGCCTGCCTGCTGCTGGTGCTGACATTAGGGCTGGTTATTGGCGGCAGCGCCGGGTTGATTGGCGGGTGCGTTGATCAGGCCACCATGCGCGTCGCCGATATGTTTATGACCTTCCCGACCTCGa       >  2:485052/1‑148 (MQ=255)
                                                                                                                                              cgcgTGTCGCTCGGTTCGGTAATGGCCTGCCTGCTGCTGGTGCTGACATTAGGGCTGGTTATTGGCGGCAGCGCCGGGTTGATTGGCGGGCGCGTTGATCAGGCCACCATGCGCGTCGCCGATATGTTTATGACCTTCCCGACCTCGAt      <  4:804790/149‑1 (MQ=255)
                                                                                                                                                cgTGTCGCTCGGTTCTGTAATGGCCTGCCTGGTGCTGGTGCTGACATTAGGGCTGGTTATTGGCGGCAGCGCCGGGTTGATTGGCGGGCGCGTTGATCAGGCCACCATGCGCGTCGCCGATATGTTTATGACCTTCCCGACCTCGAtt     <  2:338673/148‑1 (MQ=255)
                                                                                                                                                 gtgtCGCTCGGTTCGGTAATGGCCTGCCTGCTGCTGGTGCTGACATTAGGGCTGGTTATTGGCGGCAGCGCCGGGTTGATTGGCGGGCGCGTTGATCAGGCCACCATGCGCGTCGCCGATATGTTTATGACCTTCCCGACCTCGATTc    <  1:468616/148‑1 (MQ=255)
                                                                                                                                                 gtgtCGCTCGGTTCGGTAATGGCCTGCCTGCTGCTGGTGCTGACATTAGGGCTGGTTATTGGCGGCAGCGCCGGGTTGATTGGCGGGCGCGTTGATCAGGCCACCATGCGCGTCGCCGATATGTTTATGACCTTCCCGACCTCGATTCt   <  4:396545/149‑1 (MQ=255)
                                                                                                                                                  tgtCGCTCGGTTCGGTAATGGCCTGCCTGCTGCTGGTGCTGACATTAGGGCTGGTTATTGGCGGCAGCGCCGGGTTGATTGGCGGGCGCGTTGATCAGGCCACCATGCGCGTCGCCGATATGTTTATGACCTTCCCGACCTCGATTCTg  <  4:90105/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
TGCGCTCACCAGCCAGTGGTGGCTGCCGTATGACCCACAGGCGATTGATTTGCCGTCGCGCCTGCTTTCGCCGGATGCGCAGCACTGGCTGGGCACCGATCACTTAGGTCGCGATATTTTCTCGCGGCTGATGGCAGCGACCCGCGTGTCGCTCGGTTCGGTAATGGCCTGCCTGCTGCTGGTGCTGACATTAGGGCTGGTTATTGGCGGCAGCGCCGGGTTGATTGGCGGGCGCGTTGATCAGGCCACCATGCGCGTCGCCGATATGTTTATGACCTTCCCGACCTCGATTCTG  >  NC_000913/3616253‑3616547

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