Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,542,932 A→G 7.9% R199R (CGT→CGC cysP ← thiosulfate‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,542,9320AG 7.9% 316.1 / 14.2 114R199R (CGT→CGCcysPthiosulfate‑binding protein
Reads supporting (aligned to +/- strand):  ref base A (60/44);  new base G (5/4);  total (66/48)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.72e-01

Only 100 of 116 total aligned reads displayed.

CCACCGGGAATTCCGCCAGAATGTTGGTTTTCGGAATCACCACTTCAAAGCCCTGCGCTTCATACTGTTTACGGATGTTGTTCACTTCCGATTCGAAGCTAATCAGCACATCGCCCAGGCCGCGCTCGGCAAAAGTGGTGGTCGCGCCACGACCGCCAGTATCGAACACTTCAACGTTTTTCAGGAACTGGGTCATAAACTGTTCGGTTTTGCCTTTGTCACCACCGTCAGCTTTATCCGCTGCGCCCCATGCCGCCAGATAGGTATAACGCGCGTTACCCGACGTTTTCGGGTTCGG  >  NC_000913/2542784‑2543081
                                                                                                                                                    |                                                                                                                                                     
ccACCGGGAATTCCGCCAGAATGTTGGTTTTCGGAATCACCACTTCAAAGCCCTGCGCTTCATACTGTTTACGGATGTTGTTCACTTCCGATTCGAAGCTAATCAGCACATCGCCCAGGCCGCGCTCGGCAAAAGTGGTGGTCGCGCCAc                                                                                                                                                      >  2:355368/1‑150 (MQ=255)
   ccGGGAATTCCGCCAGAATGTTGGTTTTCGGAATCACCACTTCAAAGCCCTGCGCTTCATACTGTTTACGGATGTTGTTCACTTCCGATTCGAAGCTAATCAGCACATCGCCCAGGCCGCGCTCGGCAAAAGTGGTGGTCGCGCCACGAc                                                                                                                                                   >  1:255591/1‑150 (MQ=255)
       gAATTCCGCCAGAATGTTGGTTTTCGGAATCACCACTTCAAAGCCCTGCGCTTCATACTGTTTACGGATGTTGTTCACTTCCGATTCGAAGCTAATCAGCACATCGCCCAGGCCGCGCTCGGCAAAAGTGGTGGTCGCGCCACGACCGcc                                                                                                                                               <  1:773244/150‑1 (MQ=255)
       gAATTCCGCCAGAATGTTGGTGTTCGGAATCACCACTTCAAAGCCCTGCGCTTCATACTGTATACGGATGTGGTGCACTTCCGATTCGAAGCTAATCAGCACATCGCCCAGGCCGCGCTCGGCAAAAGTGGTGGTCGCGCCACGACCGcc                                                                                                                                               <  1:1011340/150‑1 (MQ=255)
         aTTCCGCCAGAATGTTGGTTTTCGGAATCACCACTTCAAAGCCCTGCGCTTCATACTGTTTACGGATGTTGTTCACTTCCGATTCGAAGCTAATCAGCACATCGCCCAGGCCGCGCTCGGCAAAAGTGGCGGTCGTGGCGCGAccaccac                                                                                                                                             <  1:686488/150‑2 (MQ=255)
          atCCGCCAGAATGTTGGATTTCGGAATCACCACTTCAAAGCCCTGCGCTTCATACTGTTTACGGATGTTGTTCACTTCCGATTCGAAGCTAATCAGCACATCGCCCAGGCCGCGCTCGGCAAAAGTGGCGGTCGTGGCGCGAccaccact                                                                                                                                            <  1:30670/149‑3 (MQ=255)
           tCCGCCAGAATGTTGGTTTTCGGAATCACCACTTCAAAGCCCTGCGCTTCATACTGTTTACGGATGTTGTTCACTTCCGATTCGAAGCTAATCAGCACATCGCCCAGGCCGCGCTCGGCAAAAGTGGTGGTCGCGCCACGACCGCCAGTa                                                                                                                                           >  1:723441/1‑150 (MQ=255)
              gCCAGAATGTTGGTTTTCGGAATCACCACTTCAAAGCCCTGCGCTTCATACTGTTTACGGATGTTGTTCACTTCCGATTCGAAGCTAATCAGCACATCGCCCAGGCCGCGCTCGGCAAAAGTGGTGGTCGCGCCACGACCGCCAGTATCg                                                                                                                                        >  2:183991/1‑150 (MQ=255)
              gCCAGAATGTTGGTTTTCGGAATCACCACTTCAAAGCCCTGCGCTTCATACTGTTTACGGATGTTGTTCACTTCCGATTCGAAGCTAATCAGCACATCGCCCAGGCCGCGCTCGGCAAAAGTGGTGGTCGCGCCACGACCGCCAGTATCg                                                                                                                                        >  2:1477193/1‑150 (MQ=255)
              gCCAGAATGTTGGTTTTCGGAATCACCACTTCAAAGCCCTGCGCTTCATACTGTTTACGGATGTTGTTCACTTCCGATTCGAAGCTAATCAGCACATCGCCCAGGCCGCGCTCGGCAAAAGTGGTGGTCGCGCCACGACCGCCAGTATCg                                                                                                                                        >  1:1681837/1‑150 (MQ=255)
                 aGAATGATGGTTTTCGGAATCACCACTTCAAAGCCCTGCGCTTCATACTGTTTACGGATGTTGTTCACTTCCGATTCGAAGCTAATCAGCACATCGCCCAGGCCGCGCTCGGCAAAAGTGGTGGTCGCGCCACGACCGCCAGTATCGAac                                                                                                                                     >  1:1404914/1‑150 (MQ=255)
                  gAATGTTGGGATTCGGAATCACCACTTCCAAGCCCTGCGCTTCATACTGTTTACGGATGTTGTTCACTTCCGATTCGAAGCTAATCAGCACATCGCCCAGGCCGCGCTCGGCAAAAGTGGTGGTCGCGCCACGACCGCCAGTATCGAaca                                                                                                                                    <  1:1540282/150‑1 (MQ=255)
                     tGTTGGTTTTCGGAATCACCACTTCAAAGCCCTGCGCTTCATACTGTTTACGGATGTTGTTCACTTCCGATTCGAAGCTAATCAGCACATCGCCCAGGCCGCGCTCGGCAAAAGTGGTGGTCGCGCCACGACCGCCAGTATCGAACACtt                                                                                                                                 >  1:1300326/1‑150 (MQ=255)
                      gTTGGTTTTCGGAATCACCACTTCAAAGCCCTGCGCTTCATACTGTTTACGGATGTTGTTCACTTCCGATTCGAAGCTAATCAGCACATCGCCCAGGCCGCGCTCGGCAAAAGTGGTGGTCGCGCCACGACCGCCAGTATCGAACACTTc                                                                                                                                <  2:1718122/150‑1 (MQ=255)
                            tttCGGAATCACCACTTCAAAGCCCTGCGCTTCATACTGTTTACGGATGTTGTTCACTTCCGATTCGCAGCTAATCAGCACATCGCCCAGGCCGCGATCGGCAAAAGTGGTGGTCGCGCCACGACCGCCAGTATCGAACACGTCAACGtt                                                                                                                          >  2:1016105/1‑150 (MQ=255)
                              tCGGAATCACCACTTCAAAGCCCTGCGCTTCATACTGTTTACGGTTGTTGTTCACTTCCGATTCGAAGCTAATCAGCACATCGCCCAGGCCGCGCTCGGCAAAAGTGGTGGTCGCGCCACGACCGCCAGTATCGAACACTTCAACGtttt                                                                                                                        <  1:767693/150‑1 (MQ=255)
                               cGGAATCACCACTTCAAAGCCCTGCGCTTCATACTGTTTACGGATGTTGTTCACTTCCGATTCGAAGCTAATCAGCACATCGCCCAGGCCGCGCTCGGCAAAAGTGGTGGTCGCGCCACGACCGCCAGTATCGAACACTTCAACGttttt                                                                                                                       >  2:1433623/1‑150 (MQ=255)
                                ggAATCACCACTTCAAAGCCCTGCGCTTCTTACTGTTTACGGATGTTGTTCACTTCCGATTCGAACCTAATCAGCACATCGCCCAGGCCGCGCTCGGCAAAAGTGGTGGTCGCGCCACGACCGCCAGTATCGAACACTTCAACGTTTTTc                                                                                                                      <  2:238894/150‑1 (MQ=255)
                                ggAATCACCACTTCAAAGCCCTGCGCTTCATACTGTTTACGGATGTTGTTCACTTCCGATTCGAAGCTAATCAGCACATCGCCCAGGCCGCGCTCGGCAAAAGTGGTGGTCGCGCCACGACCGCCAGTATCGAACACTTCAACGTTTTTc                                                                                                                      >  2:1024235/1‑150 (MQ=255)
                                 gAATCACCACTTCAAAGCCCTGCGCTTCATACTGTTTACGGATGTTGTTCACTTCCGATTCGAAGCTAATCAGCACATCGCCCAGGCCGCGCTCGGCAAAAGTGGTGGTCGCGCCACGACCGCCAGTATCGAACACTTCAACGTTTTTCa                                                                                                                     <  1:867481/150‑1 (MQ=255)
                                 gAATCACCACTTCAAAGCCCTGCGCTTCATACTGTTTACGGATGTTGTTCACTTCCGATTCGAAGCTAATCAGCACATCGCCCAGGCCGCGCTCGGCAAAAGTGGTGGTCGCGCCACGACCGCCAGTATCGAACACTTCAACGTTTTTCa                                                                                                                     >  1:878953/1‑150 (MQ=255)
                                 gAATCACCACTTCAAAGCCCTGCGCTTCATACTGTTTACGGATGTTGTTCACTTCCGATTCGAAGCTAATCAGCACATCGCCCAGGCCGCGCTCGGCAAAAGTGGTGGTCGCGCCACGACCGCCAGTATCGAACACTTCAACGTTTTTCa                                                                                                                     <  1:50286/150‑1 (MQ=255)
                                    tCACCACTTCAAAGCCCTGCGCTTCATACTGTTTACGGATGTTGTTCACTTCCGATTCGAAGCTAATCAGCACATCGCCCAGGCCGCGCTCGGCAAAAGTGGTGGTCGCGCCACGACCGCCAGTATCGAACACTTCAACGTTTTTCAGGa                                                                                                                  <  2:558206/150‑1 (MQ=255)
                                      accacTTCAAAGCCCTGCGCTTCATACTGTTTACGGATGTTGTTCACTTCCGATTCGAAGCTAATCAGCACATCGCCCAGGCCGCGCTCGGCAAAAGTGGTGGTCGCGCCACGACCGCCAGTATCGAACACTTCAACGTTTTTCAGGAAc                                                                                                                >  1:1450271/1‑150 (MQ=255)
                                      accacTTCAAAGCCCTGCGCTTCATACTGTTTACGGATGTTGTTCACTTCCGATTCGAAGCTAATCAGCACATCGCCCAGGCCGCGCTCGGCAAAAGTGGTGGTCGCGCCACGACCGCCAGTATCGAACACTTCAACGTTTTTCAGGAAc                                                                                                                >  1:1379925/1‑150 (MQ=255)
                                      accacTTCAAAGCCCTGCGCTTCATACTGTTTACGGATGTTGTTCACTTCCGATTCGAAGCAAATCAGCACATCGCCCAGGCCGCGCTCGGCAAAAGTGGTGGTCGCGCCACGACCGCCAGTATCGAACACTTCAACGTTTTTCAGGAAc                                                                                                                >  1:229609/1‑150 (MQ=255)
                                       ccacTTCAAAGCCCTGAGCTTCATACTGTTTACGGATGTTGTTCACTTCCGATTCGAAGCTAATCAGAACATCGCCCAGGCCGCGCTCGGCAAAAGTGGTGGTCGCGCCAAGACCGCCAGTATAGAACACTTCAACGTTTTTCAGGAACt                                                                                                               >  2:804240/1‑150 (MQ=255)
                                        cacTTCAAAGCCCTGCGCTTCATACTGTTTACGGATGTTGTTCACTTCCGATTCGAAGCTAATCAGCACATCGCCCAGGCCGCGCTCGGCAAAAGTGGTGGTCGCGCCACGCCCGCCAGTATCGAACACTTCAACGTTTTTCAGGAACTg                                                                                                              >  2:1035743/1‑150 (MQ=255)
                                         acTTCAAAGCCCTGCGCTTCATACTGTTTACGGATGTTGTTCACTTCCGATTCGAAGCTAATCAGCACATCTCCCAGGACGCGCTCGGCAAAAGTGGTGGTCGCGACACGACCGCCAGTATCGAACACTTCAACGGTTTTAAGGAACTgg                                                                                                             >  2:1556931/1‑150 (MQ=255)
                                          cTTCAAAGCCCTGCGCTTCATACTGTTTACGGTTGTTGTTCACTTCCGATTCGAAGCTAATCAGCACATCTCCCAGGCCGCGCTCGGCATAAGTTGTGGTCGCGGCACGACCGCCAGTATCGAACACTTCAACGTTTTTCAGGATCTggg                                                                                                            <  1:1197547/150‑1 (MQ=255)
                                          cTTCAAAGCCCTGCGCTTCATACTGTTTACGGATGTTGTTCACTTCCGATTCGAAGCTAATCAGCACATCGCCCAGGCCGCGCTCGGCAAAAGTGGTGGTCGCGCCACGACCGCCAGTATCGAACACTTCAACGTTTTTCAGGAACTggg                                                                                                            <  1:1488660/150‑1 (MQ=255)
                                          cTTCAAAGCCCTGCGCTTCATACTGTTTACGGATGTTGTTCACTTCCGATTCGAAGCTAATCAGCACATCGCCCAGGCCGCGCTCGGCAAAAGTGGTGGTCGCGCCACGACCGCCAGTATCGAACACTTCAACGTTTTTCAGGAACTggg                                                                                                            >  2:1447944/1‑150 (MQ=255)
                                            tCAAAGCCCTGCTCTTCATACTGTTTACGGATGTTGTTCACTTCCGATTCGAAGCTAATCAGCACATCGCCCAGGCCGCTCTCTGCAAAAGTTGTGGTCGCGCCACGACCGCCAGTATCGAACACTTCAACGTTTTTCAGGAACTGGGTc                                                                                                          <  2:312858/150‑1 (MQ=255)
                                             cAAAGCCCTGCGCTTCATACTGTTTACGGATGTTGTTCACTTCCGATTCGAAGCTAATCAGCACATCGCCCAGGCCGCGCTCGGCAAAAGTGGTGGTCGCGCCACGCCCGCCAGTATCGAACACTTCAACGTTTTTCAGGAACTGGGTCa                                                                                                         >  2:833591/1‑150 (MQ=255)
                                             cAAAGCCCTGCGCTTCATACTGTTTACGGATGTTGTTCACTTCCGATTCGAAGCTAATCAGCACATCGCCCAGGCCGCGCTCGGCAAAAGTGGTGGTCGCGCCACGACCGCCAGTATCGAACACTTCAACGTTTTTCAGGAACTGGGTCa                                                                                                         >  1:1584204/1‑150 (MQ=255)
                                              aaaGCCCTGCGCTTCATACTGTTTACGGATGTTGTTCACTTCCGATTCGAAGCTAATCAGCACATCGCCCAGGCCGCGCTCGGCAAAAGTGGTGGTCGCGCCACGACCGCCAGTATCGAACACTTCAACGTTTTTCAGGAACTGGGTCAt                                                                                                        >  2:1258362/1‑150 (MQ=255)
                                                aGCCCTGCGCTTCATACTGTTTACGGATGTTGTTCACTTCCGATTCGAAGCTAATCAGCACATCGCCCAGGCCGCGCTCGGCAAAAGTGGTGGTCGCGCCACGACCGCCAGTATCGAACACTTCAACGTTTTTCAGGAACTGGGTCATaa                                                                                                      <  1:164408/150‑1 (MQ=255)
                                                  cccTGCGCTTCATACTGTTTACGGATGTTGTTCACTTCCGATTCGAAGCTAATCAGCACATCGCCCAGGCCGCGCTCGGCAAAAGTGGTGGTCGCGCCACGACCGCCAGTATCGAACACTTCAACGTTTTTCAGGAACTGGGTCATAAAc                                                                                                    >  2:759703/1‑150 (MQ=255)
                                                        gcTTCATACTGTTTACGGATGTTGTTCACTTCCGATTCGAAGCTAATCAGCACATCGCCCAGGCCGCGCTCGGCAAAAGTGGTGGTCGCGCCACGACCGCCAGTATCGAACACTTCAACGTTTTTCAGGAACTGGGTCATAAACTGTTCg                                                                                              <  2:1713039/150‑1 (MQ=255)
                                                         cTTCATACTGTTTACGGATGTTGTTCACTTCCGATTCTAAGCTAATCAGCACATCGCCCAGGCAGCGCTCGGCAAAAGTGGTGGTCGCGCCCCGACCGCCAGTATCGAACACTTCACCGTTTTTCAGGAACTGGGTAATAACCTGTTCgg                                                                                             >  1:114499/1‑150 (MQ=255)
                                                         cTTCATACTGTTTACGGATGTTGTTCACTTCCGATTCGAAGCTAATCAGCACATCGCCCAGGCCGCGCTCGGCAAAAGTGGTGGTCGCGCCACGACCGCCAGTATCGAACACTTCAACGTTTTTCAGGAACTGGGTCATAAACTGTTCgg                                                                                             >  1:1423281/1‑150 (MQ=255)
                                                             aTACTGTTTACGGATGTTGTTCACTTCCGATTCGAAGCTAATCAGCACATCGCCCAGGCCGCGCTCGGCAAAAGTGGTGGTCGCGCCGCGACCGCCAGTATCGAACACTTCAACGTTTTTCAGGAACTGGGTCATAAACTGTTCGGtttt                                                                                         >  1:510058/1‑150 (MQ=255)
                                                             aTACTGTTTACGGATGTTGTTCACTTCCGATTCGAAGCTAATCAGCACATCGCCCAGGCCGCGCTCGGCAAAAGTGGTGGTCGCGCCACGACCGCCAGTATCGAACACTTCAACGTTTTTCAGGAACTGGGTCATAAACTGTTCGGtttt                                                                                         >  2:1072521/1‑150 (MQ=255)
                                                              tACTGTTTACGGATGTTGTTCACTTCCGATTCGAAGCTAATCAGCACATCGCCCAGGCCGCGCTCGGCAAAAGTGGTGGTCGCGCCACGACCGCCAGTATCGAACACTTCAACGTTTTTCAGGAACTGGGTCATAAACTGTTCGGTTTTg                                                                                        <  1:719866/150‑1 (MQ=255)
                                                                  gTTTCCGGATGTTGTTCACTTCCGATTCGAAGCTAATCTGCACATCGCCCAGGCCGCGCTCGGCAAAAGTGGTGGTCGCGCCACGACCGCCAGTATCGAACACTTCAACGTTTTTCAGGAACTGGGTCATAAACTGTTCGGTTTTGCCtt                                                                                    <  1:227091/150‑1 (MQ=255)
                                                                  gTTTACGGATGTTGTTCACTTCCGATTCGAAGCTAATCAGCACATCGCCCAGGCCGCGCTCGGCAAAAGTGGTGGTCGCGCCACGACCGCCAGTATCGAACACTTCAACGTTTTTCAGGAACTGGGTCATAAACTGTTCGGTTTTGCCtt                                                                                    >  2:1160324/1‑150 (MQ=255)
                                                                  gTTTACGGATGGTGTTCAATTCCGATTCGAAGCTAATCAGCACATCGCCCAGGCCGCGCTCGGCAAAAGTGGTGTTCGCGCCACGACCGCCAGTATCGAACACTGCAACGTTTTTCAGGAACTGGGTCATAAACTGTTCGGTTTTGCCtt                                                                                    <  1:222012/150‑1 (MQ=255)
                                                                   tttACGTATGTTGTTCACTTCCGATTCGAAGCTAATCAGCACATCGCCCAGGCCGCGCTCGGCAAAAGTGGTGGTCGCGCCACGACCGCCAGTATCGAACACTTCAACGTTTTTCAGGAACTGGGTCATAAACTGTTCGGTTTTGCCttt                                                                                   <  1:181819/150‑1 (MQ=255)
                                                                        ggATGTTGTTCACTTCCGATTCGAAGCTAATCAGCACATCGCCCAGGCCGCGCTCGGCAAAAGTGGTGGTCGCGCCACGACCGCCAGTATCGAACACTTCAACGTTTTTCAGGAACTGGGTCATAAACTGTTCGGTTTTGCCTTTGTcac                                                                              >  1:394227/1‑150 (MQ=255)
                                                                         gATGTTGTTCACTTCCGATTCGAAGCTAATCAGCACATCGCCCAGGCCGCGCTCGGCAAAAGTGGTGGTCGCGCCACGACCGCCAGTATCGAACACTTCAACGTTTTTCAGGAACTGGGTCATAAACTGTTCGGTTTTGCCTTTGTcacc                                                                             <  2:1300326/150‑1 (MQ=255)
                                                                          aTGTTGTTCACTTCCGATTCGAAGCTAATCAGCACATCGCCCAGGCCGCGCTCGGCAAAAGTGGTGGTCGCGCCACGACCGCCAGTATCGAACACTTCAACGTTTTTCAGGAACTGGGGCATAAACTGTTCGGTTTTGCCTTTGTcacca                                                                            <  1:965544/150‑1 (MQ=255)
                                                                           tgttgtTCACTTCCGATTCGAAGCTAATCAGCACATCGCCCAGGCCGCGCTCGGCAAAAGTGGTGGTCGCGCCACGACCGCCAGTATCGAACACTTCAACGTTTTTCAGGAACTGGGTCATAAACTGTTCGGTTTTGCCTTTGTcaccac                                                                           <  1:543042/150‑1 (MQ=255)
                                                                                  cACTTCCGATTCGAAGCTAATCAGCACATCGCCCAGGCCGCGCTCGGCAAAAGTGGTGGTCGCGCCACGACCGCCAGTATCGAACACTTCAACGTTTTTCAGGAACTGGGTCATAAACTGTTCGGTTTTGCCTTTGTCACCACCGTCAGc                                                                    <  1:1038257/150‑1 (MQ=255)
                                                                                     ttCCGATTCGAAGCTAATCAGCACATCGCCCAGGCCGCGCTCGGCAAAAGTGGTGGTCGCGCCACGACCGCCAGTATCGAACACTTCAACGTTTTTCAGGAACTGGGTCATAAACTGTTCGGTTTTGCCTTTGTCACCACCGTCAGCttt                                                                 <  2:1450271/150‑1 (MQ=255)
                                                                                     ttCCGATTCGAAGCTAATCAGCACATCGCCCAGGACGCGCTCGGCAAAAGTGGTGGTCGCGCCACGACCGCCAGTATCGAACACTTCAACGTTTTTCAGGAACTGGGTCATAAACTGTTCGGTTTTGCCTTTGTCTCCACCGTCAGCttt                                                                 <  2:1050718/150‑1 (MQ=255)
                                                                                       ccGATTCGAAGCTAATCAGCACATCGCCCAGGCCGCGCTCGGCAAAAGTGGTGGTCGCGCCACGACCGCCAGTATCGAACACTTCAACGTTTTTCAGGAACTGGGTCATAAACTGTTCGGTTTTGCCTTTGTCACCACCGTCAGCTTTAt                                                               >  2:540037/1‑150 (MQ=255)
                                                                                        cGATTCGAAGCTAATCAGCACATCGCCCAGGCCGCGCTCGGCAAAAGTGGTGGTCGCGCCACGCCCGCCAGTATCGAACACTTCAACGTTTTTCAGGAACTGGGTCATAGACTGTTCGGTGTTGCCTTTGTCACCACCGTCAGCTTTATc                                                              >  2:919065/1‑150 (MQ=255)
                                                                                         gATTCGAAGCTAATCAGCACATCGCCCAGGCCGCGCTCGGCAAAAGTGGTGGTCGCGCCGCGACCGCCAGTATCGAACACTTCAACGTTTTTCAGGAACTGGGTCATAAACTGTTCGGTTTTGCCTTTGTCACCACCGTCAGCTTTATcc                                                             <  2:510058/150‑1 (MQ=255)
                                                                                          aTTCGAAGCTAATCAGCACATCGCCCAGGCCGCGCTCGGCAAAAGTGGTGGTCGCGCCACGACCGCCAGTATCGAACACTTCAACGTTTTTCAGGAACTGGGTCATAAACTGTTCGGTTTTGCCTTTGTCACCACCGTCAGCTTTATCCg                                                            >  1:1084360/1‑150 (MQ=255)
                                                                                             cGAAGCTAATCAGCACATCGCCCAGGCCGCGCTCGGCAAAAGTGGTGGTCGCGCCACGACCGCCAGTATCGAACACTTCAACGTTTTTCAGGCACTGGGTCATAAACTGTTCGGTTTTGCCTTTGTCACCACCGTCAGCTTTATCCGCTg                                                         >  1:698233/1‑150 (MQ=255)
                                                                                              gAAGCTAATCAGCACATCGCCCAGGCCGCGCTCGGCAAAAGTGGTGGTCGCGCCACGACCGCCAGTATCGAACACTTCAACGTTTTTCAGGAACTGGGTCATAAACTGTTCGGTTTTGCCTTTGTCACCACCGTCAGCTTTATCCGCTgc                                                        >  1:585490/1‑150 (MQ=255)
                                                                                              gAAGCTAATCAGCACATCGCCCAGGCCGCGCTCGGCAAAAGTGGTGGTCGCGCCACGACCGCCAGTATCGAACACTTCAACGTTTTACAGGAACTGGGTCATAAACTGTTCGGTTTTGCCTTTGTCACCACCGTCAGCTTTATCCGCTgc                                                        >  1:1689526/1‑150 (MQ=255)
                                                                                               aaGCTAATCAGCACATCGCCCAGGCCGCGCTCGGCAAAAGTGGTGGTCGCGCCACGACCGCCAGTATCGAACACTTCAACGTTTTTCAGGAACTGGGTCATAAACTGTTCGGTTTTGCCTTTGTCACCACCGTCAGCTTTATCCGCTgcg                                                       <  2:526013/150‑1 (MQ=255)
                                                                                                 gCTAATCAGCACATCGCCCAGGCCGCGCTCGGCAAAAGTGGTGGTCGCGCCACGACCGCCAGTATCGAACACTTCAACGTTTTTCAGGAACTGGGTCATAAACTGTTCGGTTTTGCCTTTGTCACCACCGTCAGCTTTATCCGCTGCGcc                                                     <  1:355368/150‑1 (MQ=255)
                                                                                                 gCTAATCAGCACATCGCCCAGGCCGCGCTCGGCAAAAGTGGTGGTCGCGCCACGACCGCCAGTATCGAACACATCAACGTTTTTCAGGAACTGGGTCATAAACTGTTCGGTTTTGCCTTTGTCACCACCGTCAGCTTTATCCGCTGCGcc                                                     <  2:255591/150‑1 (MQ=255)
                                                                                                    aaTCAGCACATCGCCCAGGCCGCGCTCGGCAAAAGTGGTGGTCGCGCCACGACCGCCAGTATCGAACACTTCAACGTTTTTCAGGAACTGGGTCATAAACTGTTCGGTTTTGCCTTTGTCACCACCGTCAGCTTTATCCGCTGCGCCCCa                                                  <  1:1356448/150‑1 (MQ=255)
                                                                                                        aGCACATCGCCCAGGCCGCGCTCGGCAAAAGTGGTGGTCGCGCCACGACCGCCAGTATCGAACACTTCAACGTTTTTCAGGAACTGGGTCATAAACTGTTCGGTTTTGCCTTTGTCACCACCGTCAGCTTTATCCGCTGCGCCCCATgcc                                              >  2:1103836/1‑150 (MQ=255)
                                                                                                          cacaTCGCCCAGGCCGCGCTCGGCAAAAGTGGTGGTCGCGCCACGACCGCCAGTATCGAACACTTCAACGTTTTTCAGGAACTGGGTCATAAACTGTTCGGTTTTGCCTTTGTCACCACCGTCAGCTTTATCCGCTGCGCCCCATgccgc                                            <  1:53541/150‑1 (MQ=255)
                                                                                                           acaTCGCCCAGGCCGCGCTCGGCAAAAGTGGTGGTCGCGCCACGACCGCCAGTATCGAACACTTCAACGTTTTTCAGGAACTGGGTCATAAACTGTTCGGTTTTGCCTTTGTCACCACCGTCAGCTTTATCCGCTGCGCCCCATgccgcc                                           <  1:1477193/150‑1 (MQ=255)
                                                                                                            caTCGCCCAGGCCGCGCTCGGCAAAAGTGGTGGTCGCGCCACGACCGCCAGTATCGAACACTTCAACGTTTTTCAGGAACTGGGTCATAAACTGTTCGGTTTTGCCTTTGTCACCACCGTCAGCTTTATCCGCTGCGCCCCCTGCCGCCa                                          >  2:1221811/1‑150 (MQ=255)
                                                                                                                gCCCAGGCCGCGCTCGGCAAAAGTGGTGGTCGCGCCACGACCGCCAGTATCGAACACTTCAACGTTTTTCAGGAACTGGGTCATAAACTGTTCGGTTTTGCCTTTGTCACCACCGTCAGCTTTATCCGCTGCGCCCCATGCCGCCAGATa                                      >  1:1000660/1‑150 (MQ=255)
                                                                                                                 cccAGGCCGCGCTCGGCAAAAGTGGTGGTCGCGCCACGACCGCCAGTATCGAACACTTCAACGTTTTTCAGGAACTGGGTCATAAACTGTTCGGTTTTGCCTTTGTCACCACCGTCAGCTTTATCCGCTGCGCCCCATGCCGCCAGATAg                                     >  1:1049365/1‑150 (MQ=255)
                                                                                                                  ccAGGCCGCGCTCGGCAAAAGTGGTGGTCGCGCCACGACCGCCAGTATCGAACACTTCAACGTTTTTCAGGAACTGGGTCATAAACTGTTCGGTTTTGCCTTTGTCACCACCGTCAGCTTTATCCGCTGCGCCCCATGCCGCCAGATAgg                                    <  1:591842/150‑1 (MQ=255)
                                                                                                                          cgcTCGGCAAAAGTGGTGGTCGCGCCACGACCGCCAGTATCGAACACTTCAACGTTTTTCAGGAACTGGGTCATAAACTGTTCGGTTTTGCCTTTGTCACCACCGTCAGCTTTATCCGCTGCGCCCCATGCCGCCAGATAGGTATAAcgc                            >  2:415349/1‑150 (MQ=255)
                                                                                                                           gcTCGGCAAAAGTGGTGGTCGCGCCACGACCGCCAGTATCGAACACTTCAACGTTTTTCAGGAACTGGGTCATAAACTGTTCGGTTTTGCCTTTGTCACCACCGTCAGCTTTATCCGCTGCGCCCCATGCCGCCAGATAGGTATAAcgcg                           <  1:1447944/150‑1 (MQ=255)
                                                                                                                            cTCGGCAAAAGTGGTGGTCGCGCCACGACCGCCAGTATCGAACACTTCAACGTTTTTCAGGAACTGGGTCATAAACTGTTCGGTTTTGCCTTTGTCACCACCGTCAGCTTTATCCGCTGCGCCCCATGCCGCCAGATAGGTATAAcgcgc                          >  2:432555/1‑150 (MQ=255)
                                                                                                                             tCGGCAAAAGTGGTGGTCGCGCCACGACCGCCAGTATCGAACACTTCAACGTTTTTCAGGAACTGGGTCATAAACTGTTCGGTTTTGCCTTTGTCACCACCGTCAGCTTTATCCGCTGCGCCCCATGCCGCCAGATAGGTATAAcgcgcg                         <  1:1160324/150‑1 (MQ=255)
                                                                                                                               ggCAAAAGTGGTGGTCGCGCCACGACCGCCAGTATCGAACACTTCAACGTTTTTCAGGAACTGGGTCATAAACTGTTCGGTTTTGCCTTTGTCACCACCGTCAGCTTTATCCGCTGCGCCCCATGCCGCCAGATAGGTATAACGCGCGtt                       >  1:1298768/1‑150 (MQ=255)
                                                                                                                                 cAAAAGTGGTGGTAGCGCGACGACCGCCAGTATCGAACACTTCAACGTTTTTCAGGAACTGGGTCATAAACTGTTCGGTTTTGCCTTTGTCACCACCGTCAGCTTTATCCGCTGCGCCCCATGCCGCCAGATAGGTATAACGCGCGTTAc                     >  2:64344/1‑150 (MQ=255)
                                                                                                                                  aaaaGTGGTGGTCGCGCCACGACCGCCAGTATCGAACACTTCAACGTTTTTCAGGAACTGGGTCATAAACTGTTCGGTTTTGCCTTTGTCACCACCGTCAGCTTTATCCGCTGCGCCCCATGCCGCCAGATAGGTATAACGCGCGTTAcc                    >  1:362344/1‑150 (MQ=255)
                                                                                                                                  aaaaGTGGTGGTCGCGCCACGACCGCCAGTATCGAACACTTCAACGTTTTTCAGGAACTGGGTCATAAACTGTTCGGTTTTGCCTTTGTCACCACCGTCAGCTGTATCCGCTGCGCCCCATGCCGCCAGATAGGTATAACGCGCGTTAcc                    <  1:1072521/150‑1 (MQ=255)
                                                                                                                                  aaaaGTGGTGGGCGCGCCACGACCGCCAGTATCGAACACTTCAACGTTTTTCAGGAACTGGGTCATAAACTGTTCGGTTTTGCCTTTGTCACCACCGTCAGCTTTATCCGCTGCGCCCCATGCCGCCAGATAGGTATAACGCGCGTTAcc                    <  1:1673058/150‑1 (MQ=255)
                                                                                                                                   aaaGTGGTGGTCGCGCCACGACCGCCAGTATCGAACACTTCAACGTTTTTCAGGAACTGGGTCATAAACTGTTCGGTTTTGCCTTTGTCACCACCGTCAGCTTTATCCGCTGCGCCCCATGCCGCCAGATAGGTATAACGCGCGTTAccc                   >  2:325462/1‑150 (MQ=255)
                                                                                                                                   aaaGTGGTGGTCGCGCCACGACCGCCAGTATCGAACACTTCAACGTTTTTCAGGAACTGGGTCATAAACTGTTCGGTTTTGCCTTTGTCACCACCGTCAGCTTTATCCGCTGCGCCCCATGCCGCCAGATAGGTATAACGCGCGTTAccc                   <  1:1038066/150‑1 (MQ=255)
                                                                                                                                      ctggcggTCGTGGCGCGACCGCCAGTATCGAACACTTCAACGTTTTTCAGGAACTGGGTCATAAACTGTTCGGTTTTGCCTTTGTCACCACCGTCAGCTTTATCCGCTGCGCGCCATGCCGCCAGATAGGTATAACGCGCGTTACCCGAc                >  1:298663/2‑150 (MQ=255)
                                                                                                                                      ctggcggTCGTGGCGCGACCACCAGTATCGAACACTTCAACGTTTTTCAGGAACTGGGTCATAAACTGTTCGGTTTTGCCTTTGTCACCACCGTCAGCTTTATCCGCTGCGCCCCATGCCGCCAGATAGGTATACCGCGCGTTACCCGAc                >  1:1118406/2‑150 (MQ=255)
                                                                                                                                      ctggcggTCGTGGCGCGACCACCAGTATCGAACACTTCAACGTTTTTCAGGAACTGGGTCATAAACTGTTCGGTTTTGCCTTTGTCACCACCGTCAGCTTTATCCGCTGCGCCCCATGCCGCCAGATAGGTATAACGCGCGTTACCCgcc                >  2:1448122/2‑148 (MQ=255)
                                                                                                                                        ggcggtCGTGGCGCGACCACCAGTATCGAACACTTCAACGTTTTTCAGGAACTGGGTCATAAACTGTTCGGTTTTGCCTTTGTCACCACCGTCAGCTTTATCCGCTGCGCCCCATGCCGCCAGATAGGTATAAcgc                            >  2:1634807/4‑136 (MQ=255)
                                                                                                                                        ggcggtCGTGGCGCGACCACCAGTATCGAACACTTCAACGTTTTTCAGGAACTGGGTCATAAACTGTTCGGTTTTGCCTTTGTCACCACCGTCAGCTTTATCCGCTGCGCCCCATGCCGCCAGATAGGTATAAcgc                            <  1:1634807/133‑1 (MQ=255)
                                                                                                                                           ggtCGCGCCACGACCGCCAGTATCGAACACTTCACCGTTTTTCAGGAACTGGGTCATAAACTGTTCGGTTTTGCCTTTGTCACCACCGTCAGCTTTATCCGCTGCGCCCCATGCCGCCAGATAGGTATAACGCGCGTTACCCGACGtttt           <  2:1690010/150‑1 (MQ=255)
                                                                                                                                            gtCGCGCCACGACCGCCAGTATCGAACACTTCAACGTTTTTCAGGAACTGGGTCATAAACTGTTCGGTTTTGCCTTTGTCACCACCGTCAGCTTTATCCGCTGCGCCCCa                                                  >  1:909471/1‑110 (MQ=255)
                                                                                                                                            gtCGCGCCACGACCGCCAGTATCGAACACTTCAACGTTTTTCAGGAACTGGGTCATAAACTGTTCGGTTTTGCCTTTGTCACCACCGTCAGCTTTATCCGCTGCGCCCCa                                                  <  2:909471/110‑1 (MQ=255)
                                                                                                                                              cgcgCCACGACCGCCAGTATCGAACACTTCCACGTTTTTCAGGAACTGGGTCATAAACTGTTCGGTTTTGCCTTTGTCACCACCGTCAGCTTTATCCGCTGCGCCCCATGCCGCCAGATAGGTATAACGCGCGTTACCCGACGTTttcgg        >  2:1231550/1‑150 (MQ=255)
                                                                                                                                              cgcgCCACGACCGCCAGTATCGAACACTTCAACGTTTTTCAGGAACTGGGTCATAAACTGTTCGGTTTTGCCTTTGTCACCACCGTCAGCTTTATCCGCTGCGCCCCATGCCGCCAGATAGGTATAACGCGCGTTACCCGACGTTttcgg        <  2:1326407/150‑1 (MQ=255)
                                                                                                                                              cgcgCCACGACCGCCAGTATCGAACACTTCAACGTTTTACAGGAACTGGGTCATAAACTGTTCGGTTTTGCCTTTGTCACCACCGTCAGCTTTATCCGCTGCGCCCCATGCCGCCAGATAGGTATAACGCGCGTTACCCGACGTTttcgg        <  2:1689526/150‑1 (MQ=255)
                                                                                                                                               gcgcCACGACCGCCAGTATCGAACACTTCAACGTTTTTCAGGAACTGGGTCATAAACTGTTCGGTTTTGCCTTTGTCACCACCGTCAGCTTTATCCGCTGCGCCCCATGCCGCCAGATAGGTATAACGCGCGTTACCCGACGTTttcggg       <  2:1328858/150‑1 (MQ=255)
                                                                                                                                               gcgcCACGACCGCCAGTATCGAACACTTCAACGTTTTTCAGGAACTGGGTCATAAACTGTTCGGTTTTGCCTTTGTCACCACCGTCAGCTTTATCCGCTGCGCCCCATGCCGCCAGATAGGTATAACGCGCGTTACCCGACGTTttcggg       <  2:1049365/150‑1 (MQ=255)
                                                                                                                                                 gcCACGACCGCCAGTATCGAACACTTCAACGTTTTTCAGGAACTGGGTCATAAACTGTTCGGTTTTGCCTTTGTCCCCACCGTCAGCTTTATCCGCTGCGCCCCATGCCGCCAGATAGGTATAACGCGCGTTACCCGACGTTttcgggtg     >  2:635630/1‑149 (MQ=255)
                                                                                                                                                    aCGACCGCCAGTATCGAACACTTCAACGTTTTTCAGGAACTGGGTCATAAACTGTTCGGTTTTGCCTTTGTCACCACCGTCAGCTTTATCCGCTGCGCCCCATGCCGCCAGATAGGTATAACGCGCGTTACCCGACGTTttcgggttcgg  <  1:703081/150‑1 (MQ=255)
                                                                                                                                                    aCGACCGCCAGTATCGAACACTTCAACGTTTTTCAGGAACTGGGTCATAAACTGTTCGGTTTTGCCTTTGTCACCACCGTCAGCTTTATCCGCTGCGCCCCATGCCGCCAGATAGGTATAACGCGCGTTACCCGACGTTttcgggttcgg  <  2:878953/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
CCACCGGGAATTCCGCCAGAATGTTGGTTTTCGGAATCACCACTTCAAAGCCCTGCGCTTCATACTGTTTACGGATGTTGTTCACTTCCGATTCGAAGCTAATCAGCACATCGCCCAGGCCGCGCTCGGCAAAAGTGGTGGTCGCGCCACGACCGCCAGTATCGAACACTTCAACGTTTTTCAGGAACTGGGTCATAAACTGTTCGGTTTTGCCTTTGTCACCACCGTCAGCTTTATCCGCTGCGCCCCATGCCGCCAGATAGGTATAACGCGCGTTACCCGACGTTTTCGGGTTCGG  >  NC_000913/2542784‑2543081

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: