Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,542,930 C→G 7.0% G200A (GGC→GCC)  cysP ← thiosulfate‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,542,9300CG 7.0% 347.4 / 12.9 111G200A (GGC→GCC) cysPthiosulfate‑binding protein
Reads supporting (aligned to +/- strand):  ref base C (60/42);  new base G (4/4);  total (65/46)
Fisher's exact test for biased strand distribution p-value = 7.18e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.03e-01

Only 100 of 115 total aligned reads displayed.

CGCCACCGGGAATTCCGCCAGAATGTTGGTTTTCGGAATCACCACTTCAAAGCCCTGCGCTTCATACTGTTTACGGATGTTGTTCACTTCCGATTCGAAGCTAATCAGCACATCGCCCAGGCCGCGCTCGGCAAAAGTGGTGGTCGCGCCACGACCGCCAGTATCGAACACTTCAACGTTTTTCAGGAACTGGGTCATAAACTGTTCGGTTTTGCCTTTGTCACCACCGTCAGCTTTATCCGCTGCGCCCCATGCCGCCAGATAGGTATAACGCGCGTTACCCGACGTTTTCGGGTT  >  NC_000913/2542782‑2543078
                                                                                                                                                    |                                                                                                                                                    
  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:1477193/1‑150 (MQ=255)
                gCCAGAATGTTGGTTTTCGGAATCACCACTTCAAAGCCCTGCGCTTCATACTGTTTACGGATGTTGTTCACTTCCGATTCGAAGCTAATCAGCACATCGCCCAGGCCGCGCTCGGCAAAAGTGGTGGTCGCGCCACGACCGCCAGTATCg                                                                                                                                     >  2:183991/1‑150 (MQ=255)
                   aGAATGTTGGTTTTCGGAATCACCACTTCAAAGCCCTGCGCTTCATACTGTTTACGGATGTTGTTCACTTCCGATTCGAAGCTAATCAGCACATCGCCCAGGCCGCGCTCGGCAAAAGTGGTGGTCGCGCCACGACCGCCAGTATCGAac                                                                                                                                  >  2:1569066/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)
                              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:1379925/1‑150 (MQ=255)
                                        accacTTCAAAGCCCTGCGCTTCATACTGTTTACGGATGTTGTTCACTTCCGATTCGAAGCTAATCAGCACATCGCCCAGGCCGCGCTCGGCAAAAGTGGTGGTCGCGCCACGACCGCCAGTATCGAACACTTCAACGTTTTTCAGGAAc                                                                                                             >  1:1450271/1‑150 (MQ=255)
                                        accacTTCAAAGCCCTGCGCTTCATACTGTTTACGGATGTTGTTCACTTCCGATTCGAAGCTAATCAGCACATCGCCCAGGCCGCGCTCGGCAAAAGTGGTGGTCGAGCCACGACCGCCAGTATCGAACACTTCAACGTTTTTCAGGAAc                                                                                                             >  2:1356448/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                                                                                                         >  2:1447944/1‑150 (MQ=255)
                                            cTTCAAAGCCCTGCGCTTCATACTGTTTACGGATGTTGTTCACTTCCGATTCGAAGCTAATCAGCACATCGCCCAGGCCGCGCTCGGCAAAAGTGGTGGTCGCGCCACGACCGCCAGTATCGAACACTTCAACGTTTTTCAGGAACTggg                                                                                                         <  1:1488660/150‑1 (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)
                                                       tGCGCTTCATACTGTTTACGGATGTTGTTCACTTCCGATTCGAAGCTAATCAGCACATCGCCCAGGCCGCGCTCGGCAAAAGTGGTGGTCGCGCCACGACCGCCAGTATCGAACACTTCAACGTTTTTCAGGAACTGGGTCATAAACTGt                                                                                              >  1:1457948/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)
                                                                    gTTTCCGGATGTTGTTCACTTCCGATTCGAAGCTAATCTGCACATCGCCCAGGCCGCGCTCGGCAAAAGTGGTGGTCGCGCCACGACCGCCAGTATCGAACACTTCAACGTTTTTCAGGAACTGGGTCATAAACTGTTCGGTTTTGCCtt                                                                                 <  1:227091/150‑1 (MQ=255)
                                                                    gTTTACGGATGTTGTTCACTTCCGATTCGAAGCTAATCAGCACATCGCCCAGGCCGCGCTCGGCAAAAGTGGTGGTCGCGCCACGACCGCCAGTATCGAACACTTCAACGTTTTTCAGGAACTGGGTCATAAACTGTTCGGTTTTGCCtt                                                                                 >  2:1160324/1‑150 (MQ=255)
                                                                    gTTTACGGATGTTGTTCACATCCGATTCGAAGCTAATCAGCACATCGCCCAGGCCGCGCTCGGCAAAAGTGGTGGTCGCGCCACGACCCCCAGTATCGAACACTTCAACGTTATTCAGGAACTGGGTCATAAAATCATCGGTATTGCCtt                                                                                 >  2:485284/1‑150 (MQ=255)
                                                                    gTTTACGGATGGTGTTCAATTCCGATTCGAAGCTAATCAGCACATCGCCCAGGCCGCGCTCGGCAAAAGTGGTGTTCGCGCCACGACCGCCAGTATCGAACACTGCAACGTTTTTCAGGAACTGGGTCATAAACTGTTCGGTTTTGCCtt                                                                                 <  1:222012/150‑1 (MQ=255)
                                                                     tttACGTATGTTGTTCACTTCCGATTCGAAGCTAATCAGCACATCGCCCAGGCCGCGCTCGGCAAAAGTGGTGGTCGCGCCACGACCGCCAGTATCGAACACTTCAACGTTTTTCAGGAACTGGGTCATAAACTGTTCGGTTTTGCCttt                                                                                <  1:181819/150‑1 (MQ=255)
                                                                      ttACGGATGTTGTTCACTTCCGATTCGAAGCTAATCAGCACATCGCCCAGGCCGCGCTCGGCAAAAGTGGTGGTCGCGCCACGACCGCCAGTATCGAACACTTCAACGTTTTTCAGGAACTGGGTCATAAACTGTTCGGTTTTGCCTTTg                                                                               <  1:678055/150‑1 (MQ=255)
                                                                        aCGGATGTTGTTCACTTCCGAATCGAAGCTAATCAGCACATCGCCCAGGCCGCGCTCGGCAAAAGTGGTGGTCGCGCCACGACCGCCAGTATCGAACACTTCAACGTTTTTCAGGAACTGGGTCATAAACTGTTCGGTTTTGCCTTTGTc                                                                             <  2:1181809/150‑1 (MQ=255)
                                                                          ggATGTTGTTCACTTCCGATTCGAAGCTAATCAGCACATCGCCCAGGCCGCGCTCGGCAAAAGTGGTGGTCGCGCCACGACCGCCAGTATCGAACACTTCAACGTTTTTCAGGAACTGGGTCATAAACTGTTCGGTTTTGCCTTTGTcac                                                                           >  1:394227/1‑150 (MQ=255)
                                                                           gATGTTGTTCACTTCCGATTCGAAGCTAATCAGCACATCGCCCAGGCCGCGCTCGGCAAAAGTGGTGGTCGCGCCACGACCGCCAGTATCGAACACTTCAACGTTTTTCAGGAACTGGGTCATAAACTGTTCGGTTTTGCCTTTGTcacc                                                                          <  2:1300326/150‑1 (MQ=255)
                                                                             tgttgtTCACTTCCGATTCGAAGCTAATCAGCACATCGCCCAGGCCGCGCTCGGCAAAAGTGGTGGTCGCGCCACGACCGCCAGTATCGAACACTTCAACGTTTTTCAGGAACTGGGTCATAAACTGTTCGGTTTTGCCTTTGTcaccac                                                                        <  1:543042/150‑1 (MQ=255)
                                                                                    cACTTCCGATTCGAAGCTAATCAGCACATCGCCCAGGCCGCGCTCGGCAAAAGTGGTGGTCGCGCCACGACCGCCAGTATCGAACACTTCAACGTTTTTCAGGAACTGGGTCATAAACTGTTCGGTTTTGCCTTTGTCACCACCGTCAGc                                                                 >  1:1328858/1‑150 (MQ=255)
                                                                                    cACTTCCGATTCGAAGCTAATCAGCACATCGCCCAGGCCGCGCTCGGCAAAAGTGGTGGTCGCGCCACGACCGCCAGTATCGAACACTTCAACGTTTTTCAGGAACTGGGTCATAAACTGTTCGGTTTTGCCTTTGTCACCACCGTCAGc                                                                 <  1:1038257/150‑1 (MQ=255)
                                                                                     aCTTCCGATTCGAAGCTAATCAGCACATCGCCCAGGCCGCGCTCGGCAAAAGTGGTGGTCGCGCCACGACCGCCAGTATCGAACACTTCAACGTTTTTCAGGAACTGGGTCATAAACTGTTCGGTTTTGCCTTTGTCACCACCGTCAGCt                                                                >  1:1670602/1‑150 (MQ=255)
                                                                                      cTTCCGATTCGAAGCTAATCAGCACATCGCCCAGGCCGCGCTCGGCAAAAGTGGTGGTCGCGCCACGACCGCCAGTATCGAACACTTCAACGTTTTTCAGGAACTGGGTCATAAACTGTTCGGTTTTGCCTTTGTCACCACCGTCAGCtt                                                               >  2:703081/1‑150 (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)
                                                                                          cGATTCGAAGCTAATCAGCACATCGCCCAGGCCGCGCTCGGCAAAAGTGGTGGTCGCGCCACGACCGCCAGTATCGAACACTTCAACGTTTTTCAGGAACTGGGTCATAAACTGTTCGGTTTTGCCTTTGTCACCACCGTCAGCTTTATc                                                           >  1:394460/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)
                                                                                                     tAATCAGCACATCGCCCAGGCCGCGCTCGGCAAAAGTGGTGGTCGCGCCACGACCGCCAGTATCGAACACTTCAACGTTTTTCAGGAACTGGGTCATAAACTGTTCGGTTTTGCCTTTGTCACCACCGTCAGCTTTATCCGCTGCGcccc                                                >  1:770501/1‑150 (MQ=255)
                                                                                                          aGCACATCGCCCAGGCCGCGCTCGGCAAAAGTGGTGGTCGCGCCACGACCGCCAGTATCGAACACTTCAACGTTTTTCAGGAACTGGGTCATAAACTGTTCGGTTTTGCCTTTGTCACCACCGTCAGCTTTATCCGCTGCGCCCCATgcc                                           >  2:1103836/1‑150 (MQ=255)
                                                                                                            cacaTCGCCCAGGCCGCGCTCGGCAAAAGTGGTGGTCGCGCCACGACCGCCAGTATCGAACACTTCAACGTTTTTCAGGAACTGGGTCATAAACTGTTCGGTTTTGCCTTTGTCACCACCGTCAGCTTTATCCGCTGCGCCCCATgccgc                                         <  1:53541/150‑1 (MQ=255)
                                                                                                             acaTCGCCCAGGCCGCGCTCGGCAAAAGTGGTGGTCGCGCCACGACCGCCAGTATCGAACACTTCAACGTTTTTCAGGAACTGGGTCATAAACTGTTCGGTTTTGCCTTTGTCACCACCGTCAGCTTTATCCGCTGCGCCCCATgccgcc                                        >  2:240262/1‑150 (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)
                                                                                                                             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)
                                                                                                                                    aaaaGTGGTGGTCGCGCCACGACCGCCAGTATCGAACACTTCAACGTTTTTCAGGAACTGGGTCATAAACTGTTCGGTTTTGCCTTTGTCACCACCGTCAGCTTTATCCGCTGCGCCCCATGCCGCCAGATAGGTATAACGCGCGTTAcc                 >  1:362344/1‑150 (MQ=255)
                                                                                                                                    aaaaGTGGTGGTCGCGCCACGACCGCCAGTATCGAACACTTCAACGTTTTTCAGGAACTGGGTCATAAACTGTTCGGTTTTGCCTTTGTCACCACCGTCAGCTGTATCCGCTGCGCCCCATGCCGCCAGATAGGTATAACGCGCGTTAcc                 <  1:1072521/150‑1 (MQ=255)
                                                                                                                                    aaaaGTGGTGGGCGCGCCACGACCGCCAGTATCGAACACTTCAACGTTTTTCAGGAACTGGGTCATAAACTGTTCGGTTTTGCCTTTGTCACCACCGTCAGCTTTATCCGCTGCGCCCCATGCCGCCAGATAGGTATAACGCGCGTTAcc                 <  1:1673058/150‑1 (MQ=255)
                                                                                                                                     aaaGTGGTGGTCGCGCCACGACCGCCAGTATCGAACACTTCAACGTTTTTCAGGAACTGGGTCATAAACTGTTCGGTTTTGCCTTTGTCACCACCGTCAGCTTTATCCGCTGCGCCCCATGCCGCCAGATAGGTATAACGCGCGTTAccc                <  1:1038066/150‑1 (MQ=255)
                                                                                                                                     aaaGTGGTGGTCGCGCCACGACCGCCAGTATCGAACACTTCAACGTTTTTCAGGAACTGGGTCATAAACTGTTCGGTTTTGCCTTTGTCACCACCGTCAGCTTTATCCGCTGCGCCCCATGCCGCCAGATAGGTATAACGCGCGTTAccc                >  2:325462/1‑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)
                                                                                                                                           gcggtCGCGCCACGACCGCCAGTATCGAACACTTCAACGTTTTTCAGGAACTGGGTCATAAACTGTTCGGTTTTGCCTTTGTCACCACCGTCAGCTTTATCCGCTGCGCCCCATGCCGCCAGATAGGTATAACGCGCGTTACCCGACGtt          >  1:230607/3‑150 (MQ=255)
                                                                                                                                             ggtCGCGCCACGACCGCCAGTATCGAACACTTCACCGTTTTTCAGGAACTGGGTCATAAACTGTTCGGTTTTGCCTTTGTCACCACCGTCAGCTTTATCCGCTGCGCCCCATGCCGCCAGATAGGTATAACGCGCGTTACCCGACGtttt        <  2:1690010/150‑1 (MQ=255)
                                                                                                                                             ggtCGCGCCACGACCGCCAGTATCGAACACTTCAACGTTTTTCAGGAACTGGGTCATAAACTGTTCGGTTTTGCCTTTGTCACCACCGTCAGCTTTATCCGCTGCGCCCCATGCCGCCAGATAGGTATAACGCGCGTTACCCGACGtttt        >  2:1027303/1‑150 (MQ=255)
                                                                                                                                              gtCGCGCCACGACCGCCAGTATCGAACACTTCAACGTTTTTCAGGAACTGGGTCATAAACTGTTCGGTTTTGCCTTTGTCACCACCGTCAGCTTTATCCGCTGCGCCCCa                                               >  1:909471/1‑110 (MQ=255)
                                                                                                                                              gtCGCGCCACGACCGCCAGTATCGAACACTTCAACGTTTTTCAGGAACTGGGTCATAAACTGTTCGGTTTTGCCTTTGTCACCACCGTCAGCTTTATCCGCTGCGCCCCa                                               <  2:909471/110‑1 (MQ=255)
                                                                                                                                                cgcgCCACGACCGCCAGTATCGAACACTTCCACGTTTTTCAGGAACTGGGTCATAAACTGTTCGGTTTTGCCTTTGTCACCACCGTCAGCTTTATCCGCTGCGCCCCATGCCGCCAGATAGGTATAACGCGCGTTACCCGACGTTttcgg     >  2:1231550/1‑150 (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)
                                                                                                                                                    |                                                                                                                                                    
CGCCACCGGGAATTCCGCCAGAATGTTGGTTTTCGGAATCACCACTTCAAAGCCCTGCGCTTCATACTGTTTACGGATGTTGTTCACTTCCGATTCGAAGCTAATCAGCACATCGCCCAGGCCGCGCTCGGCAAAAGTGGTGGTCGCGCCACGACCGCCAGTATCGAACACTTCAACGTTTTTCAGGAACTGGGTCATAAACTGTTCGGTTTTGCCTTTGTCACCACCGTCAGCTTTATCCGCTGCGCCCCATGCCGCCAGATAGGTATAACGCGCGTTACCCGACGTTTTCGGGTT  >  NC_000913/2542782‑2543078

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: