Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,543,413 A→C 7.3% L39R (CTG→CGG)  cysP ← thiosulfate‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,543,4130AC 7.3% 181.6 / 5.9 65L39R (CTG→CGG) cysPthiosulfate‑binding protein
Reads supporting (aligned to +/- strand):  ref base A (30/30);  new base C (3/2);  total (33/32)
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.12e-01

TATAAGTGACAACGTCGGCTTTTAAGCCCTGTAAAATCGCCAGCGCCTGTTTTGATGACCCGGCATGAGATTGTTTTATCGTCAGTTTGTCGCCGCCGTTATCTTTTGCCCATTGTTGCTCAAACGGCGGATTCAGGGCGGCAAACAGCTCGCGGGAGACGTCATAAGAACTGTTCAGCAGTTCCGTTGCCTGTACATGGCCCGCCAGCAGCAGAGAAGCGACCAGCGCGAGTGAGTTCTTTTTCAGTAAGTTAACGGCCATTGCGCACCCTTATAAATTTAATGACTTTC  >  NC_000913/2543267‑2543557
                                                                                                                                                  |                                                                                                                                                
tataAGTGACAACGTCGGCGTTTAAGCCCTGTAAAAGCGCCAGCGCCTGTTTTGATGACCCGGCATGAGATTGTTTTATCGTCAGTTTGTCGCCGCCGTTATCTTTTGCCCATTGTTGCGCAAACGGCGGATTCAGGGCGGCAAACAGc                                                                                                                                                <  2:112559/149‑1 (MQ=255)
 ataAGTGACAACGTCGGCTTTTAAGCCCTGTAAAATCGCCAGCGCCTGTTTTGATGACCCGGCATGAGATAGTATTATCGTCAGTTAGACGCCGCCGGTAAAGTTTGCCCATTGTTGCTCAAACGGCGGATACATGGAGGCAAACAGCt                                                                                                                                               >  6:144629/1‑149 (MQ=255)
  taAGTGACAACGTCGGCTTTTAAGCCCTGTAAAATCGCCAGCGCCTGTTTTGATGACCCGGCATGAGATTGTTTTATCGTCAGTTTGTCGCCGCCGTTATCTTTTGCCCATTGTTGCTCAAACGGCGGATTCAGGGCGGCAAACAGCt                                                                                                                                               <  8:42602/148‑1 (MQ=255)
  taAGTGACAACGTCGGCTTTTAAGCCCTGTAAAATCGCCAGCGCCTGTTTTGATGACCCGGCATGAGATTGTTTTATCGTCAGTTTGTCGCCGCCGTTATCTTTTGCCCATTGTTGCTCAAACGGCGGATTCAGGGCGGCAAACAGCTc                                                                                                                                              >  3:209027/1‑149 (MQ=255)
     gTGACAACGTCGGCTTTTAAGCCCTGTAAAATCGCCAGCGCCTGTTTTGATGACCCGGCATGAGATTGTTTTATCGTCAGTTTGTCGCCGCCGTTATCTTTTGCCCATTGTTGCTCAAACGGCGGATTCAGGGCGGCAAACAGCTcgcg                                                                                                                                           >  3:163864/1‑149 (MQ=255)
         cAACGTCGGCTTTTAAGCCCTGTAAAATCGCCAGCGCCTGTTTTGATGACCCGGCATGAGATTGTTTTATCGTCAGTTTGTCGCCGCCGTTATCTTTTGCCCATTGTTGCTCAAACGGCGGATTCAGGGCGGCAAACAGCTCGCgg                                                                                                                                          >  5:136070/1‑146 (MQ=255)
           aCGTCGGCGTTTCAGCCCTGTAAACTCGCCAGCGCCTGATTTGATGACCCGGCATGAGATTGTTTTAGCGTCAGTTTGTCGCCGCCGTTATCTTTTGCCCATTGTTGCTCAAACGGCGGATTCAGGGCGGCAAACAGCTCGCGGGAGAc                                                                                                                                     <  7:228656/149‑1 (MQ=255)
             gTCGGCTTTTAAGCCCTGTAAAATCGCCAGCGCCTGTTTTGATGACCCGGCATGAGATTGTTTTATCGTCAGTTTGTCGCCGCCGTTATCTTTTGCCCATTGTTGCTCAAACGGCGGATTCAGGGCGGCAAACAGCTCGCGGGAGACGt                                                                                                                                   <  3:212167/149‑1 (MQ=255)
              tCGGCTTTTAAGCCCTGTAAAATCGCCAGCGCCTGTTTTGATGACCCGGCATGAGATTGTTTTATCGTCAGTTTGTCGCCGCCGTTATCTTTTGCCCATTGTTGCTCAAACGGCGGATTCAGGGCGGCAAACAGCTCGCGGGAGACGt                                                                                                                                   >  7:220512/1‑148 (MQ=255)
                      tAAGCCCTGTAAAATCGCCAGCGCCTGTTTTGATGACCCGGCATGAGATTGTTTTATCGTCAGTTTGTCGCCGCCGTTATCTTTTGCCCATTGTTGCTCAAACGGCGGATTCAGGGCGGCAAACAGCTCGCGGGAGACGTCATAAGaa                                                                                                                           <  2:84855/148‑1 (MQ=255)
                        aGCCCTGTAAAATCGCCAGCGCCTGTTTTGATGACCCGGCATGAGATTGTTTTATCGTCAGTTTGTCGCCGCCGTTATCTTTTGCCCATTGTTGCTCAAACGGCGGATTCAGGGCGGCAAACAGCTCGCGGGAGACGTCATAAGAACTg                                                                                                                        <  3:116521/149‑1 (MQ=255)
                         ggcgtGTAATATCGCCAGCGCCTGTTTTGATGACCCGGCATGAGATTGTTTTATGGTCAGTTCGTCGCGGCCGTTATCTTTTGCCCATTGTTGCTCAAACGGCGGATTCAGGGCGGCAAACAGCTCGCGGGAGACGTCATAAGAACTGt                                                                                                                       <  5:144629/145‑1 (MQ=255)
                           ccTGTAAAATCGCCAGCGCCTGTTTTGATGACCCGGCATGAGATTGTTTTATCGTCAGTTTGTCGCCGCCGTTATCTTTTGCCCATTGTTGCTCAAACGGCGGATTCAGGGCGGCAAACAGCTCGCGGGAGACGTCATAAGAACTGtt                                                                                                                      >  3:55197/1‑148 (MQ=255)
                             tGTAAAATCGCCAGCGCCTGTTTTGATGACCCGGCATGAGATTGTTTTATCGTCAGTTTGTCGCCGCCGTTATCTTTTGCCCATTGTTGCTCAAACGGCGGATTCAGGGCGGCAAACAGCTCGCGGGAGACGTCa                                                                                                                                 >  3:121751/1‑135 (MQ=255)
                             tGTAAAATCGCCAGCGCCTGTTTTGATGACCCGGCATGAGATTGTTTTATCGTCAGTTTGTCGCCGCCGTTATCTTTTGCCCATTGTTGCTCAAACGGCGGATTCAGGGCGGCAAACAGCTCGCGGGAGACGTCa                                                                                                                                 <  4:121751/135‑1 (MQ=255)
                             tGTAAAATCGCCAGCGCCTGTTTTGATGACCCGGCATGAGATTGTTTTATCGTCAGTTTGTCACCGCCGTTATCTTTTGCCCATTGTTGCTCAAACGGCGGATTCAGGGCGGCAAACAGCTCGCGGGAGACGTCATAAGAACTGTTcag                                                                                                                   >  6:55934/1‑149 (MQ=255)
                               tAAAATCGCCAGCGCCTGTTTTGATGACCCGGCATGAGATTGTTTTATCGTCAGTTTGTCGCCGCCGTTATCTTTTGCCCATTGTTGCTCAAACGGCGGATTCAGGGCGGCAAACCGCTCGCGGGAGACGTCATAAGAACTGTTcag                                                                                                                   >  5:239019/1‑147 (MQ=255)
                                         aGCGCCTGTTTTGATGACCCGGCATGAGATTGTTTTATCGTCAGTTTGTCGCCGCCGTTATCTTTTGCCCATTGTTGCTCAAACGGCGGATTCAGGGCGGCAAACAGCTCGCGGGAGACGTCATAAGAACTGTTCAGCAGTTCCGTTGc                                                                                                       >  2:92465/1‑149 (MQ=255)
                                         aGCGCCTGTTTTGATGACCCGGCATGAGATTGTTTTATCGTCAGTTTGTCGCCGCCGTTATCTTTTGCCCATTGTTGCTCAAACGGCGGATTCAGGGCGGCAAACAGCTCGCGGGAGACGTCATAAGAACTGTTCAGCAGTTCCGTTGc                                                                                                       <  4:55197/149‑1 (MQ=255)
                                         aGCGCCTGTTTTGATGACCCGGCATGAGATTGTTTTATCGTCAGTTTGTCGCCGCCGTTATCTTTTGCCCATTGTTGCTCAAACGGAGGATTCAGGGCGGCAAACAGCTCGCGGGAGACGTCATAAGAACTCTACAGCACTTCCGTTGc                                                                                                       >  8:149440/1‑149 (MQ=255)
                                           cgcCTGTTTTGATGACCCGGCATGAGATTGTTTTATCGTCAGTTTGTCGCCGCCGTTATCTTTTGCCCATTGTTGCTCAAACGGCGGATTCAGGGCGGCAAACAGCTCGCGGGAGACGTCATAAGAACTGtt                                                                                                                      >  7:244755/1‑132 (MQ=255)
                                           cgcCTGTTTTGATGACCCGGCATGAGATTGTTTTATCGTCAGTTTGTCGCCGCCGTTATCTTTTGCCCATTGTTGCTCAAACGGCGGATTCAGGGCGGCAAACAGCTCGCGGGAGACGTCATAAGAACTGtt                                                                                                                      <  8:244755/132‑1 (MQ=255)
                                            gcCTGTTTTGATGACCCGGCATGAGATTGTTTTATCGTCAGTTTGTCGCCGCCGTTATCTTTTGCCCATTGTTGCTCAAACGGCGGATTCAGGGCGGCAAACAGCTCGCGGGAGACGTCATAAGAACTGTTCAGCAGTTCCGTTGCCTg                                                                                                    >  5:154474/1‑149 (MQ=255)
                                             ccTGTTTTGATGTCCCGGCATGAGATTGTTTTTTCGTGAGTTTGTCGCCGCCGTTATCTTTTGCCCATTGTTGCTCAAACGGCGGATTCAGGGCGGCAAACAGCTCGCGGGAGACGTCATAAGAACTGTTCAGCAGTTCCGTTGCCTGt                                                                                                   <  6:154474/149‑1 (MQ=255)
                                              ctttttgGATGTCCCGGCATGAGATTGTTTTAGCGTAAGTTTGTCGCCGCCGTTATCTTTTGGCCATTGTTGCTCAACCGGCGGTTTCAGGGCGTCAAACAGCTCGCGGGAGACGGCATAAGAACTGTTCAGCAGTTCCGTTGCCTGTa                                                                                                  <  6:203395/146‑1 (MQ=255)
                                              cTGTTTTGATGACCCGGCATGAGATTGTTTTATCGTCAGTTTGTCGCCGCCGTTATCTTTTGCCCATTGTTGCTCAAACGGCGGATTCAGGGCGGCAAACAGCTCGCGGGAGACGTCATAAGAACTGTTCAGCAGTTCCGTTGCCTGTa                                                                                                  <  6:136070/149‑1 (MQ=255)
                                                 ttttGATGACCCGGCATGAGATTGTTTTATCGTCAGTTTGTCGCCGCCGTTATCTTTTGCCCATTGTTGCTCAAACGGCGGATTCAGGGCGGCAAACAGCTCGCGGGAGACGTCATAAGAACTGTTCAGCAGTTCCGTTGCCTGTACa                                                                                                <  8:220512/148‑1 (MQ=255)
                                                   ttGATGACCCGGCATGAGATTGTTTTATCGTCAGTTTGTCGCCGCCGTTATCTTTTGCCCATTGTTGCTCAAACGGCGGATTCAGGGCGGCAAACAGCTCGCGGGAGACGTCATAAGAACTGTTCAGCAGTTCCGTTGCCTGTACATgg                                                                                             >  7:209162/1‑149 (MQ=255)
                                                      atgaCCCGGCATGAGATTGTTTTATCGTCAGTTTGTCGCCGCCGTTATCTTTTGCCCATTGTTGCTCAAACGGCGGATTCAGGGCGGCAAACAGCTCGCGGGAGACGTCATAAGAACTg                                                                                                                        >  4:174863/1‑119 (MQ=255)
                                                      atgaCCCGGCATGAGATTGTTTTATCGTCAGTTTGTCGCCGCCGTTATCTTTTGCCCATTGTTGCTCAAACGGCGGATTCAGGGCGGCAAACAGCTCGCGGGAGACGTCATAAGAACTg                                                                                                                        <  3:174863/119‑1 (MQ=255)
                                                        gaCCCGGCATGAGATTGTTTTATCGTCAGTTTGTCGCCGCCGTTATCTTTTGCCCATTGTTGCTCAAACGGCGGATTCAGGGCGGCAAACAGCTCGCGGGAGACGTCATAAGAACTGTTCAGCAGTTCCGTTGCCTGTACATGGCCCGc                                                                                        >  7:115184/1‑149 (MQ=255)
                                                              gCATGAGATTGTTTTATCGTCTGTTTTTCGCCGCCGTTATCTTTTGCCCATTGTTGCTCAAACGGCGGATTCAGGGCGGCAAACAGCTCGCGGGAGACGTCATAAGAACTGTTCAGCAGTTCCGTTGCCTGTACATGGCCCGCcagcag                                                                                  <  7:149440/149‑1 (MQ=255)
                                                                 tGAGATTGTTTTATCGTCAGTTTGTCGCCGCCGTTATCTTTTGCCCATTGTTGCTCAAACGGCGGATTCAGGTCGGCATCCAGCTCGCGGGAGACGTCATAAGAACTGTTCAGCAGGACCGTTGCCTGTACATGCCCAGCcagcag                                                                                  >  8:97089/1‑146 (MQ=255)
                                                                     aTTGTTTTATCGTCAGTTTGTCGCCGCCGTTATCTTTTGCCCATTGTTGCTCAAACGGCGGATTCAGGGCGGCAAACAGCTCGCGGGAGCCGTCATAAGAACTGTTCAGCAGTTCCGTTGCCTGTACATGGCCCGCCAGCAGCagag                                                                             >  2:90366/1‑147 (MQ=255)
                                                                           ttATCGTCAGTTTGTCGCCGCCGTTATCTTTTGCCCATTGTTGCTCAAACGGCGGATTCAGGGCGGCAAACCGCTCGCGGGAGACGTCATaa                                                                                                                              >  1:195418/1‑92 (MQ=255)
                                                                           ttATCGTCAGTTTGTCGCCGCCGTTATCTTTTGCCCATTGTTGCTCAAACGGCGGATTCAGGGCGGCAAACCGCTCGCGGGAGACGTCATaa                                                                                                                              <  2:198600/92‑1 (MQ=255)
                                                                           ttATCGTCAGTTTGTCGCCGCCGTTATCTTTTGCCCATTGTTGCTCAAACGGCGGATTCAGGGCGGCAAACCGCTCGCGGGAGACGTCATaa                                                                                                                              <  2:195418/92‑1 (MQ=255)
                                                                           ttATCGTCAGTTTGTCGCCGCCGTTATCTTTTGCCCATTGTTGCTCAAACGGCGGATTCAGGGCGGCAAACCGCTCGCGGGAGACGTCATaa                                                                                                                              >  1:198600/1‑92 (MQ=255)
                                                                             aTCGTCAGTTTGTCGCCGCCGTTATCTTTTGCCCATTGTTGCTCAAACGGCGGATTCAGGGCGGCAAACAGCTCGCGGGAGACGTCATAAGAACTGTTCAGCAGTTCCGTTGCCTGTACATGGCCCGCCAGCAGCAGAGAAGCGACCa                                                                    <  1:92465/148‑1 (MQ=255)
                                                                                gTCAGTTTGTCGCCGCCGTTATCTTTTGCCCATTGTTGCTCAAACGGCGGATTCAGGGCGGCAAACAGCTCGCGGGAGACGTCATAAGAACTGTTCAGCAGTTCCGTTGCCTGTACATGGCCCGCCAGCAGCAGAGAAGCGACCAgcgc                                                                >  3:83890/1‑149 (MQ=255)
                                                                                 tCAGTTTGTCGCCGCCGTTATCTTTTGCCCATTGTTGCTCAAACGGCGGCTTCAGGGCGGCAAACAGCTCGCGGGAGCCGTCATAAGAACTGTTCAGCAGTTCCGTTGCCTGTACATGGCCCGCCCGCAGCAGAGCAGCGACCAgcgc                                                                >  7:124363/1‑148 (MQ=255)
                                                                                             cgccgTTATCTTTTGCCCATTGTTGCTCAAACGGCGGATTCAGGGCGGCAAACAGCTCGCGGGAGACGTCATAAGAACTGTTCAGCAGTTCCGTTGCCTGTACATGGCCCGCCAGCAGCAGAGAAGCGACCAGCGCGAGTGAGTTCttt                                                   <  8:115184/149‑1 (MQ=255)
                                                                                              gccgTTATCTTTTGCCCATTGTTGCTCAAACGGCGGATTCAGGGCGGCAAACAGCTCGCGGGAGACGTCATAAGAACTGTTCAGCAGTTCCGTTGCCTGTACATGGCCCGCCAGCAGCAGAGAAGCGACCAGCGCGAGTGAGTTCtttt                                                  >  7:38453/1‑149 (MQ=255)
                                                                                               ccgTTATCTTTTGCCCATTGTTGCTCAAACGGCGGATTCAGGGCGGCAAACAGCTCGCGGGAGACGTCATAAGAACTGTTCAGCAGTTCCGTTGCCTGTACATGGCCCGCCAGCAGCAGAGAAGCGACCAGCGCGAGTGAGTTCttttt                                                 <  4:209027/149‑1 (MQ=255)
                                                                                               ccgTTATCTTTTGCCCATTGTTGCTCAAACGGCGGATTCAGGGCGGCAAACAGCTCGCGGGAGACGTCATAAGAACTGTTCAGCAGTTCCGTTGCCTGTACATGGCCCGCCAGCAGCAGAGAAGCGACCAGCGCGAGTGAGTTCttttt                                                 <  7:168024/149‑1 (MQ=255)
                                                                                                cgTTATCTTTTGCCCATTGTTGCTCAAACGGCGGATTCAGGGCGGCAAACAGCTCGCGGGAGACGTCATAAGAACTGTTCAGCAGTTCCGTTGCCTGTACATGGCCCGCCAGCAGCAGAGAAGCGACCAGCGCGAGTGAGTTCTTTTTc                                                <  8:209162/149‑1 (MQ=255)
                                                                                                 gTTATCTTTTGCCCATTGTTGCTCAAACGGCGGATTCAGGGCGGCAAACAGCTCGCGGGAGACGCCATAAGAACTGTTCAGCAGTTCCGTTGCCTGTACATGGCCCGCCAGCAGCAGAGAAGCGACCAGCGCGAGTGAGTTCTTTTTCa                                               <  4:163864/149‑1 (MQ=255)
                                                                                                  ttATCTTTTGCCCATTGTTGCTCAAACGGCGGATTCAGGGCGGCAAACAGCTCGCGGGAGACGTCATAAGAACTGTTCAGCAGTTCCGTTGCCTGTACATGGCCCGCCAGCAGCAGAGAAGCGACCAGCGCga                                                              >  6:198413/1‑133 (MQ=255)
                                                                                                  ttATCTTTTGCCCATTGTTGCTCAAACGGCGGATTCAGGGCGGCAAACAGCTCGCGGGAGACGTCATAAGAACTGTTCAGCAGTTCCGTTGCCTGTACATGGCCCGCCAGCAGCAGAGAAGCGACCAGCGCga                                                              <  5:198413/133‑1 (MQ=255)
                                                                                                       ttttGCCCATTGTTGCTCAAACGGCGGATTCAGGGCGGCAAACAGCTCGCGGGAGACGTCATAAGAACTGTTCAGCAGTTCCGTTGCCTGTACATGGCCCGCCAGCAGCAGAGAAGCGACCAGCGCGAGTGAGTTCTTTTTCAGTAAGt                                         >  4:145179/1‑149 (MQ=255)
                                                                                                             ccATTGTTGCTCAAACGGCGGATTCAGGGCGGCAAACAGCTCGCGGGAGACGTCATAAGAACTGTTCAGCAGTTCCGTTGCCTGTACATGGCCCGCCAGCAGCAGAGAAGCGACCAGCGCGAGTGAGTTCTTTTTCa                                               >  7:100842/1‑137 (MQ=255)
                                                                                                             ccATTGTTGCTCAAACGGCGGATTCAGGGCGGCAAACAGCTCGCGGGAGACGTCATAAGAACTGTTCAGCAGTTCCGTTGCCTGTACATGGCCCGCCAGCAGCAGAGAAGCGACCAGCGCGAGTGAGTTCTTTTTCa                                               <  8:100842/137‑1 (MQ=255)
                                                                                                                       tCAAACGGCGGATTCAGGGCGGCAAACAGCTCGCGGGAGACGTCATAAGAACTGTTCAGCAGTTCCGTTGCCTGTACATGGCCCGCCAGCAGCAGAGAAGCGACCAGCGCGAGTGAGTTCTTTTTCAGTAAGTTAACGGCCATTGCGCa                         >  5:214996/1‑149 (MQ=255)
                                                                                                                       tCAAACGGCGGATTCAGGGCGGCAAACAGCTCGCGGGAGACGTCATAAGAACTGTTCAGCAGTTCCGTTGCCTGTACATGGCCCGCCAGCAGCAGAGAAGCGACCAGCGCGAGTGAGTTCTTTTTCAGTAAGTTAACGGCCATTGCGCa                         >  6:13168/1‑149 (MQ=255)
                                                                                                                         aaaCGGCGGATTCAGGGCGGCAAACAGCTCGCGGGAGACGTCATAAGAACTGTTCAGCAGTTCCGTTGCCTGTACATGGCCCGCCAGCAGCAGAGCAGCGACCAGCGCGAGTGAGTTCTTTTTCAGTAAGTTAACGGCCATTGCGCAcc                       >  3:176605/1‑149 (MQ=255)
                                                                                                                         aaaCGGCGGATTCAGGGCGGCAAACAGCTCGCGGGAGACGTCATAAGAACTGTTCAGCAGTTCCGTTGCCTGTACATGGCCCGCCAGCAGCAGAGAAGCGACCAGCGCGAGTGAGTTCTTTTTCAGTAAGTTAACGGCCATTGCGCAcc                       >  6:16949/1‑149 (MQ=255)
                                                                                                                          aaCGGCGGATTCAGGGCGGCAAACAGCTCGCGGGAGACGTCATAAGAACTGTTCAGCAGTTCCGTTGCCTGTACATGGCCCGCCAGCAGCAGAGAAGCGACCAGCGCGAGTGAGTTCTTTTTCAGTAAGTTAACGGCCATTGCGCAccc                      <  1:231428/149‑1 (MQ=255)
                                                                                                                                ggATTCAGGGCGGCAAACAGCTCGCGGGAGACGTCATAAGAACTGTTCAGCAGTTCCGTTGCCTGTACATGGCCCGCCAGCAGCAGAGAAGCGCCCAGCGCGAGTGAGTTCTTTTTCAGTAAGTTAACGGCCATTGCGCACCCTtata                 >  2:71736/1‑148 (MQ=255)
                                                                                                                                ggATTCAGGGCGGCAAACAGCTCGCGGGAGACGTCATAAGAACTGTTCAGCAGTTCCGTTGCCTGTACATGGCCCGCCAGCAGCAGAGAAGCGACCa                                                                    <  7:98219/97‑1 (MQ=255)
                                                                                                                                ggATTCAGGGCGGCAAACAGCTCGCGGGAGACGTCATAAGAACTGTTCAGCAGTTCCGTTGCCTGTACATGGCCCGCCAGCAGCAGAGAAGCGACCa                                                                    >  8:98219/1‑97 (MQ=255)
                                                                                                                                  aTTCAGGGCCGCAAACAGCTCGCGGGAGACGTCATAAGAACTGTTCAGCAGTTCCGTTGCCTGTACATGGCCCGCCAGCAGCAGAGAAGCGACCAGCGCGAGTGAGTTCTTTTTCAGTAAGTTAACGGCCATTGCGCACCCTTATAAAt              <  5:13168/149‑1 (MQ=255)
                                                                                                                                      aGGGCGGCAAACAGCTCGCGGGAGACGTCATAAGAACTGTTCAGCAGTTACGTTGCCTGTACATGGCCCGCCAGCAGCAGAGAAGAGACCAGCGCGAGTGATTTCTTTTTCAGTAAGTTAACGGCCATTGCTCACACTTATAAATTTaa          >  6:145264/1‑149 (MQ=255)
                                                                                                                                       gggCGGCAAACAGCTCGCGGGAGACGTCATAAGAACTGTTCAGCAGTTCCGTTGCCTGTACATGGCCCGCCAGCAGCAGAGAAGCGACCAGCGCGAGTGAGTTCTTTTTCAGTAAGTTAACGGCCATTGCGCACCCTTATAAATTTAAt         <  5:55934/149‑1 (MQ=255)
                                                                                                                                       gggCGGCAAACAGCTCGCGGGAGACGTCATAAGAACTGTTCAGCAGTTCCGTTGCCTGTACATGGCCCGCCAGCAGCAGAGAAGCGACCAGCGCGAGTGAGTTCTTTTTCAGTAAGTTAACGGCCATTGCGCACCCTTATAAATTTAAt         <  2:103600/149‑1 (MQ=255)
                                                                                                                                              aaaCAGCTCGCGGGAGACGTCATAAGAACTGTTCAGCAGTTCCGTTGCCTGTACATGGCCCGCCAGCAGCAGAGAAGCGACCAGCGCGAGTGAGTTCTTTTTCAGTAAGTTAACGGCCATTGCGCACCCTTATAAATTTAATGACTTTc  <  8:38453/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                
TATAAGTGACAACGTCGGCTTTTAAGCCCTGTAAAATCGCCAGCGCCTGTTTTGATGACCCGGCATGAGATTGTTTTATCGTCAGTTTGTCGCCGCCGTTATCTTTTGCCCATTGTTGCTCAAACGGCGGATTCAGGGCGGCAAACAGCTCGCGGGAGACGTCATAAGAACTGTTCAGCAGTTCCGTTGCCTGTACATGGCCCGCCAGCAGCAGAGAAGCGACCAGCGCGAGTGAGTTCTTTTTCAGTAAGTTAACGGCCATTGCGCACCCTTATAAATTTAATGACTTTC  >  NC_000913/2543267‑2543557

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: ‑