Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,970,813 Δ1 bp 8.2% coding (681/1131 nt) wecB → UDP‑N‑acetyl glucosamine‑2‑epimerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,970,8130C. 8.2% 431.8 / 32.2 122coding (681/1131 nt)wecBUDP‑N‑acetyl glucosamine‑2‑epimerase
Reads supporting (aligned to +/- strand):  ref base C (59/53);  new base . (5/5);  total (64/58)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.32e-01

Only 100 of 128 total aligned reads displayed.

GGTGCGTGACCAGGTGATGAGCAGCGACAAGCTGCGTTCAGAACTGGCGGCAAATTACCCGTTTATCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGCTGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGAAAAATGTCATTCTGATCGATCCCCAGGAGTATTTACCGTTT  >  NC_000913/3970669‑3970960
                                                                                                                                                |                                                                                                                                                    
ggTGCGTGACCAGGTGATGAGCAGCGACAAGCTGCGTTCAGAACTGGCGGCAAATTACCCGTTTATCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCAcgcg                                                                                                                                                   >  3:119078/1‑148 (MQ=255)
    cGTGACCAGGTGATGAGCAGCGACAAGCTGCGTTCAGAACTGGCGGCAAATTACCCGTTTATCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGCTgg                                                                                                                                              >  2:763627/1‑148 (MQ=255)
             gTGATGAGCAGCGACAAGCTGCGTTCAGAACTGGCGGCAAATTACCCGTTTATCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGCTGGCAGACATCg                                                                                                                                     <  1:763627/148‑1 (MQ=255)
               gatgaGCAGCGACAAGCTGCGTTCAGAACTGGCGGCAAATTACCCGTTTATCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGCTGGCAGACATCGcc                                                                                                                                   <  4:119078/148‑1 (MQ=255)
                               cTGCGTTCAGAACTGGCGGCAAATTACCCGTTTATCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGCTGGCAGACATCGCCACCACGCACCAgg                                                                                                                      >  3:1009642/1‑145 (MQ=255)
                               cTGCGTTCAGAACTGGCGGCAAATTACCCGTTTATCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGCTGGCAGACATCGCCACCACGCACCAgg                                                                                                                      <  4:1009642/145‑1 (MQ=255)
                                 gCGTTCAGAACTGGCGGCAAATTACCCGTTTATCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGCTGGCAGAc                                                                                                                                         <  2:183943/124‑1 (MQ=255)
                                 gCGTTCAGAACTGGCGGCAAATTACCCGTTTATCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGCTGGCAGAc                                                                                                                                         >  1:183943/1‑124 (MQ=255)
                                  cGTTCAGAACTGGCGGCAAATTACCCGTTTATCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGCTGGCAGACATCGCCACCACGCACCAGGACATcc                                                                                                                <  1:1074412/148‑1 (MQ=255)
                                        gAACTGGCGGCAAATTACCCGTTTATCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGCTGGCAGACATCGCCACCACGCACCAGGACATCCAGATTg                                                                                                          <  3:411070/148‑1 (MQ=255)
                                          aCTGGCGGCAAATTACCCGTTTATCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGCTGGCAGACAt                                                                                                                                       <  1:375194/117‑1 (MQ=255)
                                             ggcggcAAATTACCCGTTTATCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGCTg                                                                                                                                               <  2:963606/106‑1 (MQ=255)
                                             ggcggcAAATTACCCGTTTATCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGCTg                                                                                                                                               >  1:963606/1‑106 (MQ=255)
                                              gcggcAAATTACCCGTTTATCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGCTGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATc                                                                                                    >  4:804854/1‑148 (MQ=255)
                                                 gcAAATTACCCGTTTATCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGCTGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCgg                                                                                                 <  3:804854/148‑1 (MQ=255)
                                                 gcAAATTACCCGTTTATCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCAGCGCGTGGCAGa                                                                                                                                          <  4:317038/107‑1 (MQ=255)
                                                 gcAAATTACCCGTTTATCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCAGCGCGTGGCAGa                                                                                                                                          >  3:317038/1‑107 (MQ=255)
                                                    aaTTACCCGTTTATCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGCTGGCAGACATCGCCACCACGCa                                                                                                                           >  2:454142/1‑119 (MQ=255)
                                                    aaTTACCCGTTTATCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGCTGGCAGACATCGCCACCACGCa                                                                                                                           <  1:454142/119‑1 (MQ=255)
                                                         cccGTTTATCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGCTGGCAGACATCGCCACCACGCACCa                                                                                                                        >  1:301678/1‑117 (MQ=255)
                                                           cGTTTATCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGCTGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTg                                                                                               >  2:606509/1‑140 (MQ=255)
                                                           cGTTTATCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGCTGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCaa                                                                                       >  4:1782/1‑148 (MQ=255)
                                                           cGTTTATCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCAGCGCGTGGCAGa                                                                                                                                          <  2:234661/97‑1 (MQ=255)
                                                             tttATCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGCTGGCAGACATCGcca                                                                                                                                  >  1:636703/1‑103 (MQ=255)
                                                             tttATCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGCTGGCAGACATCGcca                                                                                                                                  <  2:636703/103‑1 (MQ=255)
                                                             tttATCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCAGCGCGTGGCAGa                                                                                                                                          <  1:482740/95‑1 (MQ=255)
                                                              ttATCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGCTGGCAGACATCGCCACCACGCACCa                                                                                                                        >  3:803321/1‑112 (MQ=255)
                                                              ttATCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGCTGGCAGACATCGCCACCACGCACCa                                                                                                                        >  3:403931/1‑112 (MQ=255)
                                                              ttATCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGCTGGCAGACATCGCCACCACGCACCa                                                                                                                        <  4:803321/112‑1 (MQ=255)
                                                              ttATCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGCTGGCAGACATCGCCACCACGCACCa                                                                                                                        <  4:403931/112‑1 (MQ=255)
                                                                   gACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGCTGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACg                                                                               >  1:347439/1‑148 (MQ=255)
                                                                    aCCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGCTGGCAGACATCGccacc                                                                                                                                <  3:56312/98‑1 (MQ=255)
                                                                       ccGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGCTGGCAGACATCGCCACCACGCACCAGGACa                                                                                                                   >  4:1036493/1‑108 (MQ=255)
                                                                       ccGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGCTGGCAGACATCGCCACCACGCACCAGGACATCCAGATTg                                                                                                          >  4:946646/1‑117 (MQ=255)
                                                                       ccGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGCTGGCAGACATCGCCACCACGCACCAGGACATCCAGATTg                                                                                                          <  3:946646/117‑1 (MQ=255)
                                                                       ccGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGCTGGCAGACATCGCCACCACCCACCAGGACa                                                                                                                   <  3:1036493/108‑1 (MQ=255)
                                                                           tAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGCTGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCa                                                                                             >  3:350939/1‑126 (MQ=255)
                                                                            aaaaaGATGATTCTGGTGACCGGTCACAGGCGTGAGAGATTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGCTGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAAc                                                                      >  2:739201/1‑148 (MQ=255)
                                                                             aaaaGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGCTGGCAGACAt                                                                                                                                       <  2:381735/82‑1 (MQ=255)
                                                                             aaaaGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCAGCGCGTGGCAGa                                                                                                                                          <  4:1038087/79‑1 (MQ=255)
                                                                               aaGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGCTGGCAGACATCGCCACCACGCACCAGGACATCCAGAt                                                                                                            <  1:297348/107‑1 (MQ=255)
                                                                               aaGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGCTGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCa                                                                                        <  3:108778/127‑1 (MQ=255)
                                                                               aaGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGCTGGCAGACATAGCCACCACGCACCAGGACATCCAGAt                                                                                                            >  2:297348/1‑107 (MQ=255)
                                                                                  atgatTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGCTGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATc                                                                                                    >  3:185349/1‑112 (MQ=255)
                                                                                  atgatTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGCTGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCa                                                                >  3:952005/1‑148 (MQ=255)
                                                                                  atgatTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGCTGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCa                                                                >  2:319124/1‑148 (MQ=255)
                                                                                      ttCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGCTGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGaa                                                                       >  4:522219/1‑137 (MQ=255)
                                                                                      ttCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGCTGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGaa                                                                       <  3:522219/137‑1 (MQ=255)
                                                                                       tCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGCTGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGc                                                           >  3:114782/1‑148 (MQ=255)
                                                                                             gACCGGTCACAGGCGTGAGAGTTTCGGTCGTAGCTTTGAAGAAATCTGCCACGCGCTGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCt                                                      <  1:1059260/147‑1 (MQ=255)
                                                                                             gACCGGTCACAGGCGTGAGAGTTTCGGTCGTAGCTTTGAAGAAATCTGCCACGCGCTGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCt                                                      >  2:1059260/1‑147 (MQ=255)
                                                                                                cGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGCTGGCAGACATCGCCACCACGCACCAGGACa                                                                                                                   <  2:805522/83‑1 (MQ=255)
                                                                                                cGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGCTGGCAGACATCGCCACCACGCACCAGGACa                                                                                                                   >  1:805522/1‑83 (MQ=255)
                                                                                                cGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGCTGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTggg                                                   >  4:64249/1‑147 (MQ=255)
                                                                                                   tCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGCTGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGt                                                                              >  2:631509/1‑117 (MQ=255)
                                                                                                   tCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGCTGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGt                                                                              <  1:631509/117‑1 (MQ=255)
                                                                                                            tGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGCTGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCAtgt                                             >  1:729713/1‑141 (MQ=255)
                                                                                                              agagTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGCTGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGAAAAATGt                                    >  1:619826/1‑148 (MQ=255)
                                                                                                               gagTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGCTGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGAAAAATGTc                                   >  1:157930/1‑148 (MQ=255)
                                                                                                                 gTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGCTGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGaaaaa                                       <  1:234245/142‑1 (MQ=255)
                                                                                                                 gTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGCTGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGAAAAATGTCAt                                 <  1:319124/148‑1 (MQ=255)
                                                                                                                   ttCGGTCGTGGCTTTGAAGAAATCTGCCACGCGCTGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTc                                                                             >  1:403261/1‑102 (MQ=255)
                                                                                                                   ttCGGTCGTGGCTTTGAAGAAATCTGCCACGCGCTGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTc                                                                             <  2:403261/102‑1 (MQ=255)
                                                                                                                      ggTCGTGGCTTTGAAGAAATCTGCCACGCGCTGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGAAAAATGTCATTCTga                            <  2:347439/148‑1 (MQ=255)
                                                                                                                          gTGGCTTTGAAGAAATCTGCCACGCGCTGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGAAAAATGTCATTCTgatcga                        <  3:746244/148‑1 (MQ=255)
                                                                                                                                ttGAAGAAATCTGCCACGCGCTGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGAAAAATGTCATTCTGATCGATCTCCa                  >  4:647307/1‑148 (MQ=255)
                                                                                                                                ttGAAGAAATCTGCCACGCGCTGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGAAAAATGTCATTCTGATCGATCCCCa                  <  4:114782/148‑1 (MQ=255)
                                                                                                                                ttGAAGAAATCTGCCACGCGCTGGCAGACATCGCCACCACGCACCAGGACATCCAGATAGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGAAAAATGTCATTCTGATCGATCCCCa                  >  4:235753/1‑148 (MQ=255)
                                                                                                                                ttGAAGAAATCTGCCACGCGCTGGCAGACATCGCCACCACGCACCAGGACATCCAGAATGTCTATCCgg                                                                                                 >  1:220711/1‑69 (MQ=255)
                                                                                                                                ttGAAGAAATCTGCCACGCGCTGGCAGACATCGCCACCACGCACCAGGACATCCAGAATGTCTATCCgg                                                                                                 <  2:220711/69‑1 (MQ=255)
                                                                                                                                   aagaaATCTGCCACGCGCTGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAAc                                                                                      <  4:209656/77‑1 (MQ=255)
                                                                                                                                   aagaaATCTGCCACGCGCTGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAAc                                                                                      >  3:209656/1‑77 (MQ=255)
                                                                                                                                   aagaaATCTGCCACGCGCTGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGAAAAATGTCATTCTGATCGATCCCCAGGa               >  4:911012/1‑148 (MQ=255)
                                                                                                                                    agaaATCTGCCACGCGCTGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCa                                                <  2:1126709/114‑1 (MQ=255)
                                                                                                                                    agaaATCTGCCACGCGCTGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCa                                                >  1:1126709/1‑114 (MQ=255)
                                                                                                                                    agaaATCTGCCACGCGCTGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGAAAAATGTCATTCTgat                           >  2:297872/1‑135 (MQ=255)
                                                                                                                                    agaaATCTGCCACGCGCTGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGAAAAATGTCATTCTgat                           <  1:297872/135‑1 (MQ=255)
                                                                                                                                    agaaATCTGCCACGCGCTGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGAAAAATGTCATTCTGATCGATCCCCAGGAg              >  4:173334/1‑148 (MQ=255)
                                                                                                                                    agaaATCTGCCACGCGCTGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGAAAAATGTCATTCTGATCGATCCCCAGGAg              >  4:827532/1‑148 (MQ=255)
                                                                                                                                     gaaATCTGCCACGCGCTGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGc                                                           >  4:510192/1‑102 (MQ=255)
                                                                                                                                     gaaATCTGCCACGCGCTGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGAAAAATGTCATTCTGATCGATCCCCa                  >  4:1021048/1‑143 (MQ=255)
                                                                                                                                     gaaATCTGCCACGCGCTGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGAAAAATGTCATTCTGATCGATCCCCa                  <  3:1021048/143‑1 (MQ=255)
                                                                                                                                     gaaATCTGCCACGCGCTGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGAAAAATGTCATTCTGATCGATCCCCAGGAGt             >  4:403755/1‑148 (MQ=255)
                                                                                                                                       tgtCTGCCAGCGCGTGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAAc                                                                      >  3:592322/3‑89 (MQ=255)
                                                                                                                                       tgtCTGCCAGCGCGTGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAAc                                                                      <  4:592322/87‑1 (MQ=255)
                                                                                                                                        aTCTGCCACGCGCTGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCAtg                                              <  1:925456/112‑1 (MQ=255)
                                                                                                                                        aTCTGCCACGCGCTGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGAAAAATGTCATTCTGATCGATCCCCAGGAGTAtt          <  3:911012/148‑1 (MQ=255)
                                                                                                                                         tCTGCCACGCGCTGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGAAAAATGTCATTCTGATCGATCCCCAGGAGTAttt         <  2:619826/148‑1 (MQ=255)
                                                                                                                                         tCTGCCACGCGCTGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGAAAAATGTCATTCTGATCGATCCCCAGGAGTAttt         >  3:192/1‑148 (MQ=255)
                                                                                                                                          cTGCCACGCGCTGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGAAAAATGTCATTCTGATCGATCCCCAGGAGTATTTa        <  4:952005/148‑1 (MQ=255)
                                                                                                                                             ccACGCGCTGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGAAAAATGTCATTCTGATCGATCCCCAGGAGTATTTACCg     <  2:1127746/148‑1 (MQ=255)
                                                                                                                                              cACGCGCTGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGAAAAATGTCATTCTGATCGATCCCCAGGAGTATTTACCGt    >  1:784876/1‑148 (MQ=255)
                                                                                                                                              cACGCGCTGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGAAAAATGTCATTCTGATCGATCCCCAGGAGTATTTACCGt    >  3:1076010/1‑148 (MQ=255)
                                                                                                                                              cACGCGCTGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCACCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATCCTGGGGCATGTGAAAAATGTCATTCTGATCGATCCCCAGGAGTATTTACCGt    <  3:647307/148‑1 (MQ=255)
                                                                                                                                               aCGCGCTGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGaaaaa                                       >  2:316913/1‑112 (MQ=255)
                                                                                                                                                cgcgCTGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCaa                                                               <  1:359774/87‑1 (MQ=255)
                                                                                                                                                cgcgCTGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCaa                                                               >  2:359774/1‑87 (MQ=255)
                                                                                                                                                cgcgCTGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGAAAAATGTCATTCTGATCGATCCCCAGGAGTAttt         >  2:1025974/1‑141 (MQ=255)
                                                                                                                                                cgcgCTGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGAAAAATGTCATTCTGATCGATCCCCAGGAGTAttt         <  1:1025974/141‑1 (MQ=255)
                                                                                                                                                cgcgCTGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGAAAAATGTCATTCTGATCGATCCCCAGGAGTATTTACCGttt  >  4:90102/1‑148 (MQ=255)
                                                                                                                                                |                                                                                                                                                    
GGTGCGTGACCAGGTGATGAGCAGCGACAAGCTGCGTTCAGAACTGGCGGCAAATTACCCGTTTATCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGCTGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGAAAAATGTCATTCTGATCGATCCCCAGGAGTATTTACCGTTT  >  NC_000913/3970669‑3970960

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