Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 516,487 C→T 6.2% A190V (GCC→GTC)  fetA → iron export ABC transporter ATPase; peroxide resistance protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913516,4870CT 6.2% 206.6 / 3.0 64A190V (GCC→GTC) fetAiron export ABC transporter ATPase; peroxide resistance protein
Reads supporting (aligned to +/- strand):  ref base C (28/32);  new base T (2/2);  total (30/34)
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.92e-01

GTGAAAAACAACGCATCTCATTGATTCGTAACCTGCAATTTATGCCGAAGGTTTTATTGCTGGATGAAATAACCAGTGCGCTGGATGAAAGTAATAAACATAACGTCAATGAGATGATCCATCGTTATGTGCGCGAGCAAAATATTGCCGTGCTGTGGGTGACACACGATAAAGACGAAATTAATCATGCGGATAAAGTGATTACACTGCAACCGCATGCCGGAGAAATGCAGGAAGCACGCTATGAACTCGCATAATATTACTAACGAATCATTAGCACTGGCATTAATGCTGGT  >  NC_000913/516340‑516635
                                                                                                                                                   |                                                                                                                                                    
gtgAAAAACAACGCATCTCATTGATTCGTAACCTGCAATTTATGCCGAAGGTTTTATTGCTGGATGAAATAACCAGTGCGCTGGATGAAAGTAATAAACATAACGTCAATGAGATGATCCATCGTTATGTGCGCGAGCAAAATATTGcc                                                                                                                                                     <  3:473036/149‑1 (MQ=255)
gtgAAAAACAACGCATCTCATTGATTCGTAACCTGCAATTTATGCCGAAGGTTTTATTGCTGGATGAAATAACCAGTGCGCTGGATGAAAGTAATAAACATAACGTCAATGAGATGATCCAACGTTATGTGCGCGAGCAAAATATTGcc                                                                                                                                                     <  1:144698/149‑1 (MQ=255)
  gAAAAACAACGCATCTCATTGATTCGTAACCTGCAATTTATGCCGAAGGTTTTATTGCTGGATGAAATAACCAGTGCGCTGGATGAAAGTAATAAACATAACGTCAATGAGATGATCCATCGTTATGTGCGCGAGCAAAATATTGCCGt                                                                                                                                                   >  4:9680/1‑149 (MQ=255)
   aaaaaCAACGCATCTCATTGATTCGTAACCTGCAATTTATGCCGAAGGTTTTATTGCTGGATGAAATAACCAGTGCGCTGGATGAAAGTAATAAACATAACGTCAATGAGATGATCCATCGTTATGTGCGCGAGCAAAATATTGCCGTg                                                                                                                                                  >  1:255152/1‑149 (MQ=255)
     aaaCAACGCATCTCATTGATTCGTAACCTGCAATTTATGCCGAAGGTTTTATTGCTGGATGAAATAACCAGTGCGCTGGATGAAAGTAATAAACATAACGTCAATGAGATGATCCATCGTTATGTGCGCGAGCAAAATATTGCCGTGCt                                                                                                                                                >  2:410373/1‑149 (MQ=255)
     aaaCAACGCATCTCATTGATTCGTAACCTGCAATTTATGCCGAAGGTTTTATTGCTGGATGAAATAACCAGTGCGCTGGATGAAAGTAATAAACATAACGTCAATGAGATGATCCATCGTTATGTGCGCGAGCAAAATATTGCCGTGCt                                                                                                                                                <  4:7619/149‑1 (MQ=255)
                ctcATTGATTCGTAACCTGCAATTTATGCCGAAGGTTTTATTGCTGGATGAAATAACCAGTGCGCTGGATGAAAGTAATAAACATAACGTCAATGAGATGATCCATCGTTATGTGCGCGAGCAAAATATTGCCGTGCTGTGGGTGacac                                                                                                                                     <  3:505909/149‑1 (MQ=255)
                 tcATTGATTCGTAACCTGCAATTTATGCCGAAGGTTTTATTGCTGGATGAAATAACCAGTGCGCTGGATGAAAGTAATAAACATAACGTCAATGAGATGATCCATCGTTATGTGCGCGAGCAAAATATTGCCGTGCTGTGGGTGacaca                                                                                                                                    >  3:385951/1‑149 (MQ=255)
                   aTTGATTCGTAACCTGCAATTTATGCCGAAGGTTTTATTGCTGGATGAAATAACCAGTGCGCTGGATGAAAGTAATAAACATAACGTCAATGAGATGATCCATCGTTATGTGCGCGAGCAAAATATTGCCGTGCTGTGGGTGACACACg                                                                                                                                  >  3:254852/1‑149 (MQ=255)
                   aTTGATTCGTAACCTGCAATTTATGCCGAAGGTTTTATTGCTGGATGAAATAACCAGTGCGCTGGATGAAAGTAATAAACATAACGTCAATGAGATGATCCATCGTTATGTGCGCGAGCAAAATATTGCCGTGCTGTGGGTGACACACg                                                                                                                                  >  3:412843/1‑149 (MQ=255)
                    ttGATTCGTAACCTGCAATTTATGCCGAAGGTTTTATTGCTGGATGAAATAACCAGTGCGCTGGATGAAAGTAATAAACATAACGTCAATGAGATGATCCATCGTTATGTGCGCGAGCAAAATATTGCCGTGCTGTGGGTGACACACGa                                                                                                                                 <  2:202030/149‑1 (MQ=255)
                     tGATTCGTAACCTGCAATTTATGCCGAAGGTTTTATTGCTGGATGAAATAACCAGTGCGCTGGATGAAAGTAATAAACATAACGTCAATGAGATGATCCATCGTTATGTGCGCGAGCAAAATATTGCCGTGCTGTGGGTGACACACGAt                                                                                                                                <  2:20504/149‑1 (MQ=255)
                        ttCGTAACCTGCAATTTATGCCGAAGGTTTTATTGCTGGATGAAATAACCAGTGCGCTGGATGAAAGTAATAAACATAACGTCAATGAGATGATCCATCGTTATGTGCGCGAGCAAAATATTGCCGTGCTGTGGGTGACACACGATaaa                                                                                                                             >  2:472624/1‑149 (MQ=255)
                          cGTAACCTGCAATTTATGCCGAAGGTTTTATTGCTGGATGAAATAACCAGTGCGCTGGATGAAAGTAATAAACATAACGTCAATGAGATGATCCATCGTTATGTGCGCGAGCAAAATATTGCCGTGCTGTGGGTGACACACGATAAAGa                                                                                                                           <  4:392680/149‑1 (MQ=255)
                               ccTGCAATTTATGCCGAAGGTTTTATTGCTGGATGAAATAACCAGTGCGCTGGATGAAAGTAATAAACATAACGTCAATGAGATGATCCATCGTTATGTGCGCGAGCAAAATATTGCCGTGCTGTGGGTGACACACGATAAAGACGaaa                                                                                                                      <  2:382194/149‑1 (MQ=255)
                                   cAATTTATGCCGAAGGTTTTATTGCTGGATGAAATAACCAGTGCGCTGGATGAAAGTAATAAACATAACGTCAATGAGATGATCCATCGTTATGTGCGCGAGCAAAATATTGCCGTGCTGTGGGTGACACACGATAAAGACGAAATTaa                                                                                                                  <  1:420035/149‑1 (MQ=255)
                                           gCCGAAGGTTTTATTGCTGGATGAAATAACCAGTGCGCTGGATGAAAGTAATAAACATAACGTCAATGAGATGATCCATCGTTATGTGCGCGAGCAAAATATTGCCGTGCTGTGGGTGACACACGATAAAGACGAAATTAATCATGCgg                                                                                                          >  3:308025/1‑149 (MQ=255)
                                              gAAGGTTTTATTGCTGGATGAAATAACCAGTGCGCTGGATGAAAGTAATAAACATAACGTCAATGAGATGGTCCATCGTTATGTGCGCGAGCAAAATATTGCCGTGCTGTGGGTGACACACGATAAAGACGAAATTAATCATGCGGATa                                                                                                       >  4:548707/1‑149 (MQ=255)
                                                 ggTTTTATTGCTGGATGAAATAACCAGTGCGCTGGATGAAAGTAATAAACATAACGTCAATGAGATGATCCATCGTTATGTGCGCGAGCAAAATATTGTCGTGCTGTGGGTGACACACGATAAAGACGAAATTAATCATGCGGATAAAg                                                                                                    >  2:575076/1‑149 (MQ=255)
                                                 ggTTTTATTGCTGGATGAAATAACCAGTGCGCTGGATGAAAGTAATAAACATAACGTCAATGAGATGATCCATCGTTATGTGCGCGAGCAAAATATTGCCGTGCTGTGGGTGACACACGATAAAGACGAAATTAATCATGCGGATAAAg                                                                                                    <  2:22775/149‑1 (MQ=255)
                                                  gTTTTATTGCTGGATGAAATAACCAGTGCGCTGGATGAAAGTAATAAACATAACGTCAATGAGATGATCCATCGTTATGTGCGCGAGCAAAATATTGCCGTGCTGTGGGTGACACACGATAAAGACGAAATTAATCATGCGGATAAAGt                                                                                                   <  1:523618/149‑1 (MQ=255)
                                                     ttATTGCTGGATGAAATAACCAGTGCGCTGGATGAAAGTAATAAACATAACGTCAATGAGATGATCCATCGTTATGTGCGCGAGCAAAATATTGCCGTGCTGTGGGTGACACACGATAAAGACGAAATTAATCATGCGGATAAAGTGAt                                                                                                <  3:308307/149‑1 (MQ=255)
                                                      tATTGCTGGATGAAATAACCAGTGCGCTGGATGAAAGTAATAAACATAACGTCAATGAGATGATCCATCGTTATGTGCGCGAGCAAAATATTGCCGTGCTGTGGGTGACACACGATAAAGACGAAATTAATCATGCGGATAAAGTGAtt                                                                                               <  2:585852/149‑1 (MQ=255)
                                                          gCTGGATGAAATAACCAGTGCGCTGGATGAAAGTAATAAACATAACGTCAATGAGATGATCCATCGTTATGTGCGCGAGCAAAATATTGCCGTGCTGTGGGTGACAAACGATAAAGACGAAATTAATCATGCGGATAAAGTGATTacac                                                                                           <  3:548707/149‑1 (MQ=255)
                                                           cTGGATGAAATAACCAGTGCGCTGGATGAAAGTAATAAACATAACGTCAATGAGATGATCCATCGTTATGTGCGCGAGCAAAATATTGCCGTGCTGTGGGTGACACACGATAAAGACGAAATTAATCATGCGGATAAAGTGATTACACt                                                                                          <  3:9680/149‑1 (MQ=255)
                                                            tGGATGAAATAACCAGTGCGCTGGATGAAAGTAATAAACATAACGTCAATGAGATGATCCATCGTTATGTGCGCGAGCAAAATATTGCCGTGCTGTGGGTGACACACGATAAAGACGAAATTAATCATGCGGATAAAGTGATTACACTg                                                                                         >  3:442996/1‑149 (MQ=255)
                                                              gATGAAATAACCAGTGCGCTGGATGAAAGTAATAAACATAACGTCAATGAGATGATCCATCGTTATGTGCGCGAGCAAAATATTGCCGTGCTGTGGGTGACACACGATAAAGACGAAATTAATCATGCGGATAAAGTGATTACACTGCa                                                                                       >  3:233543/1‑149 (MQ=255)
                                                                tGAAATAACCAGTGCGCTGGATGAAAGTAATAAACATAACGTCAATGAGATGATCCATCGTTATGTGCGCGAGCAAAATATTGCCGTGCTGTGGGTGACACACGATAAAGACGAAATTAATCATGCGGATAAAGTGATTACACTGCAAc                                                                                     >  3:359517/1‑149 (MQ=255)
                                                                tGAAATAACCAGTGCGCTGGATGAAAGTAATAAACATAACGTCAATGAGATGATCCATCGTTATGTGCGCGAGCAAAATATTGCCGTGCTGTGGGTGACACACGATAAAGACGAAATTAATCATGCGGATAAAGTGATTACACTGCAAc                                                                                     >  4:171183/1‑149 (MQ=255)
                                                                    aTAACCAGTGCGCTGGATGAAAGTAATAAACATAACGTCAATGAGATGATCCATCGTTATGTGCGCGAGCAAAATATTGCCGTGCTGTGGGTGACACACGATAAAGACGAAATTAATCATGCGGATAAAGTGATTACACTGCAACCGCa                                                                                 >  2:342281/1‑149 (MQ=255)
                                                                          aGTGCGCTGGATGAAAGTAATAAACATAACGTCAATGAGATGATCCATCGTTATGTGCGCGAGCAAAATATTGCCGTGCTGTGGGTGACACACGATAAAGACGAAATTAATCATGCGGATAAAGTGATTACACTGCAACCGCATGCCgg                                                                           >  4:539046/1‑149 (MQ=255)
                                                                              cgcTGGATGAAAGTAATAAACATAACGTCAATGAGATGATCCATCGTTATGTGCGCGAGCAAAATATTGCCGTGCTGTGGGTGACACACGATAAAGACGAAATTAATCATGCGGATAAAGTGATTACACTGCAACCGCATGCCGGAGaa                                                                       >  1:404197/1‑149 (MQ=255)
                                                                                cTGGATGAAAGTAATAAACATAACGTCAATGAGATGATCCATCGTTATGTGCGCGAGCAAAATATTGTCGTGCTGTGGGTGACACACGATAAAGACGAAATTAATCATGCGGATAAAGTGATTACACTGCAACCGCATGCCGGAGAAAt                                                                     <  1:287109/149‑1 (MQ=255)
                                                                                cTGGATGAAAGTAATAAACATAACGTCAATGAGATGATCCATCGTTATGTGCGCGAGCAAAATATTGCCGTGCTGTGGGTGACACACGATAAAGACGAAATTAATCATGCGGATAAAGTGATTACACTGCAACCGCATGCCGGAGAAAt                                                                     <  1:148049/149‑1 (MQ=255)
                                                                                  ggATGAAAGTAATAACCATAACGACAAAGAGATGATCCATCGTTATGTGCGCGAGCAAAATACGGCCGTGCTGTTGGTGACAAACGATAAAGACGAAATTAATCATGCGGATAAAGTGATTACACTGCAACCGCATGCCGGAGAAACGc                                                                   <  2:523102/149‑1 (MQ=255)
                                                                                  ggATGAAAGTAATAAACATAACGTCAATGAGATGATCCATCGTTATGTGCGCGAGCAAAATATTGCCGTGCTGTGGGTGACACACGATAAAGACGAAATTAATCATGCGGATAAAGTGATTACACTGCAACCGCATGCCGGAGAAATGc                                                                   <  1:202248/149‑1 (MQ=255)
                                                                                   gATGAAAGTAATAAACATAACGTCAATGAGATGATCCATCGTTATGTGCGCGAGCAAAATATTGCCGTGCTGTGGGTGACACACGATAAAGACGAAATTAATCATGCGGATAAAGTGATTACACTGCAACCGCATGCCGGAGAAATGCa                                                                  <  1:342281/149‑1 (MQ=255)
                                                                                   gATGAAAGTAATAAACATAACGTCAATGAGATGATCCATCGTTATGTGCGCGAGCAAAATATTGCCGTGCTGTGGGTGACACACGATAAAGACGAAATTAATCATGCGGATAAAGTGATTACACTGCAACCGCATGCCGGAGAAATGCa                                                                  <  4:275147/149‑1 (MQ=255)
                                                                                   gATGAAAGTAATAAACATAACGTCAATGAGATGATCCATCGTTATGTGCGCGAGCAAAATATTGCCGTGCTGTGGGTGACACACGATAAAGACGAAATTAATCATGCGGATAAAGTGATTACACTGCAACCGCATGCCGGAGAAATGCa                                                                  >  2:306370/1‑149 (MQ=255)
                                                                                   gATGAAAGTAATAAACATAACGTCAATGAGATGATCCATCGTTATGTGCGCGAGCAAAATATTGCCGTGCTGTGGGTGACACACGATAAAGACGAAATTAATCATGCGGATAAAGTGATTACACTGCAACCGCATGCCGGAGAAATGCa                                                                  <  3:276126/149‑1 (MQ=255)
                                                                                    aTGAAAGTAATAAACATAACGTCAATGAGATGATCCATCGTTATGTGCGCGAGCAAAATATTGCCGTGCTGTGGGTGACACACGATAAAGACGAAATTAATCATGCGGATAAAGTGATTACACTGCAACCGCATGCCGGAGAAATGCAg                                                                 <  1:149374/149‑1 (MQ=255)
                                                                                     tGAAAGTAATAAACATAACGTCAATGAGATGATCCATCGTTATGTGCGCGAGCAAAATATTGCCGTGCTGTGGGTGACACACGATAAAGACGAAATTAATCATGCGGATAAAGTGATTACACTGCAACCGCATGCCGGAGAAATGCAgg                                                                <  3:137922/149‑1 (MQ=255)
                                                                                       aaaGTAATAAACATAACGTCAATGAGATGATCCATCGTTATGTGCGCGAGCAAAATATTGTCGTGCTGTGGGTGACACACGATAAAGACGAAATTAATCATGCGGATAAAGTGATTACACTGCAACCGCATGCCGGAGAAATGCAGGaa                                                              <  1:575076/149‑1 (MQ=255)
                                                                                         aGTAATAAACATAACGTCAATGAGATGATCCATCGTTATGTGCGCGAGCAAAATATTGCCGTGCTGTGGGTGACACACGATAAAGACGAAATTAATCATGCGGATAAAGTGATTACACTGCAACCGCATGCCGGAGAAATGCAGGAAGc                                                            <  4:515095/149‑1 (MQ=255)
                                                                                           taataaACATAACGTCAATGAGATGATCCATCGTTATGTGCGCGAGCAAAATATTGCCGTGCTGTGGGTGACACACGATAAAGACGAAATTAATCATGCGGATAAAGTGATTACACTGCAACCGCATGCCGGAGAAATGCAGGAAGCAc                                                          <  1:92860/149‑1 (MQ=255)
                                                                                            aataaACATAACGTCAATGAGATGATCCATCGTTATGTGCGCGAGCAAAATATTGCCGTGCTGTGGGTGACACACGATAAAGACGAAATTAATCATGCGGATAAAGTGATTACACTGCAACCGCATGCCGGATAAATGCAGGAAGCACg                                                         >  4:381712/1‑149 (MQ=255)
                                                                                               aaaCATAACGTCAATGAGATGATCCATCGTTATGTGCGCGAGCAAAATATTGCCGTGCTGTGGGTGACACACGATAAAGACGAAATTAATCATGCGGATAAAGTGATTACACTGCAACCGCATGCCGGAGAAATGCAGGAAGCACGCTa                                                      <  4:442996/149‑1 (MQ=255)
                                                                                                  cATAACGTCAATGAGATGATCCATCGTTATGTGCGCGAGCAAAATATTGCCGTGCTGTGGGTGACACACGATAAAGACGAAATTAATCATGCGGATAAAGTGATTACACTGCAACCGCATGCCGGAGAAATGCAGGAAGCACGCTATGa                                                   <  3:252715/149‑1 (MQ=255)
                                                                                                      aCGTCAATGAGATGATCCATCGTTATGTGCGCGAGCAAAATATTGCCGTGCTGTGGGTGACACACGATAAAGACGAAATTAATCATGCGGATAAAGTGATTACACTGCAACCGCATGCCGGAGAAATGCAGGAAGCACGCTATGAACTc                                               <  4:19151/149‑1 (MQ=255)
                                                                                                       cGTCAATGAGATGATCCATCGTTATGTGCGCGAGCAAAATATTGCCGTGCTGTGGGTGACACACGATAAAGACGAAATTAATCATGCGGATaaa                                                                                                     >  2:449589/1‑94 (MQ=255)
                                                                                                       cGTCAATGAGATGATCCATCGTTATGTGCGCGAGCAAAATATTGCCGTGCTGTGGGTGACACACGATAAAGACGAAATTAATCATGCGGATaaa                                                                                                     <  1:449589/94‑1 (MQ=255)
                                                                                                         tCAATGAGATGATCCATCGTTATGTGCGCGAGCAAAATATTGCCGTGCTGTGGGTGACACACGATAAAGACGAAATTAATCATGCGGATAAAGTGATTACACTGCAACCGCATGCCGGAGAAATGCAGGAAGCACGCTATGAACTCGCa                                            <  4:385951/149‑1 (MQ=255)
                                                                                                             tGAGATGATCCATCGTTATGTGCGCGAGCAAAATATTGCCGTGCTGTGGGTGACACACGATAAAGACGAAATTAATCATGCGGATAAAGTGATTACACTGCAACCGCATGCCGGAGAAATGCAGGAAGCACGCTATGAACTCGCataat                                        >  4:99767/1‑149 (MQ=255)
                                                                                                                    atCCATCGTTATGTGCGCGAGCAAAATATTGCCGTGCTGTGGGTGACACACGATAAAGACGAAATTAATCATGCGGATAAAGTGATTACACTGCAACCGCATGCCGGAGAAATGCAGGAAGCACGCTATGAACTCGCataat                                        <  2:284511/142‑1 (MQ=255)
                                                                                                                    atCCATCGTTATGTGCGCGAGCAAAATATTGCCGTGCTGTGGGTGACACACGATAAAGACGAAATTAATCATGCGGATAAAGTGATTACACTGCAACCGCATGCCGGAGAAATGCAGGAAGCACGCTATGAACTCGCataat                                        >  1:284511/1‑142 (MQ=255)
                                                                                                                     tCCATCGTTATGTGCGCGAGCAAAATATTGTCGTGCTGTGGGTGACACACGATAAAGACGAAATTAATCATGCGGATAAAGTGATTACACTGCAACCGCATGCCGGAGAAATGCAGGAAGCACGCTATGAACTCGCATAATATTACTaa                                >  2:484664/1‑149 (MQ=255)
                                                                                                                     tCCATCGTTATGTGCGCGAGCAAAATATTGCCGTGCTGTGGGTGACACACGATAAAGACGAAATTAATCATGCGGATAAAGTGATTACACTGCAACCGCATGCCGGAGAAATGCAGGAAGCACGCTATGAACTCGCATAATATTACTaa                                >  1:35110/1‑149 (MQ=255)
                                                                                                                     tCCATCGTTATGTGCGCGAGCAAAATATTGCCGTGCTGTGGGTGACACACGATAAAGACGAAATTAATCATGCGGATAAAGTGATTACACTGCAACCGCATGCCGGAGAAATGCAGGAAGCACGCTATGAACTCGCATAATATTACTaa                                >  4:415072/1‑149 (MQ=255)
                                                                                                                       cATCGTTATGTGCGCGAGCAAAATATTGCCGTGCTGTGGGTGACACACGATAAAGACGAAATTAATCATGCGGATAAAGTGATTACACTGCAACCGCATGCCGGAGAAATGCAGGAAGCACGCTATGAACTCGCATAATATTACTAACg                              >  1:8756/1‑149 (MQ=255)
                                                                                                                        aTCGTTATGTGCGCGAGCAAAATATTGCCGTGCTGTGGGTGACACACGATAAAGACGAAATTAATCATGCGGATAAAGTGATTACACTGCAACCGCATGCCGGAGAAATGCAGGAAGCACGCTATGAACTCGCATAATATTACTAACGa                             >  2:164249/1‑149 (MQ=255)
                                                                                                                            ttATGTGCGCGAGCAAAATATTGCCGTGCTGTGGGTGACACACGATAAAGACGAAATTAATCATGCGGATAAAGTGATTACACTGCAACCGCATGCCGGAGAAATGCAGGAAGCACGCTATGAACTCGCATAATATTACTAACGAATCa                         <  1:472624/149‑1 (MQ=255)
                                                                                                                             tATGTGCGCGAGCAAAATATTGCCGTGCTGTGGGTGACACACGATAAAGACGAAATTAATCATGCGGATAAAGTGATTACACTGCAACCGCATGCCGGAGAAATGCAGGAAGCACGCTATGAACTCGCATAATATTACTAACGAATCAt                        <  4:254852/149‑1 (MQ=255)
                                                                                                                               tgtgCGCGAGCAAAATATTGCCGTGCTGTGGGTGACACACGATAAAGACGAAATTAATCATGCGGATAAAGTGATTacac                                                                                           <  1:99234/80‑1 (MQ=255)
                                                                                                                               tgtgCGCGAGCAAAATATTGCCGTGCTGTGGGTGACACACGATAAAGACGAAATTAATCATGCGGATAAAGTGATTacac                                                                                           >  2:99234/1‑80 (MQ=255)
                                                                                                                                   cgcgAGCAAAATATTGCCGTGCTGTGGGTGACACACGATAAAGACGAAATTAATCATGCGGATAAAGTGATTACACTGCAACCGCATGCCGGAGAAATGCAGGAAGCACGCTATGAACTCGCATAATATTACTAACGAATCATTAGCAc                  >  2:119757/1‑149 (MQ=255)
                                                                                                                                            aaTATTGCCGTGCTGTGGGTGACACACGATAAAGACGAAATTAATCATGCGGATAAAGTGATTACACTGCAACCGCATGCCGGAGAAATGCAGGAAGCACGCTATGAACTCGCATAATATTACTAACGAATCATTAGCACTGGCATTaa         >  1:79317/1‑149 (MQ=255)
                                                                                                                                                   ccGTGCTGTGGGTGACACACGATAAAGACGAAATTAATCATGCGGATAAAGTGATTACACTGCAACCGCATGCCGGAGAAATGCAGGAAGCACGCTATGAACTCGCATAATATTACTAACGAATCATTAGCACTGGCATTAATGCtggt  <  3:171183/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
GTGAAAAACAACGCATCTCATTGATTCGTAACCTGCAATTTATGCCGAAGGTTTTATTGCTGGATGAAATAACCAGTGCGCTGGATGAAAGTAATAAACATAACGTCAATGAGATGATCCATCGTTATGTGCGCGAGCAAAATATTGCCGTGCTGTGGGTGACACACGATAAAGACGAAATTAATCATGCGGATAAAGTGATTACACTGCAACCGCATGCCGGAGAAATGCAGGAAGCACGCTATGAACTCGCATAATATTACTAACGAATCATTAGCACTGGCATTAATGCTGGT  >  NC_000913/516340‑516635

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