Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,660,442 C→G 6.2% 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 6.2% 394.1 / 13.3 113L792V (CTG→GTG) ynfEputative selenate reductase, periplasmic
Reads supporting (aligned to +/- strand):  ref base C (55/51);  new base G (5/2);  total (60/53)
Fisher's exact test for biased strand distribution p-value = 4.45e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.10e-01

Only 100 of 114 total aligned reads displayed.

TTCGGCTACCAGCCAAAGTAACACCACGTATTCTCCCTGGAGTTAGCGCTATGGGCCAGGGAGCCTGGCACGAGGCCAATATGTCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATAAATTCTGGCATGACATTTTGTTTGAAAAGCAATAAGTGAGTAATGATGAAAATCCATACCACAGAGGCGCTCATGAA  >  NC_000913/1660296‑1660587
                                                                                                                                                  |                                                                                                                                                 
ttCGGCTACCAGCCAAAGTAACACCACGTATTCTCCCTGGAGTTAGCGCTATGGGCCAGGGAGCCTGGCACGAGGCCAATATGTCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACt                                                                                                                                                  >  1:298476/1‑148 (MQ=255)
   ggCTACCAGCCAAAGTAACACCACGTATTCTCCCTGGAGTTAGCGCTATGGGCCAGGGAGCCTGGCACGAGGCCAATATGTCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGc                                                                                                                                               >  1:110468/1‑148 (MQ=255)
      tACCAGCCAAAGTAACACCACGTATTCTCCCTGGAGTTAGCGCTATGGGCCAGGGAGCCTGGCACGAGGCCAATATGTCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAg                                                                                                                                           <  4:866412/149‑1 (MQ=255)
       aCCAGCCAAAGTAACACCACGTATTCTCCCTGGAGTTAGCGCTATGGGCCAGGGAGCCTGGCACGAGGCCAATATGTCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAgg                                                                                                                                          >  4:258387/1‑149 (MQ=255)
        ccAGCCAAAGTAACACCACGTATTCTCCCTGGAGTTAGCGCTATGGGCCAGGGAGCCTGGCACGAGGCCAATATGTCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAggg                                                                                                                                         >  4:642546/1‑149 (MQ=255)
           gCCAAAGTAACACCACGTATTCTCCCTGGAGTTAGCGCTATGGGCCAGGGAGCCTGGCACGAGGCCAATATGTCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGaaa                                                                                                                                      >  3:99575/1‑149 (MQ=255)
             cAAAGTAACACCACGTATTCTCCCTGGAGTTAGCGCTATGGGCCAGGGAGCCTGGCACGAGGCCAATATGTCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTGTGCGCCCTTCACCACTGGCGAAGGGAAAcc                                                                                                                                    >  3:776447/1‑149 (MQ=255)
              aaaGTAACACCACGTATTCTCCCTGGAGTTAGCGCTATGGGCCAGGGAGCCTGGCACGAGGCCAATATGTCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCGCCAGTGGTGAAgggcgcaga                                                                                                                                   <  3:329156/149‑7 (MQ=255)
                     caccacGTATTCTCCCTGGAGTTAGCGCTATGGGCCAGGGAGCCTGGCACGAGGCCAATATGTCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCgcagca                                                                                                                             <  2:738722/148‑1 (MQ=255)
                      accacGTATTCTCCCTGGAGTTAGCGCTATGGGCCAGGGAGCCTGGCACGAGGCCAATATGTCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGcaca                                                                                                                           <  3:866163/149‑1 (MQ=255)
                      accacGTATTCTCCCTGGAGTTAGCGCTATGGGCCAGGGAGCCTGGCACGAGGCCAATATGTCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGcac                                                                                                                            <  1:674170/148‑1 (MQ=255)
                         acGTATTCTCCCTGGAGTTAGCGCTATGGGCCAGGGAGCCTGGCACGAGGCCAATATGTCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACt                                                                                                                         >  2:889134/1‑148 (MQ=255)
                                ctcCCTGGAGTTAGCGCTATGGGCCAGGGAGCCTGGCACGAGGCCAATATGTCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGt                                                                                                                 >  3:732624/1‑149 (MQ=255)
                                 tcCCTGGAGTTAGCGCTATGGGCCAGGGAGCCTGGCACGAGGCCAATATGTCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTc                                                                                                                >  4:380823/1‑149 (MQ=255)
                                 tcCCTGGAGGTAGCGCTATGGGCCAGGGAGCCTGGCACGAGGCCAATATGTCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTc                                                                                                                <  3:753139/149‑1 (MQ=255)
                                  cccTGGAGTTAGCGCTATGGGCCAGGGAGCCTGGCACGAGGCCAATATGTCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGa                                                                                                                                             >  2:289844/1‑119 (MQ=255)
                                  cccTGGAGTTAGCGCTATGGGCCAGGGAGCCTGGCACGAGGCCAATATGTCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGa                                                                                                                                             <  1:289844/119‑1 (MQ=255)
                                   ccTGGAGTTAGCGCTATGGGCCAGGGAGCCTGGCACGAGGCCAATATGTCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCga                                                                                                              <  3:285487/149‑1 (MQ=255)
                                   ccTGGAGTTAGCGCTATGGGCCAGGGAGCCTGGCACGAGGCCAATATGTCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCga                                                                                                              <  4:838736/149‑1 (MQ=255)
                                    cTGGAGTTAGCGCTATGGGCCAGGGAGCCTGGCACGAGGCCAATATGTCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCga                                                                                                              >  1:655349/1‑148 (MQ=255)
                                     tGGAGTTAGCGCTATGGGCCAGGGAGCCTGGCACGAGGCCAATATGTCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCgag                                                                                                             >  1:65366/1‑148 (MQ=255)
                                     tGGAGTTAGCGCTATGGGCCAGGGAGCCTGGCACGAGGCCAATATGTCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCgag                                                                                                             >  1:953084/1‑148 (MQ=255)
                                     tGGAGTTAGCGCTATGGGCCAGGGAGCCTGGCACGAGGCCAATATGTCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCgag                                                                                                             >  2:413331/1‑148 (MQ=255)
                                      ggAGTTAGCGCTATGGGCCAGGGAGCCTGGCACGAGGCCAATATGTCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGAt                                                                                                           <  4:507387/149‑1 (MQ=255)
                                       gAGTTAGCGCTATGGGCCAGGGAGCCTGGCACGAGGCCAATATGTCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGAt                                                                                                           <  1:611907/148‑1 (MQ=255)
                                         gTTAGCGCTATGGGCCAGGGAGCCTGGCACGAGGCCAATATGTCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGa                                                                                                        >  4:121504/1‑149 (MQ=255)
                                             gcgcTATGGGCCAGGGAGCCTGGCACGAGGCCAATATGTCTGGCGATAAGATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGaaaa                                                                                                     >  2:474301/1‑148 (MQ=255)
                                              cgcTATGGGCCAGGGAGCCTGGCACGAGGCCAATATGTCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGaaaaa                                                                                                    >  1:614729/1‑148 (MQ=255)
                                                cTATGGGCCAGGGAGCCTGGCACAAGGCCAATATGTCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAat                                                                                                  <  2:896692/148‑1 (MQ=255)
                                                   tGGGCCAGGGAGCCTGGCACGAGGCCAATATGTCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATaa                                                                                              <  4:470958/149‑1 (MQ=255)
                                                     ggCCAGGGAGCCTGGCACGAGGCCAATATGTCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAAcc                                                                                            >  4:774829/1‑149 (MQ=255)
                                                      gCCAGGGAGCCTGGCACGAGGCCAATATGTCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAAccc                                                                                           <  3:704475/149‑1 (MQ=255)
                                                         aGGGAGCCTGGCACGAGGCCAATATGTCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACg                                                                                        >  3:462776/1‑149 (MQ=255)
                                                         aGGGAGCCTGGCACGAGGCCAATATGTCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACg                                                                                        >  3:515909/1‑149 (MQ=255)
                                                          gggAGCCTGGCACGAGGCCAATATGTCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACg                                                                                        <  1:889134/148‑1 (MQ=255)
                                                          gggAGCCTGGCACGAGGCCAATATGTCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGa                                                                                       <  4:732624/149‑1 (MQ=255)
                                                           ggAGCCTGGCACGAGGCCAATATGTCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGa                                                                                       >  2:544422/1‑148 (MQ=255)
                                                             aGCCTGGCACGAGGCCAATATGTCTGACGATACACTCGCCCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTTCTCTGGTCGAGATTGAAAAAATATAACCCACGACaa                                                                                    <  4:701301/149‑1 (MQ=255)
                                                                cTGGCACGAGGCCAATATGTCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCa                                                                                 >  4:543083/1‑149 (MQ=255)
                                                                 tGGCACGAGGCCAATATGTCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCa                                                                                 >  2:792660/1‑148 (MQ=255)
                                                                        aGGCCAATATGTCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATAAATTCt                                                                         <  4:769433/149‑1 (MQ=255)
                                                                         ggCCAATATGTCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATAAATTCt                                                                         <  1:470470/148‑1 (MQ=255)
                                                                          gCCAATATGTCTGGTGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAGTCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATAAATTCTg                                                                        <  2:968401/148‑1 (MQ=255)
                                                                          gCCAATATGTCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTAGCGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATAAATTCTgg                                                                       <  4:734477/149‑1 (MQ=255)
                                                                           ccAATATGTCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATAAATTCTGGc                                                                      <  4:322569/149‑1 (MQ=255)
                                                                            cAATATGTCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATAAATTCTGGc                                                                      >  2:287506/1‑148 (MQ=255)
                                                                               tatGTCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATAAATTCTGGCATg                                                                   >  1:563680/1‑148 (MQ=255)
                                                                                atGTCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATAAATTCTGGCATGa                                                                  >  2:404464/1‑148 (MQ=255)
                                                                                    cTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATAAATTCTGGCATGACAtt                                                              <  1:663497/148‑1 (MQ=255)
                                                                                     tGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATAAATTCTg                                                                        <  3:123578/137‑1 (MQ=255)
                                                                                     tGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATAAATTCTg                                                                        >  4:123578/1‑137 (MQ=255)
                                                                                      ggCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATAAATTCTGGCATGACAtttt                                                            >  1:571740/1‑148 (MQ=255)
                                                                                         gacaagaTCGCCGATGCCGGCTGTGTCAATATGCTTACCACTATGCGCCCTCCACCACTGCCGACGGGAAACCCGCAGCACACGAATCTGATCGAGATCGATAAAATATAACCCACAATACCCGTTAATTCTGGCATGACATtttgttt                                                        >  4:599393/7‑149 (MQ=255)
                                                                                         gATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATAAATTCTGGCATGACATtttgttt                                                        >  3:607498/1‑149 (MQ=255)
                                                                                         gATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATAAATTCTGGCATGACATtttgtt                                                         <  1:493594/148‑1 (MQ=255)
                                                                                             aaaTCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATAAATTCTGGCATGACATtttg                                                           <  1:846532/142‑1 (MQ=255)
                                                                                             aaaTCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATAAATTCTGGCATGACATtttg                                                           >  2:846532/1‑142 (MQ=255)
                                                                                              aaTCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATAAATTCTGGCATGACATTTTGTTTGaaaa                                                   <  4:930824/149‑1 (MQ=255)
                                                                                               aTCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATAAATTCTGGCATGACATTTTGTTTGAAAAg                                                  <  3:361184/149‑1 (MQ=255)
                                                                                                tCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATAAATTCTGGCATGACATTTTGTTTGAAAAg                                                  <  1:819996/148‑1 (MQ=255)
                                                                                                       tGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATAAATTCTGGCATGACATTTTGTTTGAAAAGCAATAAg                                           <  2:864409/148‑1 (MQ=255)
                                                                                                       tGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATAAATTCTGGCATGACATTTTGTTTGAAAAGCAATAAGt                                          >  3:167531/1‑149 (MQ=255)
                                                                                                       tGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATAAATTCTGGCATGACATTTTGTTTGAAAAGCAATAAGt                                          <  4:201155/149‑1 (MQ=255)
                                                                                                         gcggcTGTGTGAATACGCTAACCTCTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATAAATTCTGGCATGACATTTTGTTTGAAAAGCAATAAGTGa                                        <  3:2283/149‑1 (MQ=255)
                                                                                                         gcggcTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATAAATTCTGGCATGACATTTTGTTTGAAAAGCAATAAGTg                                         <  2:563680/148‑1 (MQ=255)
                                                                                                         gcggcTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATAAATTCTGGCATGACATTTTGTTTGAAAAGCAATAAGTg                                         <  1:571218/148‑1 (MQ=255)
                                                                                                         gcggcTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATAAATTCTGGCATGACATTTTGTTTGAAAAGCAATAAGTg                                         <  2:571740/148‑1 (MQ=255)
                                                                                                               gtgtgAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATAAATTCTGGCATGACATTTTGTTTGAAAAGCAATAAGTGAGTAat                                   >  1:497668/1‑148 (MQ=255)
                                                                                                                tgtgAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATAAATTCTGGCATGACATTTTGTTTGAAAAGCAATAAGTGAGTAatga                                 <  3:923554/149‑1 (MQ=255)
                                                                                                                tgtgAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATAAATTCTGGCATGACATTTTGTTTGAAAAGCAATAAGTGAGTAatg                                  >  1:373073/1‑148 (MQ=255)
                                                                                                                  tgAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATAAATTCTGGCATGACATTTTGTTTGAAAAGCAATAAGTGAGTAatgatg                               >  4:802262/1‑149 (MQ=255)
                                                                                                                  tgAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATAAATTCTGGCATGACATTTTGTTTGAAAAGCAATAAGTGAGTAatgatg                               <  3:121504/149‑1 (MQ=255)
                                                                                                                   gAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATAAATTCTGGCATGACATTTTGTTTGAAAAGCAATAAGTGAGTAatgatg                               >  2:14137/1‑148 (MQ=255)
                                                                                                                       aCGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATAAATTCTGGCATGACATTTTGTTTGAAAAGCAATAAGTg                                         <  1:280015/134‑1 (MQ=255)
                                                                                                                       aCGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATAAATTCTGGCATGACATTTTGTTTGAAAAGCAATAAGTGAGTAATGATGaaaa                           >  1:863008/1‑148 (MQ=255)
                                                                                                                        cGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATAAATTCTGGCATGACATTTTGTTTGAAAAGCAATAAGTGAGTAATGATGAAAAt                          >  1:441895/1‑148 (MQ=255)
                                                                                                                          cTAACCACTCTGCGCCCTTCACCACTGGCAAAGGGAAACCCGCAGCACACTAAGCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATAAATTCTGGCATGACATTTTGTTTGAAAAGCAATAAGTGAGTAATGATGAAAATCCa                       <  4:856357/149‑1 (MQ=255)
                                                                                                                           tAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATAAATTCTGGCATGACATTTTGTTTGAAAAGCAATAAGTGAGTAATGATGAAAATcta                       >  2:652818/1‑146 (MQ=255)
                                                                                                                           tAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATAAATTCTGGCATGACATTTTGTTTGAAAAGCAATAAGTGAGTAATGATGAAAATCCa                       >  1:298195/1‑148 (MQ=255)
                                                                                                                             aCCACTCTGCGCCCTTCACCACTGGCGAGGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATAAATTCTGGCATGACATTTTGTTTGAAAAGCAATAAGTGAGTAATGATGAAAATCCATAc                    >  3:72581/1‑149 (MQ=255)
                                                                                                                             aCCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATAAATTCTGGCATGACATTTTGTTTGAAAAGCAATAAGTGAGTAATGATGAAAATCCATa                     >  2:244678/1‑148 (MQ=255)
                                                                                                                             aCCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATAAATTCTGGCATGACATTTTGTTTGAAAAGCAATAAGTGAGTAATGATGAAAATCCATAc                    <  3:642546/149‑1 (MQ=255)
                                                                                                                                aCTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGt                                                                                                                 >  2:463066/1‑53 (MQ=255)
                                                                                                                                aCTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATAAATTCTGGCATGACATTTTGTTTGAAAAGCAATAAGTGAGTAATGATGAAAATCCATACca                  >  1:126401/1‑148 (MQ=255)
                                                                                                                                   ctGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATAAATTCTGGCATGATATTTTGTTTGAAAAGCAATAAGTGAGTAATGATGAAAATCCATACCACag               <  2:829591/148‑1 (MQ=255)
                                                                                                                                   ctGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATAAATTCTGGCATGACATTTTGTTTGAAAAGCAATAAGTGAGTAATGATGAAAATCCATACCACag               <  2:373073/148‑1 (MQ=255)
                                                                                                                                    tGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATAAATTCTGGCATGACATTTTGTTTGAAAAGCAATAAGTGAGTAATGATGAAAATCCATACCACaga              <  2:186849/148‑1 (MQ=255)
                                                                                                                                     gcgcCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATAAATTCTGGCATGACATTTTGTTTGAAAAGCAATAAGTGAGTAATGATGAAAATCCATACCACagag             <  1:558356/148‑1 (MQ=255)
                                                                                                                                        cccTTCGCCAGTGGTGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATAAATTCTGGCATGACATTTTGTTTGAAAAGCAATAAGTGAGTAATGATGAAAATCCATACCACAGAGgcgc         >  3:165979/1‑149 (MQ=255)
                                                                                                                                        cccTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATAAATTCTGGCATGACATTTTGTTTGAAAAGCAATAAGTGAGTAATGATGAAAATCCATACCACAGAGgcgc         >  3:464814/1‑149 (MQ=255)
                                                                                                                                         ccTTCGCCAGTGGTGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGCTCGAAAAAATATAACCCACGACAACCATAAATTCTGGCCTGACATTTTGTTTGAAAAGCAATAAGTGAGTAATGATGAAAATCCATACCACAGAGgcgc         >  2:299439/1‑148 (MQ=255)
                                                                                                                                         ccTTCGCCAGTGGTGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATAAATTCTGGCATGACATtttgt                                                          >  3:847222/1‑99 (MQ=255)
                                                                                                                                         ccTTCGCCAGTGGTGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATAAATTCTGGCATGACATtttgt                                                          <  4:847222/99‑1 (MQ=255)
                                                                                                                                         ccTTCGCCAGTGGTGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATAAATTCTGGCATGACATTTTGTTTGAAAAGCAATAAGTGAGTAATGATGAAAATCCATACCACAGAGGCGCt        >  4:500904/1‑149 (MQ=255)
                                                                                                                                         ccTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATAAATTCTGGCATGACATTTTGTTTGAAAAGCAATAAGTGAGTAATGATGAAAATCCATACCACAGAGGCGCt        >  3:141317/1‑149 (MQ=255)
                                                                                                                                          cTTCGCCAGTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATAAATTCTGGCATGACATTTTGTTTGAAAAGCAATAAGTTAGTAATGATGAAAATCCATACCACAGAGGCGCt        >  2:654746/1‑148 (MQ=255)
                                                                                                                                          cTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATAAATTCTGGCATGACATTTTGTTTGAAAAGCAATAAGTGAGTAATGATGAAAATCCATACCGCAGAGGCGCTc       <  3:543083/149‑1 (MQ=255)
                                                                                                                                          cTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATAAATTCTGGCATGACATTTTGTTTGAAAAGCAATAAGTGAGTAATGATGAAAATCCATACCACAGAGGCGCTc       >  3:481233/1‑149 (MQ=255)
                                                                                                                                            tCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATAAATTCTGGCATGACATTTTGTTTGAAAAGCAATAAGTGAGTAATGATGAAAATCCATACCACAGAGGCGCTCa      <  2:945017/148‑1 (MQ=255)
                                                                                                                                                cacTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATAAATTCTGGCATGACATTTTGTTTGAAAAGCAATAAGTGAGTAATGATGAAAATCCATACCACAGAGGCGCTCATGaa  <  2:858464/148‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                 
TTCGGCTACCAGCCAAAGTAACACCACGTATTCTCCCTGGAGTTAGCGCTATGGGCCAGGGAGCCTGGCACGAGGCCAATATGTCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATAAATTCTGGCATGACATTTTGTTTGAAAAGCAATAAGTGAGTAATGATGAAAATCCATACCACAGAGGCGCTCATGAA  >  NC_000913/1660296‑1660587

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