Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,502,093 A→G 5.6% T178A (ACC→GCC)  tehB → tellurite methyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,502,0930AG 5.6% 266.7 / 3.6 88T178A (ACC→GCC) tehBtellurite methyltransferase
Reads supporting (aligned to +/- strand):  ref base A (43/40);  new base G (3/2);  total (46/42)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.74e-01

GTTACAACCTGATTGTGGCGGCGATGGATACCGCTGATTATCCATGTACCGTCGGCTTCCCGTTTGCCTTCAAAGAGGGAGAATTACGTCGATATTACGAAGGCTGGGAGAGGGTGAAATACAATGAAGACGTCGGCGAGCTGCACCGCACCGACGCCAACGGTAATCGTATTAAACTGCGTTTCGCCACGATGCTGGCACGTAAAAAATGACCCGGTAAGCACAAAACGCGTGAAAATTCCCCACGCTGAGATGATTTACTGTTCTTCTTTTCGGTAAGCATATTTTTTATCGAAGG  >  NC_000913/1501944‑1502241
                                                                                                                                                     |                                                                                                                                                    
gtTACAACCTGATTGTGGCGGCGATGGATACCGCTGATTATCCATGTACCGTCGGCTTCCCGTTTGCCTTCAAAGAGGGAGAATTACGTCGATATTACGAAGGCTGGGAGAGGGTGAAATACAATGAAGACGTCGGCGAGCTgcaccgca                                                                                                                                                      <  1:1037100/150‑1 (MQ=255)
    cAACCTGATTGTGGCGGCGATGGATACCGCTGATTATCCATGTACCGTCGGCTTCCCGTTTGCCTTCAAAGAGGGAGAATTACGTCGATATTACGAAGGCTGGGAGAGGGTGAAATACAATGAAGACGTCGGTGCGGTGCAGCTCGCCGa                                                                                                                                                  >  1:1284656/1‑150 (MQ=255)
     aaCCTGATTGTGGCGGCGATGGATACCGCTGATTATCCATGTACCGTCGGCTTCCCGTTTGCCTTCAAAGAGGGAGAATTACGTCGATATTACGAAGGCTGGGAGAGGGTGAAATACAATGAAGACGTCGGCGAGCTGCACCGCACCGAc                                                                                                                                                 >  1:567040/1‑150 (MQ=255)
      aCCTGATTGTGGCGGCGATGGATACCGCTGATTATCCATGTACCGTCGGCTTCCCGTTTGCCTTCAAAGAGGGAGAATTACGTCGATATTACGAAGGCTGGGAGAGGGTGAAATACAATGAAGACGTCGGCGAGCTGCACCGCACCGACg                                                                                                                                                <  2:714182/150‑1 (MQ=255)
      aCCTGATTGTGGCGGCGATGGATACCGCTGATTATCCATGTACCGTCGGCTTCCCGTTTGCCTTCAAAGAGGGAGAATTACGTCGATATTACGAAGGCTGGGAGAGGGTGAAATACAATGAAGACCTCGGCGAGCTGCACCGCACCGAcc                                                                                                                                                >  1:665699/1‑149 (MQ=255)
        cTGATTGTGGCGGCGATGGATACCGCTGATTATCCATGTACCGTCGGCTTCCCGTTTGCCTTCAAAGAGGGAGAATTACGTCGATATTACGAAGGCTGGGAGAGGGTGAAATACAATGAAGACGTCGGCGAGCTGCACCGCACCGACGcc                                                                                                                                              >  1:547973/1‑150 (MQ=255)
         tGATTGTGGCGGCGATGGATACCGCTGATTATCCATGTACCGTCGGCTTCCCGTTTGCCTTCAAAGAGGGAGAATTACGTCGATATTACGAAGGCTGGGAGAGGGTGAAATACAATGAAGACGTCGGCGAGCTGCACCGCACCGACGCCa                                                                                                                                             >  1:719416/1‑150 (MQ=255)
          gATTGTGGCGGCGATGGATACCGCTGATTATCCATGTACCGTCGGCTTCCCGTTTGCCTTCAAAGAGGGAGAATTACGTCGATATTACGAAGGCTGGGAGAGGGTGAAATACAATGAAGACGTCGGTGCGGTGCAGCTCGCCGACgtctt                                                                                                                                            <  2:1284656/150‑5 (MQ=255)
          gATTGTGGCGGCGATGGATACCGCTGATTATCCATGTACCGTCGGCTTCCCGTTTGCCTTCAAAGAGGGAGAATTACGTCGATATTACGAAGGCTGGGAGAGGGTGAAATACAATGAAGACGTCGGTGCGGTGCAGCTCGCCGACgtctt                                                                                                                                            <  1:871865/150‑5 (MQ=255)
          gATTGTGGCGGCGATGGATACCGCTGATTATCCATGTACCGTCGGCTTCCCGTTTGCCTTCAAAGAGGGAGAATTACGTCGATATTACGAAGGCTGGGAGAGGGTGAAATACAATGAAGACGTCGGCGAGCTGCACCGCACCGACGCCaa                                                                                                                                            >  1:1133467/1‑150 (MQ=255)
              gtgGCGGCGATGGATACCGCTGATTATCCATGTACCGTCGGCTTCCCGTTTGCCTTCAAAGAGGGAGAATTACGTCGATATTACGAAGGCTGGGAGAGGGTGAAATACAATGAAGACGTCGGCGAGCTGCACCGCACCGACGCCAACGGt                                                                                                                                        <  2:887569/150‑1 (MQ=255)
               tgGCGGCGATGGATACCGCTGATTATCCATGTACCGTCGGCTTCCCGTTTGCCTTCAAAGAGGGAGAATTACGTCGATATTACGAAGGCTGGGAGAGGGTGAAATACAATGAAGACGTCGGCGAGCTGCACCGCACCGACGCCAACGGTa                                                                                                                                       >  1:565208/1‑150 (MQ=255)
                ggcggcGATGGATACCGCTGATTATCCATGTACCGTCGGCTTCCCGTTTGCCTTCAAAGAGGGAGAATTACGTCGATAATACGAAGGCTGGGAGAGGGTGAAATACAATGAAGACGTCGGCGAGCTGCACCGCACCGACGCCGACGGTaa                                                                                                                                      >  2:971325/1‑150 (MQ=255)
                 gcggcgATGGATACCGCTGATTATCCATGTACCGTCGGCTTCCCGTTTGCCTTCAAAGAGGGAGAATTACGTCGATATTACGAAGGCTGGGAGAGGGTGAAATACAATGAAGACGTCGGCGAGCTGCACCGCACCGACGCCAACGGTAAt                                                                                                                                     >  1:919395/1‑150 (MQ=255)
                  cggctATGGATACCGCTGATTATCCATGTACCGTCGGCTTCCCGTTTGCCTTCAAAGAGGGAGAATTACGTCGATATTACGAAGGCTGGGAGAGGGTGAAATACAATGAAGACGTCGGCGAGCTGCACCGCACCGACGCCAACGGTAATc                                                                                                                                    <  2:1055193/150‑1 (MQ=255)
                    gcgATGGATACCGCTGATTATCCATGTACCGTCGGCTTCCCGTTTGCCTTCAAAGATGGAGAATTACGTCGATATTACGAAGGCTGGGAGAGGGTGAAATACAATGAAGACGTCGGCGAGCTGCACCGCACCGACGCCAACGGTAATCGt                                                                                                                                  >  2:609251/1‑150 (MQ=255)
                    gcgATGGATACCGCTGATTATCCATGTACCGTCGGCTTCCCGTTTGCCTTCAAAGAGGGAGAATTACGTCGATATTACGAAGGCTGGGAGAGGGTGAAATACAATGAAGACGTCGGCGAGCTGCACCGCACCGACGCCAACGGTAATCGt                                                                                                                                  >  2:1586560/1‑150 (MQ=255)
                     cgATGGATACCGCTGATTATCCATGTACCGTCGGCTTCCCGTTTGCCTTCAAAGAGGGAGAATTACGTCGATATTACGAAGGCTGGGAGAGGGTGAAATACAATGAAGACGTCGGCGAGCTGCACCGCACCGACGCCAACGGTAATCGTa                                                                                                                                 <  1:966886/150‑1 (MQ=255)
                      gATGGATACCGCTGATTATCCATGTACCGTCGGCTTCCCGTTTGCCTTCAAAGAGGGAGAATTACGTCGATATTACGAAGGCTGGGAGAGGGTGAAATACAATGAAGACGTCGGCGAGCTGCACCGCACCGACGCCAACGGTAATCGTAt                                                                                                                                >  2:1196964/1‑150 (MQ=255)
                       aTGGATACCGCTGATTATCCATGTACCGTCGGCTTCCCGTTTGCCTTCAAAGAGGGAGAATTACGTCGATATTACGAAGGCTGGGAGAGGGTGAAATACAATGAAGACGTCGGCGAGCTGCACCGCACCGACGCCAACGGTAATCGTAtt                                                                                                                               >  1:1033967/1‑150 (MQ=255)
                          gATACCGCTGATTATCCATGTACCGTCGGCTTCCCGTTTGCCTTCAAAGAGGGAGAATTACGTCGATATTACGAAGGCTGGGAGAGGGTGAAATACAATGAAGACGTCGGCGAGCTGCACCGCACCGACGCCAACGGTAATCGTATTaaa                                                                                                                            >  1:1315736/1‑150 (MQ=255)
                           aTACCGCTGATTATCCATGTACCGTCGGCTTCCCGTTTGCCTTCAAAGAGGGAGAATTACGTCGATATTACGAAGGCTGGGAGAGGGTGAAATACAATGAAGACGTCGGCGAGCTGCCCCGCACCGACGCCAACCGTAATCGTATTAAAc                                                                                                                           >  1:555934/1‑150 (MQ=255)
                              ccGCTGATTATCCATGTACCGTCGGCTTCCCGTTTGCCTTCAAAGAGGGAGAATTACGTCGATATTACGAAGGCTGGGAGAGGGTGAAATACAATGAAGACGTCGGCGAGCTGCACCGCACCGACGCCAACGGTAATCGTATTAAACTGc                                                                                                                        >  2:736371/1‑150 (MQ=255)
                                gCTGATTATCCATGTACCGTCGGCTTCCCGTTTGCCTTCAAAGAGGGAGAATTACGTCGATATTACGAAGGCTGGGAGAGGGTGAAATACAATGAAGACGTCGGCGAGCTGCACCGCACCGACGCCAACGGTAATCGTATTAAACTGCGt                                                                                                                      >  2:290715/1‑150 (MQ=255)
                                   gATTATCCATGTACCGTCGGCTTCCCGTTTGCCTTCAAAGAGGGAGAATTACGTCGATATTACGAAGGCTGGGAGAGGGTGAAATACAATGAAGACGTCGGCGAGCTGCACCGCACCGACGCCAACGGTAATCGTATTAAACTGCGTTTc                                                                                                                   <  1:346884/150‑1 (MQ=255)
                                     ttATCCATGTACCGTCGGCTTCCCGTTTGCCTTCAAAGAGGGAGAATTACGTCGATATTACGAAGGCTGGGAGAGGGTGAAATACAATGAAGACGTCGGCGAGCTGCTCCGCACCGACGCCAACGGTAATCGTATTAAACTGCGTTTCGc                                                                                                                 <  2:33744/150‑1 (MQ=255)
                                      tATCCATGTACCGTCGGCTTCCCGTTTGCGTTCAAAGAGTGAGAATTACGTCGATATTACGAAGGCTGGGAGAGGGTGAAATACAATGAAGACGTCGGCTAGCTGCACCGCACCGACGCCAAGGGTAATCGTATTAAACTGCGTTTCGcc                                                                                                                <  2:722546/150‑1 (MQ=255)
                                        tCCATGTACCGTCGGCTTCCCGTTTGCCTTCAAAGAGGGAGAATTACGTCGATATTACGAAGGCTGGGAGAGGGTGAAATACAATGAAGACGTCGGCGAGCTGCACCGCACCGACGCCAACGGTAATCGTATTAAACTGCGTTTCGCCAc                                                                                                              <  2:1449563/150‑1 (MQ=255)
                                          cATGTACCGTCGGCTTCCCGTTTGCCTTCAAAGAGGGAGAATTACGTCGATATTACGAAGGCTGGGAGAGGGTGAAATACAATGAAGACGTCGGCGAGCTGCACCGCACCGACGCCAACGGTAATCGTATTAAACTGCGTTTCGCCACGa                                                                                                            >  1:870033/1‑150 (MQ=255)
                                             gTACCGTCGGCTTCCCGTTTGCCTTCAAAGAGGGAGAATTACGTCGATATTACGAAGGCTGGGAGAGGGTGAAATACAATGAAGACGTCGGCGAGCTGCACCGCACCGACGCCAACGGTAATCGTATTAAACTGCGTTTCGCCACGATGc                                                                                                         >  2:361608/1‑150 (MQ=255)
                                                ccGTCGGCTTCCCGTTTGCCTTCAAAGAGGGAGAATTACGTCGATATTACGAAGGCTGGGAGAGGGTGAAATACAATGAAGACGTCGGCGAGCTGCACCGCACCGACGCCAACGGTAATCGTATTAAACTGCGTTTCGCCACGATGCTgg                                                                                                      >  1:818695/1‑150 (MQ=255)
                                                  gTCGGCTTCCCGTTTGCCTTCAAAGAGGGAGAATTACGTCGATATTACGAAGGCTGGGAGAGGGTGAAATACAATGAAGACGTCGGCGAGCTGCACCGCACCGACGCCAACGGTAATCGTATTAAACTGCGTTTCGCCACGATGCTGGCa                                                                                                    <  2:792857/150‑1 (MQ=255)
                                                    cGGCTTCCCGTTTGCCTTCAAAGAGGGAGAATTACGTCGATATTACGAAGGCTGGGAGAGGGTGAAATACAATGAAGACGTCGGCGAGCTGCACCGCACCGACGCCAACGGTAATCGTATTAAACTGCGTTTCGCCACGATGCTGGCACg                                                                                                  <  2:385151/150‑1 (MQ=255)
                                                    cGGCTTCCCGTTTGCCTTCAAAGAGGGAGAATTACGTCGATATTACGAAGGCTGGGAGAGGGTGAAATACAATGAAGACGTCGGCGAGCTGCACCGCACCGACGCCAACGGTAATCGTATTAAACTGCGTTTCGCCACGATGCTGGCACg                                                                                                  <  1:645996/150‑1 (MQ=255)
                                                       cTTCCCGTTTGCCTTCAAAGAGGGAGAATTACGTCGATATTACGAAGGCTGGGAGAGGGTGAAATACAATGAAGACGTCGGCGAGCTGCACCGCACCGACGCCAACGGTAATCGTATTAAACTGCGTTTCGCCACGATGCTGGCACGTaa                                                                                               >  2:809898/1‑150 (MQ=255)
                                                         tCCCGTTTGCCTTCAAAGAGGGAGAATTACGTCGATATTACGAAGGCTGGGAGAGGGTGAAATACAATGAAGACGTCGGCGAGCTGCACCGCACCGACGCCAACGGTAATCGTATTAAACTGCGTTTCGCCACGATGCTGGCACGTaaaa                                                                                             >  1:1331454/1‑150 (MQ=255)
                                                         tCCCGTTTGCCTTCAAAGAGGGAGAATTACGTCGATATTACGAAGGCTGGGAGAGGGTGAAATACAATGAAGACGTCGGCGAGCTGCACCGCACCGACGCCAACGGTAATCGTATTAAACTGCGTTTCGCCACGATGCTGGCACGTaaaa                                                                                             >  1:222324/1‑150 (MQ=255)
                                                           ccGTTTGCCTTCAAAGAGGGAGAATTACGTCGATATTACGAAGGCTGGGAGAGGGTGAAATACAATGAAGACGTCGGCGAGCTGCACCGCACCGACGCCAACGGTAATCGTATTAAACTGCGTTTCGCCACGATGCTGGCACGTaaaaaa                                                                                           <  1:1196964/150‑1 (MQ=255)
                                                                tGCCTTCAAAGAGGGAGAATTACGTCGATATTACGAAGGCTGGGAGAGTGTGAAATACAATGAAGACGTCGGCGAGCTGCACCGCACCGACGCCAACGGTAATCGTATTAAACTGCGTTTCGCCACGATGCTGGCACGTAAAAAATGAcc                                                                                      <  2:936423/150‑1 (MQ=255)
                                                                tGCCTTCAAAGAGGGAGAATTACGTCGATATTACGAAGGCTGGGAGAGGGTGAAATACAATGAAGACGTCGGCGAGCTGCACCGCACCGACGCCAACGGTAATCGTATTAAACTGCGTTTCGCCACGATGCTGGCACGTAAAAAATGAcc                                                                                      >  2:1461454/1‑150 (MQ=255)
                                                                tGCCTTCAAAGAGGGAGAATTACGTCGATATTACGAAGGCTGGGAGAGGGTGAAATACAATGAAGACGTCGGCGAGCTGCACCGCACCGACGCCAACGGTAATCGTATTAAACTGCGTTTCGCCACGATGCTGGCACGTAAAAAATGAcc                                                                                      <  1:1466108/150‑1 (MQ=255)
                                                                 gCCTTCAAAGAGGGAGAATTACGTCGATATTACGAAGGCTGGGAGAGGGTGAAATACAATGAAGACGTCGGCGAGCTGCACCGCACCGACGCCAACGGTAATCGTATTAAACTGCGTTTCGCCACGATGCTGGCACGTAAAAAATGAccc                                                                                     >  1:1370776/1‑150 (MQ=255)
                                                                  ccTTCAAAGAGGGAGAATTACGTCGATATTACGAAGGCTGGGAGAGGGTGAAATACAATGAAGACGTCGGCGAGCTGCACCGCACCGACGCCAACGGTAATCGTATTAAACTGCGTTTCGCCACGATGCTGGCACGTAAAAAATGACCCg                                                                                    >  2:994069/1‑150 (MQ=255)
                                                                     tCAAAGAGGGAGAATTACGTCGATATTACGAAGGCTGGGAGAGGGTGAAATACAATGAAGACGTCGGCGAGCTGCACCGCACCGACGCCAACGGTAATCGTATTAAACTGCGTTTCGCCACGATGCTGGCACGTAAAAAATGACCCGGTa                                                                                 <  1:609251/150‑1 (MQ=255)
                                                                      cAAAGAGGGAGAATTACGTCGATATTACGAAGGCTGGGAGAGGGTGAAATACAATGAAGACGTCGGCGAGCTGCACCGCACCGACGCCAACGGTAATCGTATTAAACTGCGTTTCGCCACGATGCTGGCACGTAAAAAATGACCCGGTaa                                                                                <  2:1329258/150‑1 (MQ=255)
                                                                       aaaGAGGGAGAATTACGTCGATATTACGAAGGCTGGGAGAGGGTGAAATACAATGAAGACGTCGGCGAGCTGCACCGCACCGACGCCAACGGTAATCGTATTAAACTGCGTTTCGCCACGATGCTGGCACGTAAAAAATGACCCGGTAAg                                                                               <  2:952684/150‑1 (MQ=255)
                                                                       aaaGAGGGAGAATTACGTCGATATTACGAAGGCTGGGAGAGGGTGAAATACAATGAAGACGTCGGCGAGCTGCACCGCACCGACGCCAACGGTAATCGTATTAAACTGCGTTTCGCCACGATGCTGGCACGTAAAAAATGACCCGGTAAg                                                                               >  1:642707/1‑150 (MQ=255)
                                                                          gagGGAGAATTACGTCGATATTACGAAGGCTGGGAGAGGGTGAAATACAATGAAGACGTCGGCGAGCTGCACCGCACCGACGCCAACGGTAATCGTATTAAACTGCGTTTCGCCACGATGCTGGCACGTAAAAAATGACCCGGTAAGcac                                                                            <  2:1370776/150‑1 (MQ=255)
                                                                                 aaTTACGTCGATATTACGAAGGCTGGGAGAGGGTGAAATACAATGAAGACGTCGGCGAGCTGCACCGCACCGACGCCAACGGTAATCGTATTAAACTGCGTTTCGCCACGATGCTGGCACGTAAAAAATGACCCGGTAAGCACAAAAcgc                                                                     <  2:421430/150‑1 (MQ=255)
                                                                                  aTTACGTCGATATTACGAAGGCTGGGAGAGGGTGAAATACAATGAAGACGTCGGCGAGCTGCACCGCACCGACGCCAACGGTAATCGTATTAAACTGCGTTTCGCCACGATGCTGGCACGTAAAAAATGACCCGGTAAGCACAAAAcgcg                                                                    <  1:337779/150‑1 (MQ=255)
                                                                                    tACGTCGATATTACGAAGGCTGGGAGAGGGTGAAATACAATGAAGACGTCGGCGAGCTGCACCGCACCGACGCCAACGGTAATCGTATTAAACTGCGTTTCGCCACGATGCTGGCACGTAAAAAATGACCCGGTAAGCACAAAACGCGTg                                                                  >  2:758222/1‑150 (MQ=255)
                                                                                    tACGTCGATATTACGAAGGCTGGGAGAGGGTGAAATACAATGAAGACGTCGGCGAGCTGCACCGCACCGACGCCAACGGTAATCGTATTAAACTGCGTTTCGCCACGATGCTGGCACGTAAAAAATGACCCGGTAAGCACAAAACGCGTg                                                                  <  2:1325911/150‑1 (MQ=255)
                                                                                       gTCGATATTACGAAGGCTGGGAGAGGGTGAAATACAATGAAGACGTCGGCGAGCTGCACCGCACCGACGCCAACGGTAATCGTATTAAACTGCGTTTCGCCACGATGCTGGCACGTAAAAAATGACCCGGTAAGCACAAAACGCGTGaaa                                                               <  2:567040/150‑1 (MQ=255)
                                                                                        tCGATATTACGAAGGCTGGGAGAGGGTGAAATACAATGAAGACGTCGGCGAGCTGCACCGCACCGACGCCAACGGTAATCGTATTAAACTGCGTTTCGCCACGATGCTGGCACGTAAAAAATGACCCGGTAAGCACAAAACGCGTGaaaa                                                              <  1:809898/150‑1 (MQ=255)
                                                                                         cGATATTACGAAGGCTGGGAGAGGGTGAAATACAATGAAGACGTCGGCGAGCTGCACCGCACCGACGCCAACGGTAATCGTATTAAACTGCGTTTCGCCACGATGCTGGCACGTAAAAAATGACCCGGTAAGCACAAAACGCGTGAAAAt                                                             >  2:588555/1‑150 (MQ=255)
                                                                                          gATATTACGAAGGCTGGGAGAGGGTGAAATACAATGAAGACGTCGGCGAGCTGCACCGCACCGACGCCAACGGTAATCGTATTAAACTGCGTTTCGCCACGATGCTGGCACGTAAAAAATGACCCGGTAAGCACAAAACGCGTGAAAAtt                                                            >  2:924765/1‑150 (MQ=255)
                                                                                           atatTACGAAGGCTGGGAGAGGGTGAAATACAATGAAGACGTCGGCGAGCTGCACCGCACCGACGCCAACGGTAATCGTATTAAACTGCGTTTCGCCACGATGCTGGCACGTAAAAAATGACCCGGTAAGCACAAAACGCGTGAAAATTc                                                           >  1:888493/1‑150 (MQ=255)
                                                                                              ttACGAAGGCTGGGAGAGGGTGAAATACAATGAAGACGTCGGCGAGCTGCACCGCACCGACGCCAACGGTAATCGTTTTAAACTTTGTTTGGGCACGATGCTTGCACGTAAAAAATGTCCCGGTAAGCACAAAATGCTTGAAAATTcgcc                                                        <  2:536345/150‑1 (MQ=255)
                                                                                               tACGAAGGCTGGGAGAGGGTGAAATACAATGAAGACGTCGGCGAGCTGCACCGCACCGACGCCAACGGTAATCGTATTAAACTGCGTTTCGCCACGATGCTGGCACGTAAAAAATGACCCGGTAAGCACAAAAcgc                                                                     >  2:1287575/1‑136 (MQ=255)
                                                                                               tACGAAGGCTGGGAGAGGGTGAAATACAATGAAGACGTCGGCGAGCTGCACCGCACCGACGCCAACGGTAATCGTATTAAACTGCGTTTCGCCACGATGCTGGCACGTAAAAAATGACCCGGTAAGCACAAAAcgc                                                                     <  1:1287575/136‑1 (MQ=255)
                                                                                               tACGAAGGCTGGGAGAGGGTGAAATACAATGAAGACGTCGGCGAGCTGCACCGCACCGACGCCAACGGTAATCGTATTAAACTGCGTTTCGCCACGATGCTGGCACGTAAAAAATGACCCGGTAAGCACAAAACGCGTGAAAATTCCCCa                                                       >  1:136515/1‑150 (MQ=255)
                                                                                                aCGAAGGCTGGGAGAGGGTGAAATACAATGAAGACGTCGGCGAGCTGCACCGCACCGACGCCAACGGTAATCGTATTAAACTGCGTTTCGCCACGATGCTGGCACGTAAAAAATGACCCGGTAAGCACAAAACGCGTGAAAATTCCCCAc                                                      >  2:1248964/1‑150 (MQ=255)
                                                                                                 cGAAGGCTGGGAGAGGGTGAAATACAATGAAGACGTCGGCGAGCTGCACCGCACCGACGCCAACGGTAATCGTATTAAACTGCGTTTCGCCACGATGCTGGCACGTAAAAAATGACCCGGTAAGCACAAAACGCGTGAAAATTCCCCACg                                                     >  1:38313/1‑150 (MQ=255)
                                                                                                     ggCTGGGAGAGGGTGAAATACAATGAAGACGTCGGCGAGCTGCACCGCACCGACGCCAACGGTAATCGTATTAAACTGCGTTTCGCCACGATGCTGGCACGTAAAAAATGACCCGGTAAGCACAAAACGCGTGAAAATTCCCCACGCTga                                                 <  2:547973/150‑1 (MQ=255)
                                                                                                       cTGGGAGAGGGTGAAATACAATGAAGACGTCGGCGAGCTGCACCGCACCGACGCCAACGGTAATCGTATTAAACTGCGTTTCGCCACGATGCTGGCACGTAAAAAATGACCCGGTAAGCACAAAACGCGTGAAAATTCCCCACGCTgaga                                               <  1:1586560/150‑1 (MQ=255)
                                                                                                        tGGGAGAGGGTGAAATACAATGAAGACGTCGGCGAGCTGCACCGCACCGACGCCAACGGTAATCGTATTAAACTGCGTTTCGCCACGATGCTGGCACGTAAAAAATGACCCGGTAAGCACAAAACGCGTGAAAATTCCCCACGCTGAgat                                              <  2:1120782/150‑1 (MQ=255)
                                                                                                         gggAGAGGGTGAAATACAATGAAGACGTCGGCGAGCTGCACCGCACCGACGCCAACGGTAATCGTATTAAACTGCGTTTCGCCACGATGCTGGCACGTAAAAAATGACCCGGTAAGCACAAAACGCGTGAAAATTCCCCACGCTGAgatg                                             >  2:1356381/1‑150 (MQ=255)
                                                                                                               gggTGAAATACAATGAAGACGTCGGCGAGCTGGACCGCACCGACGCCAACGGTAATCGTATTAAACTGCGTTTCGCCACGATGCTGGCACGTAAAAAATGACCCGGTAAGCACAAAACGCGTGAAAATTCCCCACGCTGAGATGATTTAc                                       <  1:1158872/150‑1 (MQ=255)
                                                                                                               gggTGAAATACAATGAAGACGTCGGCGAGCTGCACCGCACCGACGCCAACGGTAATCGTATTAAACTGCGTTTCGCCACGATGCTGGCACGTAAAAAATGACCCGGTAAGCACAAAACGCGTGAAAATTCCCCACGCTGAGATGATTTAc                                       <  2:1295785/150‑1 (MQ=255)
                                                                                                               gggTGAAATACAATGAAGACGTCGGCGAGCTGCACCGCACCGACGCCAACGGTAATCGTATTAAACTGCGTTTCGCCACGATGCTGGCACGTAAAAAATGACCCGGTAAGCACAAAACGCGTGAAAATTCCCCACGCTGAGATGATTTAc                                       <  1:736371/150‑1 (MQ=255)
                                                                                                                  tGAAATACAATGAAGACGTCGGCGAGCTGCACCGCACCGACGCCAACGGTAATCGTATTAAACTGCGTTTCGCCACGATGCTGGCACGTAAAAAATGACCCGGTAAGCACAAAACGCGTGAAAATTCCCCACGCTGAGATGATTTACTGt                                    <  1:108645/150‑1 (MQ=255)
                                                                                                                   gAAATACAATGAAGACGTCGGCGAGCTGCACCGCACCGACGCCAACGGTAATCGTATTAAACTGCGTTTCGCCACGATGCTGGCACGTAAAAAATGACCCGGTAAGCACAAAACGCGTGAAAATTCCCCACGCTGAGATGATTTACTGtt                                   >  1:1090965/1‑150 (MQ=255)
                                                                                                                       tACAATGAAGACGTCGGCGCGCTGCACCGCACCGCCGCCAACGGTAATCGTATTAAACTGCGTTTCGCCACGATGATGGCACGTAAAAAATGACCCGGTACGCCCAAAACGCGTGAAACTTCCCGACGCTGAGATGATTTCCTGttcttc                               >  2:234922/1‑150 (MQ=255)
                                                                                                                          aaTGAAGACGTCGGCGAGCTGCACCGCACCGACGCCAACGGTAATCGTATTAAACTGCGTTTCGCCACGATGCTGGCACGTAAAAAATGACCCGGTAAGCACAAAACGCGTGAAAATTCCCCACGCTGAGATGATTTACTGTTCTTCttt                            <  2:818695/150‑1 (MQ=255)
                                                                                                                           ccgttggcGTCGGTGCGGTGCAGCTCGCCGACGCCAACGGTAATCGTATTAAACTGCGTTTCGCCACGATGCTGGCACGTAAAAAATGACACGGTAAGCACAAAACGCGTGAAAATTCCCCACGCTGAGATGATTTACTGTTCTTCtttt                           >  2:1089814/8‑150 (MQ=255)
                                                                                                                           aTGAAGACGTCGGCGAGCTGCACCGCACCGACGCCAACGGTAATCGTATTAAACTGCGTTTCGCCACGATGCTGGCACGTAAAAAATGACCCGGTAAGCACAAAACGCGTGAAAATTCCCCACGCTGAGATGATTTACTGTTCTTCtttt                           >  1:634293/1‑150 (MQ=255)
                                                                                                                           aTGAAGACGTCGGCGAGCTGCACCGCACCGACGCCAACGGTAATCGTATTAAACTGCGTTTCGCCACGATGCTGGCACGTAAAAAATGACCCGGTAAGCACAAAACGCGTGAAAATTCCCCACGCTGAGATGATTTACTGTTCTTCtttt                           >  2:1465964/1‑150 (MQ=255)
                                                                                                                            tGAAGACGTCGGCGAGCTGCACCGCACCGACGCCAACGGTAATCGTATTAAACTGCGTTTCGCCACGATGCTGGCACGTAAAAAATGACCCGGTAAGCACAAAACGCGTGAAAATTCCCCACGCTGAGATGATTTACTGTTCTTCTTTTc                          <  2:1603342/150‑1 (MQ=255)
                                                                                                                            tGAAGACGTCGGCGAGCTGCACCGCACCGACGCCAACGGTAATCGTATTAAACTGCGTTTCGCCACGATGCTGGCACGTAAAAAATGACCCGGTAAGCACAAAACGCGTGAAAATTCCCCACGCTGAGATGATTTACTGTTCTTCTTTTc                          <  2:565208/150‑1 (MQ=255)
                                                                                                                            tGAAGACGTCGGCGAGCTGCACCGCACCGACGCCAACGGTAATCGTATTAAACTGCGTTTCGCCACGATGCTGGCACGTAAAAAATGACCCGGTAAGCACAAAACGCGTGAAAATTCCCCACGCTGAGATGATTTACTGTTCTTCTTTTc                          <  2:222324/150‑1 (MQ=255)
                                                                                                                              aaGACGTCGGCGAGCTGCACCGCACCGACGCCAACGGTAATCGTATTAAACTGCGTTTCGCCACGATGCTGGCACGTAAAAAATGACCCGGTAAGCACAAAACGCGTGAAAATTCCCCACGCTGAGATGATTTACTGTTCTTCTTTTCgg                        >  2:1394057/1‑150 (MQ=255)
                                                                                                                              aaGACGTCGGCGAGCTGCACCGCACCGACGCCAACGGTAATCGTATTAAACTGCGTTTCGCCACGATGCTGGCACGTAAAAAATGACCCGGTAAGCACAAAACGCGTGAAAATTCCCCACGCTGAGATGATTTACTGTTCTTCTTTTCgg                        <  2:1033967/150‑1 (MQ=255)
                                                                                                                              aaGACGTCCGCGAGCTGCACCGCACCGACGCCAACGGTAATCGTATTAAACTGCGTTTCGCCACGATGCTGGCACGTAAAAAATGACCAGGTAAGCACAAAACGCGTGAAAATTACCCACGCTGAGATGATTTACTGTTCTTCTTTTCgg                        >  2:1245120/1‑150 (MQ=255)
                                                                                                                                 aCGTCGGCGAGCTGCACCGCACCGACGCCAACGGTAATCGTATTAAACTGCGTTTCGCCACGATGCTGGCACGTAAAAAATGACCCGGTAAGCACAAAACGCGTGAAAATTCCCCACGCTGAGATGATTTACTGTTCTTCTTTTCGGTaa                     <  1:341287/150‑1 (MQ=255)
                                                                                                                                  cGTCGGTGCGGTGCAGCTCGCCGACGCCAACGGTAATCGTATTAAACTGCGTTTCGCCACGATGCTGGCACGTAAAAAATGACCCGGTAAGCACAAAACGCGTGAAAATTCCCCACGCTGAGATGATTTACTGTTCTTCTTTTCGGTAAg                    >  1:247098/1‑150 (MQ=255)
                                                                                                                                  cGTCGGCGAGCTGCACCGCACCGACGCCAACGGTAATCGTATTAAACTGCGTTTCGCCACGATGCTGGCACGTAAAAAATGACCCGGTAAGCACAAAACGCGTGAAAATTCCCCACGCTGAGATGATTTACTGTTCTTCTTTTCGGTAAg                    >  1:1139870/1‑150 (MQ=255)
                                                                                                                                  cGTCGGCGAGCTGCACCGCACCGACGCCAACGGTAATCGTATTAAACTGCGTTTCGCCACGATGCTGGCACGTAAAAAATGACCCGGTAAGCACAAAACGCGTGAAAATTCCCCACGCTGAGATGATTTACTGTTCTTCTTTTCGGTAAg                    <  2:266995/150‑1 (MQ=255)
                                                                                                                                    tCGGCGAGCTGCACCGCACCGACGCCAACGGTAATCGTATTAAACTGCGTTTCGCCACGATGCTGGCACGTAAAAAATGACCCGGTAAGCACAAAACGCGTGAAAATTCCCCACGCTGAGATGATTTACTGTTCTTCTTTTCGGTAAGCa                  <  1:1356381/150‑1 (MQ=255)
                                                                                                                                             tGCACCGCACCGACGCCAACGGTAATCGTATTAAACTGCGTTTCGCCACGATGCTGGCACGTAAAAAATGACCCGGTAAGCACAAAACGCGTGAAAATTCCCCACGCTGAGATGATTTACTGTTCTTCTTTTCGGTAAGCATATTTTTTa         >  1:381591/1‑150 (MQ=255)
                                                                                                                                               caccgcaccgACGCCAACGGTAATCGTATTAAACTGCGTTTCGCCACGATGCTGGCACGTAAAAAATGACCCGGTAAGCACAAAACGCGTGAAAATTCCCCACGCTGAGATGATTTACTGTTCTTCTTTTCGGTAAGCATATTTTTTATc       >  1:879258/1‑150 (MQ=255)
                                                                                                                                                accgcaccgACGCCAACGGTAATCGTATTAAACTGCGTTTCGCCACGATGCTGGCACGTAAAAAATGACCCGGTAAGCACAAAACGCGTGAAAATTCCCCACGCTGAGATGATTTACTGTTCTTCTTTTCGGTAAGCATATTTTTTATCg      >  1:887793/1‑150 (MQ=255)
                                                                                                                                                 ccgcaccgACGCCAACGGTAATCGTATTAAACTGCGTTTCGCCACGATGCTGGCACGTAAAAAATGACCCGGTAAGca                                                                             <  2:228792/78‑1 (MQ=255)
                                                                                                                                                 ccgcaccgACGCCAACGGTAATCGTATTAAACTGCGTTTCGCCACGATGCTGGCACGTAAAAAATGACCCGGTAAGca                                                                             >  1:228792/1‑78 (MQ=255)
                                                                                                                                                 ccgcaccgACGCCAACGGTAATCGTATTAAACTGCGTTTCGCCACGATGCTGGCACGTAAAAAATGACCCGGTAAGCACAAAACGCGTGAAAATTCCCCACGCTGAGATGATTTACTGTTCTTCTTTTCGGTAAGCATATTTTTTATCGa     >  2:824328/1‑150 (MQ=255)
                                                                                                                                                 ccgcaccgACGCCAACGGTAATCGTATTAAAATGCGTTTCGCCACGATGCTGGCACGTAAAAAATGCCCCGGTAAGCACAAAACGCGTGAAAATTCCCCACGCTGAGATGATTTACTGTTCTTCTGTTCGGTAAGCATATTTTTTATCGa     >  2:580733/1‑150 (MQ=255)
                                                                                                                                                    caccgACGCCAACGGTAATCGTATTAAACTGCGTTTCGCCACGATGCTGGCACGTAAAAAATGACCCGGTAAGCACAAAACGCGTGAAAATTCCCCACGCTGAGATGATTTACTGTTCTTCTTTTCGGTAAGCATATTTTTTATCGAAgg  <  1:1248964/150‑1 (MQ=255)
                                                                                                                                                    caccgACGCCAACGGTAATCGTATTAAACTGCGTTTCGCCACGATGCTGGCACGTAAAAAATGACCCGGTAAGCACAAAACGCGTGAAAATTCCCCACGCTGAGATGATTTACTGTTCTTCTTTTCGGTAAGCATATTTTTTATCGAAgg  <  1:290715/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
GTTACAACCTGATTGTGGCGGCGATGGATACCGCTGATTATCCATGTACCGTCGGCTTCCCGTTTGCCTTCAAAGAGGGAGAATTACGTCGATATTACGAAGGCTGGGAGAGGGTGAAATACAATGAAGACGTCGGCGAGCTGCACCGCACCGACGCCAACGGTAATCGTATTAAACTGCGTTTCGCCACGATGCTGGCACGTAAAAAATGACCCGGTAAGCACAAAACGCGTGAAAATTCCCCACGCTGAGATGATTTACTGTTCTTCTTTTCGGTAAGCATATTTTTTATCGAAGG  >  NC_000913/1501944‑1502241

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: