Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,502,084 C→G 5.4% L175V (CTG→GTG)  tehB → tellurite methyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,502,0840CG 5.4% 315.3 / 6.5 92L175V (CTG→GTG) tehBtellurite methyltransferase
Reads supporting (aligned to +/- strand):  ref base C (44/43);  new base G (3/2);  total (47/45)
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.83e-01

CCTGGTGGTTACAACCTGATTGTGGCGGCGATGGATACCGCTGATTATCCATGTACCGTCGGCTTCCCGTTTGCCTTCAAAGAGGGAGAATTACGTCGATATTACGAAGGCTGGGAGAGGGTGAAATACAATGAAGACGTCGGCGAGCTGCACCGCACCGACGCCAACGGTAATCGTATTAAACTGCGTTTCGCCACGATGCTGGCACGTAAAAAATGACCCGGTAAGCACAAAACGCGTGAAAATTCCCCACGCTGAGATGATTTACTGTTCTTCTTTTCGGTAAGCA  >  NC_000913/1501937‑1502225
                                                                                                                                                   |                                                                                                                                             
ccTGTTGGTTACAACCTGATTGTGGCGTCGATGGATACCGCTGATTATCCATGTACCGTCGTCTTCCCGTTTGCCTTCAAAGAGGGAGAATTACGTCGATATTACGAAGGCTGGGAGAGGGTTAAATACAATGAAGACGTCGGCGAGCTg                                                                                                                                             <  1:578518/150‑1 (MQ=255)
    gtggtTACAACCTGATTGTGGCGGCGATGGATACCGCTGATTATCCATGTACCGTCGGCTTCCCGTTTGCCTTCAAAGAGGGAGAATTACGTCGATATTACGAAGGCTGGGAGAGGGTGAAATACAATGAAGACGTCGGCGAGCTgcacc                                                                                                                                         >  1:266995/1‑150 (MQ=255)
    gtggtTACAACCTGATTGTGGCGGCGATGGATACCGCTGATTATCCATGTACCGTCGGCTTCCCGTTTGCCTTCAAAGAGGGAGAATTACGTCGATATTACGAAGGCTGGGAGAGGGTGAAATACAATGAAGACGTCGGCGAGCTgcacc                                                                                                                                         <  2:1502803/150‑1 (MQ=255)
     tggtTACAACCTGATTGTGGCGGCGATGGATACCGCTGATTATCCATGTACCGTCGGCTTCCCGTTTGCCTTCAAAGAGGGAGAATTACGTCGATATTACGAAGGCTGGGAGAGGGTGAAATACAATGAAGACGTCGGCGAGCTgcaccg                                                                                                                                        >  1:556304/1‑150 (MQ=255)
       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                                                                                  <  1:645996/150‑1 (MQ=255)
                                                           cGGCTTCCCGTTTGCCTTCAAAGAGGGAGAATTACGTCGATATTACGAAGGCTGGGAGAGGGTGAAATACAATGAAGACGTCGGCGAGCTGCACCGCACCGACGCCAACGGTAATCGTATTAAACTGCGTTTCGCCACGATGCTGGCACg                                                                                  <  2:385151/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                                                               >  1:642707/1‑150 (MQ=255)
                                                                              aaaGAGGGAGAATTACGTCGATATTACGAAGGCTGGGAGAGGGTGAAATACAATGAAGACGTCGGCGAGCTGCACCGCACCGACGCCAACGGTAATCGTATTAAACTGCGTTTCGCCACGATGCTGGCACGTAAAAAATGACCCGGTAAg                                                               <  2:952684/150‑1 (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:1325911/150‑1 (MQ=255)
                                                                                           tACGTCGATATTACGAAGGCTGGGAGAGGGTGAAATACAATGAAGACGTCGGCGAGCTGCACCGCACCGACGCCAACGGTAATCGTATTAAACTGCGTTTCGCCACGATGCTGGCACGTAAAAAATGACCCGGTAAGCACAAAACGCGTg                                                  >  2:758222/1‑150 (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                       <  1:736371/150‑1 (MQ=255)
                                                                                                                      gggTGAAATACAATGAAGACGTCGGCGAGCTGCACCGCACCGACGCCAACGGTAATCGTATTAAACTGCGTTTCGCCACGATGCTGGCACGTAAAAAATGACCCGGTAAGCACAAAACGCGTGAAAATTCCCCACGCTGAGATGATTTAc                       <  2:1295785/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:565208/150‑1 (MQ=255)
                                                                                                                                   tGAAGACGTCGGCGAGCTGCACCGCACCGACGCCAACGGTAATCGTATTAAACTGCGTTTCGCCACGATGCTGGCACGTAAAAAATGACCCGGTAAGCACAAAACGCGTGAAAATTCCCCACGCTGAGATGATTTACTGTTCTTCTTTTc          <  2:222324/150‑1 (MQ=255)
                                                                                                                                   tGAAGACGTCGGCGAGCTGCACCGCACCGACGCCAACGGTAATCGTATTAAACTGCGTTTCGCCACGATGCTGGCACGTAAAAAATGACCCGGTAAGCACAAAACGCGTGAAAATTCCCCACGCTGAGATGATTTACTGTTCTTCTTTTc          <  2:1603342/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)
                                                                                                                                                   |                                                                                                                                             
CCTGGTGGTTACAACCTGATTGTGGCGGCGATGGATACCGCTGATTATCCATGTACCGTCGGCTTCCCGTTTGCCTTCAAAGAGGGAGAATTACGTCGATATTACGAAGGCTGGGAGAGGGTGAAATACAATGAAGACGTCGGCGAGCTGCACCGCACCGACGCCAACGGTAATCGTATTAAACTGCGTTTCGCCACGATGCTGGCACGTAAAAAATGACCCGGTAAGCACAAAACGCGTGAAAATTCCCCACGCTGAGATGATTTACTGTTCTTCTTTTCGGTAAGCA  >  NC_000913/1501937‑1502225

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: