Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,660,442 C→G 5.8% L792V (CTG→GTG)  ynfE → putative selenate reductase, periplasmic

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,660,4420CG 5.8% 424.3 / 12.9 133L792V (CTG→GTG) ynfEputative selenate reductase, periplasmic
Reads supporting (aligned to +/- strand):  ref base C (62/63);  new base G (4/4);  total (66/67)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.44e-01

Only 100 of 142 total aligned reads displayed.

TTCGGCTACCAGCCAAAGTAACACCACGTATTCTCCCTGGAGTTAGCGCTATGGGCCAGGGAGCCTGGCACGAGGCCAATATGTCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATAAATTCTGGCATGACATTTTGTTTGAAAAGCAATAAGTGAGTAATGATGAAAATCCATACCACAGAGGCGCTCATGAAGGCT  >  NC_000913/1660296‑1660591
                                                                                                                                                  |                                                                                                                                                     
ttCGGCTACCAGCCAAAGTAACACCACGTATTCTCCCTGGAGTTAGCGCTATGGGCCAGGGAGCCTGGCACGAGGCCAATATGTCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTgg                                                                                                                                                    <  2:1110927/150‑1 (MQ=255)
  cGGCTACCAGCCAAAGTAACACCACGTATTCTCCCTGGAGTTAGCGCTATGGGCCAGGGAGCCTGGCACGAGGCCAATATGTCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCg                                                                                                                                                  >  2:1835933/1‑150 (MQ=255)
   ggCTACCAGCCAAAGTAACACCACGTATTCTCCCTGGAGTTAGCGCTATGGGCCAGGGAGCCTGGCACGAGGCCAATATGTCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGa                                                                                                                                                 >  2:470803/1‑150 (MQ=255)
     cTACCAGCCAAAGTAACACCACGTATTCTCCCTGGAGTTAGCGCTATGGGCCAGGGCGCCTGGCACGAGGCCAATATGTCTGGCGATAAAATCGACCATGGCGGCTGTGGGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAg                                                                                                                                               >  2:1577207/1‑150 (MQ=255)
       aCCAGCCAAAGTAACACCACGTATTCTCCCTGGAGTTAGCGCTATGGGCCAGGGAGCCTGGCACGAGGCCAATATGTCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAggg                                                                                                                                             <  1:1417186/150‑1 (MQ=255)
       aCCAGCCAAAGTAACACCACGTATTCTCCCTGGAGTTAGCGCTATGGGCCAGGGAGCCGGGCACGAGGCCAATATGTCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTCTCCCCTCTGCGCCCTTCACCACTGGCGAAggg                                                                                                                                             >  1:1032058/1‑150 (MQ=255)
        ccAGCCAAAGTAACACCACGTATTCTCCCTGGAGTTAGCGCTATGGGCCAGGGAGCCTGGCACGAGGCCAATATGTCTGGCGCTAAAATCGCCCATGGCGGCGGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGATGGGa                                                                                                                                            >  1:1470172/1‑150 (MQ=255)
        ccAGCCAAAGTAACACCACGTATTCTCCCTGGAGTTAGCGCTATGGGCCAGGGAGCCTGGCACGAGGCCAATATGTCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCGCCAGTGGTGAAgggc                                                                                                                                            <  1:927535/150‑2 (MQ=255)
        ccAGCCAAAGTAACACCACGTATTCTCCCTGGAGTTAGCGCTATGGGCCAGGGAGCCTGGCACGAGGCCAATATGTCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCGCCAGTGGTGAAgggc                                                                                                                                            <  1:794845/150‑2 (MQ=255)
        ccAGCCAAAGTAACACCACGTATTCTCCCTGGAGTTAGCGCTATGGGCCAGGGAGCCTGGCACGAGGCCAATATGTCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGa                                                                                                                                            >  2:1394477/1‑150 (MQ=255)
          aGCCAAAGTAACACCACGTATTCTCCCTGGAGTTAGCGCTATGGGCCAGGGAGCCTGGCACGAGGCCAATATGTCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGaaa                                                                                                                                          >  2:1808027/1‑150 (MQ=255)
          aGCCAAAGTAACACCACGTATTCTCCCTGGAGTTAGCGCTAGGGGCCAGGGAGCCTGGCACGAGGCCAATATGTCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGaaa                                                                                                                                          <  2:1364880/150‑1 (MQ=255)
            ccAAAGTAACACCACGTATTCTCCCTGGAGTTAGCGCTATGGGCCAGGGAGCCTGGCACGAGGCCAATATGTCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAAcc                                                                                                                                        <  1:448062/150‑1 (MQ=255)
               aaGTAACACCACGTATTCTCCCTGGAGTTAGCGCTATGGGCCAGGGAGCCTGGCACGAGGCCAATATGTCTGGCGATAAAATCGGCCATGGCGGCTGTGTGAATACGCTACCCACTCTGCGCCCTGCACCCCTGGCGAAGGGACCCCCgc                                                                                                                                     >  2:381792/1‑150 (MQ=255)
                aGTAACACCACGTATTCTCCCTGGAGTTAGCGCTATGGGCCAGGGAGCCTGGCACGAGGCCAATATGTCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCgcc                                                                                                                                    >  2:163935/1‑149 (MQ=255)
                aGTAACACCACGTATTCTCCCTGGAGTTAGCGCTATGGGCCAGGGAGCCTGGCACGAGGCCAATATGTCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCgca                                                                                                                                    <  1:244987/150‑1 (MQ=255)
                 ttAACACCACGTATTCTCCCTGGAGTTAGCGCTATGGGCCAGGGAGCCTGGCACGAGGCCAATATGTCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCgcag                                                                                                                                   <  1:699698/149‑1 (MQ=255)
                    acacCACGTATTCTCCCTGGAGTTAGCGCTATGGGCCAGGGAGCCTGGCACGAGGCCAATATGTCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGcac                                                                                                                                <  2:756748/150‑1 (MQ=255)
                       ccacGTATTCTCCCTGGAGTTAGCGCTATGGGCCAGGGAGCCTGGCACGAGGCCAATATGTCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACt                                                                                                                             >  2:484140/1‑150 (MQ=255)
                       ccacGTATTCTCCCTGGAGATAGCGCTATGGGCCAGGGAGCCTGGCACGAGGCCAATATGTCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACt                                                                                                                             >  1:569145/1‑150 (MQ=255)
                        cacGTATTCTCCCTGGAGTTAGCGCTATGGGCCAGGGAGCCTGGCACGAGGCCAATATGTCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTa                                                                                                                            <  2:55693/150‑1 (MQ=255)
                             aTTCTCCCTGGAGTTAGCGCTATGGGCCAGGGAGCCTGGCACGAGGCCAATATGTCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTg                                                                                                                       >  1:78221/1‑150 (MQ=255)
                                ctcCCTGGAGTTAGCGCTATGGGCCAGGGAGCCTGGCACGAGGCCAATATGTCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTc                                                                                                                    >  1:377269/1‑150 (MQ=255)
                                ctcCCTGGAGTTAGCGCTATGGGCCAGGGAGCCTGGCACGAGGCAAATATGTCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGTGAAACCCGCAGCACACTAATCTGGTc                                                                                                                    <  2:1476974/150‑1 (MQ=255)
                                 tcCCTGGAGTTAGCGCTATGGGCCAGGGAGCCTGGCACGAGGCCAATATGTCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCg                                                                                                                   >  1:1400488/1‑150 (MQ=255)
                                 tcCCTGGAGTTAGCGCTATGGGCCAGGGAGCCTGGCACGAGGCCAATATGTCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCg                                                                                                                   >  1:81165/1‑150 (MQ=255)
                                   ccTGGAGTTAGCGCTATGGGCCAGGGAGCCTGGCACGAGGCCAATATGTCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCgag                                                                                                                 <  1:873869/150‑1 (MQ=255)
                                   ccTGGAGTTAGCGCTATGGGCCAGGGAGCCTGGCACGAGGCCAATATGTCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCgag                                                                                                                 >  1:1943728/1‑150 (MQ=255)
                                   ccTGGAGTTAGCGCTATGGGCCAGGGAGCCTGGCACGAGGCCAATATGTCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCgag                                                                                                                 <  2:914054/150‑1 (MQ=255)
                                     tGGAGTTAGCGCTATGGGCCAGGGAGCCTGGCACGAGGCCAATATGTCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGAt                                                                                                               >  2:1580157/1‑150 (MQ=255)
                                      ggAGTTAGCGCTATGGGCCAGGGAGCCTGGCACGAGGCCAATATGTCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCTCCACTGGCGAAGGGAATCCCGCAGCACACTAATCTGGTCGAGATc                                                                                                              <  1:234394/150‑1 (MQ=255)
                                      ggAGTTAGCGCTATGGGCCAGGGAGCCTGGCACGAGGCCAATATGTCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATc                                                                                                              <  1:418959/150‑1 (MQ=255)
                                       gAGTTAGCGCTATGGGCCAGGGAGCCTGGCACGATGCCAATATGTCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCg                                                                                                             <  2:44270/150‑1 (MQ=255)
                                       gAGTTAGCGCTATGGGCCAGGGAGCCTGGCACGAGGCCAATATGTCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCg                                                                                                             <  1:1056791/150‑1 (MQ=255)
                                           tAGCGCTATGGGCCAGGGAGCCTGGCACGAGGCCAATATGTCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACAACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGaaaa                                                                                                         >  1:1729418/1‑150 (MQ=255)
                                             gcgcTATGGGCGAGGGAGCCTGGCACGAGGCCAATATGTCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGaaaaaa                                                                                                       <  2:527304/150‑1 (MQ=255)
                                               gcTATGGGCCAGGGAGCCTGGCACGAGGCCAATATGTCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCCCCACTGGCGACGGGAAACCCGCAGCGCACTAATCTGGTCGAGATAGAAAAAata                                                                                                     >  1:209093/1‑150 (MQ=255)
                                                      gCCAGGGAGCCTGGCACGAGGCCAATATGTCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATATCCCa                                                                                              >  2:124745/1‑150 (MQ=255)
                                                            gAGCCTGGCACGAGGCCAATATGTCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATACCCCACGACaa                                                                                        >  2:1384362/1‑150 (MQ=255)
                                                            gAGCCTGGCACGAGGCCAATATGTCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACaa                                                                                        >  1:584382/1‑150 (MQ=255)
                                                                 tGGCACGAGGCCAATATGTCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATa                                                                                   <  2:931177/150‑1 (MQ=255)
                                                                    cACGAGGCCAATATGTCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATAAAt                                                                                >  1:697093/1‑150 (MQ=255)
                                                                       gAGGCCAATATGTCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATAAATTCt                                                                             >  1:1970168/1‑150 (MQ=255)
                                                                        acgCCAATATGTCTGGCGATAAAATCGACCATTGCGGCTGTGTGAATACTCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTATTCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATAAATTCTg                                                                            <  1:277508/148‑1 (MQ=255)
                                                                         ggCCAATATGTCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATAAATTCTgg                                                                           >  2:1714451/1‑150 (MQ=255)
                                                                          gCCAATATGTCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATAAATTCTGGc                                                                          <  1:1033807/150‑1 (MQ=255)
                                                                          gCCAATATGTCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATAAATTCTGGc                                                                          <  2:1046542/150‑1 (MQ=255)
                                                                          gCCAATATGTCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATAAATTCTGGc                                                                          <  1:484140/150‑1 (MQ=255)
                                                                            cAATATGTCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATAAATTCTGGCAt                                                                        <  1:611078/150‑1 (MQ=255)
                                                                             aaTATGTCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATAAATTCTGGCATg                                                                       <  1:1808028/150‑1 (MQ=255)
                                                                              atatGTCTGGCGATAAAATCGTCCATGGTGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATAAATTCTGGCATGa                                                                      <  2:266494/150‑1 (MQ=255)
                                                                                 tGTCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATAAATTCTGGCATGACAt                                                                   >  2:779186/1‑150 (MQ=255)
                                                                                    cTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAATCCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAGCCCACGCCAACCATAAATTCTGGAATGACAtttt                                                                >  1:1283337/1‑150 (MQ=255)
                                                                                    cTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATAAATTCTGGCATGACAtttt                                                                <  2:1892899/150‑1 (MQ=255)
                                                                                     tGGCGATAAAATCGACCATGGCGTCTGGGTGAATACTCTAACCACTCTGCTCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATAAATTCTGGCCTGACATtttg                                                               <  1:1340431/150‑1 (MQ=255)
                                                                                       gCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATAAATTCTGGCATGACATtttgtt                                                             <  1:1621754/150‑1 (MQ=255)
                                                                                        cGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCAGGGGCGAAGCGAAACCCGCAGCACACTAAGCGGGTCGAGAACGAAAAAATATAACCCACGACAACAATAATTTCGGGCAAGACATtttgttt                                                            >  2:1846068/1‑150 (MQ=255)
                                                                                        cGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATAAATTCTGGCATGACATtttgttt                                                            <  1:743821/150‑1 (MQ=255)
                                                                                         taTAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATAAATTCTGGCATGACATtttgtttg                                                           <  2:2029102/149‑1 (MQ=255)
                                                                                         taTAAAATCGACCATGGCGGCTGTGTGAATACGCAAACCTCTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACTACCATAAATTCTGGCATGACATtttgtttg                                                           <  2:1361867/149‑1 (MQ=255)
                                                                                         gATAAAATCGAACATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATAAATTCTGGCATGACATtttgtttg                                                           <  1:333670/150‑1 (MQ=255)
                                                                                          aTAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATAAATTCTGGCATGACATTTTGTTTGa                                                          >  1:1513407/1‑150 (MQ=255)
                                                                                          aTAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATAAATTCTGGCATGACATTTTGTTTGa                                                          <  1:470803/150‑1 (MQ=255)
                                                                                            aaaaTCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATAAATTCTGGCATGACATTTTGTTTGaaa                                                        <  1:779411/150‑1 (MQ=255)
                                                                                            aaaaTCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATAAATTCTGGCATGACATTTTGTTTGaaa                                                        >  1:1748117/1‑150 (MQ=255)
                                                                                               aTCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATAAATTCTGGCATGACATTTTGTTTGAAAAGc                                                     <  2:774535/150‑1 (MQ=255)
                                                                                               aTCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATAAATTCTGGCATGACATTTTGTTTGAAAAGc                                                     <  1:1611406/150‑1 (MQ=255)
                                                                                                 cGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATAAATTCTGGCATGACATTTTGTTTGAAAAGCaa                                                   >  1:1839397/1‑150 (MQ=255)
                                                                                                  gACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATAAATTCTGGCATGACATTTTGTTTGAAAAGCAAt                                                  <  2:377269/150‑1 (MQ=255)
                                                                                                  gACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATAAATTCTGGCATGACATTTTGTTTGAAAAGCAAt                                                  <  2:1771007/150‑1 (MQ=255)
                                                                                                  gACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATAAATTCTGGCATGACATTTTGTTTGAAAAGCAAt                                                  >  1:636382/1‑150 (MQ=255)
                                                                                                         gcggcTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATAAATTCTGGCATGACATTTTGTTTGAAAAGCAATAAGTGAg                                           <  2:1400487/150‑1 (MQ=255)
                                                                                                          cggcTGTGTTAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCTCACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATAAATTCTGGCATGACATTTTGTTTGAAAAGCAATAAGTGAGt                                          <  2:265463/150‑1 (MQ=255)
                                                                                                          cggcTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATAAATTCTGGCATGACATTTTGTTTGAAAAGCAATAAGTGAGt                                          >  1:648590/1‑150 (MQ=255)
                                                                                                           ggcTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATAAATTCTGGCATGACATTTTGTTTGAAAAGCAATAAGTGAGTa                                         >  1:1477565/1‑150 (MQ=255)
                                                                                                                tgtgAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATAAATTCTGGCATGACATTTTGTTTGAAAAGCAATAAGTGAGTAatgat                                    >  2:736198/1‑150 (MQ=255)
                                                                                                                tgtgAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATAAATTCTGGCATGACATTTTGTTTGAAAAGCAATAAGTGAGTAatgat                                    >  2:1037492/1‑150 (MQ=255)
                                                                                                                   gAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATAAATTCTGGCATGACATTTTGTTTGAAAAGCAATAAGTGAGTAATGATGaa                                 >  2:1653408/1‑150 (MQ=255)
                                                                                                                   gAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATAAATTCTGGCATGACATTTTGTTTGAAAAGCAATAAGTGAGTAATGATGaa                                 <  1:1107562/150‑1 (MQ=255)
                                                                                                                       aCGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTAGAGATCGAAAAAATATAACCCACGACAACCATAAATTCTGGCATGACATTTTGTTTGAAAAGCAATAAGTGAGTAATGATGAAAATc                             <  1:342215/150‑1 (MQ=255)
                                                                                                                        cGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCAACGACAACCATAAATTCTGGCTGGACATTTTGTTTGAAAAGCAATAAGTGAGTAATGATGAAAATcc                            >  1:862846/1‑150 (MQ=255)
                                                                                                                         gCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATAAATTCTGGCATGACATTTTGTTTGAAAAGCAATAAGTGAGTAATGATGAAAATCCa                           <  1:449057/150‑1 (MQ=255)
                                                                                                                            aaCCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATAAATTCTGGCATGACATTTTGTTTGAAAAGCAATAAGTGAGTAATGATGAAAATCCATAc                        >  2:957523/1‑150 (MQ=255)
                                                                                                                            aaCCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATAAATTCTGGCATGACATTTTGTTTGAAAAGCAATAAGTGAGTAATGATGAAAATCCATAc                        >  2:919574/1‑150 (MQ=255)
                                                                                                                             aCCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATAAATTCTGGCATGACATTTTGTTTGAAAAGCAATAAGTGAGTAATGATGAAAATCCATAcc                       <  1:957524/150‑1 (MQ=255)
                                                                                                                               cACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATAAATTCTGGCATGACATTTTGTTTGAAAAGCAATAAGTGAGTAATGATGAAAATCCATACcac                     <  1:628354/150‑1 (MQ=255)
                                                                                                                               cACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATAAATTCTGGCATGACATTTTGTTTGAAAAGCAATAAGTGAGTAATGATGAAAATCCATACcac                     >  1:20503/1‑150 (MQ=255)
                                                                                                                                  tctTCGCCCTTCACCACTTGCGAAGGGAAACCAGCAGCACAATAATCTGGTCGAGATCGAAATAATATAACCGACGACAACCATAAATTCTGGCATGACATTTTGTTTGAACAGCAATAAGTGAGTAATGATGAAAATCCATACCACaga                  <  1:378256/150‑1 (MQ=255)
                                                                                                                                    tGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATAAATTCTGGCATGACATTTTGTTTGAAAAGCAATAAGTGAGTAATGATGAAAATCCATACCACAGAgg                <  2:1934694/150‑1 (MQ=255)
                                                                                                                                        cccTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATAAATTCTGGCATGACATTTTGTTTGAAAAGCAATAAGTGAGTAATGATGAAAATCCATACCACAGAGGCGCt            >  2:1586764/1‑150 (MQ=255)
                                                                                                                                         ccTTCGCCAGTGGTGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACACACGACAACCATAAATTCTGGCATGACATTTTGTTTCAAAAGAAATAAGTGAGTAATGATGAAAATCAATACCACAGAGGCGCta           >  1:648979/1‑149 (MQ=255)
                                                                                                                                          cTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATAAATTCTGGCATGACATTTTGTTTGAAAAGCAATAAGTGAGTAATGATGAAAATCCATACCACAGAGGCGCTCa          >  1:1379832/1‑150 (MQ=255)
                                                                                                                                           ttCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATAAATTCTGGCATGACATTTTGTTTGAAAAGCAATAAGTGAGTAATGATGAAAATCCATACCACAGAGGCGCTCAt         <  1:2000672/150‑1 (MQ=255)
                                                                                                                                            tCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATAAATTCTGGCATGACATTTTGTTTGAAAAGCAATAAGTGAGTAATGATGAAAATCCATACCACAGAGGCGCTCATg        <  2:1970167/150‑1 (MQ=255)
                                                                                                                                              accacTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATAAATTCTGGCATGACATTTTGTTTGAAAAGCAATAAGTGTGTAATGATGAAAATCCATACCACAGAGGCGCTCATGaa      <  2:1467165/150‑1 (MQ=255)
                                                                                                                                               ccacTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATAAATTCTGGCATGACATTTTGTTTGAAAAGCAATAAGTGAGTAATGATGAAAATCCATACCACAGAGGCGCTCATGAAg     >  1:941806/1‑150 (MQ=255)
                                                                                                                                               ccacTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATAAATTCTGGCATGACATTTTGTTTGAAAAGCAATAAGTGAGTAATGATGAAAATCCATACCACAGAGGCGCTCATGAAg     >  2:1818700/1‑150 (MQ=255)
                                                                                                                                               ccacTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTAGAGATCGAAAAAATATAACCCACGACAACCATAAATTCTGGCATGACATTTTGTTTGAAAAGCAATAAGTGAGTAATGATGAAAATCCATACCACAGAGGCGCTCATGAAg     <  2:113584/150‑1 (MQ=255)
                                                                                                                                                aacTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATAAATTCTGGCATGACATTTTGTTTGAAAAGCAATAAGTGAGTAATGATGAAAATCAATACCACATATGAGCTCATGAAgg    >  1:1389350/2‑150 (MQ=255)
                                                                                                                                                  cTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATAAATTCTGGCATGACATTTTGTTTGAAAAGCAATAAGTGAGTAATGATGTAAATCCATACCACAGAGGCGCTCATGAAGGCt  <  2:78221/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                     
TTCGGCTACCAGCCAAAGTAACACCACGTATTCTCCCTGGAGTTAGCGCTATGGGCCAGGGAGCCTGGCACGAGGCCAATATGTCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATAAATTCTGGCATGACATTTTGTTTGAAAAGCAATAAGTGAGTAATGATGAAAATCCATACCACAGAGGCGCTCATGAAGGCT  >  NC_000913/1660296‑1660591

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: