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.6% 226.4 / NA 70P599S (CCC→TCC) ynfFS‑ and N‑oxide reductase, A subunit, periplasmic
Reads supporting (aligned to +/- strand):  ref base C (0/1);  major base T (35/33);  minor base G (0/1);  total (35/35)
Fisher's exact test for biased strand distribution p-value = 4.93e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CAAACCCATCTACTGGATGCTGAGTGAAGTCGCCAAACGCTTAGGACCAGACGTTTATCAAACCTTTACTGAAGGTCGCAGTCAGCATGAATGGATCAAATATCTCCATGCGAAAACGAAGGAACGTAACCCTGAGATGCCCGACTACGAGGAGATGAAAACGACCGGGATCTTTAAGAAAAAATGCCCGGAAGAACACTACGTCGCTTTCCGCGCATTCCGTGAAGATCCACAGGCAAACCCGTTGAAAACACCTTCGGGGAAAATCGAAATTTATTCTGAA  >  NC_000913/1662211‑1662493
                                                                                                                                           |                                                                                                                                               
cAAACCCATCTACTGGATGCTGAGTGAAGTCGCCAAACGCTTAGGACCAGACGTTTATCAAACCTTTACTGAAGGTCGCAGTCAGCATGAATGGATCAAATATCTCCATGCGAAAACGAAGGAACGTAACCCTGAGAAGTCCGACTAc                                                                                                                                         >  4:348971/1‑148 (MQ=255)
      cATCTACTGGATGCTGAGTGAAGTCGCCAAACGCTTAGGACCAGACGTTTATCAAACCTTTACTGAAGGTCGCAGTCAGCATGAATGGATCAAATATCTCCATGCGAAAACGAAGGAACGTAACCCTGAGATGTCCGACTACGAGgaga                                                                                                                                  <  7:202148/149‑1 (MQ=255)
            cttgatgatTAGTGAAGTCTCCAAACGCTTAGGACCACACGTTTATCAAACCTTTACTGAAGGTCGCAGTCAGCATGAATGGATCAAATATCTCCATGCGAAAACGAAGCAACGTAACCCTGAGATGTCCGACTACGATGAGATGaaaa                                                                                                                            <  1:132428/146‑1 (MQ=255)
            cTGGATGCTTAGTGAACTCGCCAAACGCTTAGGACCAGACGTTTATCAAACCTTTACTGAAGGTCGCAGTCAGCATGAATGGATCAAATATCTCCATGCGAAAACGAAGGAACGTAACCCTGAGATGTCCGACTACGAGgaga                                                                                                                                  <  6:140477/143‑1 (MQ=255)
               gATGCTGAGTGAAGTCGCCAAACGCTTAGGACCAGACGTTTATCAAACCTTTACTGAAGGTCGCAGTCAGCATGAATGGATCAAATATCTCCATGCGAAAACGAAGGAACGTAACCCTGAGATGTCCGACTACGAGGAGATGAAAAc                                                                                                                           <  8:52266/147‑1 (MQ=255)
                aTGCTGAGTGAAGTCGCCAAACGCTTAGGACCAGACGTTTATCAAACCTTTACTGAAGGTCGCAGTCAGCATGAATGGATCAAATATCTCCATGCGAAAACGAAGGAACGTAACCCTGAGATGTCCGACTACGAGGAGATGAAAACGa                                                                                                                         <  4:256410/148‑1 (MQ=255)
                 tGCTGAGTGAAGTCGCCAAACGCTTAGGACCAGACGTTTATCAAACCTTTACTGAAGGTCGCAGTCAGCATGAATGGATCAAATATCTCCATGCGAAAACGAAGGAACGTAACCCTGAGATGTCCGACTACGAGGAGATGAAAACGAcc                                                                                                                       <  1:188196/149‑1 (MQ=255)
                  gCTGAGTGAAGTCGCCAAACGCTTAGGACCAGACGTTTATCAAACCTTTACTGAAGGTCGCAGTCAGCATGAATGGATCAAATATCTCCATGCGAAAACGAAGGAACGTAACCCTGAGATGTCCGACTACGAGGAGATGAAAACGACCg                                                                                                                      >  5:56619/1‑149 (MQ=255)
                    tGAGTGAAGTCGCCAAACGCTTAGGACCAGACGTTTATCAAACCTTTACTGAAGGTCGCAGTCAGCATGAATGGATCAAATATCTCCATGCGAAAACGAAGGAACGTAACCCTGAGATGTCCGACTACGAGGAGATGAAAACGACCgg                                                                                                                     >  7:248024/1‑148 (MQ=255)
                      aGTGAAGTCGCCAAACGCTTAGGACCAGACGTTTATCAAACCTTTACTGAAGGTCGCAGTCAGCATGAATGGATCAAATATCTCCATGCGAAAACGAAGGAACGTAACCCTGAGATGTCCGACTACGAGGAGATGAAAACGACCggg                                                                                                                    <  3:243130/147‑1 (MQ=255)
                        tGAAGTCGCCAAACGCTTAGGACCAGACGTTTATCAAACCTTTACTGAAGGTCGCAGTCAGCATGAATGGATCAAATATCTCCATGCGAAAACGAAGGAACGTAACCCTGAGATGTCCGACTACGAGGAGATGAAAACGACCGGGATCt                                                                                                                <  6:214331/149‑1 (MQ=255)
                        tGAAGTCGCCAAACGCTTAGGACCAGACGTTTATCAAACCTTTACTGAAGGTCGCAGTCAGCATGAATGGATCAAATATCTCCATGCGAAAACGAAGGAACGTAACCCTGAGACGTCCGACTACGAGGAGATGAAAACGACCGGGATCt                                                                                                                <  3:348971/149‑1 (MQ=255)
                          aaGTCGCCAAACGCTTAGGACCAGACGTTTATCAAACCTTTACTGAAGGTCGCAGTCAGCATGAATGGATCAAATATCTCCATGCGAAAACGAAGGAACGTAACCCTGAGATGTCCGACTACGAGGAGATGAAAACGACCGGGATCtt                                                                                                               <  5:376833/148‑1 (MQ=255)
                           aGTCGCCAAACTCTTAGGACCAGACGTTTATCAAACCTTTACTGAAGGTCGCATTCAGCATGAATGGATCAAATATCTCCATGCGAAAACGAAGGAAAGTAACCCTGAGATGTCCGACTACGAGGAGATGAAAACGACCGGGATCttt                                                                                                              <  1:80352/148‑1 (MQ=255)
                            gTCGCCAAACTCTTAGGACCAGACGTTTATCAAACCTTTACTGAAGGTCGCAGTAAGCATGAATGGATCAAATATCTCCATGCGAAAACGAAGGAACTTAACCCTGATATGTCCGACTACGAGGATATGAAAACGACCGGGATCTTTaa                                                                                                            <  5:324718/149‑1 (MQ=255)
                            gTCGCCAAACGCTTGGGACCAGACGTTTATCAAACCTTTACTGAAGGTCGCAGTCAGCATGAATGGATCAAATATCTCCATGCGAAAACGAAGGAACGTAACCCTGAGATGTCCGACTACGAGGAGATGAAAACGa                                                                                                                         >  3:167540/1‑136 (MQ=255)
                            gTCGCCAAACGCTTGGGACCAGACGTTTATCAAACCTTTACTGAAGGTCGCAGTCAGCATGAATGGATCAAATATCTCCATGCGAAAACGAAGGAACGTAACCCTGAGATGTCCGACTACGAGGAGATGAAAACGa                                                                                                                         <  4:167540/136‑1 (MQ=255)
                               gCCAAACGCTTAGGACCAGACGTTTATCAAACCTTTACTGAAGGTCGCAGTCAGCATGAATGGATCAAATATCTCCATGCGAAAACGAAGGAACGTAACCCTGAGATGTCCGACTACGAGGAGATGAAAACGACCGGGATCTTTAAGa                                                                                                          >  6:83669/1‑148 (MQ=255)
                               gCCAAACGATTAGGACCAGACGTTTATCAAACATTTACTTAAGGTCGCAGTCAGCATTAATGTATCAAATATCTCCATGCGAAAACGAAGTAACGTAACCCTGATATGTCCGAATACTATGAGATGAAAACGACCGGGATCTTTAAGa                                                                                                          <  7:336314/148‑1 (MQ=255)
                                     ccgatctgACCAGACGTTTATCAAACCTTTACTGAAGGTCGCAGTCAGCATGAATGGATCAAATATCTCCATGCGAAAACGAAGGAACGTAACCCTGAGATGTCCGACTACGAGGAGATGAAAACGACCGGGATCTTTAAGAAAAAAt                                                                                                    <  5:229125/141‑1 (MQ=255)
                                         cctatctaGACGTTTATCAAACCTTTACTGAAGGTCGCAGTCAGCATGAATGGATCAAATATCTCCATGCGAAAACGAAGGAACGTAACCCTGAGATGTCCGACTACGAGGAGATGAAAACGACCGGGATCTTTAAGAAAAAATGccc                                                                                                <  5:13160/141‑1 (MQ=255)
                                            gACCAGACGTTTATCAAACCTTTACTGAAGGTCGCAGTCAGCATGAATGGATCAAATATCTCCATGCGAAAACGAAGGAACGTAACCCTGAGATGTCCGACTACGAGGAGATGaaaa                                                                                                                            <  1:152211/117‑1 (MQ=255)
                                            gACCAGACGTTTATCAAACCTTTACTGAAGGTCGCAGTCAGCATGAATGGATCAAATATCTCCATGCGAAAACGAAGGAACGTAACCCTGAGATGTCCGACTACGAGGAGATGaaaa                                                                                                                            >  2:152211/1‑117 (MQ=255)
                                             aCCAGACGTTTATCAAACCTTTACTGAAGGTCGCAGTCAGCATGAATGGATCAAATATCTCCATGCGAAAACGAAGGAACGTAACCCTGAGATGTCCGACTACGAGGAGATGAAAACGACCGGGATCTTTAAGAAAAAAt                                                                                                    >  6:229125/1‑140 (MQ=255)
                                             aCCAGACGTTTATCAAACCTTTACTGAAGGTCGCAGTCAGCATGAATGGATCAAATATCTCCATGCGAAAACGAAGGAACGTAACCCTGAGATGTCCGACTACGAGGAGATGAAAACGACCGGGATCTTTAAGAAAAAATGCCCGgaag                                                                                           <  6:56619/149‑1 (MQ=255)
                                               cAGACGTTTATCAAACCTTTACTGAAGGTCGCAGTCAGCATGAATGGATCAAATATCTCCATGCGAAAACGAAGGAACGTAACCCTGAGATGTCCGACTACGAGGAGATGAAAACGACCGGGATCt                                                                                                                <  1:225996/126‑1 (MQ=255)
                                               cAGACGTTTATCAAACCTTTACTGAAGGTCGCAGTCAGCATGAATGGATCAAATATCTCCATGCGAAAACGAAGGAACGTAACCCTGAGATGTCCGACTACGAGGAGATGAAAACGACCGGGATCt                                                                                                                >  2:225996/1‑126 (MQ=255)
                                                 gACGTTTATCAAACCTTTACTGAAGGTCGCAGTCAGCATGAATGGATCAAATATCTCCATGCGAAAACGAAGGAACGTACCCCTGAGATGTCCGACTACGAGGAGATGAAAACGACCGGGATCTTTAAGAAAAAATGccc                                                                                                >  6:13160/1‑140 (MQ=255)
                                                    gTTTATCAAACCTTTACTGAAGGTCGCAGTCAGCATGAATGGATCAAATATCTCCATGCGAAAACGAAGGAACGTAACCCTGAGATGTCCGACTACGAGGAGATGAAAACGACCGGGATCTTTAAGAAAAAATGCCCGGAAGAACACTa                                                                                    >  2:292115/1‑149 (MQ=255)
                                                     tttATCAAACCTTTACTGAAGGTCGCAGTCAGCATGAATGGATCAAATATCTCCATGCGAAAACGAAGGAACGTAACCCTGAGATGTCCGACTACGAGGAGATGAAAACGACCGGGATCTTTAAGAAAAAATGcc                                                                                                 <  7:293403/135‑1 (MQ=255)
                                                     tttATCAAACCTTTACTGAAGGTCGCAGTCAGCATGAATGGATCAAATATCTCCATGCGAAAACGAAGGAACGTAACCCTGAGATGTCCGACTACGAGGAGATGAAAACGACCGGGATCTTTAAGAAAAAATGcc                                                                                                 >  8:293403/1‑135 (MQ=255)
                                                     tttATCAAACCTTTACTGAAGGTCGCAGTCAGCATGAATGGATCAAATATCTCCATGCGAAAACGAAGGAACGTAACCCTGAGATGTCCGACTACGAGGAGATGAAAACGACCGGGATCTTTAAGAAAAAATGCCCGGAAGAACACTAc                                                                                   >  5:307686/1‑149 (MQ=255)
                                                          cAAACCTTTACTGAAGGTCGCAGTCAGCATGAATGGATCAAATATCTCCATGCGAAAACGAAGGAACGTAACCCTGAGATGTCCGACTACGAGGAGATGAAAACGACCGGGATCTTTAAGAAAAAATGCCCGGAAGAACACTACGTCGc                                                                              >  6:211828/1‑149 (MQ=255)
                                                            aaCCTTTACTGAAGGTCGCAGTCAGCATGAATGGATCAAATATCTCCATGCGAAAACGAAGGAACGTAACCCTGAGATGTCCGACTACGAGGAGATGAAAACGACCGGGATCTTTAAGAAAAAATGCCCGGAAGAACACTACGTCGCtt                                                                            >  4:194476/1‑149 (MQ=255)
                                                                   aCTGAAGGTCGCAGTCAGCATGAATGGATCAAATATCTCCATGCGAAAACGAAGGAACGTAACCCTGAGATGTCCGACTACGAGGAGATGAAAACGACCGGGATCTTTAAGAAAAAATGCCCGGAAGAACACTACGTCGCTTTCagcg                                                                      >  4:39587/1‑148 (MQ=255)
                                                                     tGAAGGTCGCAGTCAGCATGAATGGATCAAATATCTCCATGCGAAAACGAAGGAACGTAACCCTGAGATGTCCGACTACGAGGAGATGAAAACGACCGGGATCTTTAAGAAAAAATGCCCGGAAGAACACTACGTCGCTTTCCGCGCAt                                                                   >  7:151128/1‑149 (MQ=255)
                                                                      gAAGGTCGCAGTCAGCATGAATGGATCAAATATCTCCATGCGAAAACGAAGGAACGTAACCCTGAGATGTCCGACTACGAGGAGATGAAAACGACCGGGATCTTTAAGAAAAAATGCCCGGAAGAACACTACGTCGCTTTCCGCGCAtt                                                                  >  7:39926/1‑149 (MQ=255)
                                                                      gAAGGTCGCAGTCAGCATGAATGGATCAAATATCTCCATGCGAAAACGAAGGAACGTAACCCTGAGATGTCCGACTACGAGGAGATGAAAACGACCGGGATCTTTAAGAAAAAATGCCCGGAAGAACACTACGTCGCTTTCCGCGCAtt                                                                  >  7:42097/1‑149 (MQ=255)
                                                                       aaGGTCGCAGTCAGCATGAATGGATCAAATATCTCCATGCGAAAACGAAGGAACGTAACCCTGAGATGTCCGACTACGAGGAGATGAAAACGACCGGGATCTTTATGAAAAAATGCCCGGAAGAACACTACGTCGCTTTCCGAGCATTc                                                                 >  8:132746/1‑149 (MQ=255)
                                                                        aGGTCGCAGTCAGCATGAATGGATCAAATATCTCCATGCGAAAACGAAGGAACGTAACCCTGAGATGTCCGACTACGAGGAGATGAAAACGACCGGGATCTTTAAGAAAAAATGCCCGGAAGAACACTACGTCGCTTTCCGCGCATTcc                                                                <  3:194476/149‑1 (MQ=255)
                                                                          gTCGCAGTAAGCATGAATGCATAAAATAGCTCCATGCCAAAACGAAGGAACGTAACCCTGAGATGGCCGACTACGAGGAGATGAAAACGACCGGGATCTTTAAGAAAAAATGCCCGGAAGAACACTACGTCGCTTTCCGCGCATTCCGt                                                              <  8:248024/149‑1 (MQ=255)
                                                                                 tCAGCTTGAATGGATCTAATTTCTCCATGCTAAAACGAAGGAACTTATCCCTTAGATGTCCGCCTACGAGGTGATTAAAACGACCGGGATCTTTAAGAAAAAATGCCCGGAAGAACACTACGTCGCTTTCCGCGCATTCCGTGAAGa                                                         <  3:284616/147‑1 (MQ=255)
                                                                                      aTGAATGGATCAAATATCTCCATGCGAAAACGAAGGAACGTAACCCTGAGATGTCCGACTACGAGGAGATGAAAACGACCGGGATCTTTAAGAAAAAATGCCCGGAAGAACACTACGTCGCTTTCCGCGCATTCCGTGAAGATCCACAg                                                  <  5:83669/149‑1 (MQ=255)
                                                                                       tGAATGGATCAAATATCTCCATGCGAAAACGAAGGAACGTAACCCTGAGATGTCCGACTACGAGGAGATGAAAACGACCGGGATCTTTAAGAAAAAATGCCCGGAAGAACACTACGTCGCTTTCCGCGCATTCCGTGAAGATCCACAgg                                                 >  1:383813/1‑149 (MQ=255)
                                                                                       tGAATGGATCAAATATCTCCATGCGAAAACGAAGGAACGTAACCCTGAGATGTCCGACTACGAGGAGATGAAAACGACCGGGATCTTTAAGAAAAAATGCCCGGAAGAACACTACGTCGCTTTCCGCGCATTCCGTGAAGATCCACAgg                                                 >  4:336596/1‑149 (MQ=255)
                                                                                         aaTGGTTCAAATATCTCCATGCGAAAACGAAGGAACGTAACCCTGAGATGTCCGACTACGAGGAGATGAAAACGACCGGGATCTTTAAGAAAAAATGCCCGGAAGAACACTACGTCGCTTTCCGCGCATTCCGTGAAGATCCACAGGc                                                >  2:321106/1‑148 (MQ=255)
                                                                                              aTCAAATATCTCCATGCGAAAACGAAGGAACGTAACCCTGAGATGTCCGACTACGAGGAGATGAAAACGACCGGGATCTTTAAGAAAAAATGCCCGGAAGAACACTACGTCGCTTTCCGCGCATTCCGTGAAGATCCACAGGCAAAccc                                          >  7:105335/1‑149 (MQ=255)
                                                                                               tCAAATATCTCCATGCGAAAACGAAGGAACGTAACCCTGAGATGTCCGACTACGAGGAGATGAAAACGACCGGGATCTTTAAGAAAAAATGCCCGGAAGAACACTACGTCGCTTTCCGCGCATTCCGTGAAGATCCACAg                                                  <  6:290336/140‑1 (MQ=255)
                                                                                               tCAAATATCTCCATGCGAAAACGAAGGAACGTAACCCTGAGATGTCCGACTACGAGGAGATGAAAACGACCGGGATCTTTAAGAAAAAATGCCCGGAAGAACACTACGTCGCTTTCCGCGCATTCCGTGAAGATCCACAg                                                  >  5:290336/1‑140 (MQ=255)
                                                                                               tCAAATATCTCCATGCGAAAACGAAGGAACGTAACCCTGAGATGTCCGACTACGAGGAGATGAAAACGACCGGGATCTTTAAGAAAAAATGCCCGGAAGAACACTACGTCGCTTTCCGCGCATTCCGTGAAGATCCACAGGCAAACCCg                                         >  2:152293/1‑149 (MQ=255)
                                                                                                   atatCTCCATGCGAAAACGAAGGAACGTAACCCTGAGATGTCCGACTACGAGGAGATGAAAACGACCGGGATCTTTAAGAAAAAATGCCCGGAAGAACACTACGTCGCTTTCCGCGCATTCCGTGAAGATCCACAGGCAAACCCGTTg                                      >  7:262158/1‑148 (MQ=255)
                                                                                                   atatCTCCATGCGAAAACGAAGGAACGTAACCCTGAGATGTCCGACTACGAGGAGATGAAAACGACCGGGATCTTTAAGAAAAAATGCCCGGAAGAACACTACGTCGCTTTCCGCGCATTCCGTGAAGATCCACAGGCAAACCCGTTGa                                     <  1:94514/149‑1 (MQ=255)
                                                                                                     atCTCCATGCGAAAACGAAGGAACGTAACCCTGAGATGTCCGACTACGAGGAGATGAAAACGACCGGGATCTTTAAGAAAAAATGCCCGGAAGAACACTACGTCGCTTTCCGCGCATTCCGTGAAGATc                                                       >  1:50332/1‑129 (MQ=255)
                                                                                                     atCTCCATGCGAAAACGAAGGAACGTAACCCTGAGATGTCAGACTACGAGGAGATGAAAACGACCGGGATCTTTAAGAAAAAATGCCCGGAAGAACACTACGTCGCTTTCCGCGCATTCCGTGAAGATc                                                       <  2:50332/129‑1 (MQ=255)
                                                                                                         ccATGCGAAAACGAAGGAACGTAACCCTGAGATGTCCGACTACGAGGAGATGAAAACGACCGGGATCTTTAAGAAAAAATGCCCGGAAGAACACTACGTCGCTTTCCGCGCATTCCGTGAAGATCCACAGGCAAACCCGTTGAAAacac                               >  4:68419/1‑149 (MQ=255)
                                                                                                          cATGCGAAAACGAAGGAACGTAACCCTGAGATGTCCGACTACGAGGAGATGAAAACGACCGGGATCTTTAAGAAAAAATGCCCGGAAGAACACTACGTCGCTTTCCGCGCATTCCGTGAAGATCCACAGGCAAACCCGTTGaaaa                                  <  1:321106/145‑1 (MQ=255)
                                                                                                           aTGCGAAAACGAAGGAACGTAACCCTGAGATGTCCGACTACGAGGAGATGAAAACGACCGGGATCTTTAAGAAAAAATGCCCGGAAGAACACTACGTCGCTTTCCGCGCATTCCGTGAAGATCCACAGGCAAACCCGTTGAAACCACCt                             >  5:362358/1‑149 (MQ=255)
                                                                                                             gCGAAAACGAAGGAACGTAACCCTGAGATGTCCGACTACGAGGAGATGAAAACGACCGGGATCTTTAAGAAAAAATGCCCGGAAGAACACTACGTCGCTTTCCGCGCATTCCGTGAAGATCCACAGGCAAACCCGTTGAAAACACCTTc                           >  5:173923/1‑149 (MQ=255)
                                                                                                              cGAAAACGAAGGAACGTAACACTGAGATGTCCGACTACGAGGAGATGAAAACGACCGGGATCTTTAAGAAAAAATGCCCGGAAGAACACTACGTCGCTTTCCGCGCATTCCGTGAAGATCCACAGGCAAACCCGTTGAAAACACCTTc                           <  3:149028/148‑1 (MQ=255)
                                                                                                                  aaCGAAGGAACGTAACCCTGAGATGTCCGACTACGAGGAGATGAAAACGACCGGGATCTTTAAGAAAAAATGCCCGGAAGAACACTACGTCGCTTTCCGCGCATTCCGTGAAGATCCACAGGCAAACCCGTTGAAAACACCTTCgggg                       >  4:19314/1‑148 (MQ=255)
                                                                                                                      aaGGAACGTAACCCTGAGATGTTCGACTACGAGGAGATGAAAACGACCGGGATCTTTAAGAAAAAATGCCCGGAAGAACACTACGTCGCTTTCCGCGCATTCCGTGAAGATCCACAGGCAAACCCGTTGAAAACACCTTCGGGGAAAAt                  >  6:250671/1‑149 (MQ=255)
                                                                                                                       cggAACGTAACCCTGAGATGCCCGAATACGAGGACATGAAAACGACCGGGATCTTTAAGAAAAAATGCCCGGAAGAACACTACGTCGCTTTCCGCGCATTCCGTGAAGATCCACAGGCAAACCCGTTGAAAACACCTTCGGGGAAAATc                 <  8:58467/148‑1 (MQ=255)
                                                                                                                         gAACGTAACCCTGAGATGTCCGACTACGAGGAGATGAAAACGACCGGGATCTTTAAGAAAAAATGCCCGGAAGAACACTACGTCGCTTTCCGCGCATTCCGTGAAGATCCACAGGCAAACCCGTTGCACACACCTTCGGGGAAAATCg                >  1:217680/1‑148 (MQ=255)
                                                                                                                          aaCGTAACCCTGAGATGTCCGACTACGAGGAGATGAAAACGACCGGGATCTTTAAGAAAAAATGCCCGGAAGAACACTACGTCGCTTTCCGCGCATTCCGTGAAGATCCACAGGCAAACCCGTTGAAAACACCTTCGGGGAAAATCGaa              <  8:262158/149‑1 (MQ=255)
                                                                                                                            cGTAACCCTGAGATGTCCGACTACGAGGAGATGAAAACGACCGGGATCTTTAAGAAAAAATGCCCGGAAGAACACTACGTCGCTTTCCGCGCATTCCGTGAAGATCCACAGGCAAACCCGTTGAAAACACCTTCGGGGAAAATCGAAAt            <  8:42097/149‑1 (MQ=255)
                                                                                                                             gTAACCCTGAGATGTCCGACTACGAGGAGATGAAAACGACCGGGATCTTTAAGAAAAAATGCCCGGAAGAACACTACGTCGCTTTCCGCGCATTCCGTGAAGATCCACAGGCAAACCCGTTGAAAACACCTTCGGGGAAAATCGAAAtt           >  3:79851/1‑149 (MQ=255)
                                                                                                                                 cccTGAGATGTCCGACTACGAGGAGATGAAAACGACCGGGATCTTTAAGAAAAAATGCCCGGAAGAACACTACGTCGCTTTCCGCGCATTCCGTGAAGATCCACAGGCAAACCCGTTGAAAACACCTTCGGGGAAAATCGAAATTTAtt       <  3:39587/149‑1 (MQ=255)
                                                                                                                                    tGAGATGTCCGACTACGAGGAGATGAAAACGACCGGGATCTTTAAGAAAAAATGCCCGGAAGAACACTACGTCGCTTTCCGCGCATTCCGTGAAGATCCACAGGCAAACCCGTTGAAAACACCTTCGGGGAAAATc                 <  2:95135/136‑1 (MQ=255)
                                                                                                                                    tGAGATGTCCGACTACGAGGAGATGAAAACGACCGGGATCTTTAAGAAAAAATGCCCGGAAGAACACTACGTCGCTTTCCGCGCATTCCGTGAAGATCCACAGGCAAACCCGTTGAAAACACCTTCGGGGAAAATc                 >  1:95135/1‑136 (MQ=255)
                                                                                                                                      agaTGTCCGACTACGAGGAGATGAAAACGACCGGGATCTTTAAGAAAAAATGCCCGGAAGAACACTACGTCGCTTTCCGCGCATTCCGTGAAGATCCACAGGCAAACCCGTTGAAAACACCTTCGGGGAAAATCGAAATTTATTCTGaa  <  6:173923/149‑1 (MQ=255)
                                                                                                                                           |                                                                                                                                               
CAAACCCATCTACTGGATGCTGAGTGAAGTCGCCAAACGCTTAGGACCAGACGTTTATCAAACCTTTACTGAAGGTCGCAGTCAGCATGAATGGATCAAATATCTCCATGCGAAAACGAAGGAACGTAACCCTGAGATGCCCGACTACGAGGAGATGAAAACGACCGGGATCTTTAAGAAAAAATGCCCGGAAGAACACTACGTCGCTTTCCGCGCATTCCGTGAAGATCCACAGGCAAACCCGTTGAAAACACCTTCGGGGAAAATCGAAATTTATTCTGAA  >  NC_000913/1662211‑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: