Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,211,186 G→C 5.3% intergenic (+86/‑39) yehQ → / → yehR pseudogene/lipoprotein, DUF1307 family

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,211,1860GC 5.3% 274.2 / 5.0 75intergenic (+86/‑39)yehQ/yehRpseudogene/lipoprotein, DUF1307 family
Reads supporting (aligned to +/- strand):  ref base G (34/37);  new base C (2/2);  total (36/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.91e-01

CAATGAATAACGGTGTTTTGCTTTGCGAACAATTATTCCCCATGCTTCAGCAACAAGGATGAAATAGTGCTTTTTACTAAGAGTTCTACTCCAGTTCCGGACTGCTCACGCCACGGTATTAGGCATATCCTATATAGCCCCTGGTGAGAGTCACCAGTTCCTTGATTAAATAAAATGGAGTTTTACATGAAGGCTTTCAATAAGCTGTTTTCCCTCGTTGTTGCATCTGTTCTGGTTTTCTCTCTTGCTGGCTGCGGTGACAAAGA  >  NC_000913/2211039‑2211304
                                                                                                                                                   |                                                                                                                      
cAATGAATAACGGTGTTTTGCTTTGCGAACAATTATTCCCCATGCTTCAGCAACAAGGATGAAATAGTGCTTTTTACTAAGAGTTCTACTCCAGTTCCGGACTGCTCACGCCACGGTATTAGGCATATCCTATATAGCCCCTGGTgag                                                                                                                        <  4:1046629/148‑1 (MQ=255)
  aTGAATAACGGTGTTTTGCTTTGCGAACAATTATTCCCCATGCTTCAGCAACAAGGATGAAATAGTGCTTTTTACTAAGAGTTCTACTCCAGTTCCGGACTGCTCACGCCACGGTATTAGGCATATCCTATATAGCCCCTGGTGagag                                                                                                                      <  1:1119850/148‑1 (MQ=255)
          cGGTGTTTTGCTTTGCGAACAATTATTCCCCATGCTTCAGCAACAAGGATGAAATAGTGCTTTTTACTAAGAGTTCTACTCCAGTTCCGGACTGCTCACGCCACGGTATTAGGCATATCCTATATAGCCCCTGGTGAGAGTCACCAGt                                                                                                              <  1:471844/148‑1 (MQ=255)
          cGGTGTTTTGCTTTGCGAACAATTATTCCCCATGCTTCAGCAACAAGGATGAAATAGTGCTTTTTACTAAGAGTTCTACTCCAGTTCCGGACTGCTCACGCCACGGTATTAGGCATATCCTATATAGCCCCTGGTGAGAGTCACCAGt                                                                                                              <  3:791149/148‑1 (MQ=255)
           ggTGTTTTGCTTTGCGAACAATTATTCCCCATGCTTCAGCAACAAGGATGAAATAGTGCTTTTTACTAAGAGTTCTACTCCAGTTCCGGACTGCTCACGCCACGGTATTAGGCATATCCTATATAGCCCCTGGTGAGAGTCACCAGtt                                                                                                             <  3:357951/148‑1 (MQ=255)
                      ttgcGAACAATTATTCCCCATGCTTCAGCAACAAGGATGAAATAGTGCTTTTTACTAAGAGTTCTACTCCAGTTCCGGACTGCTCACGCCACGGTATTAGGCATATCCTATATAGCCCCTGGTgag                                                                                                                        <  3:320449/126‑1 (MQ=255)
                      ttgcGAACAATTATTCCCCATGCTTCAGCAACAAGGATGAAATAGTGCTTTTTACTAAGAGTTCTACTCCAGTTCCGGACTGCTCACGCCACGGTATTAGGCATATCCTATATAGCCCCTGGTgag                                                                                                                        >  4:320449/1‑126 (MQ=255)
                            aCAATTATTCCCCATGCTTCAGCAACAAGGATGAAATAGTGCTTTTTACTAAGAGTTCTACTCCAGTTCCGGACTGCTCACGCCACGGTATTAGGCATATCCTATATAGCCCCTGGTGAGAGTCACCAGTTCCt                                                                                                          >  3:306580/1‑134 (MQ=255)
                            aCAATTATTCCCCATGCTTCAGCAACAAGGATGAAATAGTGCTTTTTACTAAGAGTTCTACTCCAGTTCCGGACTGCTCACGCCACGGTATTAGGCATATCCTATATAGCCCCTGGTGAGAGTCACCAGTTCCt                                                                                                          <  4:306580/134‑1 (MQ=255)
                            aCAATTATTCCCCATGCTTCAGCAACAAGGATGAAATAGTGCTTTTTACTAAGAGTTCTACTCCAGTTCCGGACTGCTCACGCCACGGTATTAGGCATATCCTATATAGCCCCTGGTGAGAGTCACCAGTTCCTTGATTAAATAAAAt                                                                                            >  3:909826/1‑148 (MQ=255)
                              aaTTATTCCCCATGCTTCAGCAACAAGGATGAAATAGTGCTTTTTACTAAGAGTTCTACTCCAGTTCCGGACTGCTCACGCCACGGTATTAGGCATATCCTATATAGCCCCTGGTGAGAGTCACCAGTTCCTTGATTAAATAAAATgg                                                                                          <  4:909826/148‑1 (MQ=255)
                                    tCCCCATGCTTCAGCAACAAGGATGAAATAGTGCTTTTTACTAAGAGTTCTACTCCAGTTCCGGACTGCTCACGCCACGGTATTAGGCATATCCTATATAGCCCCTGGTGACTCTCACCAg                                                                                                               <  3:770543/121‑1 (MQ=255)
                                    tCCCCATGCTTCAGCAACAAGGATGAAATAGTGCTTTTTACTAAGAGTTCTACTCCAGTTCCGGACTGCTCACGCCACGGTATTAGGCATATCCTATATAGCCCCTGGTGACTCTCACCAg                                                                                                               >  4:770543/1‑121 (MQ=255)
                                       ccATGCTTCAGCAACAAGGATGAAATAGTGCTTTTTACTAAGAGTTCTACTCCAGTTCCGGACTGCTCACGCCACGGTATTAGGCATATCCTATATAGCCCCTGGTGAGAGTCACCAGTTcc                                                                                                           >  3:895036/1‑122 (MQ=255)
                                       ccATGCTTCAGCAACAAGGATGAAATAGTGCTTTTTACTAAGAGTTCTACTCCAGTTCCGGACTGCTCACGCCACGGTATTAGGCATATCCTATATAGCCCCTGGTGAGAGTCACCAGTTcc                                                                                                           <  4:895036/122‑1 (MQ=255)
                                       ccATGCTTCAGCAACAAGGATGAAATAGTGCTTTTTACTAAGAGTTCTACTCCAGTTCCGGACTGCTCACGCCACGGTATTAGGCATATCCTATATAGCCCCTGGTGAGAGTCACCAGTTCCt                                                                                                          >  1:339774/1‑123 (MQ=255)
                                       ccATGCTTCAGCAACAAGGATGAAATAGTGCTTTTTACTAAGAGTTCTACTCCAGTTCCGGACTGCTCACGCCACGGTATTAGGCATATCCTATATAGCCCCTGGTGAGAGTCACCAGTTCCt                                                                                                          <  2:339774/123‑1 (MQ=255)
                                       ccATGCTTCAGCAACAAGGATGAAATAGTGCTTTTTACTAAGAGTTCTACTCCAGTTCCGGACTGCTCACGCCACGGTATTAGGCATATCCTATATAGCCCCTGGTGAGAGTCACCAGTTCCTTGATTAAATAAAATGGAGtt                                                                                      >  1:107124/1‑143 (MQ=255)
                                       ccATGCTTCAGCAACAAGGATGAAATAGTGCTTTTTACTAAGAGTTCTACTCCAGTTCCGGACTGCTCACGCCACGGTATTAGGCATATCCTATATAGCCCCTGGTGAGAGTCACCAGTTCCTTGATTAAATAAAATGGAGtt                                                                                      <  2:107124/143‑1 (MQ=255)
                                         aTGCTTCAGCAACAAGGATGAAATAGTGCTTTTTACTAAGAGTTCTACTCCAGTTCCGGACTGCTCACGCCACGGTATTAGGCATATCCTATATAGCCCCTGGTGAGAGTc                                                                                                                    <  1:1034644/111‑1 (MQ=255)
                                         aTGCTTCAGCAACAAGGATGAAATAGTGCTTTTTACTAAGAGTTCTACTCCAGTTCCGGACTGCTCACGCCACGGTATTAGGCATATCCTATATAGCCCCTGGTGAGAGTc                                                                                                                    >  2:65373/1‑111 (MQ=255)
                                         aTGCTTCAGCAACAAGGATGAAATAGTGCTTTTTACTAAGAGTTCTACTCCAGTTCCGGACTGCTCACGCCACGGTATTAGGCATATCCTATATAGCCCCTGGTGAGAGTc                                                                                                                    >  2:1034644/1‑111 (MQ=255)
                                         aTGCTTCAGCAACAAGGATGAAATAGTGCTTTTTACTAAGAGTTCTACTCCAGTTCCGGACTGCTCACGCCACGGTATTAGGCATATCCTATATAGCCCCTGGTGAGAGTc                                                                                                                    <  1:65373/111‑1 (MQ=255)
                                              tCAGCAACAAGGATGAAATAGTGCTTTTTACTAAGAGTTCTACTCCAGTTCCGGACTGCTCACGCCACGGTATTAGGCATATCCTATATAGCCCCTGGTgag                                                                                                                        >  3:831790/1‑102 (MQ=255)
                                              tCAGCAACAAGGATGAAATAGTGCTTTTTACTAAGAGTTCTACTCCAGTTCCGGACTGCTCACGCCACGGTATTAGGCATATCCTATATAGCCCCTGGTgag                                                                                                                        <  4:831790/102‑1 (MQ=255)
                                                     caaGGATGAAATAGTGCTTTTTACTAAGAGTTCTACTCCAGTTCCGGACTGCTCACGCCACGGTATTAGGCATATCCTATATAGCCCCTGGTGAGAGTCACCAGTTCCTTGATTAAATAAAAt                                                                                            >  1:16794/1‑123 (MQ=255)
                                                     caaGGATGAAATAGTGCTTTTTACTAAGAGTTCTACTCCAGTTCCGGACTGCTCACGCCACGGTATTAGGCATATCCTATATAGCCCCTGGTGAGAGTCACCAGTTCCTTGATTAAATAAAAt                                                                                            <  2:16794/123‑1 (MQ=255)
                                                               aTAGTGCTTTTTACTAAGAGTTCTACTCCAGTTCCGGACTGCTCACGCCACGGTATTAGGCATATCCTATATAGCCCCTGGTGAGAGTCACCAGt                                                                                                              >  2:462883/1‑95 (MQ=255)
                                                               aTAGTGCTTTTTACTAAGAGTTCTACTCCAGTTCCGGACTGCTCACGCCACGGTATTAGGCATATCCTATATAGCCCCTGGTGAGAGTCACCAGt                                                                                                              <  1:462883/95‑1 (MQ=255)
                                                                 aGTGCTTTTTACTAAGAGTTCTACTCCAGTTCCGGACTGCTCACGCCACGGTATTAGGCATATCCTATATAGCCCCTGGTGAGAGTCa                                                                                                                   <  1:120741/88‑1 (MQ=255)
                                                                 aGTGCTTTTTACTAAGAGTTCTACTCCAGTTCCGGACTGCTCACGCCACGGTATTAGGCATATCCTATATAGCCCCTGGTGAGAGTCa                                                                                                                   >  2:120741/1‑88 (MQ=255)
                                                                    gCTTTTTACTAAGAGTTCTACTCCAGTTCCGGACTGCTCACGCCACGGTATTAGGCATATCCTATATAGCCCCTGGTGAGAGTCACCAGTTCCTTGATTAAATAAAATGGAGTTTTACATGAAGGCTTTCAATAAGCTGtttt                                                         <  4:238516/143‑1 (MQ=255)
                                                                    gCTTTTTACTAAGAGTTCTACTCCAGTTCCGGACTGCTCACGCCACGGTATTAGGCATATCCTATATAGCCCCTGGTGAGAGTCACCAGTTCCTTGATTAAATAAAATGGAGTTTTACATGAAGGCTTTCAATAAGCTGtttt                                                         >  3:238516/1‑143 (MQ=255)
                                                                       ttttACTAAGAGTTCTACTCCAGTTCCGGACTGCTCACGCCACGGTATTAGGCATATCCTATATAGCCCCTGGTGAGAGTCACCAGTTCCTTGATTAAATAAAATGGAGtttt                                                                                    <  4:751821/113‑1 (MQ=255)
                                                                       ttttACTAAGAGTTCTACTCCAGTTCCGGACTGCTCACGCCACGGTATTAGGCATATCCTATATAGCCCCTGGTGAGAGTCACCAGTTCCTTGATTAAATAAAATGGAGtttt                                                                                    >  3:751821/1‑113 (MQ=255)
                                                                              aaGAGTTCTACTCCAGTTCCGGACTGCTCACGCCACGGTATTAGGCATATCCTATATAGCCCCTGGTGAGAGTCACCAGTTCCTTGATtaaataaa                                                                                              >  2:1092202/1‑96 (MQ=255)
                                                                              aaGAGTTCTACTCCAGTTCCGGACTGCTCACGCCACGGTATTAGGCATATCCTATATAGCCCCTGGTGAGAGTCACCAGTTCCTTGATtaaataaa                                                                                              <  1:1092202/96‑1 (MQ=255)
                                                                              aaGAGTTCTACTCCAGTTCCGGACTGCTCACGCCACGGTATTAGGCATATCCTATATAGCCCCTGGTGAGAGTCACCAGTTCCTTGATTAAATAAAATGGAGTTTTACATGAAGGCtt                                                                        <  2:107456/118‑1 (MQ=255)
                                                                              aaGAGTTCTACTCCAGTTCCGGACTGCTCACGCCACGGTATTAGGCATATCCTATATAGCCCCTGGTGAGAGTCACCAGTTCCTTGATTAAATAAAATGGAGTTTTACATGAAGGCtt                                                                        >  1:107456/1‑118 (MQ=255)
                                                                              aaGAGTTCTACTCCAGTTCCGGACTGCTCACGCCACGGTATTAGGCATATCCTATATAGCCCCTGGTGACTCTCACCAg                                                                                                               <  3:488279/79‑1 (MQ=255)
                                                                              aaGAGTTCTACTCCAGTTCCGGACTGCTCACGCCACGGTATTAGGCATATCCTATATAGCCCCTGGTGACTCTCACCAg                                                                                                               >  4:488279/1‑79 (MQ=255)
                                                                                gagTTCTACTCCAGTTCCGGACTGCTCACGCCACGGTATTAGGCATATCCTATATAGCCCCTGGTGAGAGTCa                                                                                                                   <  1:102635/73‑1 (MQ=255)
                                                                                gagTTCTACTCCAGTTCCGGACTGCTCACGCCACGGTATTAGGCATATCCTATATAGCCCCTGGTGAGAGTCa                                                                                                                   >  2:102635/1‑73 (MQ=255)
                                                                                  gTTCTACTCCAGTTCCGGACTGCTCACGCCACGGTATTAGGCATATCCTATATAGCCCCTGGTGAGAGTCACCAGTTCCTTGATTAAATAAAATGGAGTTTTACATGAAGGCTTTCAATAAGCTGTTTTCCCTCGTTGTTGCATCTGt                                      >  1:788283/1‑148 (MQ=255)
                                                                                   ttCTACTCCAGTTCCGGACTGCTCACGCCACGGTATTAGGCATATCCTATATAGCCCCTGGTGAGAGTCACCAGTTCCTTGATTAAATAAAATGGAGTTTTACATGAAGGCTTTCAATAAGCTGTTTTCCCTCgt                                                  >  1:1116136/1‑135 (MQ=255)
                                                                                   ttCTACTCCAGTTCCGGACTGCTCACGCCACGGTATTAGGCATATCCTATATAGCCCCTGGTGAGAGTCACCAGTTCCTTGATTAAATAAAATGGAGTTTTACATGAAGGCTTTCAATAAGCTGTTTTCCCTCgt                                                  <  2:1116136/135‑1 (MQ=255)
                                                                                    tCTACTCCAGTTCCGGACTGCTCACGCCACGGTATTAGGCATATCCTATATAGCCCCTGGTGAGAGTCACCAGTTCCTTGATtaaataaa                                                                                              <  2:874089/90‑1 (MQ=255)
                                                                                    tCTACTCCAGTTCCGGACTGCTCACGCCACGGTATTAGGCATATCCTATATAGCCCCTGGTGAGAGTCACCAGTTCCTTGATtaaataaa                                                                                              >  1:874089/1‑90 (MQ=255)
                                                                                    tCTACTCCAGTTCCGGACTGCTCACGCCACGGTATTAGGCATATCCTATATAGCCCCTGGTGAGAGTCACCAGTTCCTTGATTAAATAAAATGGAGTTTTACATGAAGGCTTTCAATAAGCTg                                                             >  2:725403/1‑123 (MQ=255)
                                                                                    tCTACTCCAGTTCCGGACTGCTCACGCCACGGTATTAGGCATATCCTATATAGCCCCTGGTGAGAGTCACCAGTTCCTTGATTAAATAAAATGGAGTTTTACATGAAGGCTTTCAATAAGCTg                                                             <  1:725403/123‑1 (MQ=255)
                                                                                         tCCAGTTCCGGACTGCTCACGCCACGGTATTAGGCATATCCTATATAGCCCCTGGTGAGAGTCACCAGTTCCt                                                                                                          <  4:322568/73‑1 (MQ=255)
                                                                                         tCCAGTTCCGGACTGCTCACGCCACGGTATTAGGCATATCCTATATAGCCCCTGGTGAGAGTCACCAGTTCCt                                                                                                          >  3:322568/1‑73 (MQ=255)
                                                                                            aGTTCCGGACTGCTCACGCCACGGTATTAGGCATATCCTATATAGCCCCTGGTGAGAGTCACCAGTTCCTTGATTAAATAAAATGGAGTTTTACATGaa                                                                             >  2:446819/1‑99 (MQ=255)
                                                                                            aGTTCCGGACTGCTCACGCCACGGTATTAGGCATATCCTATATAGCCCCTGGTGAGAGTCACCAGTTCCTTGATTAAATAAAATGGAGTTTTACATGaa                                                                             <  1:446819/99‑1 (MQ=255)
                                                                                             gTTCCGGACTGCTCACGCCACGGTATTAGGCATATCCTATATAGCCCCTGGTGAGAGTCACCAGTTCCTTGATTAAATAAAATGGAGTTTTACATGAAGGCTTTCAATAAGCTGTTTTCCCTCGTTGTTGCATCTGTTCTGGtttt                             >  2:4434/1‑146 (MQ=255)
                                                                                             gTTCCGGACTGCTCACGCCACGGTATTAGGCATATCCTATATAGCCCCTGGTGAGAGTCACCAGTTCCTTGATTAAATAAAATGGAGTTTTACATGAAGGCTTTCAATAAGCTGTTTTCCCTCGTTGTTGCATCTGTTCTGGtttt                             <  1:4434/146‑1 (MQ=255)
                                                                                               tCCGGACTGCTCACGCCACGGTATTAGGCATATCCTATATAGCCCCTGGTGAGAGTCACCAGTTCCTTGATTAAATAAAATgg                                                                                          >  1:92641/1‑83 (MQ=255)
                                                                                               tCCGGACTGCTCACGCCACGGTATTAGGCATATCCTATATAGCCCCTGGTGAGAGTCACCAGTTCCTTGATTAAATAAAATgg                                                                                          <  2:92641/83‑1 (MQ=255)
                                                                                                   gACTGCTCACGCCACGGTATTAGGCATATCCTATATAGCCCCTGGTGAGAGTCACCAGTTCCTTGATTAAATAAAATGGAGTTTTACATGAAGGCTTTCAATAAGCTGTTTTCCCTCGTTGTTGCATCTGTTCTGGTTTTCTCTCtt                      <  4:248633/147‑1 (MQ=255)
                                                                                                   gACTGCTCACGCCACGGTATTAGGCATATCCTATATAGCCCCTGGTGAGAGTCACCAGTTCCTTGATTAAATAAAATGGAGTTTTACATGAAGGCTTTCAATAAGCTGTTTTCCCTCGTTGTTGCATCTGTTCTGGTTTTCTCTCtt                      >  3:248633/1‑147 (MQ=255)
                                                                                                         tCACGCCACGGTATTAGGCATATCCTATATAGCCCCTGGTGAGAGTCACCAGTTCCTTGATtaaataaa                                                                                              <  1:84440/69‑1 (MQ=255)
                                                                                                         tCACGCCACGGTATTAGGCATATCCTATATAGCCCCTGGTGAGAGTCACCAGTTCCTTGATtaaataaa                                                                                              >  2:84440/1‑69 (MQ=255)
                                                                                                                aCGGTATTAGGCATATCCTATATAGCCCCTGGTGAGAGTCACCAGTTCCTTGATTAAATAAAATGGAGTTTTACATGAAGGCTTTCAATAAGCTGTTTTCCCTCGTTGTTGCATCTGTTCTGGTTTTCTCTCTTGCTGGCTGCGGTGa        <  2:788283/148‑1 (MQ=255)
                                                                                                                       tAGGCATATCCTATATAGCCCCTGGTGAGAGTCACCAGTTCCTTGATTAAATAAAATGGAGTTTTACATGAAGGCtt                                                                        >  3:417702/1‑77 (MQ=255)
                                                                                                                       tAGGCATATCCTATATAGCCCCTGGTGAGAGTCACCAGTTCCTTGATTAAATAAAATGGAGTTTTACATGAAGGCtt                                                                        <  4:417702/77‑1 (MQ=255)
                                                                                                                         ggCATATCCTATATAGCCCCTGGTGAGAGTCACCAGTTCCTTGATTAAATAAAATGGAGTTTTACATGAAGGCTTTCAATAAGCTGTTTTCCCTCGTTGTTGCATCTGt                                      >  3:147625/1‑109 (MQ=255)
                                                                                                                         ggCATATCCTATATAGCCCCTGGTGAGAGTCACCAGTTCCTTGATTAAATAAAATGGAGTTTTACATGAAGGCTTTCAATAAGCTGTTTTCCCTCGTTGTTGCATCTGt                                      <  4:147625/109‑1 (MQ=255)
                                                                                                                          gCATATCCTATATAGCCCCTGGTGAGAGTCACCAGTTCCTTGATTAAATAAAATGGAGTTTTACATGAAGGCTTTCAATAAGCTGTTTTCCCTCGTTGTTGCATCTGTTCTGGTTTTCTCTCTTgc                    <  1:218917/126‑1 (MQ=255)
                                                                                                                          gCATATCCTATATAGCCCCTGGTGAGAGTCACCAGTTCCTTGATTAAATAAAATGGAGTTTTACATGAAGGCTTTCAATAAGCTGTTTTCCCTCGTTGTTGCATCTGTTCTGGTTTTCTCTCTTgc                    >  2:218917/1‑126 (MQ=255)
                                                                                                                                 cTATATAGCCCCTGGTGAGAGTCACCAGTTCCTTGATTAAATAAAATGGAGTTTTACATGAAGGCTTTCAATAAGCTGTTTTCCCTCGTTGTTGCATCTGTTCTGGTTTtctct                         >  1:153262/1‑114 (MQ=255)
                                                                                                                                 cTATATAGCCCCTGGTGAGAGTCACCAGTTCCTTGATTAAATAAAATGGAGTTTTACATGAAGGCTTTCAATAAGCTGTTTTCCCTCGTTGTTGCATCTGTTCTGGTTTtctct                         <  2:153262/114‑1 (MQ=255)
                                                                                                                                    tataGCCCCTGGTGAGAGTCACCAGTTCCTTGATTAAATAAAATGGAGTTTTACATGAAGGCTTTCAATAAGCTGTTTTCCCTCGTTGTTGCATCTGTTCTGGTTTtctct                         >  1:552745/1‑111 (MQ=255)
                                                                                                                                    tataGCCCCTGGTGAGAGTCACCAGTTCCTTGATTAAATAAAATGGAGTTTTACATGAAGGCTTTCAATAAGCTGTTTTCCCTCGTTGTTGCATCTGTTCTGGTTTtctct                         <  2:552745/111‑1 (MQ=255)
                                                                                                                                     ataGCCCCTGGTGAGAGTCACCAGTTCCTTGATTAAATAAAATGGAGTTTTACATGAAGGCTTTCAGTAAGCTGTTTTCCCTGGTTGTTGCATCTGTTCTGGTTTTCTCTCTTgctggct               >  2:758617/1‑120 (MQ=255)
                                                                                                                                     ataGCCCCTGGTGAGAGTCACCAGTTCCTTGATTAAATAAAATGGAGTTTTACATGAAGGCTTTCAATAAGCTGTTTTCCCTCGTTGTTGCATCTGTTCTGGTTTTCTCTCTTgctggct               <  1:758617/120‑1 (MQ=255)
                                                                                                                                         ccccTGGTGAGAGTCACCAGTTCCTTGATTAAATAAAATGGAGTTTTACATGAAGGCTTTCAATAAGCTGTTTTCCCTCGTTGTTGCATCTGTTCTGGTTTTCTCTCTTgctg                  <  3:611661/113‑1 (MQ=255)
                                                                                                                                         ccccTGGTGAGAGTCACCAGTTCCTTGATTAAATAAAATGGAGTTTTACATGAAGGCTTTCAATAAGCTGTTTTCCCTCGTTGTTGCATCTGTTCTGGTTTTCTCTCTTgctg                  >  4:611661/1‑113 (MQ=255)
                                                                                                                                             tGGTGAGAGTCACCAGTTCCTTGATTAAATAAAATGGAGTTTTACATGAAGGCTTTCAATAAGCTGTTTTCCCTCGttgttg                                             <  4:354498/82‑1 (MQ=255)
                                                                                                                                             tGGTGAGAGTCACCAGTTCCTTGATTAAATAAAATGGAGTTTTACATGAAGGCTTTCAATAAGCTGTTTTCCCTCGttgttg                                             >  3:354498/1‑82 (MQ=255)
                                                                                                                                               gTGAGAGTCACCAGTTCCTTGATTAAATAAAATGGAGTTTTACATGAAGGCTTTCAATAAGCTGTTTTCCCTCGTTGTTGCATCTGTTCTGGTTTTCTCTCTTGCTGGCTGCGGTGACAaaga  <  2:26000/123‑1 (MQ=255)
                                                                                                                                               gTGAGAGTCACCAGTTCCTTGATTAAATAAAATGGAGTTTTACATGAAGGCTTTCAATAAGCTGTTTTCCCTCGTTGTTGCATCTGTTCTGGTTTTCTCTCTTGCTGGCTGCGGTGACAaaga  >  1:26000/1‑123 (MQ=255)
                                                                                                                                                 gagaGTCACCAGTTCCTTGATTAAATAAAATGGAGTTTTACATGAAGGCTTTCAATAAGCTGTTTTCCCTCGTTGTTGCATCTGTTCTGGTTTtctct                         <  1:184035/98‑1 (MQ=255)
                                                                                                                                                 gagaGTCACCAGTTCCTTGATTAAATAAAATGGAGTTTTACATGAAGGCTTTCAATAAGCCGTTTTCCCTCGTTGTTGCATCTGTTCTGGTTTtctct                         >  2:184035/1‑98 (MQ=255)
                                                                                                                                                   |                                                                                                                      
CAATGAATAACGGTGTTTTGCTTTGCGAACAATTATTCCCCATGCTTCAGCAACAAGGATGAAATAGTGCTTTTTACTAAGAGTTCTACTCCAGTTCCGGACTGCTCACGCCACGGTATTAGGCATATCCTATATAGCCCCTGGTGAGAGTCACCAGTTCCTTGATTAAATAAAATGGAGTTTTACATGAAGGCTTTCAATAAGCTGTTTTCCCTCGTTGTTGCATCTGTTCTGGTTTTCTCTCTTGCTGGCTGCGGTGACAAAGA  >  NC_000913/2211039‑2211304

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: