Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,662,350 C→T P599S (CCC→TCC)  ynfF → S‑ and N‑oxide reductase, A subunit, periplasmic

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,662,3500CT98.9% 285.9 / ‑4.3 87P599S (CCC→TCC) ynfFS‑ and N‑oxide reductase, A subunit, periplasmic
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base T (37/49);  total (38/49)
Fisher's exact test for biased strand distribution p-value = 4.37e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.34e-01

TGAACGCAAACCCATCTACTGGATGCTGAGTGAAGTCGCCAAACGCTTAGGACCAGACGTTTATCAAACCTTTACTGAAGGTCGCAGTCAGCATGAATGGATCAAATATCTCCATGCGAAAACGAAGGAACGTAACCCTGAGATGCCCGACTACGAGGAGATGAAAACGACCGGGATCTTTAAGAAAAAATGCCCGGAAGAACACTACGTCGCTTTCCGCGCATTCCGTGAAGATCCACAGGCAAACCCGTTGAAAACACCTTCGGGGAAAATCGAAATTTATTCTGAA  >  NC_000913/1662205‑1662493
                                                                                                                                                 |                                                                                                                                               
tGAACGCAAACCCATCTACTGGATGCTGAGTGAAGTCGCCAAACGCTTAGGACCACACGTTTATCAAACCTTTACCGAAGGTCGCAGTCAGCATGAATGGATCAAATATCTCCATGCGAAAACGAAGGAACGTAACCCTGAGATGtcc                                                                                                                                               <  5:133364/148‑1 (MQ=255)
  aaCGCAAACCCATCTACTGGATGCTGAGTGAAGTCGCCAAACGCTTAGGACCAGACGTTTATCAAACCTTTACTGAAGGTCGCAGTCAGCATGAATGGATCAAATATCTCCATGCGAAAACGAAGGAACGTAACCCTGAGATGTCCGAc                                                                                                                                            <  5:378673/149‑1 (MQ=255)
     gCAAACCCATCTACTGGATGCTGAGTGAAGTCGCCAAACGCTTAGGACCAGACGTTTATCAAACCTTTACTGAAGGTCGCAGTCAGCATGAATGGATCAAATATCTCCATGCGAAAACGAAGGAACGTAACCCTGAGATGTCCGACTAc                                                                                                                                         >  4:202601/1‑149 (MQ=255)
     gCAAACCCATCTACTGGATGCTGAGTGAAGTCGCCAAACGCTTAGGACCAGACGTTTATCAAACCTTTACTGAAGGTCGCAGTCAGCATGAATGGATCAAATATCTCCATGCGAAAACGAAGGAACGTAACCCTGAGATGTCCGACTAc                                                                                                                                         <  7:105634/149‑1 (MQ=255)
       aaaCCCATCTACTGGATGCTGAGTGAAGTCGCCAAACGCTTAGGACCAGACGTTTATCAAACCTTTACTGAAGGTCGCAGTCAGCATGAATGGATCAAATATCTCCATGCGAAAACGAAGGAACGTAACCCTGAGATGTCCGACTAc                                                                                                                                         >  1:412263/1‑147 (MQ=255)
        aaCCCATCTACTTGATGCTGAGTGAAGTCGCCAAACGCTTAGGACCAGACGTTTATCAAACCTTTACTGAAGGTCGCAGTCAGCATGAATGGATCAAATATCTCCATGCGAAAACGAAGGAACGTAACCCTGAGATGTCCGACTACgag                                                                                                                                      <  4:210002/149‑1 (MQ=255)
        aaCCCATCTACTGGATGCTGAGTGAAGTCGCCAAACGCTTAGGACCAGACGTTTATCAAACCTTTACTGAAGGTCGCAGTTAGCATGAATGGATCAAATATCTCCATGCGAAAACGAAGGAACGTAACCCTGAGATGTCCGACTACgag                                                                                                                                      <  7:127446/149‑1 (MQ=255)
           ccATCTACTGGATGCTGAGTGAAGTCGCCAAACGCTTAGGACCAGACGTTTATCAAACCTTTACTGAAGGTCGCAGTCAGCATGAATGGATCAAATATCTCCATGCGAAAACGAAGGAACGTAACCCTGAGATGTCCGACTACga                                                                                                                                       >  2:264136/1‑145 (MQ=255)
           ccATCTACTGGATGCTGAGTGAAGTCGCCAAACGCTTAGGACCAGACGTTTATCAAACCTTTACTGAAGGTCGCAGTCAGCATGAATGGATCAAATATCTCCATGCGAAAACGAAGGAACGTAACCCTGAGATGTCCGACTACga                                                                                                                                       <  1:264136/145‑1 (MQ=255)
            cATCTACTGGATGCTGAGTGAAGTCGCCAAACGCTTAGGACCAGACGTTTATCAAACCTTTACTGAAGGTCGCAGTCAGCATGAATGGATCAAATATCTCCATGCGAAAACGAAGGAACGTAACCCTGAGATGTCCGACTACGAGgaga                                                                                                                                  <  2:352594/149‑1 (MQ=255)
            cATCTACTGGATGCTGAGTGAAGTCGCCAAACGCTTAGGACCAGACGTTTATCAAACCTTTACTGAAGGTCGCAGCCAGCATGAATGGATCAAATATCTCCATGCGAAAACGAAGGAACGTAACCCTGAGATGTCCGACTACgag                                                                                                                                      >  7:324237/1‑145 (MQ=255)
             aTCTACTGGATGCTGAGTGAAGTCGCCAAACGCTTAGGACCAGACGTTTATCAAACCTTTACTGAAGGTCGCAGTCAGCATGAATGGATCAAATATCTCCATGCGAAAACGAAGGAACGTAACCCTGAGATGTCCGACTACGAGGAGAt                                                                                                                                 >  1:237648/1‑149 (MQ=255)
              tcaaCTGGATGCTGAGTGAAGTCGCCAAACGCTTAGGACCAGACGTTTATCAAACCTTTACTGAAGGTCGCAGTCAGCATGAATGGATCAAATATCTCAATGCGAAAACGAAGGAACGTAACCCTGAGATGTCCGACTACGAGGAGAt                                                                                                                                 >  8:409760/4‑148 (MQ=255)
                      aTGCTGAGTGAAGTCGCCAAACGCTTAGGACCAGACGTTTATCAAACCTTTACTGACGGTCGCAGTCAGCATGAATGGATCAAATATCTCCATGCGAAAACGAAGGAACGTAACCCTGAGATGTCCGACTACGAGGAGATGAAAACGAc                                                                                                                        <  6:86120/149‑1 (MQ=255)
                          tGAGTGAAGTCGCCAAACGCTTAGGACCAGACGTTTATCAAACCTTTACTGAAGGTCGCAGTCAGCATGAATGGATCAAATATCTCCATGCGAAAACGAAGGAACGTAACCCTGAGATGTCCGACTACGAGGAGATGAAAACGACCgg                                                                                                                     >  3:116670/1‑148 (MQ=255)
                          tGAGTGAAGTCGCCAAACGCTTAGGACCAGACGTTTATCAAACCTTTACTGAAGGTCGCAGTCAGCATGAATGGATCAAATATCTCCATGCGAAAACGAAGGAACGTAACCCTGAGATGTCCGACTACGAGGAGATGAAAACGACCgg                                                                                                                     >  4:315541/1‑148 (MQ=255)
                            aGTGAAGTCGCCAAACGCTTAGGACCAGACGTTTATCAAACCTTTACTGAAGGTCGCAGTCAGCATGAATGGATCAAATATCTCCATGCGAAAACGAAGGAACGTAACCCTGAGATGTCCGACTACGAGGAGATGAAAACGACCGGGa                                                                                                                   <  3:202601/148‑1 (MQ=255)
                             gTGAAGTCGCCAAACGCTTAGGACCAGACGTTTATCAAACCTTTACTGAAGGTCGCAGTCAGCATGAATGGATCAAATATCTCCATGCGAAAACGAAGGAACGTAACCCTGAGATGTCCGACTACGAGGAGATGAAAACGACCGGGATc                                                                                                                 <  4:131410/149‑1 (MQ=255)
                                aaGTCGCCAAACGCTTAGGACCAGACGTTTATCAAACCTTTACTGAAGGTCGCAGTCAGCATGAATGGATCAAATATCTCCATGCGAAAACGAAGGATCTTAACCCTGAGATGTCCGACTACGAGGAGATGAAAACGACCGGGATCtt                                                                                                               <  5:394531/148‑1 (MQ=255)
                                  gTCGCCAAACGCTTAGGACCAGACGTTTATCAAACCTTTACTGAAGGTCGCAGTCAGCATGAATGGATCAAATATCTCCATGCGAAAACGAAGGAACGTAACCCTGAGATGTCCGACTACGAGGAGATGAAAACGACCGGGATCttt                                                                                                              >  1:28181/1‑147 (MQ=255)
                                       cAAACGCTTAGGACCAGACGTTTATCAAACCTTTACTGAAGGTCGCAGTCAGCATGAATGGATCAAATATCTCCATGCGAAAACGAAGGAACGTAACCCTGAGATGTCCGACTACGAGGAGATGAAAACGACCGGGATCTTTa                                                                                                             <  7:337107/143‑1 (MQ=255)
                                       cAAACGCTTAGGACCAGACGTTTATCAAACCTTTACTGAAGGTCGCAGTCAGCATGAATGGATCAAATATCTCCATGCGAAAACGAAGGAACGTAACCCTGAGATGTCCGACTACGAGGAGATGAAAACGACCGGGATCTTTa                                                                                                             >  8:337107/1‑143 (MQ=255)
                                       cAAACGCTTAGGACCAGACGTTTATCAAACCTTTACTGAAGGTCGCAGTCAGCATGAATGGATCAAATATCTCCATGCGAAAACGAAGGAACGTAACCCTGAGATGTCCGACTACGAGGAGATGAAAACGACCGGGATCTTTAAGaaa                                                                                                        >  6:244738/1‑148 (MQ=255)
                                        aaaCGCTTAGTACCAGACGTTTATCAAACCTTTACTGAAGGTCGCAGTCAGCATGAATGGATCAAATATCTCCATGCGAAAACGATGGAACGTAACCCTGAGATGTCCGACTACGAGGAGATGAAAACGACCGGGATCTTTAAGaaaaa                                                                                                      <  3:312959/149‑1 (MQ=255)
                                            gCTTAGGACCAGACGTTTATCAAACCTTTACTGAAGGTCGCAGTCAGCATGAATGGATCAAATATCTCCCTGCGAAAACGAAGGAACGTAACCCTGAGATGTCCGACTACGAGGAGATGAAAACGACCGGGATCTTTAAGAAAAAATGc                                                                                                  <  4:192747/149‑1 (MQ=255)
                                                aGGACCAGACGTTTATCAAACCTTTACTGAAGGTCGCAGTCAGCATGAATGGATCAAATATCTCCATGCGAAAACGAAGGAACGTAACCCTGAGATGTCCGACTACGAGGAGATGAAAACGACCGGGAGCTTTAAGAAAAAATGccc                                                                                                <  2:237648/147‑1 (MQ=255)
                                                aGGACCAGACGTTTATCAAACCTTTACCGAAGGGCGCAGTCAGCATGAATGGATCAAATATCTCCATGCGAAAACGAAGGAACGTAACCCTGAGATGTCCGACTACGAGGAGATGAAAACGACCGGGATCTTTAAGAAAAAATGCCCg                                                                                               >  4:13286/1‑148 (MQ=255)
                                                  gACCAGACGTTTATCAAACCTTTACTGAAGGTCGCAGTCAGCATGAATGGATCAAATATCTCCATGCGAAAACGAAGGAACGTAACCCTGAGATGTCCGACTACGAGGAGATGAAAACGACCGGGATCTTTAAGAAAAAATGCCCgg                                                                                              >  1:15084/1‑147 (MQ=255)
                                                    ccAGACGTTTATCAAACCTTTACTGAAGGTCGCAGTCAGCATGAATGGATCAAATATCTCCATGCGAAAACGAAGGAACGTAACCCTGAGATGTCCGACt                                                                                                                                           <  1:141741/100‑1 (MQ=255)
                                                    ccAGACGTTTATCAAACCTTTACTGAAGGTCGCAGTCAGCATGAATGGATCAAATATCTCCATGCGAAAACGAAGGAACGTAACCCTGAGATGTCCGACt                                                                                                                                           >  2:141741/1‑100 (MQ=255)
                                                     cAGACGTTTATCAAACCTTTACTGAAGGTCGCAGCCAGCATGAATGGATCAAATATCTCCATGCGAAAACGAAGGAACGTAACCCTGAGATGTCCGACTACGAGGAGATGAAAACGACCGGGATCTTTAAGAAAAAATGCCCGgaagaa                                                                                         <  8:324237/149‑1 (MQ=255)
                                                      aGACGTTTATCAAACCTTTACTGAAGGTCGCAGTCAGCATGAATGGATCAAATATCTCCATGCGAAAACGAAGGAACGTAACCCTGAGATGTCCGACTACGAGGAGATGAAAACGACCGGGATCTTTAAGAAAAAATGCCCGGAAGAac                                                                                        <  4:77890/149‑1 (MQ=255)
                                                       gACGTTTATCAAACCTTTACTGAAGGTCGCAGTCAGCAAGAATGGATCAAATATCACCATGCGAAAACTAACGAACGTAACCCTGAGATGTCCGACTACGAGGAGATGAAAACGACCGGGATCTTTAAGAAAAAATGCCCGGAAGAaca                                                                                       <  8:215624/149‑1 (MQ=255)
                                                       gACGTTTATCAAACCTTTACTGAAGGTCGCAGACAGCATGAATGGATCAAATATCTCCATGCGAAAACGAAGGAACGTAACCCTGAGATGTCCGACTACGAGGAGATGAAAACGACCGGGATCTTTAAGAAAAAATGCCCGGAAGAaca                                                                                       <  8:373248/149‑1 (MQ=255)
                                                          gTTTATCAAACCTTTACTGAAGGTCGCAGTCAGCATGAATGGATCAAATATCTCCATGCGAAAACGAAGGAACGTAACCCTGAGATGTCCGACTACGAGGAGATGAAAACGACCGGGATCTTTAAGAAAAAATGCCCGGAAGAACACTa                                                                                    >  7:149082/1‑149 (MQ=255)
                                                              aTCAAACCTTTACTGAAGGTCGCAGTCAGCATGAATGGATCAAATATCTCCATGCGAAAACGAAGGAACGTAACCCTGAGATGTCCGACTACGAGGAGATGAAAACGACCGGGATCTTTAAGAAAAAATGCCCGGAAGAACACTACGTc                                                                                >  1:425135/1‑149 (MQ=255)
                                                              aTCAAACATTTACTGAAGGTCGCAGTCAGCATGAATGGATCAAATATCTCCATGCTAAAACGAAGGAACGTAACCCTGAGATGTCCGACTACGAGGAGATGAAAACGACCGGGATCTTTAAGAAAAAATGCCCGGAAGAACACTACGt                                                                                 <  6:195259/148‑1 (MQ=255)
                                                                   aCCTTTACTGAAGGTCGCAGTCAGCATGAATGGATCAAATATCTCCATGCGAAAACGAAGGAACGTAACCCTGAGATGTCCGACTACGAGGAGATGAAAACGACCGGGATCTTTAAGAAAAAATGCCCGGAAGAACACTACGTCGCttt                                                                           >  1:271195/1‑149 (MQ=255)
                                                                   aCCTTTACTGAAGGTCGCAGTCAGCATGAATGGATCAAATATCTCCATGCGAAAACGAAGCAACGTAACCATGAGATGTCCGACTACGAGGAGTTGAAAACGACCGGGATCTTTAAGAAAAAATGCCCGGAAGAACACTACGTCGCttt                                                                           <  1:423241/149‑1 (MQ=255)
                                                                   aCCTTTACTGAAGGTCGCAGTCAGCATGAATGGATCAAAAATCTCCATGCGAAAACGAAGGAACGTAACCCTGAGATGTCCGACTACGAGGAGATGAAAACGACCGGGATCTTTAAGAAAAAATGCCCGGAAGAACACTACGTCGCttt                                                                           >  3:143485/1‑149 (MQ=255)
                                                                     cTTTACTGAAGGTCGCAGTCAGCATGAATGGATCAAATATCTCCATGCGAAAACGAAGGAACGTAACCCTGAGATGtcc                                                                                                                                               >  5:334906/1‑79 (MQ=255)
                                                                      tttACTGAAGGTCGCAGTCAGCATGAATGGATCAAATATCTCCATGCGAAAACGAAGGAACGTAACCCTGAGATGTCCGACTACGAGGAGATGAAAACGACCGGGATc                                                                                                                 >  6:79810/1‑108 (MQ=255)
                                                                      tttACTGAAGGTCGCAGTCAGCATGAATGGATCAAATATCTCCATGCGAAAACGAAGGAACGTAACCCTGAGATGTCCGACTACGAGGAGATGAAAACGACCGGGATc                                                                                                                 <  5:79810/108‑1 (MQ=255)
                                                                        tACTGAAGGTCGCAGTCAGCATGAATGGATCAAATATCTCCATGCGAAAACGAAGGAACGTAACCCTGAGATGTCCGACTACGAGGAGATGAAAACGACCGGGATCTTTAAGAAAAAATGCCCGGAAGAACACTACGTCGCTTTCcgcg                                                                      >  8:11467/1‑149 (MQ=255)
                                                                        tACTGAAGGTCGCAGTCAGCATGAATGGATCAAATATCTCCATGCGAAAACGAAGGAACGTAACCCTGAGATGTCCGACTACGAGGAGATGAAAACGACCGGGATCTTTAAGAAAAAATGCCCGGAAGAACACTACGTCGCTTTCcgc                                                                       <  2:271195/148‑1 (MQ=255)
                                                                            gAAGGTCGCAGTCAGCATGAATGGATCAAATATCTCCATGCGAAAACGAAGGAACGTAACCCTGAGATGtcc                                                                                                                                               >  1:333429/1‑72 (MQ=255)
                                                                            gAAGGTCGCAGTCAGCATGAATGGATCAAATATCTCCATGCGAAAACGAAGGAACGTAACCCTGAGATGtcc                                                                                                                                               <  2:333429/72‑1 (MQ=255)
                                                                               ggTCGCAGTCAGCATGAATGGATCAAATATCTCCATGCGAAAACGAAGGAACGTAACCCTGAGATGTCCGACTACGAGGAGATGAAAACGACCGGGATCTTTAAGAAAAAATGCCCGGAAGAACACTACGTCGCTTTCCGCGCATTCCg                                                               <  4:66679/149‑1 (MQ=255)
                                                                                gTCGCAGTCAGCATGAATGGATCAAATATCTCCATGCGAAAACGAAGGAACGTAACCCTGAGATGCCCGACTACGAGGAGATGAAAACGACCGGGATCTTTAAGAAAAAATGCCCGGAAGAACACTACGTCGCTTTCCGCGCATTCAGt                                                              >  3:109919/1‑149 (MQ=255)
                                                                                 tCGCAGTCAGCATGAATGGATCAAATATCTCCATGCGAAAACGAAGGAACGTAACCCTGAGATGTCCGACTACGAGGAGATGAAAACGACCGGGATCTTTAAGAAAAAATGCCCGGAAGAACACTACGTCGCTTTCCGCGCATTCCGTg                                                             <  2:203488/149‑1 (MQ=255)
                                                                                   gCAGTCAGCATGAATGGATCAAATATCTCCATGCGAAAACGAAGGAACGTAACCCTGAGATGTCCGACTACGAGGAGATGAAAACGACCGGGATCTTTAAGAAAAAATGCCCGGAAGAACACTACGTCGCTTTCCGCGCATTCCGTGa                                                            <  2:28181/148‑1 (MQ=255)
                                                                                    cAGTCAGCATGAATGGATCAAATATCTCCATGCGAAAACGAAGGAACGTAACCCTGAGATGTCCGACTACGAGGAGATGAAAACGACCGGGATCTTTAAGAAAAAATGCCCGGAAGAACACTACGTAGCATTCAGCGCATTCCGTg                                                             >  2:193718/1‑146 (MQ=255)
                                                                                       tCAGCATGAATGGATCAAATATCTCCATGCGAAAACGAAGGAACGTAACCCTGAGATGTCCGACTACGAGGAGATGAAAACGACCGGGATCTTTAAGAAAAAATGCCCGGAAGAACACTACGTCGCTTTCCGCGCATTCCGTGAAGATc                                                       <  8:149082/149‑1 (MQ=255)
                                                                                         aGCATGAATGGATCAAATATCTCCATGCGAAAACGAAGGAACGTAACCCTGAGATGTCCGACTACGAGGAGATGAAAACGACCGGGATCTTTAAGAAAAAATGCCCGGAAGAACACTACGTCGCTTTCCGCGCATTCCGTGAAGATCca                                                     <  1:147388/149‑1 (MQ=255)
                                                                                         aGCATGAATGGATCAAATATCTCCATGCGAAAACGAAGGAACGTAACCCTGAGATGTCCGACTACGAGGAGATGAAAACGACCGGGATCTTTAAGAAAAAATGCCCGGAAGAACACTACGTCGCTTTCCGCGCATTCCGTGAAGATCca                                                     <  3:315541/149‑1 (MQ=255)
                                                                                          gCATGAATGGATCAAATATCTCCATGCGAAAACGAAGGAACGTAACCCTGAGATGTCCGACTACGAGGAGATGAAAACGACCGGGATCTTTAAGAAAAAATGCCCGGAAGAACACTACGTCGCTTTCCGCGCATTCCGTGAAGATCca                                                     >  3:325971/1‑148 (MQ=255)
                                                                                            aTGAATGGATCAAATATCTCCATGCGAAAACGAAGGAACGTAACCCTGAGATGTCCGACTACGAGGAGATGAAAACGACCGGGATCTTTAAGAAAAAATGCCCGGAAGAACACTACGTCGCTTTCCGCGCATTCCGTGAAGATCCACAg                                                  <  2:425135/149‑1 (MQ=255)
                                                                                             tGAATGGATCAAATATCTCCATGCGAAAACGAAGGAACGTAACCCTGAGATGTCCGACTACGAGGAGATGAAAACGACCGGGATCTTTAAGAAAAAATGcc                                                                                                 >  4:259936/1‑101 (MQ=255)
                                                                                             tGAATGGATCAAATATCTCCATGCGAAAACGAAGGAACGTAACCCTGAGATGTCCGACTACGAGGAGATGAAAACGACCGGGATCTTTAAGAAAAAATGcc                                                                                                 <  3:259936/101‑1 (MQ=255)
                                                                                               aaTGGATCAAATATCTCCATGCGAAAACGAAGGAAAGTAACCCTGAGATGTCCGACTACGAGGAGATGAAAACGACCGGGATCTTTAAGAAAAAATGCCCGGAAGAACACTACGTCGCTTTCCGCGCATTCCGTGAAGATCCACAGGc                                                >  6:310832/1‑148 (MQ=255)
                                                                                                       aaaTATCTCCATGCGAAAACGAAGGAACGTAACCCTGAGATGTCCGACTACGAGGAGATCAAAACGACCGGGATCTTTAAGAAAAAATGCCCGGAAGAACACTACGTCGCTTTCCGCGCATTCCGTGAAGATCCACAGGCAAACCCGtt                                       <  8:285682/149‑1 (MQ=255)
                                                                                                        tatatCTCCATGCGAAAACGAAGGAACGTAACCCTGAGATGTCCGACTACGAGGAGATGAAAACGACCGGGATCTTTAAGAAAAAATGCCCGGAAGAACACTACGTCGCTTTCCGCGCATTCCGTGAAGa                                                         <  4:67687/129‑1 (MQ=255)
                                                                                                        tatatCTCCATGCGAAAACGAAGGAACGTAACCCTGAGATGTCCGACTACGAGGAGATGAAAACGACCGGGATCTTTAAGAAAAAATGCCCGGAAGAACACTACGTCGCTTTCCGCGCATTCCGTGAAGa                                                         >  3:67687/2‑130 (MQ=255)
                                                                                                        aaTATCTCCATGCGAAAACGAAGGAACGTAACCCTGAGATGTCCGACTACGAGGAGATGAAAACGACCGGGATCTTTAAGAAAAAATGCCCGGAAGAACACTACGTCGCTTTCCGCGCATTCCGTGAAGATCCACAGGCAAACCCGtt                                       >  4:187969/1‑148 (MQ=255)
                                                                                                           atCTCCATGCGAAAACGAAGGAACGTAACCCTGAGATGTCCGACTACGAGGAGATGAAAACGACCGGGATCTTTAAGAAAAAATGCCCGGAAGAACACTACGTCGCTTTCCGCGCATTCCGTGAAGATCCACAGGCAAACCCGTTGaaa                                   >  8:423698/1‑149 (MQ=255)
                                                                                                            tctcCATGCGAAAACGAAGGAACGTAACCCTGAGATGTCCGACTACGAGGAGATGAAAACGACCGGGATCTTTAAGAAAAAATGCCCGGAAGAACACTACGTCGCTTTCCGCGCATTCCGTGAAGATCCACAGGCAAACCCGTTGaaa                                   <  4:325971/148‑1 (MQ=255)
                                                                                                             ctcCATGCGAAAACGAAGGAACGTAACCCTGAGATGTCCGACTACGAGGAGATGAAAACGACCGGGATCTTTAAGAAAAAATGCCCGGAAGAACACTACGTCGCTTTCCGCGCATTCCGTGAAGATCcac                                                    <  6:233236/130‑1 (MQ=255)
                                                                                                             ctcCATGCGAAAACGAAGGAACGTAACCCTGAGATGTCCGACTACGAGGAGATGAAAACGACCGGGATCTTTAAGAAAAAATGCCCGGAAGAACACTACGTCGCTTTCCGCGCATTCCGTGAAGATCcac                                                    >  5:233236/1‑130 (MQ=255)
                                                                                                              tcCATGCGAAAACGAAGGAACGTAACCCTGAGATGTCCGACTACGAGGAGATGAAAACGACCGGGATCTTTAAGAAAAAATGCCCGGAAGAACACTACGTCGCTTTCCGCGCATTCCGTGAAGATCcac                                                    >  4:249313/1‑129 (MQ=255)
                                                                                                                    cGAAAACGAAGGAACGTAACCCTGAGATGTCCGACTACGAGGAGATGAAAACGACCGGGATCTTTAAGAAAAAATGCCCGGAAGAACACTACGTCGCTTTCCGCGCATTCCGTGAAGATCCACAGGCAAACCCGTTGAAAACACCTTCg                          >  1:232230/1‑149 (MQ=255)
                                                                                                                        aaCGAAGGAACGTAACCCTGAGATGTCCGACTACGAGGAGATGAAAACGACCGGGATCTTTAAGAAAAAATGCCCGGAAGAACACTACGTCGCTTTCCGCGCATTCCGTGAAGATCCACAGGCAAACCCGTTGAAAACACCTTCgggg                       >  2:385765/1‑148 (MQ=255)
                                                                                                                           gacggAACGTAACCCTGAGATGTCCGACTACGAGGAGATGAAAACGACCGGGATCTTTAAGAAAAAATGCCCGGAAGAACACTACGTCGCTTTcc                                                                         <  6:291717/92‑1 (MQ=255)
                                                                                                                           gAAGGAACGTAACCCTGAGATGTCCGACTACGAGGAGATGAAAACGACCGGGATCTTTAAGAAAAAATGCCCGGAAGAACACTACGTCGCTTTCCGCGCATTCCGTGAAGATCCACAGGCAAACCCGTTGAAAACACCTTCGGGGaaaa                   <  7:423698/149‑1 (MQ=255)
                                                                                                                             aGGAACGTAACCCTGAGATGTCCGACTACGAGGAGATGAAAACGACCGGGATCTTTAAGAAAAAATGCCCGGAAGAACACTACGTCGCTTTCCGCGCATTCCGTGAAGATCCACAGGCAAACCCGTTGAAAACACCTTCGGGGAAAATc                 <  4:256700/149‑1 (MQ=255)
                                                                                                                               gAACGTAACCCTGAGATGTCCGACTACGAGTATATGAAAACGACCGTTATCTTTACGAAAAAATGCCCGGAAGAACACTACGTCGCTTTCCGCGCATTCCGTGAAGATCCACAGGCAAACCCGTTGaaaa                                  <  7:297285/130‑1 (MQ=255)
                                                                                                                               gAACGTAACCCTGAGATGTCCGACTACGAGGAGATGAAAACGCCCGGGCTCTTTAAGAAAAAATGCCCGGAAGAACACTCCGTCGCTTTCCGCGCATTCCGTGAAGCTCCACAGGCAACCCCGTTGaaac                                  >  8:297285/1‑129 (MQ=255)
                                                                                                                               gAACGTAACCCTGAGATGTCCGACTACGAGGAGATGAAAACGACCGGGATCTTTAAGAAAAAATGCCCGGAAGAACACTACGTCGCTTTCCGCGCATTCCGTGAAGATCCACAGGCAAACCCGTTGAAAACACCTTCGGGGAAAATCg                >  7:70223/1‑148 (MQ=255)
                                                                                                                                aaCGTAACCCTGAGATGTCCGACTACGAGGAGATGAAAACGACCGGGATCTTTAAGAAAAAATGcc                                                                                                 >  6:9506/1‑66 (MQ=255)
                                                                                                                                aaCGTAACCCTGAGATGTCCGACTACGAGGAGATGAAAACGACCGGGATCTTTAAGAAAAAATGcc                                                                                                 <  5:9506/66‑1 (MQ=255)
                                                                                                                                aaCGTAACCCTGAGATGTCCGACTACGAGGAGATGAAAACGACCGGGATCTTTAAGAAAAAATGCCCGGAAGAACACTACGTCGCTTTCCGCGCATTCCGTGAAGATCCACATGCAAACCCGTTGAAAACACCTTCGGGGAAAATCGaa              <  3:13286/149‑1 (MQ=255)
                                                                                                                                aaCGTAACCCTGAGATGTCCGACTACGAGGAGATGAAAACGACCGGGATCTTTAAGAAAAAATGCCCGGAAGAACACTACGTCGCTTTCCGCGCATTCCGTGAAGATCCACAGGCAAACCCGTTGAAAACACCTTCGGGGAAAATCGaa              <  1:107186/149‑1 (MQ=255)
                                                                                                                                 aCGTAACCCTGAGATGTCCGACTACGAGGAGATGAAAACGACCGGGATCTTTAAGAAAAAATGCCCGGAAGAACACTACGTCGCTTTCCGCGCATTCCGTGAAGATCCACAGGCAAACCCGTTGAAAACACCTTCGGGGAAAATCGaaa             >  6:223954/1‑149 (MQ=255)
                                                                                                                                 aCGTAACACTGAGATGTCCGTCTACGAGGAGATTAAAACGACCGGGATCTTTAAGAAAAAATGCCCGTAAGAACACTACGTCGCTTTCCGCGCATTCCGTGAAGATCCACAGGCAAAGCCGTTGAAAACACCTTCGGGGAAAATCGaaa             <  1:193718/149‑1 (MQ=255)
                                                                                                                                  cGTAACCCTGAGATGTCCGACTACGAGGAGATGAAAACGACCGGGATCTTTAAGAAAAAATGCCCGGAAGAACACTACGTCGCTTTCCGCGCATTCCGTGAAGATCCACAGGCAAACCCGTTGAAAACACCTTCGGGGAAAATCGAAAt            <  4:116670/149‑1 (MQ=255)
                                                                                                                                    tAACCCTGAGATGTCCGACTACGAGGAGATGAAAACGACCGGGATCTTTAAGAAAAAATGCCCGGAAGAACACTACGTCGCTTTCCGCGCATTCCGTGAAGATCCACAGGCAAACCCGTTGAAAACACCTTCGGGGAAAATCGAAAttt          <  2:223296/149‑1 (MQ=255)
                                                                                                                                       cccTGAGATGTCCGACTACGAGGAGATGAAAACGACCGGGATCTTTAAGAAAAAATGCCCGGAAGAACACTACGTCGCTTTCCGCGCATTCCGTGAAGATCCACAGGCAAACCCGTTGAAAACACCTTCGGGGAAAATCGAAATTTAtt       <  7:11467/149‑1 (MQ=255)
                                                                                                                                        ccTGAGATGTCCGACTACGAGGAGATGAAAACGACCGGGATCTTTAAGAAAAAATGCCCGGAAGAACACTACGTCGCTTTCCGCGCATTCCGTGAAGATCCACAGGCAAACCCGTTGAAAACACCTTCGGGGAAAATCGAAAttt          <  1:398357/145‑1 (MQ=255)
                                                                                                                                        ccTGAGATGTCCGACTACGAGGAGATGAAAACGACCGGGATCTTTAAGAAAAAATGCCCGGAAGAACACTACGTCGCTTTCCGCGCATTCCGTGAAGATCCACAGGCAAACCCGTTGAAAACACCTTCGGGGAAAATCGAAAttt          >  2:398357/1‑145 (MQ=255)
                                                                                                                                        ccTGACAAGTCCGACTACGAGGAGATGAAAACGACCGGGATCTTTAAGAAAAAATTCCCGGAAGAACACTACGTCGCATTCCGCGCATTCCGTGAAGATCCACAGGCAAACCCGTTGAAAACACCTTCGGGGAAAATCGAAATTTATTc      <  1:184826/149‑1 (MQ=255)
                                                                                                                                          tGAGATGTCCGACTACGAGGAGATGAAAACGACCGGGATCTTTAAGAAAAAATGCCCGGAAGAACACTACGTCGCTTTCCGCGCATTCCGTGAAGATCCACAGGCAAACCCGTTGAAAACACCTTCGGGGAAAATCGAAATTTATTCt     >  7:86529/1‑148 (MQ=255)
                                                                                                                                            agaTGTCCGACTACGAGGAGATGAAAACGACCGGGATCTTTAAGAAAAAATGCCCGGAAGAACACTACGTCGCTTTCCGCGCATTCCGTGAAGATCCACAGGCAAACCCGTTGAAAACACCTTCGGGGAAAATCGAAATTTATTCTGaa  <  2:412263/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                               
TGAACGCAAACCCATCTACTGGATGCTGAGTGAAGTCGCCAAACGCTTAGGACCAGACGTTTATCAAACCTTTACTGAAGGTCGCAGTCAGCATGAATGGATCAAATATCTCCATGCGAAAACGAAGGAACGTAACCCTGAGATGCCCGACTACGAGGAGATGAAAACGACCGGGATCTTTAAGAAAAAATGCCCGGAAGAACACTACGTCGCTTTCCGCGCATTCCGTGAAGATCCACAGGCAAACCCGTTGAAAACACCTTCGGGGAAAATCGAAATTTATTCTGAA  >  NC_000913/1662205‑1662493

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