Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,660,438 A→G 6.5% S790S (TCA→TCG ynfE → putative selenate reductase, periplasmic

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,660,4380AG 6.5% 346.1 / 10.9 108S790S (TCA→TCGynfEputative selenate reductase, periplasmic
Reads supporting (aligned to +/- strand):  ref base A (46/55);  new base G (4/3);  total (50/58)
Fisher's exact test for biased strand distribution p-value = 7.02e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.68e-01

Only 100 of 113 total aligned reads displayed.

CGAAGTTCGGCTACCAGCCAAAGTAACACCACGTATTCTCCCTGGAGTTAGCGCTATGGGCCAGGGAGCCTGGCACGAGGCCAATATGTCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATAAATTCTGGCATGACATTTTGTTTGAAAAGCAATAAGTGAGTAATGATGAAAATCCATACCACAGAGGCGCTCAT  >  NC_000913/1660291‑1660584
                                                                                                                                                   |                                                                                                                                                  
cGAAGTTCGGCTACCAGCCAAAGTAACACCACGTATTCTCCCTGGAGTTAGCGCTATGGGCCAGGGAGCCTGGCACGAGGCCAATATGTCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTcac                                                                                                                                                   <  2:562544/149‑1 (MQ=255)
 gAAGTTCGGCTACCAGCCAAAGTAACACCACGTATTCTCCCTGGAGTTAGCGCTATGGGCCAGGGAGCCTGGCACGAGGCCAATATGTCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTcacc                                                                                                                                                  >  3:301889/1‑149 (MQ=255)
 gAAGTTCGGCTACCAGCCAAAGTAACACCACGTATTCTCCCTGGAGTTAGCGCTATGGGCCAGGGAGCCTGGCACGAGGCCAATATGTCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTcacc                                                                                                                                                  >  2:510933/1‑149 (MQ=255)
     ttCGGCTAGCAGCCAAAGTAGCGCGACGTATTCTCCCTGGAGTTAGCGCTATGGGCCAGGGAGCCTCGCACGAGGCCAATATGTCTGGGGATAAAATCGACCATGGCGGCTGTGTGAATACGCTATCCACTCTGCGATCTTCACCACTg                                                                                                                                              <  4:718474/149‑1 (MQ=255)
     ttCGGCTACCAGCCAAAGTAACACCACGTATTCTCCCTGGAGTTAGCGCTATGGGCAAGGGAGCCTGGCACGAGGCCAATATGTCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTCCACCACTg                                                                                                                                              >  3:710038/1‑149 (MQ=255)
      tCGGCTACCAGCCAAAGTAACACCACGTATTCTCCCTGGAGTTAGCGCTATGGGCCAGGGAGCCTGGCACGAGGCCAATATGTCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTgg                                                                                                                                             <  1:731395/149‑1 (MQ=255)
      tCGGCTACCAGCCAAAGTAACACCACGTATTCTCCCTGGAGTTAGCGCTATGGGCCAGGGAGCCTGGCACGAGGCCAATATGTCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTgg                                                                                                                                             <  2:38641/149‑1 (MQ=255)
        ggCTACCGGCCAAAGTACCACCACGTCTTCTCCCTGGAGTTAGCGCTATGGGCCAGGGAGCCTGGCACGAGGCCAATATGTCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCg                                                                                                                                           <  2:463932/149‑1 (MQ=255)
           tACCAGCCAAAGTAACACCACGTATTCTCCCTGGAGTTAGCGCTATGGGCCAGGGAGCCTGGCACGAGGCCAATATGTCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAg                                                                                                                                        <  1:611509/149‑1 (MQ=255)
           tACCAGCCAAAGTAACACCACGTATTCTCCCTGGAGTTAGCGCTATGGGCCAGGGAGCCTGGCACGAGGCCAATATGTCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAg                                                                                                                                        <  1:628748/149‑1 (MQ=255)
           tACCAGCCAAAGTAACACCACGTATTCTCCCTGGAGTTAGCGCTATGGGCCAGGGAGCCTGGCACGAGGCCAATATGTCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAg                                                                                                                                        <  4:327799/149‑1 (MQ=255)
            aCCAGCCAAAGTAACACCACGTATTCTCCCTGGAGTTAGCGCTATGGGCCAGGGAGCCTGGCACGAGGCCAATATGTCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAgg                                                                                                                                       <  1:590280/149‑1 (MQ=255)
               acccAAAGTAACACCACGTATTCTCCCTGGAGTTAGCGCTATGGGCCAGGGAGCCTGGCACGAGGCCAATATGTCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGaa                                                                                                                                    <  2:332694/147‑1 (MQ=255)
                gCCAAAGTAACACCACGTATTCTCCCTGGAGTTAGCGCTATGGGCCAGGGAGCCTGGCACGAGGCCAATATGTCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGaaa                                                                                                                                   >  3:321225/1‑149 (MQ=255)
                 ccAAAGTAACACCACGTATTCTCCCTGGAGTTAGCGCTATGGGCCAGGGAGCCTGGCACGAGGCCAATATGTCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAAc                                                                                                                                  <  4:788328/149‑1 (MQ=255)
                   aaaGTAACACCACGTATTCTCCCTGGAGTTAGCGCTATGGGCCAGGGAGCCTGGCACGAGGCCAATATGTCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAAccc                                                                                                                                >  4:181491/1‑149 (MQ=255)
                      gTAACACCACGTATTCTCCCTGGAGTTAGCGCTATGGGCCAGGGAGCCTGGCACGAGGCCAATATGTCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCgca                                                                                                                             >  4:149052/1‑149 (MQ=255)
                          caccacGTATTCTCCCTGGAGTTAGCGCTATGGGCCAGGGAGCCTGGCACGAGGCCAATATGTCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGcac                                                                                                                         <  3:181491/149‑1 (MQ=255)
                           accacGTATTCTCCCTGGAGTTAGCGCTATGGGCCAGGGAGCCTGGCACGAGGCCAATATGTCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGcaca                                                                                                                        >  1:1074239/1‑149 (MQ=255)
                           accacGTATTCTCCCTGGAGTTAGCGCTATGGGCCAGGGAGCCTGGCACGAGGCCAATATGTCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGcaca                                                                                                                        <  1:1050226/149‑1 (MQ=255)
                           accacGTATTCTCCCTGGAGTTAGCGCTATGGGCCAGGGAGCCTGGCACGAGGCCAATATGTCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGcaca                                                                                                                        >  2:899073/1‑149 (MQ=255)
                             cacGTATTCTCCCTGGAGTTAGCGCTATGGGCCAGGGAGCCTGGCACGAGGCCAATATGTCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACt                                                                                                                      >  1:261833/1‑149 (MQ=255)
                              acGTATTCTCCCTGGAGTTAGCGCTATGGGCCAGGGAGCCTGGCACGAGGCCAATATGTCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTa                                                                                                                     >  2:14357/1‑149 (MQ=255)
                               cGTATTCTCCCTGGAGTTAGCGCTATGGGCCAGGGAGCCTGGCACGAGGCCAATATGTCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTaa                                                                                                                    <  2:1074239/149‑1 (MQ=255)
                                     ctcCCTGGAGTTAGCGCTATGGGCCAGGGAGCCTGGTACGAGGCCAATATGTCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGt                                                                                                              >  2:590656/1‑149 (MQ=255)
                                      tcCCTGGAGTTAGCGCTATGGGCCAGGGAGCCTGGCACGTGGCCAATATGTCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTc                                                                                                             >  3:331556/1‑149 (MQ=255)
                                      tcCCTGGAGTTAGCGCTATGGGCCAGGGAGCCTGGCACGAGGCCAATATGTCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTc                                                                                                             >  3:997832/1‑149 (MQ=255)
                                      tcCCTGGAGTTAGCGCTATGGGCCAGGGAGCCTGGCACGAGGCCAATATGTCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTc                                                                                                             >  4:1118100/1‑149 (MQ=255)
                                        ccTGGAGTTAGCGCTATGGGCCAGGGAGCCTGGCACGAGGCCAATATGTCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCga                                                                                                           <  1:38008/149‑1 (MQ=255)
                                        ccTGGAGTTAGCGCTATGGGCCAGGGAGCCTGGCACGAGGCCAATATGTCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCga                                                                                                           <  1:14357/149‑1 (MQ=255)
                                        ccTGGAGTTAGCGCTATGGGCCAGGGAGCCTGGCACGAGGCCAATATGTCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACACTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCga                                                                                                           <  2:519748/149‑1 (MQ=255)
                                         cTGGAGTTAGCGCTATGGGCCAGGGAGCCTGGCACGAGGCCAATATGTCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCgag                                                                                                          >  1:172797/1‑149 (MQ=255)
                                           ggAGTTAGCGCTATGGGCCAGGGAGCCTGGCACGAGGCCAATATGTCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGAt                                                                                                        <  3:103531/149‑1 (MQ=255)
                                              gTTAGCGCTATGGGCCAGGGAGCCTGGCACGAGGCCAATATGTCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGa                                                                                                     >  1:205738/1‑149 (MQ=255)
                                                 aGCGCTATGGGCCAGGGAGCCTGGCACGAGGCCAATATGTCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTcac                                                                                                                                                   <  1:460374/100‑1 (MQ=255)
                                                 aGCGCTATGGGCCAGGGAGCCTGGCACGAGGCCAATATGTCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGaaaa                                                                                                  >  1:1033308/1‑149 (MQ=255)
                                                   cgcTATGGGCCAGGGAGCCTGGCACGAGGCCAATATGTCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGaaaaaa                                                                                                <  4:873412/149‑1 (MQ=255)
                                                       aTGGGCCAGGGAGCCTGGCACGAGGCCAATATGTCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAAtata                                                                                            <  2:1033308/149‑1 (MQ=255)
                                                          tgCCAGGGAGCCTGGCACGAGGCCAATATGTCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGaaa                                                                                                                                   >  4:309342/2‑107 (MQ=255)
                                                          ggCCAGGGAGCCTGGCACGAGGCCAATATGTCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAAcc                                                                                         >  3:681182/1‑149 (MQ=255)
                                                          ggCCAGGGAGCCTGGCACGAGGCCAATATGGCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGaaa                                                                                                                                   <  3:309342/107‑1 (MQ=255)
                                                            ccAGGGAGCCTGGCACGAGGCCAATATGTCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCa                                                                                       <  4:321225/149‑1 (MQ=255)
                                                            ccAGGGAGCCTGGCACGAGGCCAATATGTCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCa                                                                                       >  2:212352/1‑149 (MQ=255)
                                                             cAGGGAGCCTGGCACGAGGCCAATATGTCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCAc                                                                                      <  3:953169/149‑1 (MQ=255)
                                                                          aCGAGGCCAATATGTCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATAAAt                                                                         >  1:268675/1‑149 (MQ=255)
                                                                            gAGGCCAATATGTCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATc                                                                                                                  <  1:315971/106‑1 (MQ=255)
                                                                            gAGGCCAATATGTCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATc                                                                                                                  >  2:315971/1‑106 (MQ=255)
                                                                            gAGGCCAATATGTCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATAAATTc                                                                       <  4:681182/149‑1 (MQ=255)
                                                                               gCCAATATGTCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATAAATTCTgg                                                                    <  2:205738/149‑1 (MQ=255)
                                                                               gCCAATATGTCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATAAATTCTgg                                                                    >  3:606514/1‑149 (MQ=255)
                                                                                  aaTATGTCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCg                                                                                                            >  2:141977/1‑106 (MQ=255)
                                                                                  aaTATGTCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAata                                                                                              >  4:551912/1‑120 (MQ=255)
                                                                                  aaTATGTCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAata                                                                                              <  3:551912/120‑1 (MQ=255)
                                                                                  aaTATGTCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATaa                                                                                           >  2:344477/1‑123 (MQ=255)
                                                                                  aaTATGTCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATaa                                                                                           <  1:344477/123‑1 (MQ=255)
                                                                                        tCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTCACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATAAATTCTGGCATGACAtt                                                           >  2:697053/1‑149 (MQ=255)
                                                                                          tGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATAAATTCTGGCATGACAtttt                                                         >  2:294655/1‑149 (MQ=255)
                                                                                           ggCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATAAATTCTGGCATGACATtttg                                                        >  2:188910/1‑149 (MQ=255)
                                                                                             cGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATAAATTCTGGCATGACATtttgtt                                                      <  2:172797/149‑1 (MQ=255)
                                                                                             cGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATAAATTCTGGCATGACATtttgtt                                                      <  1:908041/149‑1 (MQ=255)
                                                                                              gATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATAAATTCTGGCATGACATtttgttt                                                     <  4:301889/149‑1 (MQ=255)
                                                                                               aTAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATAAATTCTGGCATGACATtttgtttg                                                    >  2:1111069/1‑149 (MQ=255)
                                                                                               aTAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATAAATTCTGGCATGACATtttgtt                                                      <  2:998711/147‑1 (MQ=255)
                                                                                               aTAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATAAATTCTGGCATGACATtttgtt                                                      >  1:998711/1‑147 (MQ=255)
                                                                                                tAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATAAATTCTGGCATGACATTTTGTTTGa                                                   >  1:1030590/1‑149 (MQ=255)
                                                                                                tAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATAAATTCTGGCATGACATTTTGTTTGa                                                   >  1:451417/1‑149 (MQ=255)
                                                                                                   aaTCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATAAATTCTGGCATGACATTTTGTTTGaaaa                                                <  1:697053/149‑1 (MQ=255)
                                                                                                   aaTCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATAAATTCTGGCATGACATTTTGTTTGaaaa                                                <  2:930585/149‑1 (MQ=255)
                                                                                                     tCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATAAATTCTGGCATGACATTTTGTTTGAAAAGc                                              <  3:149052/149‑1 (MQ=255)
                                                                                                       gACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCTCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATAAATTCTGGCATGACATTTTGTTTGAAAAGCaa                                            >  4:662299/1‑149 (MQ=255)
                                                                                                         ccATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATAAATTCTGGCATGACATTTTGTTTGAAAAGCAATa                                          >  1:680802/1‑149 (MQ=255)
                                                                                                            tGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATAAATTCTGGCATGACATTTTGTTTGAAAAGCAATAAGt                                       <  1:188910/149‑1 (MQ=255)
                                                                                                                    gtgtgAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATAAATTCTGGCATGACATTTTGTTTGAAAAGCAATAAGTGAGTAat                                <  1:309181/148‑1 (MQ=255)
                                                                                                                    gtgtgAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATAAATTCTGGCATGACATTTTGTTTGAAAAGCAATAAGTGAGTAat                                >  2:309181/1‑148 (MQ=255)
                                                                                                                         aaTACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATAAATTCTGGCATGACATtttgt                                                       >  1:121456/1‑120 (MQ=255)
                                                                                                                         aaTACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATAAATTCTGGCATGACATtttgt                                                       <  2:121456/120‑1 (MQ=255)
                                                                                                                          aTACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAAGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATAAATTCTGGCATGACATTTTGTTTGAAAAGCAATAAGTGAGTAATGATGaaa                         <  4:331556/149‑1 (MQ=255)
                                                                                                                           tACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATAAATTCTGGCATGACATTTTGTTTGAAAAGCAATAAGTGAGTAATGATGaaaa                        >  1:3064/1‑149 (MQ=255)
                                                                                                                              gCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATAAATTCTGGCATGACATtttgt                                                       <  2:610456/115‑1 (MQ=255)
                                                                                                                              gCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATAAATTCTGGCATGACATtttgt                                                       >  1:610456/1‑115 (MQ=255)
                                                                                                                              gCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATAAATTCTGGCATGACATTTTGTTTGAAAAGCAATAAGTGa                                     >  2:466381/1‑133 (MQ=255)
                                                                                                                              gCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATAAATTCTGGCATGACATTTTGTTTGAAAAGCAATAAGTGa                                     <  1:466381/133‑1 (MQ=255)
                                                                                                                              gCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATAAATTCTGGCATGACATTTTGTTTGAAAAGCAATAAGTGAGTAATGATGAAAATcc                     <  3:526879/149‑1 (MQ=255)
                                                                                                                                tAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCATCACACGAATCTGGTCGAGATCTAACCAATATAACCCACGACAACCATAAATTCTGGCATGACATTTTGTTTGAAAAGCAATAAGTGAGTAATGATGAAAATCCAt                   >  1:138984/1‑149 (MQ=255)
                                                                                                                                 acccACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATAAATTCTGGCATGACATTTTGTTTGAAAAGCAATTAGTGAGTAATGATGAAATTCCATa                  <  2:680802/147‑1 (MQ=255)
                                                                                                                                         tGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATAAATTCTGGCATGACATTTTGTTTGAAAAGCAATAAGTGAGTAATGATGAAAATCCATACCACagag          <  2:3064/149‑1 (MQ=255)
                                                                                                                                         tGCGCCCTTCACCACTGCCGAAGGGAAACCCGAAGCACACTACTCTGGTCGAGATAGAAACAATATAAACCACGACAACCATAAATTCTGGCATGACATTTTGTTTGAAAAGCAATAAGTGAGTAATGATGAAAAACCATACCACaaag          <  2:138984/149‑1 (MQ=255)
                                                                                                                                          tttcccTTCGCCAGTGGTGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATAAATTCTGGCATGACATTTTGTTTGAAAAGCAATAAGTGAg                                    <  2:1042856/119‑1 (MQ=255)
                                                                                                                                          tttcccTTCGCCAGTGGTGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATAAATTCTGGCATGACATTTTGTTTGAAAAGCAATAAGTGAg                                    >  1:1042856/4‑122 (MQ=255)
                                                                                                                                            gcCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATAAATTCTGGCATGACATTTTGTTTGAAAAGCAATAAGTGAGTAATGATGAAAATCCATACCACAGAGgcg       >  1:902431/1‑149 (MQ=255)
                                                                                                                                              ccTTCGCCAGTGGTGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATAAATTCTGGCATGACATTTTGTTTGAAAAGCAATAAGTGa                                     <  4:739713/117‑1 (MQ=255)
                                                                                                                                              ccTTCGCCAGTGGTGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATAAATTCTGGCATGACATTTTGTTTGAAAAGCAATAAGTGa                                     >  3:739713/1‑117 (MQ=255)
                                                                                                                                              ccTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATAAATTCTGGCATGACATTTTGTTTGAAAAGCAATAAGTGAGTAATGATGAAAATCCATAcc                >  1:633156/1‑138 (MQ=255)
                                                                                                                                              ccTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATAAATTCTGGCATGACATTTTGTTTGAAAAGCAATAAGTGAGTAATGATGAAAATCCATAcc                <  2:633156/138‑1 (MQ=255)
                                                                                                                                               cTTCGCCAGTGGTGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATAAATTCTGGCATGACAtttt                                                         <  1:174928/96‑1 (MQ=255)
                                                                                                                                               cTTCGCCAGTGGTGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATAAATTCTGGCATGACAtttt                                                         >  2:174928/1‑96 (MQ=255)
                                                                                                                                               cTTCGCCAGTGGTGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATAAATTCTGGCATGACATTTTGTTTGAAAAGCAATAAGTGAGTAATGATGAAAATCCATACCACAGAGGCGCTc    >  3:1144933/1‑149 (MQ=255)
                                                                                                                                               cTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATAAATTCTGGCCTGACATTTTGTTTGAAAAGCAATAAGTGAGTAATGATGAAAATCCATACCACCGAGGCGCTc    >  3:13322/1‑149 (MQ=255)
                                                                                                                                               cTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATAAATTCTGGCATGACATTTTGTTTGAAAAGCAATAAGTGAGTAATGATGAAAATCCATACCACAGAGGCGCTc    <  1:1111069/149‑1 (MQ=255)
                                                                                                                                                 tCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATAAATTCTGGCATGACATTTTGTTTGAAAAGCAATAAGTGAGTAATGATGAAAATCCATACCACAGAGGCGCTCAt  <  3:1069076/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
CGAAGTTCGGCTACCAGCCAAAGTAACACCACGTATTCTCCCTGGAGTTAGCGCTATGGGCCAGGGAGCCTGGCACGAGGCCAATATGTCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATAAATTCTGGCATGACATTTTGTTTGAAAAGCAATAAGTGAGTAATGATGAAAATCCATACCACAGAGGCGCTCAT  >  NC_000913/1660291‑1660584

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: