Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,502,089 C→G 5.7% H176Q (CAC→CAG tehB → tellurite methyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,502,0890CG 5.7% 300.3 / 6.6 88H176Q (CAC→CAGtehBtellurite methyltransferase
Reads supporting (aligned to +/- strand):  ref base C (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 = 6.10e-01

GTGGTTACAACCTGATTGTGGCGGCGATGGATACCGCTGATTATCCATGTACCGTCGGCTTCCCGTTTGCCTTCAAAGAGGGAGAATTACGTCGATATTACGAAGGCTGGGAGAGGGTGAAATACAATGAAGACGTCGGCGAGCTGCACCGCACCGACGCCAACGGTAATCGTATTAAACTGCGTTTCGCCACGATGCTGGCACGTAAAAAATGACCCGGTAAGCACAAAACGCGTGAAAATTCCCCACGCTGAGATGATTTACTGTTCTTCTTTTCGGTAAGCATATTTTTTATCGA  >  NC_000913/1501941‑1502238
                                                                                                                                                    |                                                                                                                                                     
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                                                                            <  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: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                        >  2:1465964/1‑150 (MQ=255)
                                                                                                                              aTGAAGACGTCGGCGAGCTGCACCGCACCGACGCCAACGGTAATCGTATTAAACTGCGTTTCGCCACGATGCTGGCACGTAAAAAATGACCCGGTAAGCACAAAACGCGTGAAAATTCCCCACGCTGAGATGATTTACTGTTCTTCtttt                        >  1:634293/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:1033967/150‑1 (MQ=255)
                                                                                                                                 aaGACGTCGGCGAGCTGCACCGCACCGACGCCAACGGTAATCGTATTAAACTGCGTTTCGCCACGATGCTGGCACGTAAAAAATGACCCGGTAAGCACAAAACGCGTGAAAATTCCCCACGCTGAGATGATTTACTGTTCTTCTTTTCgg                     >  2:1394057/1‑150 (MQ=255)
                                                                                                                                 aaGACGTCCGCGAGCTGCACCGCACCGACGCCAACGGTAATCGTATTAAACTGCGTTTCGCCACGATGCTGGCACGTAAAAAATGACCAGGTAAGCACAAAACGCGTGAAAATTACCCACGCTGAGATGATTTACTGTTCTTCTTTTCgg                     >  2:1245120/1‑150 (MQ=255)
                                                                                                                                    aCGTCGGCGAGCTGCACCGCACCGACGCCAACGGTAATCGTATTAAACTGCGTTTCGCCACGATGCTGGCACGTAAAAAATGACCCGGTAAGCACAAAACGCGTGAAAATTCCCCACGCTGAGATGATTTACTGTTCTTCTTTTCGGTaa                  <  1:341287/150‑1 (MQ=255)
                                                                                                                                     cGTCGGTGCGGTGCAGCTCGCCGACGCCAACGGTAATCGTATTAAACTGCGTTTCGCCACGATGCTGGCACGTAAAAAATGACCCGGTAAGCACAAAACGCGTGAAAATTCCCCACGCTGAGATGATTTACTGTTCTTCTTTTCGGTAAg                 >  1:247098/1‑150 (MQ=255)
                                                                                                                                     cGTCGGCGAGCTGCACCGCACCGACGCCAACGGTAATCGTATTAAACTGCGTTTCGCCACGATGCTGGCACGTAAAAAATGACCCGGTAAGCACAAAACGCGTGAAAATTCCCCACGCTGAGATGATTTACTGTTCTTCTTTTCGGTAAg                 <  2:266995/150‑1 (MQ=255)
                                                                                                                                     cGTCGGCGAGCTGCACCGCACCGACGCCAACGGTAATCGTATTAAACTGCGTTTCGCCACGATGCTGGCACGTAAAAAATGACCCGGTAAGCACAAAACGCGTGAAAATTCCCCACGCTGAGATGATTTACTGTTCTTCTTTTCGGTAAg                 >  1:1139870/1‑150 (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                                                                          >  1:228792/1‑78 (MQ=255)
                                                                                                                                                    ccgcaccgACGCCAACGGTAATCGTATTAAACTGCGTTTCGCCACGATGCTGGCACGTAAAAAATGACCCGGTAAGca                                                                          <  2:228792/78‑1 (MQ=255)
                                                                                                                                                    ccgcaccgACGCCAACGGTAATCGTATTAAACTGCGTTTCGCCACGATGCTGGCACGTAAAAAATGACCCGGTAAGCACAAAACGCGTGAAAATTCCCCACGCTGAGATGATTTACTGTTCTTCTTTTCGGTAAGCATATTTTTTATCGa  >  2:824328/1‑150 (MQ=255)
                                                                                                                                                    ccgcaccgACGCCAACGGTAATCGTATTAAAATGCGTTTCGCCACGATGCTGGCACGTAAAAAATGCCCCGGTAAGCACAAAACGCGTGAAAATTCCCCACGCTGAGATGATTTACTGTTCTTCTGTTCGGTAAGCATATTTTTTATCGa  >  2:580733/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
GTGGTTACAACCTGATTGTGGCGGCGATGGATACCGCTGATTATCCATGTACCGTCGGCTTCCCGTTTGCCTTCAAAGAGGGAGAATTACGTCGATATTACGAAGGCTGGGAGAGGGTGAAATACAATGAAGACGTCGGCGAGCTGCACCGCACCGACGCCAACGGTAATCGTATTAAACTGCGTTTCGCCACGATGCTGGCACGTAAAAAATGACCCGGTAAGCACAAAACGCGTGAAAATTCCCCACGCTGAGATGATTTACTGTTCTTCTTTTCGGTAAGCATATTTTTTATCGA  >  NC_000913/1501941‑1502238

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: