Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,542,932 A→G 12.5% R199R (CGT→CGC cysP ← thiosulfate/sulfate ABC transporter periplasmic binding protein CysP

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,542,9320AG12.5% 183.2 / 15.1 72R199R (CGT→CGCcysPthiosulfate/sulfate ABC transporter periplasmic binding protein CysP
Reads supporting (aligned to +/- strand):  ref base A (40/23);  new base G (3/6);  total (43/29)
Fisher's exact test for biased strand distribution p-value = 1.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.65e-01

GCCACCGGGAATTCCGCCAGAATGTTGGTTTTCGGAATCACCACTTCAAAGCCCTGCGCTTCATACTGTTTACGGATGTTGTTCACTTCCGATTCGAAGCTAATCAGCACATCGCCCAGGCCGCGCTCGGCAAAAGTGGTGGTCGCGCCACGACCGCCAGTATCGAACACTTCAACGTTTTTCAGGAACTGGGTCATAAACTGTTCGGTTTTGCCTTTGTCACCACCGTCAGCTTTATCCGCTGCGCCCCATGCCGCCAGATAGGTATAACGCGCGTTACCCGACGTTTTCGG  >  NC_000913/2542783‑2543075
                                                                                                                                                     |                                                                                                                                               
gCCACCGGGAATTCCGCCAGAATGTTGGTTTTCGGAATCACCACTTCAAAGCCCTGCGCTTCATACTGTTTACGGATGTTGTTCACTTCCGATTCGAAGCTAATCAGCACATCGCCCAGGCCGCGCTCGGCAAAAGTGGTGGTCGCGCCa                                                                                                                                                 <  1:342243/150‑1 (MQ=255)
    ccGGGAATTCCGCCAGAATGTTGGTTTTCGGAATCACCACTTCAAAGCCCTGCGCTTCATACTGTTTACGGATGTTGTTCACTTCCGATTCGAAGCTAATCAGCACATCGCCCAGGCCGCGCTCGGCAAAAGTGGTGGTCGCGCCACGAc                                                                                                                                             <  2:40539/150‑1 (MQ=255)
          aTTCCGCCAGAATGTTGGTTTTCGGAATCACCACTTCAAAGCCCTGCGCTTCATACTGTTTACGGATGTTGTTCACTTCCGATTCGGAGCTAATCAGCACATCGCCCAGGCCGCGCTCGGCAAAAGTGGCGGTCGTGGCGCGAccaccac                                                                                                                                       <  2:1098558/150‑2 (MQ=255)
           ttCCGCCAGAATGTTGGTTTTCGGAATCACCACTTCAAAGCCCTGCGCTTCATACTGTTTACGGATGTTGTTCACTTCCGATTCGAAGCTAATCAGCACATCGCCCAGGCCGCGCTCGGCAAAAGTGGTGGTCGTGGCGCGAccaccact                                                                                                                                      <  1:859187/150‑3 (MQ=255)
           ttCCGCCAGAATGTTGGTTTTCGGAATCACCACTTCAAAGCCCTGCGCTTCATACTGTTTACGGATGTTGTTCACTTCCGATTCGAAGCTAATCAGCACATCGCCCAGGCCGCGCTCGGCAAAAGTGGTGGTCGCGCCACGACCGCCAGt                                                                                                                                      >  1:866217/1‑150 (MQ=255)
            tCCGCCAGAATGTTGGTTTTCGGAATCACCACTTCAAAGCCCTGCGCTTCATACTGTTTACGGATGTTGTTCACTTCCGATTCGAAGCTAATCAGCACATCGCCCAGGCCGCGCTCGGCAAAAGTGGCGGTCGTGGCGCGAccaccactt                                                                                                                                     <  2:471657/150‑4 (MQ=255)
             ccGCCATAATGTTGGTTTTCGGAATCACCACTTCAAAGCCCTGCGCTTCATACTGTTTACGGATGTTGTTCACTTCCGATTCGAAGCTAATCAGCACATCGCCCAGGCCGCGCTCGGCAAAAGTGGCGGTCGTGGCGCGAccaccacttt                                                                                                                                    <  2:1092499/150‑5 (MQ=255)
              cGCCAGAATGTTGGTTTTCGGAATCACCACTTCAAAGCCCTGCGCTTCATACTGTTTACGGATGTTGTTCACTTCCGATTCGAAGCTAATCAGCACATCGCCCAGGCCGCGCTCGGCAAAAGTGGTGGTCGCGCCACGACCGCCAGTATc                                                                                                                                   <  1:601505/150‑1 (MQ=255)
               gCCAGAATGTTGGTTTTCGGAATCACCACTTCAAAGCCCTGCGCTTCATACTGTTTACGGATGTTGTTCACTTCCGATTCGAAGCTAATCAGCACATCGCCCAGGCCGCGCTCGGCAAAAGTGGTGGTCGCGCCACGACCGCCAGTATCg                                                                                                                                  <  2:1090811/150‑1 (MQ=255)
                ccAGAATGTTGGTTTTCGGAATCACCACTTCAAAGCCCTGCGCTTCATACTGTTTACGGATGTTGTTCACTTCCGATTCGAAGCTAATCAGCACATCGCCCAGGCCGCGCTCGGCAAAAGTGGCGGTCGTGGCGCGAccaccacttttgc                                                                                                                                 <  2:1424137/150‑8 (MQ=255)
                  aGAATGTTGGTTTTCGGAATCACCACTTCAAAGCCCTGCGCTTCATACTGTTTACGGATGTTGTTCACTTCCGATTCGAAGCTAATCAGCACATCGCCCAGGCCGCGCTCGGCAAAAGTGGTGGTCGCGCCACGACCGCCAGTATCGAac                                                                                                                               >  1:514732/1‑150 (MQ=255)
                          ggTTTTCGGAATCACCACTTCAAAGCCCTGCGCTTCATACTGTTTACGGATGTTGTTCACTTCCGATTCGAAGCTAATCAGCACATCGCCCAGGCCGCGCTCGGCAAAAGTGGTGGTCGCGCCACGACCGCCAGTATCGAACACTTCAAc                                                                                                                       <  2:1296429/150‑1 (MQ=255)
                                cGGAATCACCACTTCAAAGCCCTGCGCTTCATACTGTTTACGGATGTTGTTCACTTCCGATTCGAAGCTAATCAGCACATCGCCCAGGCCGCGCTCGGCAAAAGTGGTGGTCGCGCCACGACCGCCAGTATCGAACACTTCAACGttttt                                                                                                                 >  2:186531/1‑150 (MQ=255)
                                  gAATCACCACTTCAAAGCCCTGCGCTTCATACTGTTTACGGATGTTGTTCACTTCCGATTCGAAGCTAATCAGCACATCGCCCAGGCCGCGCTCGGCAAAAGTGGTGGTCGCGCCACGACCGCCAGTATCGAACACTTCAACGTTTTTCa                                                                                                               <  1:351690/150‑1 (MQ=255)
                                    aTCACCACTTCAAAGCCCTGCGCTTCATACTGTTTACGGATGTTGTTCACTTCCGATTCGAAGCTAATCAGCACATCGCCCAGGCCGCGCTCGGCAAAAGTGGTGGTCGCGCCACGACCGCCAGTATCGAACACTTCAACGTTTTTCAgg                                                                                                             >  1:754741/1‑150 (MQ=255)
                                    aTCACCACTTCAAAGCCCTGCGCTTCATACTGTTTACGGATGTTGTTCACTTCCGATTCGAAGCTAATCAGCACATCGCCCAGGCCGCGCTCGGCAAAAGTGGTGGTCGCGCCACGACCGCCAGTATCGAACACTTCAACGTTTTTCAgg                                                                                                             <  2:588426/150‑1 (MQ=255)
                                      caccacTTCAAAGCCCTGCGCTTCATACTGTTTACGGATGTTGTTCACTTCCGATTCGAAGCTAATCAGCACATCGCCCAGGCCGCGCTCGGCAAAAGTGGTGGTCGCGCCACGACCGCCAGTATCGCACACTTCAACGTTTTTCAGGaa                                                                                                           >  1:1119396/1‑150 (MQ=255)
                                      caccacTTCAAAGCCCTGCGCTTCATACTGTTTACGGATGTTGTTCACTTCCGATTCGAAGCTAATCAGCACATCGCCCAGGCCGCGCTCGGCAAAAGTGGTGGTCGCGCCACGACCGCCAGTATCGAACACTTCAACGTTTTTCAGGaa                                                                                                           >  1:505650/1‑150 (MQ=255)
                                       accacTTCAAAGCCCTGCGCTTCATACTGTTTACGGATGTTGTTCACTTCCGATTCGAAGCTAATCAGCACATCGCCCAGGCCGCGCTCGGCAAAAGTGGTGGTCGCGCCACGACCGCCAGTATCGAACACTTCAACGTTTTTCAGGAAc                                                                                                          >  2:429103/1‑150 (MQ=255)
                                       accacTTCAAAGCCCTGCGCTTCATACTGTTTACGGATGTTGTTCACTTCCGATTCGAAGCTAATCAGCACATCGCCCAGGCCGCGCTCGGCAAAAGTGGTGGTCGCGCCACGACCGCCAGTATCGAACACTTCAACGTTTTTCAGGAAc                                                                                                          >  2:741855/1‑150 (MQ=255)
                                                     cTGCGCTTCATACTGTTTACGGATGTTGTTCACTTCCGATTCGAAGCTAATCAGCACATCGCCCAGGCCGCGCTCGGCAAAAGTGGTGGTCGCGCCACGACCGCCAGTATCGAACACTTCAACGTTTTTCAGGAACTGGGTCATAAACt                                                                                             <  1:761663/149‑1 (MQ=255)
                                                     cTGCGCTTCATACTGTTTACGGATGTTGTTCACTTCCGATTCGAAGCTAATCAGCACATCGCCCAGGCCGCGCTCGGCAAAAGTGGTGGTCGCGCCACGACCGCCAGTATCGAACACTTCAACGTTTTTCAGGAACTGGGTCATAAACt                                                                                             >  2:761663/1‑149 (MQ=255)
                                                          cTTCATACTGTTTACGGATGTTGTTCACTTCCGATTCGAAGCTAATCAGCACATCGCCCAGGCCGCGCTCGGCAAAAGTGGTGGTCGCGCCACGACCGCCAGTATCGAACACTTCAACGTTTTTCAGGAACTGGGTCATAAACTGTTCgg                                                                                       >  1:23240/1‑150 (MQ=255)
                                                           ttCATACTGTTTACGGATGTTGTTCACTTCCGATTCGAAGCTAATCAGCACATCGCCCAGGCCGCGCTCGGCAAAAGTGGTGGTCGCGCCACGACCGCCAGTATCGAACACTTCAACGTTTTTCAGGAACTGGGTCATAAACTGTTCGGt                                                                                      <  1:581338/150‑1 (MQ=255)
                                                               tACTGTTTACGGATGTTGTTCACTTCCGATTCGAAGCTAATCAGCACATCGCCCAGGCCGCGCTCGGCAAAAGTGGTGGTCGCGCCACGACCGCCAGTATCGAACACTTCAACGTTTTTCAGGAACTGGGTCATAAACTGTTCGGTTTTg                                                                                  <  2:514732/150‑1 (MQ=255)
                                                               tACTGTTTACGGATGTTGTTCACTTCCGATTCGAAGCTAATCAGCACATCGCCCAGGCCGCGCTCGGCAAAAGTGGTGGTCGCGCCACGACCGCCAGTATCGAACACTTCAACGTGTTTCAGGAACTGGGTCATAAACTGTTCGGTTTTg                                                                                  >  2:481601/1‑150 (MQ=255)
                                                                         ggATGTTGTTCACTTCCGATTCGAAGCTAATCAGCACATCGCCCAGGCCGCGCTCGTCAAAAGTGGTGGTCGCGCCACGACCGCCAGTATCGAACACTTCAACGTTTTTCAGGAACTGGGTCATAAACTGTTCGGTTTTGCCTTTGTcac                                                                        <  1:304751/150‑1 (MQ=255)
                                                                         ggATGTTGTTCACTTCCGATTCGAAGCTAATCAGCACATCGCCCAGGCCGCGCTCGGCAAAAGTGGTGGTCGCGCCACGACCGCCAGTATCGAACACTTCAACGTTTTTCAGGAACTGGGTCATAAACTGTTCGGTTTTGCCTTTGTcac                                                                        >  1:253207/1‑150 (MQ=255)
                                                                          gATGTTGTTCACTTCCGATTCGAAGCTAATCAGCACATCGCCCAGGCCGCGCTCGGCAAAAGTGGTGGTCGCGCCACGACCGCCAGTATCGAACACTTCAACGTTTTTCAGGAACTGGGTCATAAACTGTTCGGTTTTGCCTTTGTcacc                                                                       >  2:504303/1‑150 (MQ=255)
                                                                          gATGTTGTTCACTTCCGATTCGAAGCTAATCAGCACATCGCCCAGGCCGCGCTCGGCAAAAGTGGTGGTCGCGCCACGACCGCCAGTATCGAACACTTCAACGTTTTTCAGGAACTGGGTCATAAACTGTTCGGTTTTGCCTTTGTcacc                                                                       >  1:625883/1‑150 (MQ=255)
                                                                            tgttgtTCACTTCCGATTCGAAGCTAATCAGCACATCGCCCAGGCCGCGCTCGGCAAAAGTGGTGGTCGCGCCACGACCGCCAGTATCGAACACTTCAACGTTTTTCAGGAACTGGGTCATAAACTGTTCGGTTTTGCCTTTGTcaccac                                                                     <  2:112021/150‑1 (MQ=255)
                                                                                 gtCACTTCCGATTCGAAGCTAATCAGCACATCGCCCAGGCCGCGCTCGGCAAAAGTGGTGGTCGCGCCACGACCGCCAGTATCGAACACTTCAACGTTTTTCAGGAACTGGGTCATAAACTGTTCGGTTTTGCCTTTGTCACCACCGTCa                                                                <  2:320005/149‑1 (MQ=255)
                                                                                    aCTTCCGATTCGAAGCTAATCAGCACATCGCCCAGGCCGCGCTCGGCAAAAGTGGTGGTCGCGCCACGACCGCCAGTATCGAACACTTCAACGTTTTTCAGGAACTGGGTCATAAACTGTTCGGTTTTGCCTTTGTCACCACCGTCAGCt                                                             <  2:147283/150‑1 (MQ=255)
                                                                                     cTTCCGATTCGAAGCTAATCAGCACATCGCCCAGGCCGCGCTCGGCAAAAGTGGTGGTCGCGCCACGACCGCCAGTATCGAACACTTCAACGTTTTTCAGGAACTGGGTCATAAACTGTTCGGTTTTGCCTTTGTCACCACCGTCAGCtt                                                            >  1:1297590/1‑150 (MQ=255)
                                                                                       tCCGATTCGAAGCTAATCAGCACATCGCCCAGGCCGCGCTCGGCAAAAGTGGTGGTCGCGCCACGACCGCCAGTATCGAACACTTCAACGTTTTTCAGGAACTGGGTCATAAACTGTTCGGTTTTGCCTTTGTCACCACCGTCAGCTTTa                                                          >  2:839113/1‑150 (MQ=255)
                                                                                          gATTCGAAGCTAATCAGCACATCGCCCAGGCCGCGCTCGGCAAAAGTGGTGGTCGCGCCACGACCGCCAGTATCGAACACTTCAACGTTTTTCAGGAACTGGGTCATAAACTGTTCGGTTTTGCCTTTGTCACCACCGTCAGCTTTATcc                                                       >  1:1041936/1‑150 (MQ=255)
                                                                                          gATTCGAAGCTAATCAGCACATCGCCCAGGCCGCGCTCGGCAAAAGTGGTGGTCGCGCCACGACCGCCAGTATCGAACACTTCAACGTTTTTCAGGAACTGGGTCATAAACTGTTCGGTTTTGCCTTTGTCACCACCGTCAGCTTTATcc                                                       >  1:1036984/1‑150 (MQ=255)
                                                                                             tCGAAGCTAATCAGCACATCTCCCAGGCCGCGCTCGGCAAAAGTGGTGGTCGCGCCACGACCGCCAGTATCGAACACTTCAACGTTTTTCAGGAACTGGGTCATAAACTGTTCGGTTTTGCCTTTGTCACCACCGTCAGCTTTATCCGCt                                                    <  1:741855/150‑1 (MQ=255)
                                                                                               gAAGCTAATCAGCACATCGCCCAGGCCGCGCTCGGCAAAAGTGGTGGTCGCGCCACGACCGCCAGTATCGAACACTTCAACGTTTTTCAGGAACTGGGTCATAAACTGTTCGGTTTTGCCTTTGTCACCACCGTCAGCTTTATCCGCTgc                                                  >  1:1412589/1‑150 (MQ=255)
                                                                                                  gCTAATCAGCACATCGCCCAGGCCGCGCTCGGCAAAAGTGGTGGTCGCGCCACGACCGCCAGTATCGAACACTTCAACGTTTTTCAGGAACTGGGTCATAAACTGTTCGGTTTTGCCTTTGTCACCACCGTCAGCTTTATCCGCTGCGcc                                               >  2:737041/1‑150 (MQ=255)
                                                                                                    tAATCAGCACATCGCCCAGGCCGCGCTCGGCAAAAGTGGTGGTCGCGCCACGACCGCCAGTATCGAACACTTCAACGTTTTTCAGGAACTGGGTCATAAACTGTTCGGTTTTGCCTTTGTCACCACCGTCAGCTTTATCCGCTGCGcccc                                             >  1:352911/1‑150 (MQ=255)
                                                                                                       tCAGCACATCGCCCAGGCCGCGCTCGGCAAAAGTGGTGGTCGCGCCACGACCGCCAGTATCGAACACTTCAACGTTTTTCAGGAACTGGGTCATAAACTGTTCGGTTTTGCCTTTGTCACCACCGTCAGCTTTATCCGCTGCGCCCCATg                                          <  2:23240/150‑1 (MQ=255)
                                                                                                         aGCACATCGCCCAGGCCGCGCTCGGCAAAAGTGGTGGTCGCGCCACGACCGCCAGTATCGAACACTTCAACGTTTTTCAGGAACTGGGTCATAAACTGTTCGGTTTTGCCTTTGTCACCACCGTCAGCTTTATCCGCTGCGCCCCATgcc                                        >  2:410233/1‑150 (MQ=255)
                                                                                                          gCACATCGCCCAGGCCGCGCTCGGCAAAAGTGGTGGTCGCGCCACGACCGCCAGTATCGAACACTTCAACGTTTTTCAGGAACTGGGTCATAAACTGTTCGGTTTTGCCTTTGTCACCACCGTCAGCTTTATCCGCTGCGCCCCATgccg                                       >  2:26475/1‑150 (MQ=255)
                                                                                                           cacaTCGCCCAGGCCGCGCTCGGCAAAAGTGGTGGTCGCGCCACGACCGCCAGTATCGAACACTTCAACGTTTTTCAGGAACTGGGTCATAAACTGTTCGGTTTTGCCTTTGTCACCATCGTCAGCTTTATCCGCTGCGCCCCATgccgc                                      >  1:164335/1‑150 (MQ=255)
                                                                                                           cacaTCGCCCAGGCCGCGCTCGGCAAAAGTGGTGGTCGCGCCACGACCGCCAGTATCGAACACTTCAACGTTTTTCAGGAACTGGGTCATAAACTGTTCGGTTTTGCCTTTGTCACCACCGTCAGCTTTATCCGCTGCGCCCCATgccgc                                      >  2:1005335/1‑150 (MQ=255)
                                                                                                             caTCGCCCATGCCGCGATCGGCAAAAGTGGTGGTCGCGCCACGACCGCCAGTAGCGAACACTTCAACGTTTTTCAGGAACTGGGTCATAAACTGTTCGGTTTTTCGTTTGTCACCACCGTCAGCTTTATCCGCTGCGCCCCATGCCGCCa                                    <  2:1036528/150‑1 (MQ=255)
                                                                                                             caTCGCCCAGGCCGCGCTCGGCAAAAGTGGTGGTCGCGCAACGACCGCCAGTATCGAACACTTCAACGTTTTTCAGGAACTGGGTCATAAACTGTTCGGTTTTGCCTTTGTCACCACCGTCAGCTTTATCCGCTGCGCCCCATGCCGCCa                                    >  1:519412/1‑150 (MQ=255)
                                                                                                               tCGCCCAGGCCGCGCTCGGCAAAAGTGGTGGTCGCGCCACGACCGCCAGTATCGAACACTTCAACGTTTTTCAGGAACTGGGTCATAAACTGTTCGGTTTTGCCTTTGTCACCACCGTCAGCTTTATCCGCTGCGCCCCATGCCGCCAGa                                  <  1:265067/150‑1 (MQ=255)
                                                                                                                cGCCCAGGCCGCGCTCGGCAAAAGTGGTGGTCGCGCCACGACCGCCAGTATCGAACACTTCAACGTTTTTCAGGAACTGGGTCATAAACTGTTCGGTTTTGCCTTTGTCACCACCGTCAGCTTTATCCGCTGCGCCCCATGCCGCCAGAt                                 >  1:967341/1‑150 (MQ=255)
                                                                                                                 gCCCAGGCCGCGCTCGGCAAAAGTGGTGGTCGCGCCACGACCGCCAGTATCGAACACTTCAACGTTTTTCAGGAACTGGGTCATAAACTGTTCGGTTTTGCCTTTGTCACCACCGTCAGCTTTATCCGCTGCGCCCCATGCCGCCAGATa                                >  1:627991/1‑150 (MQ=255)
                                                                                                                   ccAGGCCGCGCTCGGCAAAAGTGGTGGTCGCGCCACGACCGCCAGTATCGAACACTTCAACGTTTTTCAGGAACTGGGTCATAAACTGTTCGGTTTTGCCTTTGTCACCACCGTCAGCTTTa                                                          <  1:933246/122‑1 (MQ=255)
                                                                                                                   ccAGGCCGCGCTCGGCAAAAGTGGTGGTCGCGCCACGACCGCCAGTATCGAACACTTCAACGTTTTTCAGGAACTGGGTCATAAACTGTTCGGTTTTGCCTTTGTCACCACCGTCAGCTTTa                                                          >  2:933246/1‑122 (MQ=255)
                                                                                                                    cAGGCCGCGCTCGGCAAAAGTGGTGGTCGCGCCACGACCGCCAGTATCGAACACTTCAACGTTTTTCAGGAACTGGGTCATAAACTGTTCGGTTTTGCCTTTGTCACCACCGTCAGCTTTATCCGCTGCGCCCCATGCCGCCAGATAGGt                             >  1:1136206/1‑150 (MQ=255)
                                                                                                                      ggCCGCGCTCGGCAAAAGTGGTGGTCGCGCCACGACCGCCAGTATCGAACACTTCAACGTTTTTCAGGAACTGGGTCATAAACTGTTCGGTTTTGCCTTTGTCACCACCGTCAGCTTTATCCGCTGCGCCCCATGCCGCCAGATAGGtat                           >  2:683794/1‑150 (MQ=255)
                                                                                                                         cgcgCTCGGCAAAAGTGGTGGTCGCGCCACGACCGCCAGTATCGAACACTTCAACGTTTTTCAGGAACTGGGTCATAAACTGTTTGGTTTTGCCTTTGTCACCACCGTCAGCTTTATCCGCTGCGCCCCATGCCGCCAGATAGGTATAAc                        <  2:754741/150‑1 (MQ=255)
                                                                                                                            gcTCGGCAAAAGTGGTGGTCGCGCCACGACCGCCAGTATCGAACACTTCAACGTTTTTCAGGAACTGGGTCATAAACTGTTCGGTTTTGCCTTTGTCACCACCGTCAGCTTTATCCGCTGCGCCCCATGCCGCCAGATAGGTATAAcgcg                     >  1:999100/1‑150 (MQ=255)
                                                                                                                              tCGGCAAAAGTGGTGGTCGCGCCACGACCGCCAGTATCGAACACTTCAACGTTTTTCAGGAACTGGGTCATAAACTGTTCGGTTTTGCCTTTGTCACCACCGTCAGCTTTATCCGCTGCGCCCCATGCCGCCAGATAGGTATAAcgcgcg                   <  2:352911/150‑1 (MQ=255)
                                                                                                                                ggCAAAAGTGGTGGTCGCGCCACGACCGCCAGTATCGAACACTTCAACGTTTTTCAGGAACTGGGTCATAAACTGTTCGGTTTTGCCTTTGTCACCACCGTCAGCTTTATCCGCTGCGCCCCATGCCGCCAGATAGGTATAACGCGCGtt                 >  2:925397/1‑150 (MQ=255)
                                                                                                                                ggCAAAAGTGGTGGTCGCGCCACGACCGCCAGTATCGAACACTTCAACGTTTTTCAGGAACTGGGTCATAAACTGTTCGGTTTTGCCTTTGTCACCACCGTCAGCTTTATCCGCTGCGCCCCATGCCGCCAGATAGGTATAACGCGCGtt                 >  2:986913/1‑150 (MQ=255)
                                                                                                                                    aaaGTGGTGGTCGCGCCACGACCGCCAGTATCGAACACTTCAACGTTTTTCAGGAACTGGGTCATAAACTGTTCGGTTTTGCCTTTGTCACCACCGTCAGCTTTATCCGCTGCGCCCCATGCCGCCAGATAGGTATAACGCGCGTTAccc             <  2:1136206/150‑1 (MQ=255)
                                                                                                                                     aaGTGGTGGTCGCGCCACGACCGCCAGTATCGAACACTTCAACGTTTTTCAGGAACTGGGTCATAAACTGTTCGGTTTTGCCTTTGTCACCACCGTCAGCTTTATCCGCTGCGCCCCATGCCGCCAGATAGGTATAACGCGCGTTACCCg            >  1:1249930/1‑150 (MQ=255)
                                                                                                                                     aaGTGGTGGTCGCGCCACGACCGCCAGTATCGAACACTTCAACGTTTTTCAGGAACTGGGTCATAAACTGTTCGGTTTTGCCTTTGTCACCACCGTCAGCTTTATCCGCTGCGCCCCATGCCGCCAGATAGGTATAACGCGCGTTACCCg            >  1:1245169/1‑150 (MQ=255)
                                                                                                                                       ctggcggTCGTGGCGCGACCGCCAGTATCGAACACTTCAACGTTTTTCAGGAACTGGGTCATAAACTGTTCGGTTTTGCCTTTGTCACCACCGTCAGCTTTATCCGCTGCGCCCCATGCCGCCAGATAGGTATAACGCGCGTTACCCGAc          >  2:95282/2‑150 (MQ=255)
                                                                                                                                         ggtggtCGCGCCACGACCGCCAGTATCGAACACTTCAACGTTTTTCAGGAACTGGGTCATAAACTGTTCGGTTTTGCCTTTGTCACCACCGTCAGCTTTATCCGCTGCGCCCCATGCCGCCAGATAGGTATAACGCGCGTTACCCGACGt        >  1:1325854/1‑150 (MQ=255)
                                                                                                                                         ggcggtCGTGGCGCGACCGCCAGTATCGAACACTTCAACGTTTTTCAGGAACTGGGTCATAAACTGTTCGGTTTTGCCTTTGTCACCACCGTCAGCTTTATCCGCTGCGCCCCATGCCGCCAGATAGGTATAACGCGCGTTACCCGACGt        >  1:373760/4‑150 (MQ=255)
                                                                                                                                         ggcggtCGTGGCGCGACCACCAGTATCGAACACTTCAACGTTTTTCAGGAACTGGGTCATAAACTGTTCGGTTTTGCCTTTGTCACCACCGTCTGCTTTATCCGCTGCGCCCCATGCCGCCAGATAGGtat                           <  2:499085/128‑1 (MQ=255)
                                                                                                                                         ggcggtCGTGGCGCGACCACCAGTATCGAACACTTCAACGTTTTTCAGGAACTGGGTCATAAACTGTTCGGTTTTGCCTTTGTCACCACCGTCTGCTTTATCCGCTGCGCCCCATGCCGCCAGATAGGtat                           >  1:499085/4‑131 (MQ=255)
                                                                                                                                             gtCGCGCCACGACCGCCAGTATCGAACACTTCAACGTTTTTCAGGAACTGGGTCATAAACTGTTCGGTTTTGCCTTTGTCACCACCGTCAGCTTTATCCGCTGCGCCCCATGCCGCCAGATAGGTATAACGCGCGTTACCCGACGTTttc    <  1:737041/150‑1 (MQ=255)
                                                                                                                                              tCGCGCCACGACCGCCAGTATCGAACACTTCAACGTTTTTCAGGAACTGGGTCATAAACTGTTCGGTTTTGCCTTTGTCACCACCGTCAGCTTTATCCGCTGCGCCCCATGCCGCCAGATGGGTATAACGCGCGTTACCCGACGTTttcg   >  1:1330859/1‑150 (MQ=255)
                                                                                                                                               cgcgCCACGACCGCCAGTATCGAACACTTCAACGTTTTTCAGGAACTGGGTCATAAACTGTTCGGTTTTGCCTTTGTCACCACCGTCAGCTTTATCCGCTGCGCCCCATGCCGCCAGATAGGTATAACGCGCGTTACCCGACGTTttcgg  >  2:105874/1‑150 (MQ=255)
                                                                                                                                               cgcaCCACGACCGCCAGTATCGAACACTTCAACGTTTTTCAGGAACTGGGTCATAAACTGTTCGGTTTTGCCTTTGTCACCACCGTCAGCTTTATCCGCTGCGCCCCATGCCGCCAGATAGGta                            >  2:346800/1‑124 (MQ=255)
                                                                                                                                               cgcaCCACGACCGCCAGTATCGAACACTTCAACGTTTTTCAGGAACTGGGTCATAAACTGTTCGGTTTTGCCTTTGTCACCACCGTCAGCTTTATCCGCTGCGCCCCATGCCGCCAGATAGGta                            <  1:346800/124‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                               
GCCACCGGGAATTCCGCCAGAATGTTGGTTTTCGGAATCACCACTTCAAAGCCCTGCGCTTCATACTGTTTACGGATGTTGTTCACTTCCGATTCGAAGCTAATCAGCACATCGCCCAGGCCGCGCTCGGCAAAAGTGGTGGTCGCGCCACGACCGCCAGTATCGAACACTTCAACGTTTTTCAGGAACTGGGTCATAAACTGTTCGGTTTTGCCTTTGTCACCACCGTCAGCTTTATCCGCTGCGCCCCATGCCGCCAGATAGGTATAACGCGCGTTACCCGACGTTTTCGG  >  NC_000913/2542783‑2543075

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