Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,481,131 G→A L11L (CTG→TTG)  kefG ← glutathione‑regulated potassium‑efflux system ancillary protein KefG

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,481,1310GA97.7% 121.8 / ‑6.3 44L11L (CTG→TTG) kefGglutathione‑regulated potassium‑efflux system ancillary protein KefG
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base A (20/23);  total (20/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.35e-01

CAAAACCACGACTTAATACCCGGTCCAGCCACTCTTTCAGTAGCGCCGGGCAGCTATAGGTATAAAGAGGATGCTGAAAGACAATCACCTCGTGCTCGCGCAGTAATGCCTGCTCACGGGGGATATCAATAAAAAAATCGGGAT‑AGTGCGCGTAAAGGTCGTGCACGGTAACATTGCTGAGCTGCGTGGCCGGTTTAAGCAG‑TACCCGGTTTGCCACCGAGTCCTGAGATTCCGGATGGGCATACAGCAGCAAAACTTTCGCTGGCTGAGACATCATCCCCCTCCCGGTATGGTTTTTGTGTATAGTCGCTGTTTTGGGCTACCATTGCGCCCGGTGCGGCAGCTCGCCCATACATTACATTATCATAATGATAAGTTAACATAGTCTGAACATACGGCACCTTATGATTGTTTTCTCCTCGTTACAAATTCGTCGCGGCGTGCGCGTCCTGCTGGATAATGCCACCGCCACCATCAACCCTGGGCAGA  >  NC_000913/3480885‑3481373
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                  
cAAAACCACGACTTAATACCCGGTCCAGCCACTCTTTCAGTAGCGCCGGGCAGCTATAGGTATAAAGAGGATGCTGAAAGACAATCACCTCGTGCTCGCGCAGTAATGCCTGCTCACGGGGGATATCAATAAAAAAATCGGGAT‑AGTGCGCGTAAAGGTCGTGCACGGTAACATTGCTGAGCTGCGTGGCCGGTTTAAGCAG‑TACCCGGTTTGCCACCGAGGCCTGCACTTCCGGCTGGGCATACaacagc                                                                                                                                                                                                                                                >  2:113340/1‑251 (MQ=255)
   aaCCACGACTTAATACCCGGTCCAGCCACTCTTTCAGTAGCGCCGGGCAGCTATAGGTATAAAGAGGATGCTGAAAGACAATCACCTCGTGCTCGCGCAGTAATGCCTGCTCACGGGGGATATCAATAAAAAAATCGGGAT‑AGTGCGCGTAAAGGTCGTGCACGGTAACATTGCTGAGCTGCGTGGCCGGTTTAAGCAG‑TACCCGGTTTGCCACCGAGTCCTGAGATTCCGGATGGGCATACAACAGCaaa                                                                                                                                                                                                                                             <  1:478488/251‑1 (MQ=255)
            ttAATACCCGGTCCAGCCACTCTTTCAGTAGCGCCGGGCAGCTATAGGTATAAAGAGGATGCTGAAAGACAATCACCTCGTGCTCGCGCAGTAATGCCTGCTCACGGGGGATATCAATAAAAAAATCGGGAT‑AGTGCGCGTAAAGGTCGTGCACGGTAACATTGCTGAGCTGCGTGGCCGGTTTAAGCAG‑TACCCGGTTTGCCACCGAGTCCTGAGATTCCGGATGGGCATACAACAGCAAACCTTTCgct                                                                                                                                                                                                                                    >  2:91858/1‑251 (MQ=255)
              aaTACCCGGTCCAGCCACTCTTTCAGTAGCGCCGGGCAGCTATAGGTATAAAGAGGATGCTGAAAGACAATCACCTCGTGCTCGCGCAGTAATGCCTGCTCACGGGGGATATCAATAAAAAAATCGGGAT‑AGTGCGCGTAAAGGTCGTGCACGGTAACATTGCTGAGCTGCGTGGCCGGTTTAAGCAG‑TACCCGGTTTGCCACCGAGTCCTGAGATTCCGGATGGGCATACAACAGCAAAACTTTCgctgg                                                                                                                                                                                                                                  <  1:193207/251‑1 (MQ=255)
              aaTACCCGGTCCAGCCACTCTTTCAGTAGCGCCGGGCAGCTATAGGTATAAAGAGGATGCTGAAAGACAATCACCTCGTGCTCGCGCAGTAATGCCTGCTCACGGGGGATATCAATAAAAAAATCGGGAT‑AGTGCGCGTAAAGGTCGTGCACGGTAACATTGCTGAGCTGCGTGGCCGGTTTAAGCAG‑TACCCGGTTTGCCACCGAGTCCTGAGATTCCGGATGGGCATACAACAGCAAAACTTTCgctgg                                                                                                                                                                                                                                  >  1:333407/1‑251 (MQ=255)
                      gTCCAGCCACTCTTTCAGTAGCGCCGGGCAGCTATAGGTATAAAGAGGATGCTGAAAGACAATCACCTCGTGCTCGCGCAGTAATGCCTGCTCACGGGGGATATCAATAAAAAAATCGGGAT‑AGTGCGCGTAAAGGTCGTGCACGGTAACATTGCTGAGCTGCGTGGCCGGTTTAAGCAG‑TACCCGGTTTGCCACCGAGTCCTGAGATTCCGGATGGGCATACAACAGCAAAACTTTCGCTGGCTGAGAca                                                                                                                                                                                                                          <  1:131330/251‑1 (MQ=255)
                                 ctTTCAGTAGCGCCGGGCAGCTATAGGTATAAAGAGGATGCTGAAAGACAATCACCTCGTGCTCGCGCAGTAATGCCTGCTCACGGGGGATATCAATAAAAAAATCGGGAT‑AGTGCGCGTAAAGGTCGTGCACGGTAACATTGCTGAGCTGCGTGGCCGGTTTAAGCAG‑TACCCGGTTTGCCACCGAGTCCTGAGATTCCGGATGGGCATACAACAGCAAAACTTTCGCTGGCTGAGACATCATccccc                                                                                                                                                                                                                 <  1:381094/249‑1 (MQ=255)
                                  tttCAGTAGCGCCGGGCAGCTATAGGTATAACGCGGATGCTGAAAGACAATCACCTCGTGCTCGCGCAGTAATGCCTGCTCACGGGGGATATCAATAAACAAATCGGGAT‑AGTGCGCGTAAAGGTCGGGCACGGTAACATTGCTGAGCTGCGTGGCCGGTTTAAGCAG‑TACCCGGTTTGCCACCGAGTCCTAAGATTCCGGATGGGCATCCAACAGCAAAACTTTCGCTGGCTGAGACATCCTCCCCCTcc                                                                                                                                                                                                              >  2:32911/1‑251 (MQ=255)
                                        tAGCGCCGGGCAGCTATAGGTATAAAGAGGATGCTGAAAGACAATCACCTCGTGCTCGCGCAGTAATGCCTGCTCACGGGGGATATCAATAAAAAAATCGGGAT‑AGTGCGCGTAAAGGTCGTGCACGGTAACATTGCTGAGCTGCGTGGCCGGTTTAAGCAG‑TACCCGGTTTGCCACCGAGTCCTGAGATTCCGGATGGGCATACAACAGCAAAACTTTCGCTGGCTGAGACATCATCCCCCTCCCGGTAt                                                                                                                                                                                                        >  1:36885/1‑251 (MQ=255)
                                              cGGGCAGCTATAGGTATCAAGAGGATGCTGAAAGACAATCACCTCGTGCTCGCGCAGTAATGCCTGCTCACGGGGGATATCAATAAAAAAATCGGGAT‑AGTGCGCGTAAAGGTCGTGCACGGTAACATTGCTGAGCTGCGTGGCCGGTTTAAGCAG‑TACCCGGTTTGCCACCGAGTCCTGAGATTCCGGATGGGCATACAACAGCAAAACTTTCGCTGGCTGAGACATCATCCCCCTCCCGGTATGGtttt                                                                                                                                                                                                  >  1:365055/1‑251 (MQ=255)
                                               gggCAGCTATAGGTATAAAGAGGATGCTGAAAGACAATCACCTCGTGCTCGCGCAGTAATGCCTGCTCACGGGGGATATCAATAAAAAAATCGGGAT‑AGTGCGCGTAAAGGTCGTGCACGGTAACATTGCTGAGCTGCGTGGCCGGTTTAAGCAG‑TACCCGGTTTGCCACCGAGTCCTGAGATTCCGGATGGGCATACAACAGCAAAACTTTCGCTGGCTGAGACATCATCCCCCTCCCGGTATGGttttt                                                                                                                                                                                                 <  1:17078/251‑1 (MQ=255)
                                               gggCAGCTATAGGTATAAAGAGGATGCTGAAAGACAATCACCTCGTGCTCGCGCAGTAATGCCTGCTCACGGGGGATATCAATAAAAAAATCGGGAT‑AGTGCGCGTAAAGGTCGTGCACGGTAACATTGCTGAGCTGCGTGGCCGGTTTAAGCAG‑TACCCGGTTTGCCACCACGTCCTGAGATTCCGGCTGGGCATACAACAGCAAAACTTTCGCTGGCTGAGACATCATCCCCCTCCCGGTATGGatttt                                                                                                                                                                                                 >  2:350251/1‑251 (MQ=255)
                                                    gCTATAGGTATAAAGAGTATGCTGAAAGACAATCACCTCGTGCTCGCGCAGTAAGGCCTGCGACCGGGGGATATCAATAAAAAAATCGTGAT‑AGGGCGCGTAAAGGTCGTGCACGGTAACATTGCTGAGCTGCGTGGCCGGTTTAAGCAG‑TACCCGGTTGGCCACCGAGTCCTGAGATTCCGGATGGGCATACAGCAGCAAAACTTGCGCTGGTTGAGACATCATCCCCCTCCCGGTATGGTTGTTGTGta                                                                                                                                                                                            <  2:260685/251‑1 (MQ=255)
                                                         aGGTATAAAGAGGATGCTGAAAGACAATCACCTCGTGCTCGCGCAGTAATGCCTGCTCACGGGGGATATCAATAAAAAAATCGGGAT‑AGTGCGCGTAAAGGTCGTGCACGGTAACATTGCTGAGCTGCGTGGCCGGTTTAAGCAG‑TACCCGGTTTGCCACCGAGTCCTGAGATTCCGGATGGGCATACAACAGCAAAACTTTCGCTGGCTGAGACATCATCCCCCTCCCGGTATGGTTTTTGTGTATAGt                                                                                                                                                                                        <  1:91858/250‑1 (MQ=255)
                                                                aaGAGGATGCGGAAAGACAATCACCTCGTTCTCGCGCAGTAATGCCTGCGCACGGGGGATATCAATAAAAAAAGCGGGAT‑AGTGCGCGTAAAGGTCGTGCACGGTAACATTGCTGAGCTGCGTGGCCGGTTTAAGCAG‑TACCCGGTTTGCCACCGAGTCCTGAGATTCCGGATGGGCATACAACAGCAAAACTTTCGCTGGCTGAGACATCATCCCCCTCCCGGTATGGTTTTTGTGTATAGTCGCTGttt                                                                                                                                                                                <  2:506156/251‑1 (MQ=255)
                                                                                         ccGTGCTCGCGCAGGAATGCCTGCTCACGGGGGATATCAATAAAAAAATCGGGAT‑AGTGCGCGTAAAGGTCGTGCACGGTAACATTGCTGAGCTGCGTGGCCGGTTTAAGCAG‑TACCCGGTTTGCCACCGAGTCCTGAGATTCCGGATGGGCATACAACAGCAAAACTTTCGCTGGCTGAGACATCATCCCCCTCCCGGTATGGTTTTTGTGTATAGTCGCTGTTTTGGGCTACCATTGCGCCCGGTGCgg                                                                                                                                                       <  2:365055/250‑1 (MQ=255)
                                                                                             gCTCGCGCAGTAATGCCTGCTCACGGGGGATATCAATAAAAAAATCGGGAT‑AGTGCGCGTAAAGGTCGTGCACGGTAACATTGCTGAGCTGCGTGGCCGGTTTAAGCAG‑TACCCGGTTTGCCACCGAGTCCTGAGATTCCGGATGGGCATACAACAGCAAAACTTTCGCTGGCTGAGACATCCTCCCCCTCCCGGTCTGGTTTTTGTGTATAGTCGCTGTTTTGGGCTACCATTGCGCCCGGTGCGGCAGc                                                                                                                                                   >  2:447226/1‑251 (MQ=255)
                                                                                                     agtaattcccgctctcgtgggagtacAATAAAAAAAGCGGGATCAG‑GTGCGTAAAGGTATTGTACGGTAACAGTGCTGGGCTGCGTTGCCGGGTTAAGCAT‑TTCCCGTTTTGCCACCGAGACCTGTGATTCCGGATGGGCATACAACAGCAAAAATGTCGCTGGCGGAGACAGCATCCCCCTCCCGGTATGGTTGGTGTGTATAGTCGCTGTTTTGGGCTACCATTGCGCCCAGTTCGGCAGCTCGCCCAt                                                                                                                                           <  2:377486/226‑1 (MQ=255)
                                                                                                        aaTGCCTGCTCACGGGGGATATCAATAAAAAAATCGGGAT‑AGTGCGCGTAAAGGTCGTGCACGGTAACATTGCTGAGCGGCGTGGCCGGTTTAAGCAG‑TACCCGGTTTGCCACCGAGTCCTGAGATTCCGGATGGGCATACAACAGCAAAACTTTCGCTGGCTGAGACATCATCCCCCTCCCGGTATGGTTTTTGTGTATAGTCGCTGTTTTGGGCTACCATTGCGCCCGGTGCGGCAGCTCGCCCAtac                                                                                                                                         <  2:41843/250‑1 (MQ=255)
                                                                                                          tGCCTGCTCACGGGGGATATCAATAAAAAAGTCTGGAT‑AGTTAGCGTCAAGGTCGTGCACGGTAACATTGCTGAGCTGCGTGGCCGGTTTAAGCAG‑TACCCGGTTTGCCACCGAGTCCTGCGATTCCGGAGGGGCATACAACAGCAAAACTTTCGCTGGCTGAGAAATCATCCCCCTCCCGGTATGAGTTTTGTGGATAGTCGCTGTTTTGGGCTACCCTTGCTCCCGGTGCGGCAGCACGGCCAtacatt                                                                                                                                      >  2:132157/1‑251 (MQ=255)
                                                                                                           gCCTGCTCACGGGGGATATCAATAAAAAAAGCGGGAT‑AGTGCGCGTAAAGGTCGTGCACGGTAACATTGCTGAGCTGCGTGGCCGGTTTAAGCAG‑TACCCGGTTTGCCACCGAGTCCTGAGATTCCGGATGGGCATACAACAGCAAAACTTTCGCTGGCTGAGACATCATCCCCCTCCCGGTATGGTTTTTGTGTATAGTCGCTGTTTTGGGCTACCATTGCTCCCGGTGCGGCAGCTCGCCCAtacatt                                                                                                                                      <  1:429501/250‑1 (MQ=255)
                                                                                                                     gggggAGATCAAGAAAAAAATCGGGAT‑AGGGCGCGTAAATGTCGTGCACGGTAACATTGCTGAGCTGCGTGGCCGGTTTAAGCAG‑TACCCGGTTTGCCACCGAGTCCTGAGATTCCGGATGGGCATACAACAGCAAAACTTTCGCTGGCTGAGACATCATCCCCCTCCCGGTATGGTTTTTGTGTATAGGCGCTGTTTTGGGCTACCATTGCGCCCGGTGCGGCAGCTCGCCCATACATTACATTATCATa                                                                                                                           <  2:165788/251‑1 (MQ=255)
                                                                                                                        gtatatCAATAAAAAAATCGGGAT‑AGTGCGCGTAAAGGTCGTGCACGGTAACATTGCTGAGCTGCGTGGCCGGTTTAAGCAG‑TACCCGGTTTGCCACCGAGTCCTGAGATTCCGGATGGGCATACAACAGCAAAACTTTCGCTGGCTGAGACATCATCCCCCTCCCGGTATGGTTTTTGTGTATAGTCGCTGTTTTGGGCTACCATTGCGCCCGGTGCGGCAGCTCGCCCATACATTACATTATCATAATg                                                                                                                        <  1:483723/249‑1 (MQ=255)
                                                                                                                         gATATCAATAAAAAAATCGGGAT‑AGTGCGCGTAAAGGTCGTGCACGGTAACATTGCTGAGCTGCGTGGCCGGTTTAAGCAG‑TACCCGGTTTGCCACCGAGTCCTGAGATTCCGGATGGGCATACAACAGCAAAACTTTCGCTGGCTGAGACATCATCCCCCTCCCGGTATGGTTTTTGTGTATAGTCGCTGTTTTGGGCTACCCTTGCGCCCGGTGCGGCAGCGCGCCCATACATTACATTATCATAATg                                                                                                                        >  2:119282/1‑250 (MQ=255)
                                                                                                                               aaTAAAAAAATCGGGAT‑AGTGCGCGTAAAGGTCGTGCACGGTAACATTGCTGAGCTGCGTGGCCGGTTTAAGCAG‑TACCCGGTTTGCCACCGAGTCCTGAGATTCCGGATGGGCATACAACAGCAAAACTTTCGCTGGCTGAGACATCATCCCCCTCCCGGTATGGTTTTTGTGTATAGTCGCTGTTTTGGGCTACCATTGCGCCCGGTGCGGCAGCTCGCCCATACATTACATTATCATAATGATAAGtt                                                                                                                 <  1:85259/251‑1 (MQ=255)
                                                                                                                                   aaaaaaTCGGGAT‑AGTGCGCGTAAAGGTCGTGCACGGTAACATTGCTGAGCTGCGTGGCCGGTTTAAGCAG‑TACCCGGTTTGACACCGAGTCCTGAGATTCCGGATGGGCATACAACAGCAAAACTTTCGCTGGCTGAGACATCATCCCCCTCCCGGTATGGTTTTTGTGTAAAGTCGCTGTTTTGGGCTACCATTGCGCCCGGTGCGGCAGCTCGCCCATACATTACATTATCATAATGATAAGTTAAc                                                                                                              >  2:65111/1‑250 (MQ=255)
                                                                                                                                       aaTCGGGAT‑AGGGCGCGTAAAGGCCGTGCATGGTAACATTGCTGAGCTGCGTGGCCGGTGTAAGCAG‑TACCCGGATTGCCACCGAGTCCTGAGATGCCGGATGGGCATACAACAGCAAAACTTTCGCTGGCTGAGACATCATCCCCCTCCCGGTATGGTTTTTGTGTATAGTCGCTGTTTTGGGCTACCATTGCGCCCGGTGCGGCAGCACGCCCATACATTACATTATCATAATGATAAGTTAACATAGt                                                                                                         <  2:139181/251‑1 (MQ=255)
                                                                                                                                                 ggggtgcgTAAAGGTCGTCCACGGAAACATTGCTGGGCTCCGTGGCCGGTTTAAGCAT‑TACCCGGTTTCCCCCCGA‑TCCTGAGATTCCGGATGGGCATACAACAGCAAAACTTTCGCTGGCTGAGACATCACCCCCCTCCCGGTATGGTTTTTGTGTATAGTCGCTGTTTTGGGCTACCATTGCGCCCGGTGCGGCAGCTCGCCCATACATTACATTATCATAATGATAAGTTAACATAGTCTGAACATAc                                                                                               <  2:389841/246‑1 (MQ=255)
                                                                                                                                                   tGCGCGTAAAGGTCGTGCACGGTAACATTGCTGAGCTGCGTGGCCGGTTTAAGCAG‑TACCCGGTTTGCCACCGAGTCCTGAGATTCCGGATGGGCATACAACAGCAAAACTTTCGCTGGCTGAGACATCATCCCCCTCCCGGTATGGTTTTTGTGTATAGTCGCTGTTTTGGGCTACCATTGCGCCCGGTGCGGCAGCTCGCCCATACATTACATTCTCATAATGATAAGTTAACATAGTCTGAACAtccg                                                                                              >  2:251454/1‑248 (MQ=255)
                                                                                                                                                      gcgTAAAGGTCGTGCACGGTAACATTGCTGAGCTGCGTGGCCGGTTTAAGCAG‑TACCCGGTTTGCCACCGAGTCCTGAGATTCCGGATGGGCATACAACAGCAAAACTTTCGCTGGCTGAGACATCATCCCCCTCCCGGTATGGTTTTTGTGTATAGTCGCTGTTTTGGGCTACCATTGCGCCCGGTGCGGCAGCTCGCCCATACATTACATTATCATAATGAGAAGTTAACATAGTCTGAACATACGGCa                                                                                           <  1:434634/251‑1 (MQ=255)
                                                                                                                                                                      cGGTAACATTGCTGAGCTGCGTGGCCGGTTTAAGCAG‑TACCCGGTTTGCCACCGAGTCCTGAGATTCCGGATGGGCATACAACAGCAAAACTTTCGCTGGCTGAGACATCATCCCCCTCCCGGTATGGTTTTTGTGTATAGTCGCTGTTTTGGGCTACCATTGCGCCCGGTGCGGCAGCTCGCCCATACATTACATTATCATAATGATAAGTTAACATAGTCTGAACATACGGCACCTTATGATTGTTTtc                                                                           >  1:66287/1‑251 (MQ=255)
                                                                                                                                                                        gTAACAGTGCTGTGCTGCGTGGCCGGTTTAAGCAG‑TACCCGGTTTGCCCCCGAGTCCTGAGAGTTCGTATGGGCATACAACAGCAAAACTTTCGCTGGCTGAGACATCATCCCCCTCCCGGTATGGTTTTTGTGTATAGTCGCTTTTTTGGGCTACCATTGCGCCCGGTGCGGCAGCTCGTCCATACATTACATTATCATAATGATAAGTTAACATAGTCTGAACATACGGCACCTTATGATTGTTTtctc                                                                         <  1:314129/251‑1 (MQ=255)
                                                                                                                                                                            aaTTGCGGAGCTGCGGGGCCGGTTTAAGCAGGTACC‑GGTTTGCACCCGAGTCCTGAGATTCCGGATGGGAATACAACAACAAAACGTGCGCTGGCTGAGACATCATACCCCTACCGGTATGGTTTTTGTGTATAGTCGCTGTTTTGGGCTAAAATTGCGCCCGGTGCGGCAGCTCGCCCATACATTACATTATCATAATGATAAGTTAACATAGTATGAACATACGGCCACTTATGATGGTTTTctccgcg                                                                     <  2:174286/250‑4 (MQ=255)
                                                                                                                                                                                    aGCTTCGTGGCCGGTTTAAGCAG‑TACCCGGTTTGCCACCGAGTCCTGAGATTCCGGATGGGCATACAACAGCAAAACTTTCGCTGGCTGAGTCATCATCCCCCTCCCGGTATGGTTTTTGTGTATAGTCGCTGTTTTGGGCTACCATTGCGCCCGGTGCGGCAGCTCGCCCATACATTACATTATCATAATGATAAGTTAACATAGTCTGAACATACGGCTCCTTATGATTTTTTTCTCCTCGTTACAAAt                                                             <  1:396181/251‑1 (MQ=255)
                                                                                                                                                                                      cggCGTTGCCGGTTTAAGCAG‑TACCCGGTTTGCCACCGAGTCCTGAGATTCCGGATGGGCATACAACAGCAAAACTTTCGCTGGCTGAGACATCATCCCCCTCCCGGTATGGTTTTTGTGTATAGTCGCTGTTTTGGGCTACCATTGCGCCCGGTGCGGCAGCTCGCCCATACATTACATTATCATAATTATAAGTTAACATAGTCTGAACATACGGCACCTTATGATTGTGTTCTCCTCTTTACAAATtc                                                           <  2:331187/249‑1 (MQ=255)
                                                                                                                                                                                         cGTGGCCGGTTTAAGCAG‑TACCCGGTTTGCCACCGAGTCCTGAGATTCCGGATGGGCATACAACAGCAAAACTTTCGCTGGCTGAGACATCATCCCCCTCCCGGTATGGTTTTTGTGTATAGTCGCTGTTTTGGGCTACCATTGCGCCCGGTGCGGCAGCTCGCCCATACATTACATTATCATAATGATAAGTTAACATAGTCTGAACATACGGCACCTTATGATTGTTTTCTCCTCGTTACAAATtcgtc                                                        >  1:144420/1‑251 (MQ=255)
                                                                                                                                                                                         cGTGGCCGGTTTAAGCAG‑TACCCGGTTTGACACCGAGTCCTGAGATTCCGGATGGGCATACAACAGCAAAACTTTCGCTGGCTGAGACATCATCCCCCTCCCGGTATGGTTTTTGTGTAAAGTCGCTGTTTTGGGCTACCATTGCGCCCGGTGCGGCAGCTCGCCCATACATTACATTATCATAATGATAAGTTAACATAGTCTGAACATACGGCACCTTATGATTGTTTTCTCCTCGTTACAAATtcgtc                                                        <  1:65111/251‑1 (MQ=255)
                                                                                                                                                                                                 gTTTAAGCAG‑TACCCGGTTTGCCACCGAGTCCTGAGATTCCGGATGGGCATACAACAGCAAAACTTTCGCTGGCTGAGACATCATCCCCCTCCCGGTATGGTTTTTGTGTATAGTCGCTGTTTTGGGCTACCATTGCGCCCGGTGCGGCAGCTCGCCCATACATTACATTATCATAATGATAAGTTAACATAGTCTGAACATACGGCACCTTATGATTGTTTTCTCCTCGTTACAAATTCGTCGCGGCGTg                                                >  1:77015/1‑251 (MQ=255)
                                                                                                                                                                                                       gCAG‑TACCCGGTTTGCCACCGAGTCCTGAGATTCCGGATGGGCATACAACAGCAAAACTTTCGCTGGCTGAGACATCATCCCCCTCCCGGTATGGTTTTTGTGTATAGTCGCTGTTTTGGGCTACCATTGCGCCCGGTGCGGCAGCTCGCCCATACATTACATTATCATAATGATAAGTTAACATAGTCTGAACATACGGCCCCTTATGATTGTTTTCTCCTCGTTACCAATTCGTCTCGGCGTGCGCGTc                                          >  2:138936/1‑251 (MQ=255)
                                                                                                                                                                                                                   tttGCCACCGAGTCCTGAGATTCCGGATGGGCATACAACAGCAAAACTTTCGCTGGCTGAGACATCATCCCCCTCCCGGTATGGTTTTTGTGTATAGTCGCTGTTTTGGGCTACCATTGCGCCCGGTGCGGCAGCTCGCCCATACATTACATTATCATAATGATAAGTTAACATAGTCTGAACATACGGCACCTTATGATTGTTTTCTCCTCGTTACAAATTAGTCGCGGCGTGCGCGTCCTGCTGGATa                                >  1:478331/1‑250 (MQ=255)
                                                                                                                                                                                                                           cGAGTCCTGAGATTCCGGATGGGCATTCAACAGCAAAACTTTCGCTGGCTGAGACATCATCCCCCTCCCGGTATTGTTTTTGTGTATAGTCGCTGTTTTGGGCTACCATTGCGCCCGGTGCGGCAGCTCGCCCATACATTACATTATCATAATGATAAGTTAACATAGTCTGAACATACGGCACCTTATGATTGTTTTCTCCCCGTTACAAATTCGTCGCGGCGTGCGCGTCCTGCTGGATAATgccaccg                       <  2:252652/251‑1 (MQ=255)
                                                                                                                                                                                                                               tCCTGAGATTCCGGATGGGCATACAACAGCAAAACTTTCGCTGGCTGAGACATCATCCCCCTCCCGGTATGGTTTTTGTGTATAGTCGCTGTTTTGGGCTACCATTGCGCCCGGTGCGGCAGCTCGCCCATACATTACATTATCATAATGATAAGTTAACATAGTCTGAACATACGGCACCTTATGATTGTTTTCTCCTCGTTACAAATTCGTCGCGACGTGCGCGTCCTGCTGGATAATgccaccgccac                   >  1:487252/1‑251 (MQ=255)
                                                                                                                                                                                                                                           ggATGGGCATACAACAGCAAAACTTTCGCTGGCTGAGACATCATCCCCCTCCCGGTATGGTTTTTGTGTATAGTCGCTGTTTTGGGCTACCATTGCGCCCGGTGCGGCAGCTCGCCCATACATTACATTATCATAATGATAAGTTAACATAGTCTGAACATACGGCACCTTATGATTGTTTTCTCCTCGTTACAAATTCGTCGCGGCGTGCGCGTCCTGCTGGATAATGCCACCGCCACCATCAACCCTgg       >  1:88113/1‑251 (MQ=255)
                                                                                                                                                                                                                                             aTGGGCATACAACAGCAAAACTTTCGCTGGCTGAGACATCATCCCCCTCCCGGTATGGTTTTTGTGTATAGTCGCTGTTTTGGGCTACCATTGCGCCCGGTGCGGCAGCTCGCCCATACATTACATTATCATAATGATAAGTTAACATAGTCTGAACATACGGCACCTTATGATTGTTTTCTCCTCGTTACAAATTCGTCGCGGCGTGCGCGTCCTGCTGGATAATGCCACCGCCACCATCAACCCTGGGc     >  1:221646/1‑251 (MQ=255)
                                                                                                                                                                                                                                                 gCATACAACAGCAAAACTTTCGCTGGCTGAGACATCATCCCCCTCCCGGTATGGTTTTTGTGTATAGTCGCTGTTTTGGGCTACCATTGCGCCCGGTGCGGCAGCTCGCCCATACATTACATTATCATAATGATAAGTTAACATAGTCTGAACATACGGCACCTTATGATTGTTTTCTCCTCGTTACAAATTCGTCGCGGCGTGCGCGTCCTGCTGGATAATGCCACCGCCACCATCAACCCTGGGCAGa  >  1:147125/1‑250 (MQ=255)
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                  
CAAAACCACGACTTAATACCCGGTCCAGCCACTCTTTCAGTAGCGCCGGGCAGCTATAGGTATAAAGAGGATGCTGAAAGACAATCACCTCGTGCTCGCGCAGTAATGCCTGCTCACGGGGGATATCAATAAAAAAATCGGGAT‑AGTGCGCGTAAAGGTCGTGCACGGTAACATTGCTGAGCTGCGTGGCCGGTTTAAGCAG‑TACCCGGTTTGCCACCGAGTCCTGAGATTCCGGATGGGCATACAGCAGCAAAACTTTCGCTGGCTGAGACATCATCCCCCTCCCGGTATGGTTTTTGTGTATAGTCGCTGTTTTGGGCTACCATTGCGCCCGGTGCGGCAGCTCGCCCATACATTACATTATCATAATGATAAGTTAACATAGTCTGAACATACGGCACCTTATGATTGTTTTCTCCTCGTTACAAATTCGTCGCGGCGTGCGCGTCCTGCTGGATAATGCCACCGCCACCATCAACCCTGGGCAGA  >  NC_000913/3480885‑3481373

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑