Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,660,442 C→G 11.4% 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,4420CG11.4% 466.4 / 40.8 140L792V (CTG→GTG) ynfEputative selenate reductase, periplasmic
Reads supporting (aligned to +/- strand):  ref base C (55/69);  new base G (9/7);  total (64/76)
Fisher's exact test for biased strand distribution p-value = 4.30e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01

Only 100 of 143 total aligned reads displayed.

GTTCGGCTACCAGCCAAAGTAACACCACGTATTCTCCCTGGAGTTAGCGCTATGGGCCAGGGAGCCTGGCACGAGGCCAATATGTCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATAAATTCTGGCATGACATTTTGTTTGAAAAGCAATAAGTGAGTAATGATGAAAATCCATACCACAGAGGCGCTCATGAA  >  NC_000913/1660295‑1660587
                                                                                                                                                   |                                                                                                                                                 
gTTCGGCTACCAGCCAAAGTAACACCACGTATTCTCCCTGGAGTTAGCGCTATGGGCCAGGGAGCCTGGCACGAGGCCAATATGTCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACt                                                                                                                                                  >  2:500547/1‑149 (MQ=255)
 ttCGGCTACCAGCCAAAGTAACACCACGTATTCTCCCTGGAGTTAGCGCTATGGGCCAGGGAGCCTGGCACGAGGCCAATATGTCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTg                                                                                                                                                 <  1:608991/149‑1 (MQ=255)
 ttCGGCTACCAGCCAAAGTAACACCACGTATTCTCCCTGGAGTTAGCGCTATGGGCCAGGGAGCCTGGCACGAGGCCAATATGTCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTg                                                                                                                                                 >  1:408619/1‑149 (MQ=255)
   cGGCTACCAGCCAAAGTAACACCACGTATTCTCCCTGGAGTTAGCGCTATGGGCCAGGGAGCCTGGCACGAGGCCAATATGTCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGc                                                                                                                                               <  3:707887/149‑1 (MQ=255)
    ggCTACCAGCCAAAGTAACACCACGTATTCTCCCTGGAGTTAGCGCTATGGGCCAGGGAGCCTGGCACGAGGCCAATATGTCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCg                                                                                                                                              <  4:975311/149‑1 (MQ=255)
        aCCAGCCAAAGTAACACCACGTATTCTCCCTGGAGTTAGCGCTATGGGCCAGGGAGCCTGGCACGAGGCCAATATGTCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAgg                                                                                                                                          <  4:131970/149‑1 (MQ=255)
        aCCAGCCAAAGTAACACCACGTATTCTCCCTGGAGTTAGCGCTATGGGCCAGGGAGCCTGGCACGAGGCCAATATGTCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAgg                                                                                                                                          <  1:882462/149‑1 (MQ=255)
        aCCAGCCAAAGTAACACCACGTATTCTCCCTGGAGTTAGCGCTATGGGCCAGGGAGCCTGGCACGAGGCCAATATGTCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAgg                                                                                                                                          <  4:9119/149‑1 (MQ=255)
          cAGCCAAAGTAACACCACGTATTCTCCCTGGAGTTAGCGCTATGGGCCAGGGAGCCTGGCACGAGGCCAATATGTCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGa                                                                                                                                        <  4:471036/149‑1 (MQ=255)
            gCCAAAGTAACACCACGTATTCTCCCTGGAGTTAGCGCTATGGGCCAGGGAGCCTGGCACGAGGCCAATATGTCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGaaa                                                                                                                                      >  1:1084937/1‑149 (MQ=255)
               aaaGTAACACCACGTATTCTCCCTGGAGTTAGCGCTATGGGCCAGGGAGCCTGGCACGAGGCCAATATGTCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAAccc                                                                                                                                   >  3:1162767/1‑149 (MQ=255)
                aaGTAACACCACGTATTCTCCCTGGAGTTAGCGCTATGGGCCAGGGAGCCTGGCACGAGGCCAATATGTCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCg                                                                                                                                  <  3:453118/149‑1 (MQ=255)
                aaGTAACACCACGTATTCTCCCTGGAGTTAGCGCTATGGGCCAGGGAGCCTGGCACGAGGCCAATATGTCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCg                                                                                                                                  <  1:1231653/149‑1 (MQ=255)
                   tAACACCACGTATTCTCCCTGGAGTTAGCGCTATGGGCCAGGGAGCCTGGCACGAGGCCAATATGTCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCgcag                                                                                                                               <  4:547619/149‑1 (MQ=255)
                   tAACACCACGTATTCTCCCTGGAGTTAGCGCTATGGGCCAGGGAGCCTGGCACGAGGCCAATATGTCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCgcag                                                                                                                               <  1:649556/149‑1 (MQ=255)
                       accacGTATTCTCCCTGGAGTTAGCGCTATGGGCCAGGGAGCCTGGCACGAGGCCAATATGTCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCGCCAGTGGTGAAgggc                                                                                                                                        >  3:292596/1‑135 (MQ=255)
                       accacGTATTCTCCCTGGAGTTAGCGCTATGGGCCAGGGAGCCTGGCACGAGGCCAATATGTCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGcaca                                                                                                                           <  1:817905/149‑1 (MQ=255)
                         cacGTATTCTCCCTGGAGTTAGCGCTATGGGCCAGGGAGCCTGGCACGAGGCCAATATGTCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACt                                                                                                                         >  4:881894/1‑149 (MQ=255)
                          acGTATTCTCCCTGGAGTTAGCGCTATGGGCCAGGGAGCCTGGCACGAGGCCAATATGTCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAAccc                                                                                                                                   >  4:759852/1‑138 (MQ=255)
                          acGTATTCTCCCTGGAGTTAGCGCTATGGGCCAGGGAGCCTGGCACGAGGCCAATATGTCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAAccc                                                                                                                                   <  3:759852/138‑1 (MQ=255)
                              aTTCTCCCTGGAGTTAGCGCTATGGGCCAGGGAGCCTGGCACGAGGCCAATATGTCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAgg                                                                                                                                          <  3:759162/127‑1 (MQ=255)
                              aTTCTCCCTGGAGTTAGCGCTATGGGCCAGGGAGCCTGGCACGAGGCCAATATGTCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCt                                                                                                                    >  1:1397114/1‑149 (MQ=255)
                              aTTCTCCCTGGAGTTAGCGCTATGGGCCAGGGAGCCTGGCACGAGGCCAATATGTCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCt                                                                                                                    >  4:612096/1‑149 (MQ=255)
                                  tcCCTGGAGTTAGCGCTATGGGCCAGGGAGCCTGGCACGAGGCCAATATGTCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGTACACTAATCTGGTc                                                                                                                >  2:1187161/1‑149 (MQ=255)
                                    ccTGGAGTTAGCGCTATGGGCCAGGGAGCCTGGCACGAGGCCAATATGTCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCga                                                                                                              <  1:806070/149‑1 (MQ=255)
                                               cgcTATGGGCCAGGGAGCCTGGCACGAGGCCAATATGTCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGaaaaaa                                                                                                   <  4:1162767/149‑1 (MQ=255)
                                                      ggCCAGGGAGCCTGGCACGAGGCCAATATGTCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCgaga                                                                                                            >  3:1000159/1‑133 (MQ=255)
                                                      ggCCAGGGAGCCTGGCACGAGGCCAATATGTCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGATAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCgaga                                                                                                            <  4:1000159/133‑1 (MQ=255)
                                                       gCCAGGGAGCCTGGCACGAGGCCAATATGTCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCga                                                                                                              <  2:519572/130‑1 (MQ=255)
                                                       gCCAGGGAGCCTGGCACGAGGCCAATATGTCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCga                                                                                                              >  1:519572/1‑130 (MQ=255)
                                                       gCCAGGGAGCCTGGCACGAGGCCAATATGTCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAAccc                                                                                           <  4:475535/149‑1 (MQ=255)
                                                           gggAGCCTGGCACGAGGCCAATATGTCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGa                                                                                       <  2:853634/149‑1 (MQ=255)
                                                              aGCCTGGCACGAGGCCAATATGTCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACaa                                                                                    >  4:526996/1‑149 (MQ=255)
                                                                      aCGAGGCCAATATGTCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATaa                                                                                              <  4:316289/131‑1 (MQ=255)
                                                                      aCGAGGCCAATATGTCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATAAAt                                                                            >  3:1317640/1‑149 (MQ=255)
                                                                        gAGGCCAATATGTCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAata                                                                                                 <  4:1001980/126‑1 (MQ=255)
                                                                         aGGCCAATATGTCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAAc                                                                                                                                     >  3:237682/1‑89 (MQ=255)
                                                                         aGGCCAATATGTCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAAc                                                                                                                                     <  4:237682/89‑1 (MQ=255)
                                                                           gCCAATATGTCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGaa                                                                                                       >  4:255875/1‑117 (MQ=255)
                                                                           gCCAATATGTCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAata                                                                                                 <  2:1152826/123‑1 (MQ=255)
                                                                           gCCAATATGTCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAata                                                                                                 >  1:1152826/1‑123 (MQ=255)
                                                                           gCCAATATGTCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATAAATTCTgg                                                                       >  1:1086682/1‑149 (MQ=255)
                                                                           gCCAATATATCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATAAATTCTgg                                                                       >  2:244377/1‑149 (MQ=255)
                                                                                      tGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAAc                                                                                   >  2:357470/1‑126 (MQ=255)
                                                                                         cGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATAAATTCTGGCATGACATtttgtt                                                         <  4:277604/149‑1 (MQ=255)
                                                                                              aaaTCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATAAATTCTGGCATGACATTTTGTTTGaaa                                                    >  4:1357081/1‑149 (MQ=255)
                                                                                               aaTCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAAccc                                                                                           >  4:718213/1‑109 (MQ=255)
                                                                                               aaTCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATAAATTCTGGCATGACATTTTGTTTGaaaa                                                   <  3:881894/149‑1 (MQ=255)
                                                                                                aTCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCGCCAGTGGTGAAgg                                                                                                                                          <  4:1410906/61‑1 (MQ=38)
                                                                                                aTCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCGCCAGTGGTGAAgg                                                                                                                                          >  3:1410906/1‑61 (MQ=38)
                                                                                                aTCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATAAATTCTGGCATGACATtttg                                                           <  1:1133728/140‑1 (MQ=255)
                                                                                                 tCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCgca                                                                                                                                >  3:1429775/1‑70 (MQ=255)
                                                                                                    aCCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATAAATTCTGGCATGACATTTTGTTTGAAAAGCa                                                <  3:257802/147‑1 (MQ=255)
                                                                                                    aCCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATAAATTCTGGCATGACATTTTGTTTGAAAAGCAAt                                              <  4:521721/149‑1 (MQ=255)
                                                                                                      cATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATAAATTCTGGCATGACATTTTGTTTGAAAAGCAATaa                                            <  4:1317640/149‑1 (MQ=255)
                                                                                                      cATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATAAATTCTGGCATGACATTTTGTTTGAAAAGCAATaa                                            <  3:1357081/149‑1 (MQ=255)
                                                                                                          gcggcTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATAAATTCTGGCATGACATTTTGTTTGAAAAGCAATAAGTGa                                        <  2:1086682/149‑1 (MQ=255)
                                                                                                          gcggcTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATAAATTCTGGCATGACATTTTGTTTGAAAAGCAATAAGTGa                                        <  3:526996/149‑1 (MQ=255)
                                                                                                            ggcTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGAt                                                                                                           >  2:1099915/1‑80 (MQ=255)
                                                                                                               tgtgtgAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATAAATTCTgg                                                                       <  3:26837/113‑1 (MQ=255)
                                                                                                               tgtgtgAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATAAATTCTgg                                                                       >  4:26837/1‑113 (MQ=255)
                                                                                                                gtgtgAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATAAATTCTGGCATGACATtttgtt                                                         >  4:863206/1‑126 (MQ=255)
                                                                                                                gtgtgAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATAAATTCTGGCATGACATtttgtt                                                         <  3:863206/126‑1 (MQ=255)
                                                                                                                  gtgAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATAAATTCTGGCATGACATTTTGTTTGAAAAGCAATAAGTGAGTAatgat                                >  2:935371/1‑149 (MQ=255)
                                                                                                                     aaTACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAgg                                                                                                                                          <  2:1273303/40‑1 (MQ=255)
                                                                                                                     aaTACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAgg                                                                                                                                          >  1:1273303/1‑40 (MQ=255)
                                                                                                                     aaTACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATAAATTCTGGCATGACATtttgtt                                                         <  1:1370674/121‑1 (MQ=255)
                                                                                                                     aaTACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATAAATTCTGGCATGACATtttgtt                                                         >  2:1370674/1‑121 (MQ=255)
                                                                                                                     aaTACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATAAATTCTGGCATGACATTTTGTTTGAAAAGCAATAAGTGAGTAATGATGaa                             >  4:1367195/1‑149 (MQ=255)
                                                                                                                       tACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATAAATTCTGGCATGACATTTTGTTTGAAAAGCAATAAGTGAGTAATGATGaaaa                           >  4:1017498/1‑149 (MQ=255)
                                                                                                                          gCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATAAATTCTGGCATGACATTTTGTTTGAAAAGCAATAAGTGAg                                       <  1:240647/134‑1 (MQ=255)
                                                                                                                          gCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATAAATTCTGGCATGACATTTTGTTTGAAAAGCAATAAGTGAg                                       >  2:240647/1‑134 (MQ=255)
                                                                                                                             aaCCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATAAATTCTGGCATGACATTTTGTTTGAAAAGCAATAAGTGAGTAATGATGAAAATCCATa                     <  1:259365/149‑1 (MQ=255)
                                                                                                                              aCCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATAAATTCTGGCATGACATTTTGTTTGAAAAGCAATAAGTGa                                        >  2:1315159/1‑129 (MQ=255)
                                                                                                                              aCCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATAAATTCTGGCATGACATTTTGTTTGAAAAGCAATAAGTGAGTAATGATGAAAATCCa                       >  2:346126/1‑146 (MQ=255)
                                                                                                                              aCCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATAAATTCTGGCATGACATTTTGTTTGAAAAGCAATAAGTGAGTAATGATGAAAATCCa                       <  1:346126/146‑1 (MQ=255)
                                                                                                                                cACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTc                                                                                                                <  4:1107065/55‑1 (MQ=255)
                                                                                                                                cACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTc                                                                                                                >  3:1107065/1‑55 (MQ=255)
                                                                                                                                  ctctGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATAAATTCTGGCATGACATTTTGTTTGAAAAGCAATAAGTGAGTAATGATGAAAATCCATACcaca                <  3:544953/149‑1 (MQ=255)
                                                                                                                                        tcccTTCGCCAGTGGTGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATAAATTCTGGCATGACATTTTGTTTGAAAAGCAATAAGTGAGTAATGATGAAAATCCATACCACAGAGgcg          >  4:1025306/2‑149 (MQ=255)
                                                                                                                                         cccTTCGCCAGTGGTGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATAAATTCTGGc                                                                      <  2:1005711/88‑1 (MQ=255)
                                                                                                                                         cccTTCGCCAGTGGTGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATAAATTCTGGc                                                                      >  1:1005711/1‑88 (MQ=255)
                                                                                                                                         cccTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCTCGACAACCATAAATTCTGGCATGACATTTTGTTTGAAAAGCAATAAGTGAGTaa                                    <  3:1021442/122‑1 (MQ=255)
                                                                                                                                         cccTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATAAATTCTGGCATGACATTTTGTTTGAAAAGCAATAAGTGAGTaa                                    >  4:1021442/1‑122 (MQ=255)
                                                                                                                                         cccTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATAAATTCTGGCATGACATTTTGTTTGAAAAGCAATAAGTGAGTAATGATGAAAATCCATACCACAGAGgcgc         <  1:935371/149‑1 (MQ=255)
                                                                                                                                          ccTTCGCCAGTGGTGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATAAATTCTGGCATGACATtttgtt                                                         >  4:365612/1‑100 (MQ=255)
                                                                                                                                          ccTTCGCCAGTGGTGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATAAATTCTGGCATGACATtttgtt                                                         <  3:365612/100‑1 (MQ=255)
                                                                                                                                          ccTTCGCCAGTGGTGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATAAATTCTGGCATGACATTTTGTTTGAAAAGc                                                 <  1:873963/108‑1 (MQ=255)
                                                                                                                                          ccTTCGCCAGTGGTGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATAAATTCTGGCATGACATTTTGTTTGAAAAGc                                                 >  2:873963/1‑108 (MQ=255)
                                                                                                                                          ccTTCGCCAGTGGTGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATAAATTCTGGCATGACATTTTGTTTGAAAAGCaa                                               >  4:306521/1‑110 (MQ=255)
                                                                                                                                          ccTTCGCCAGTGGTGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATAAATTCTGGCATGACATTTTGTTTGAAAAGCaa                                               <  3:306521/110‑1 (MQ=255)
                                                                                                                                          ccTTCGCCAGTGGTGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATAAATTCTGGCATGACATTTTGTTTGAAAAGCAATAAGTGAGTAatga                                 >  3:381798/1‑124 (MQ=255)
                                                                                                                                          ccTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATAAATTCTGGCATGAc                                                                 >  4:995473/1‑92 (MQ=255)
                                                                                                                                          ccTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATAAATTCTGGCATGAc                                                                 <  3:995473/92‑1 (MQ=255)
                                                                                                                                          ccTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATAAATTCTGGCATGACATTTTGTTTGAAAAGCAATAAGTGAGTAATGATGAAAATCCATACCACAGAGGCGCt        >  3:709185/1‑149 (MQ=255)
                                                                                                                                          ccTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATAAATTCTGGCATGACATTTTGTTTGAAAAGCAATAAGTGAGTAATGATGAAAATCCATACCACAGAGGCGCt        >  3:1058991/1‑149 (MQ=255)
                                                                                                                                           cTTCGCCAGTGGTGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATAAATTCTGGCATGACATTTTGTTTGAAAAGCAATAAGTGAGTAATGATGAAAATCCATACCACAGAGGCGCTc       >  2:1084917/1‑149 (MQ=255)
                                                                                                                                                ccagTGGTGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATAAATTCTGGCATGACATTTTGTTTGAAAAGCAATAAGTGAGTAATGATGAAAATCCATACCACAGAGGCGCTCATGaa  <  3:1025306/149‑1 (MQ=255)
                                                                                                                                                 cacTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATAAATTCTGGCATGACATTTTGTTTGAAAAGCAATAAGTGAGTAATGATGAAAATCCATACcaca                >  1:41418/1‑134 (MQ=255)
                                                                                                                                                 cacTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATAAATTCTGGCATGACATTTTGTTTGAAAAGCAATAAGTGAGTAATGATGAAAATCCATACcaca                <  2:41418/134‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
GTTCGGCTACCAGCCAAAGTAACACCACGTATTCTCCCTGGAGTTAGCGCTATGGGCCAGGGAGCCTGGCACGAGGCCAATATGTCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATAAATTCTGGCATGACATTTTGTTTGAAAAGCAATAAGTGAGTAATGATGAAAATCCATACCACAGAGGCGCTCATGAA  >  NC_000913/1660295‑1660587

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