Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,616,394 C→A 5.6% T71T (ACC→ACA nikC → nickel ABC transporter permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,616,3940CA 5.6% 375.0 / 9.8 107T71T (ACC→ACAnikCnickel ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base C (49/52);  new base A (2/4);  total (51/56)
Fisher's exact test for biased strand distribution p-value = 6.81e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.72e-01

Only 100 of 110 total aligned reads displayed.

GCTGATTGCGCTCACCAGCCAGTGGTGGCTGCCGTATGACCCACAGGCGATTGATTTGCCGTCGCGCCTGCTTTCGCCGGATGCGCAGCACTGGCTGGGCACCGATCACTTAGGTCGCGATATTTTCTCGCGGCTGATGGCAGCGACCCGCGTGTCGCTCGGTTCGGTAATGGCCTGCCTGCTGCTGGTGCTGACATTAGGGCTGGTTATTGGCGGCAGCGCCGGGTTGATTGGCGGGCGCGTTGATCAGGCCACCATGCGCGTCGCCGATATGTTTATGACCTTCCCGACCT  >  NC_000913/3616247‑3616539
                                                                                                                                                   |                                                                                                                                                 
gcTGATTGCGCTCACCAGCCAGTGGTGGCTGCCGTATGACCCACAGGCGATTGATTTGCCGTCGCGCCTGCTTTCGCCGGATGCGCAGCACTGGCTGGGCACCGATCACTTAGGTCGCGATATTTTCTCGCGGCTGATGGCAGCGAccc                                                                                                                                                  >  3:360531/1‑149 (MQ=255)
gcTGATTGCGCTCACCAGCCAGTGGTGGCTGCCGTATGACCCACAGGCGATTGATTTGCCGTCGCGCCTGCTTTCGCCGGATGCGCAGCACTGGCTGGGCACCGATCACTTAGGTCGCGATATTTTCTCGCGGCTGATGGCAGCGAcc                                                                                                                                                   <  2:947520/148‑1 (MQ=255)
 cTGATTGCGCTCACCAGCCAGTGGTGGCTGCCGTATGACCCACAGGCGATTGATTTGCCGTCGCGCCTGCTTTCGCCGGATGCGCAGCACTGGCTGGGCACCGATCACTTAGGTCGCGATATTTTCTCGCGGCTGATGGCAGCGAccc                                                                                                                                                  >  2:773811/1‑148 (MQ=255)
    aTTGCGCTCACCAGCCAGTGGTGGCTGCCGTATGACCCACAGGCGATTGATTTGCCGTCGCGCCTGCTTTCGCCGGATGCGCAGCACTGGCTGGGCACCGATCACTTAGGTCGCGATATTTTCTCGCGGCTGATGGCAGCGACCcgcg                                                                                                                                               >  2:809093/1‑148 (MQ=255)
    aTTGCGCTCACCAGCCAGTGGTGGCTGCCGTATGACCCACAGGCGATTGATTTGCCGTCGCGCCTGCTTTCGCCGGATGCGCAGCACTGGCTGGGCACCGATCACTTAGGTCGCGATATTTTCTCGCGGCTGATGGCAGCGACCcgcg                                                                                                                                               <  2:550642/148‑1 (MQ=255)
      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)
              ccagccagTGGTGGCTGCCGTATGACCCACAGGCGATTGATTTGCCGTCGCGCCTGCTTTCGCCGGATGCGCAGCACTGGCTGGGCACCGATCACTTAGGTCGCGATATTTTCTCGCGGCTGATGGCAGCGACACGCGGGTCGCtgcca                                                                                                                                    <  4:266217/149‑5 (MQ=255)
                agccagTGGTGGCTGCCGTATGACCCACAGGCGATTGATTTGCCGTCGCGCCTGCTTTCGCCGGATGCGCAGCACTGGCTGGGCACCGATCACTTAGGTCGCGATATTTTCTCGCGGCTGATGGCAGCGACACGCGGGTCGCtgc                                                                                                                                      <  4:354063/145‑3 (MQ=255)
                agccagTGGTGGCTGCCGTATGACCCACAGGCGATTGATTTGCCGTCGCGCCTGCTTTCGCCGGATGCGCAGCACTGGCTGGGCACCGATCACTTAGGTCGCGATATTTTCTCGCGGCTGATGGCAGCGACACGCGGGTCGCtgc                                                                                                                                      >  3:354063/1‑143 (MQ=255)
                 gccagTGGTGGCTGCCGTATGACCCACAGGCGATTGATTTGCCGTCGCGCCTGCTTTCGCCGGATGCGCAGCACTGGCTGGGCACCGATCACTTAGGTCGCGATATTTTCTCGCGGCTGATGGCAGCGACCCGCGTGTCGCtcggttc                                                                                                                                  >  1:825264/1‑148 (MQ=255)
                  ccagTGGTGGCTGCCGTATGACCCACAGGCGATTGATTTGCCGTCGCGCCTGCTTTCGCCGGATGCGCAGCACTGGCTGGGCACCGATCACTTAGGTCGCGATATTTTCTCGCGGCTGATGGCAGCGACCCGCGTGTCGCtcggttcg                                                                                                                                 <  2:850490/148‑1 (MQ=255)
                      tggtggCTGCCGTATGACCCACAGGCGATTGATTTGCCGTCGCGCCTGCTTTCGCCGGATGCGCAGCACTGGCTGGGCACCGATCACTTAGGTCGCGATATTTTCTCGCGGCTGATGGCAGCGACCCGCGTGTCGCTCGGTTCGGTaa                                                                                                                             <  1:530893/148‑1 (MQ=255)
                       ggtggCTGTCGTATGACCCACAGGCGATTGATTTGCCGTCGCGTCTGCTTTCGCCGGATGCGCAGCACTGGCTGGGCAGCGATCGCTTATGTCGCGATATTTTCACGCGGCTGATGGCAGCGACCCGCGTGTCGCTCGGTTCGGTAATg                                                                                                                           <  4:565122/149‑1 (MQ=255)
                       ggtggCTGCCGTATGACCCACAGGCGATTGATTTGCCGTCGCGCCTGCTTTCGCCGGATGCGCAGCACTGGCTGGGCACCGATCACTTAGGTCGCGATATTTTCTCGCGGCTGATGGCAGCGACCCGCGTGTCGCtcggt                                                                                                                                    <  2:416566/140‑1 (MQ=255)
                       ggtggCTGCCGTATGACCCACAGGCGATTGATTTGCCGTCGCGCCTGCTTTCGCCGGATGCGCAGCACTGGCTGGGCACCGATCACTTAGGTCGCGATATTTTCTCGCGGCTGATGGCAGCGACCCGCGTGTCGCtcggt                                                                                                                                    >  1:416566/1‑140 (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)
                                 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                                                                                                            >  2:2717/1‑148 (MQ=255)
                                       cccACAGGCGATTGATTTGCCGTCGCGCCTGCTTTCGCCGGATGCGCAGCACTGGCTGGGCACCGATCACTTAGGTCGCGATATTTTCTCGCGGCTGATGGCAGCGACCCGCGTGTCGCTCGGTTCGGTAATGGCCTGCctgctgctg                                                                                                            <  1:809093/148‑1 (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)
                                                        tGCCGTCGCGCCTGCTTTCGCCGGATGCGCAGCACTGGCTGGGCACCGATCACTTAGGTCGCGATATTTTCTCGCGGCTGATGGCAGCGACCCGCGTGTCGCTCGGTTCGGTAATGGCCTGCCTGCTGCTGGTGCTGACATTAGGGCt                                                                                           <  1:773811/148‑1 (MQ=255)
                                                         gCCGTCGCGCCTGCTTTCGCCGGATGCGCAGCACTGGCTGGGCACCGATCACTTAGGTCGCGATATTTTCTCGCGGCTGATGGCAGCGACCCGCGTGTCGCTCGGTTCGGTAATGGCCTGCCTGCTGCTGGTGCTGACATTAGGGCTgg                                                                                         <  4:45494/149‑1 (MQ=255)
                                                         gCCGTCGCGCCTGCTTTCGCCGGATGCGCAGCACTGGCTGGGCACCGATCACTTAGGTCGCGATATTTTCTCGCGGCTGATGGCAGCGACCCGCGTGTCGCTCGGTTCGGTAATGGCCTGCCTGCTGCTGGTGCTGACATTAGGGCTg                                                                                          <  2:12956/148‑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)
                                                                              ggATGCGCAGCACTGGCTGGGCACCGATCACTTAGGTCGCGATATTTTCTCGCGGCTGATGGCAGCGACCCGCGTGTCGCTCGGTTCGGTAATGGCCTGCCTGCTGCTGGTGCTGACATTAGGGCTGGTTATTGGCGGCAGCGCCgg                                                                      <  4:756821/147‑1 (MQ=255)
                                                                              ggATGCGCAGCACTGGCTGGGCACCGATCACTTAGGTCGCGATATTTTCTCGCGGCTGATGGCAGCGACCCGCGTGTCGCTCGGTTCGGTAATGGCCTGCCTGCTGCTGGTGCTGACATTAGGGCTGGTTATTGGCGGCAGCGCCgg                                                                      >  3:756821/1‑147 (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:290784/148‑1 (MQ=255)
                                                                                      agcaCTGGCTGGGCACCGATCACTTAGGTCGCGATATTTTCTCGCGGCTGATGGCAGCGACCCGCGTGTCGCTCGGTTCGGTAATGGCCTGCCTGCTGCTGGTGCTGACATTAGGGCTGGTTATTGGCGGCAGCGCCGGGTTGATTgg                                                             <  2:879323/148‑1 (MQ=255)
                                                                                          ctggctggGCACCGATCACTTAGGTCGCGATATTTTCTCGCGGCTGATGGCAGCGACCCGCGTGTCGCTCGGTTCGGTAATGGCCTGCCTGCTGCTGGTGCTGACATTAGGGCTGGTTATTGGCGGCAGCGCCGGGTTGATTggcgggc                                                        <  3:126694/149‑1 (MQ=255)
                                                                                          ctggctggGCACCGATCACTTAGGTCGCGATATTTTCTCGCGGCTGATGGCAGCGACCCGCGTGTCGCTCGGTTCGGTAATGGCCTGCCTGCTGCTGGTGCTGACATTAGGGCTGGTTATTGGCGGCAGCGCCGGGTTGATTggcgggc                                                        >  3:211946/1‑149 (MQ=255)
                                                                                            ggctggGCACCGATCACTTAGGTCGCGATATTTTCTCGCGGCTGATGGCAGCGACACGCGGGTCGCtgcc                                                                                                                                     <  4:324630/70‑4 (MQ=38)
                                                                                            ggctggGCACCGATCACTTAGGTCGCGATATTTTCTCGCGGCTGATGGCAGCGACACGCGGGTCGCtgcc                                                                                                                                     >  3:324630/1‑67 (MQ=38)
                                                                                             gctggGCACCGATCACTTAGGTCGCGATATTTTCTCGCGGCTGATGGCAGCGACCCGCGTGTCGCTCGGTTCGGTAATGGCCTGCctgctgc                                                                                                              >  4:616054/1‑92 (MQ=255)
                                                                                               tggGCACCGATCACTCAGGTCGCGATATTTTCTCGCGGCTGATGGCAGCGACCCGCGTGTCGCTCGGTTCGGTAATGGCCTGCCTGCTGCTGGCGCTGACATTAGGGCTGGTTATTGGCGGCAGCGCCGGGTTGATTGGCGGGCGCGt                                                    <  2:935722/148‑1 (MQ=255)
                                                                                                gggCACCGATCACTTAGGTCGCGATATTTTCTCGCGGCTGATGGCAGCGACCCGCGTGTCGCTCGGTTCGGTAATGGCCTGCCTGCTGCTGGTGCTGACATTAGGGCTGGTTATTGGCGGCAGCGCCGGGTTGATTGGCGGGCGCGtt                                                   >  2:921430/1‑148 (MQ=255)
                                                                                                gggCACCGATCACTTAGGTCGCGATATTTTCTCGCGGCTGATGGCAGCGACCCGCGTGTCGCTCGGTTCGGTAATGGCCTGCCTGCTGCTGGTGCTGACATTAGGGCTGGGTATTGGCGGCAGCGCCGGGTTGATTGGCGGGCGCGTTg                                                  >  3:7455/1‑149 (MQ=255)
                                                                                                    aCCGATCACTTAGGTCGCGATATTTTCTCGCGGCTGATGGCAGCGACCCGCGTGTCGCTCGGTTCGGTAATGGCCTGCCTGCTGCTGGTGCTGACATTAGGGCTGGTTATTGGCGGCAGCGCCGGGTTGATTGGCGGGCGCGTTGATCa                                              <  4:325496/149‑1 (MQ=255)
                                                                                                     ccGATCACTTAGGTCGCGATATTTTCTCGCGGCTGATGGCAGCGACCCGCGTGTCGCTCGGTTCGGTAATGGCCTGCCTGCTGCTGGTGCTGACATTAGGGCTGGTTATTGGCGGCAGCGCCGGGTTGATTGGCGGGCGCGTTGATCa                                              >  1:838853/1‑148 (MQ=255)
                                                                                                     ccGATCACTTAGGTCGCGATATTTTCTCGCGGCTGATGGCAGCGACCCGCGTGTCGCTCGGTTCGGTAATGGCCTGCCTGCTGCTGGTGCTGACATTAGGGCTGGTTATTGGCGGCAGCGCCGGGTTGATTGGCGGGCGCGTTGATCa                                              >  1:322764/1‑148 (MQ=255)
                                                                                                     ccGATCACTTAGGTCGCGATATTTTCTCGCGGCTGATGGCAGCGACCCGCGTGTCGCTCGGTTCGGTAATGGCCTGCCTGCTGCTGGTGCTGACATTAGGGCTGGTTATTGGCGGCAGCGCCGGGTTGATTGGCGGGCGCGTTGATCAg                                             <  3:386571/149‑1 (MQ=255)
                                                                                                          cACTTAGGTCGCGATATTTTCTCGCGGCTGATGGCAGCGACCCGCGTGTCGCTCGGTTCGGTAATGGCCTGCCTGCTGCTGGTGCTGACATTAGGGCTGGTTATTGGCGGCAGCGCCGGGTTGATTGGCGGGcgcg                                                     >  2:390107/1‑136 (MQ=255)
                                                                                                          cACTTAGGTCGCGATATTTTCTCGCGGCTGATGGCAGCGACCCGCGTGTCGCTCGGTTCGGTAATGGCCTGCCTGCTGCTGGTGCTGACATTAGGGCTGGTTATTGGCGGCAGCGCCGGGTTGATTGGCGGGcgcg                                                     <  1:390107/136‑1 (MQ=255)
                                                                                                              tAGGTCGCGATATTTTCTCGCGGCTGATGGCAGCGACCCGCGTGTCGCTCGGTTCGGTAATGGCCTGCCTGCTGCTGGTGCTGACATTAGGGCTGGTTATTGGCGGCAGCGCCGGGTTGATTGGCGGGCGCGTTGATCAGGCCACCAt                                     >  1:245641/1‑148 (MQ=255)
                                                                                                                ggTCGCGATATTTTCTCGCGGCTGATGGCAGCGACCCGCGTGTCGCTCGGTTCGGTAATGGCCTGCCTGCTGCTGGTGCTGACATTGGGGCTGGTTATTGGCGGCAGCGCCGGGTTGCTTGGCGGGCGCGTTGATCAGGCCACCATgc                                   >  2:635940/1‑148 (MQ=255)
                                                                                                                 gTCGCGATATTTTCTCGCGGCTGATGGCAGCGACCCGCGTGTCGCTCGGTTCGGTAATGGCCTGCCTGCTGCTGGTGCTGACATTAGGGCTGGTTATTGGCGGCAGCGCCGGGTTGATTGGCGGGCGCGTTGATCAGGCCACCATgcgc                                 >  3:873038/1‑149 (MQ=255)
                                                                                                                 gTCGCGATATTTTCTCGCGGCTGATGGCAGCGACCCGCGTGTCGCTCGGTTCGGTAATGGCCTGCCTGCTGCTGGTGCTGACATTAGGGCTGGTTATTGGCGGCAGCGCCGGGTTGATTGGCGGGCGCGTTGATCAGGCCACCATgcgc                                 >  4:758949/1‑149 (MQ=255)
                                                                                                                      gATATTTTCTCGCGGCTGATGGCAGCGACCCGCGTGTCGCTCGGTTCGGTAATGGCCTGCCTGCTGCTGGTGCTGACATTAGGGCTGGTTATTGGCGGCAGCGCCGGGTTGATTGGCGGGCGCGTTGATCAGGCCACCATGCGCGTCg                             >  2:379177/1‑148 (MQ=255)
                                                                                                                      gATATTTTCTCGCGGCTGATGGCAGCGACCCGCGTGTCGCTCGGTTCGGTAATGGCCTGCCTGCTGCTGGTGCTGACATTAGGGCTGGTTATTGGCGGCAGCGCCGGGTTGATTGGCGGGCGCGTTGATCAGGCCACCATGCGCGTCg                             <  2:682635/148‑1 (MQ=255)
                                                                                                                         aatccctcGCGGCTGATGGCAGCGACCCGCGTGTCGCTCGGTTCGGTAATGGCCTGCCTGCTGCTGGTGCTGACATTAGGGCTGGTTATTGGCGGCAGCGCCGGGTTGATTGGCGGGCGCGTTGATCAGGCCACCATGCGCGTCGCCg                          <  1:445630/143‑1 (MQ=255)
                                                                                                                          ttttCTCGCGGCTGATGGCAGCGACCCGCGTGTCGCTCGGTTCGGTAATGGCCTGCCTGCTGCTGGTGCTGACATTAGGGCTGGTTATTGGCGGCAGCGCCGGGTTGATTGGCGGGCGCGTTGATCAGGCCACCATGCGCGTCGCCGat                        >  3:742704/1‑149 (MQ=255)
                                                                                                                             tctcGCGGCTGATGGCAGCGACCCGCGTGTCGCTCGGTTCGGTAATGGCCTGCCTGCTGCTGGTGCTGACATTAGGGCTGGTTATTGGCGGCAGCGCCGGGTTGATTGGCGGGCGCGTTGATCAGGCCACCATGCGCGTCGCCGATATg                     <  4:540773/149‑1 (MQ=255)
                                                                                                                               tcGCGGCTGATGGCAGCGACCCGCGTGTCGCTCGGTTCGGTAATGGCCTGCctgctgctg                                                                                                            >  4:234653/1‑60 (MQ=255)
                                                                                                                               tcGCGGCTGATGGCAGCGACCCGCGTGTCGCTCGGTTCGGTAATGGCCTGCctgctgctg                                                                                                            <  3:234653/60‑1 (MQ=255)
                                                                                                                                  cgGCTGATGGCAGCGACCCGCGTGTCGCTCGGTTCGGTAATGGCCTGCCTGCTGCTGGTGCTGACATTAGGGCTGGTTATTGGCGGCAGCGCCGGGTTGATTGGCGGGCGCGTTGATCAGGCCACCATGCGCGTCGCCGATATGTTTa                 <  1:2717/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)
                                                                                                                                                   |                                                                                                                                                 
GCTGATTGCGCTCACCAGCCAGTGGTGGCTGCCGTATGACCCACAGGCGATTGATTTGCCGTCGCGCCTGCTTTCGCCGGATGCGCAGCACTGGCTGGGCACCGATCACTTAGGTCGCGATATTTTCTCGCGGCTGATGGCAGCGACCCGCGTGTCGCTCGGTTCGGTAATGGCCTGCCTGCTGCTGGTGCTGACATTAGGGCTGGTTATTGGCGGCAGCGCCGGGTTGATTGGCGGGCGCGTTGATCAGGCCACCATGCGCGTCGCCGATATGTTTATGACCTTCCCGACCT  >  NC_000913/3616247‑3616539

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