Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 114,511 A→C 6.3% intergenic (+24/+11) guaC → / ← hofC GMP reductase/assembly protein in type IV pilin biogenesis, transmembrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913114,5110AC 6.3% 294.6 / 10.7 96intergenic (+24/+11)guaC/hofCGMP reductase/assembly protein in type IV pilin biogenesis, transmembrane protein
Reads supporting (aligned to +/- strand):  ref base A (46/44);  new base C (3/3);  total (49/47)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.94e-01

Only 100 of 103 total aligned reads displayed.

TATTTTGGGCGGCCTGCGTTCAGCTTGTACATACGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCACTCATCGCATCGCCTAAATGGAAAATTGGCAGATACATTGCCACCACCAGCGTACCAATAATTCCTCCCGTTATGATCAGCAACGCCGGTTCCAGTAAGGCTGCGAGGTTATCCGC  >  NC_000913/114364‑114650
                                                                                                                                                   |                                                                                                                                           
tatTTTGGGCGGCCTGCGTTCAGCTTGTACATACGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCaa                                                                                                                                             >  6:388443/1‑148 (MQ=255)
  ttttGGGCGGCCTGCGTTCAGCTTGTACATACGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAacac                                                                                                                                          >  4:347259/1‑149 (MQ=255)
  ttttGGGCGGCCTGCGTTCAGCTTGTACATACGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAacac                                                                                                                                          >  2:198162/1‑149 (MQ=255)
       ggcggcCTGCGTTCAGCTTGTACATACGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCa                                                                                                                                      <  5:2717/148‑1 (MQ=255)
       ggcggcCTGCGTTCAGCTTGTACATACGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCa                                                                                                                                      <  1:37272/148‑1 (MQ=255)
       ggcggcCTGCGGTCAGCTTGTACATACGTTGGGGCTTCACGCCTGAAAGAGCGGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCa                                                                                                                                      <  1:240618/148‑1 (MQ=255)
             cTGCGTTCAGCTTGTACATACGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGTAGCAACACGCCAGCGttgg                                                                                                                               <  5:186142/149‑3 (MQ=255)
             cTGCGTTCAGCTTGTACATACGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTat                                                                                                                               >  4:97096/1‑149 (MQ=255)
             cTGCGTTCAGCTTGTACATACGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTat                                                                                                                               <  3:392451/149‑1 (MQ=255)
                cGTTCAGCTTGTACATACGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTatc                                                                                                                              <  5:158166/147‑1 (MQ=255)
                cGTTCAGCTTGTACATACGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTatc                                                                                                                              <  7:368692/147‑1 (MQ=255)
                cGTTCAGCTTGTACATACGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGGGGAGCAACACGCCAGCGTTatccc                                                                                                                            <  2:34257/149‑1 (MQ=255)
                    cAGCTTGTACATACGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTatcccatc                                                                                                                         <  5:388443/148‑1 (MQ=255)
                     aGCTTGTACATACGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCCTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTatcccatcc                                                                                                                        <  3:97096/148‑1 (MQ=255)
                     aGCTTGTACAGACGTTGGGGCGTCACGCCTGAAAGAGCGGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGGGGAGCAACACGCCAGCGTTatcccatcc                                                                                                                        <  6:168906/148‑1 (MQ=255)
                        ttGTACATACGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCACt                                                                                                                    >  7:23935/1‑149 (MQ=255)
                          gTACATACGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGTTGCTCCACGCCa                                                                                                                                      <  1:288015/129‑1 (MQ=255)
                          gTACATACGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGTTGCTCCACGCCa                                                                                                                                      >  2:288015/1‑129 (MQ=255)
                            aCATACGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTAtcccatccca                                                                                                                      <  5:224223/143‑1 (MQ=255)
                            aCATACGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTAtcccatccca                                                                                                                      >  6:224223/1‑143 (MQ=255)
                             cATACGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCACTcatc                                                                                                                <  1:198162/148‑1 (MQ=255)
                               tACGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCCCGcatcgca                                                                                                             >  5:374982/1‑149 (MQ=255)
                                aCGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGc                                                                                                                                        >  1:58420/1‑121 (MQ=255)
                                aCGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGc                                                                                                                                        <  2:58420/121‑1 (MQ=255)
                                  gTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCACTcatcgcatcg                                                                                                          >  2:139517/1‑149 (MQ=255)
                                   ttGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCACTCatcgcatcgc                                                                                                         >  3:151100/1‑149 (MQ=255)
                                     ggggCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCACTcatcgcat                                                                                                            >  4:17229/1‑144 (MQ=255)
                                     ggggCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCACTcatcgcat                                                                                                            <  3:17229/144‑1 (MQ=255)
                                        gCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGc                                                                                                                                    >  6:61441/1‑117 (MQ=255)
                                        gCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGc                                                                                                                                    <  5:61441/117‑1 (MQ=255)
                                         cTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCACTCATCGCATCGCCTAAAt                                                                                                   <  4:151100/149‑1 (MQ=255)
                                             aCGCCTGAAAGAGCTGACCAAGCCCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCACTCATCGCATCGCCTAAATGGa                                                                                                >  3:354447/1‑148 (MQ=255)
                                               gCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCACTCATCGCATCGCCTAAAt                                                                                                   >  3:117610/1‑143 (MQ=255)
                                               gCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCACTCATCGCATCGCCTAAAt                                                                                                   <  4:117610/143‑1 (MQ=255)
                                                 cTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCACTCATCGCATCGCCTAAATGGAAAAtt                                                                                           >  5:324099/1‑149 (MQ=255)
                                                   gAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCACTCATCGCATCGCCTAAATGGAAAATTgg                                                                                         >  2:308474/1‑149 (MQ=255)
                                                    aaaGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCACTCATCGCATCGCCTAAATGGAAAATTgg                                                                                         >  2:246870/1‑148 (MQ=255)
                                                    aaaGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCACTCATCGCATCGCCTAAATGGAAAATTGGc                                                                                        <  4:147056/149‑1 (MQ=255)
                                                     aaGAGCTGACCAAGCGCACCATGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACTCTGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCACTCATCGCATCGCCTAAATGGAAAATTGGCa                                                                                       <  3:251466/149‑1 (MQ=255)
                                                      agagCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCACTCATCGCATCGCCTAAATGGAAAATTgg                                                                                         >  6:242679/1‑146 (MQ=255)
                                                        agCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCACTCATCGCATCGCCTAAATGGAAAATTGGCAGAt                                                                                    <  6:405207/149‑1 (MQ=255)
                                                              ccAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCACTCATCGCATCGCCTAAATGGAAAATTGGCAGATACATTg                                                                              <  2:97745/149‑1 (MQ=255)
                                                                aaGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCaa                                                                                                                                             >  2:207153/1‑84 (MQ=255)
                                                                aaGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCaa                                                                                                                                             <  1:207153/84‑1 (MQ=255)
                                                                aaGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCACTCATCGCATCGCCTAAATGGAAAATTGGCAGATACATTGc                                                                             >  6:203098/1‑148 (MQ=255)
                                                                      accacGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCACTCATCGCATCGCCTAAATGGAAAATTGGCAGATACATTGccaccac                                                                       >  2:28509/1‑148 (MQ=255)
                                                                        cacGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCACTCATCGCATCGCCTAAATGGaa                                                                                               >  3:165381/1‑122 (MQ=255)
                                                                        cacGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCACTCATCGCATCGCCTAAATGGaa                                                                                               <  4:165381/122‑1 (MQ=255)
                                                                        cacGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCACTCATCGCATCGCCTAAATGGAAAATTGGCAGATACATTGCCACCACCAg                                                                    >  5:371174/1‑149 (MQ=255)
                                                                           gTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCACTCATCGCATCGCCTAAATGGAAAATTGGCAGATACATTGCCACCACCAGCg                                                                  >  3:68321/1‑148 (MQ=255)
                                                                              tATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCACTCATCGCATCGCCTAAATGGAAAATTGGCAGATACATTGCCACCACCAGCGTAc                                                               <  1:139517/148‑1 (MQ=255)
                                                                                   gtgtGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCACTCATCGCATCGCCTAAATGGAAAATTGGCAGATACATTGCCACCACCAGCGTACCaataa                                                         <  7:49930/149‑1 (MQ=255)
                                                                                       gCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCACTCATCGCATCGCCTAAATGGAAAATTGGCAGATACATTGCCACCACCAGCGTACCAATAATtcc                                                     >  6:251708/1‑149 (MQ=255)
                                                                                        cAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCTTGGAGCAACACGCCAGCGTTATCCCATCCCCCTCATCGCATCGCCTAAATGGAAAATTGGCAGATACATTGCCACCACCAGCGTACCAATAATtcc                                                     <  8:322321/148‑1 (MQ=255)
                                                                                         aGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTatcccatcc                                                                                                                        >  6:199263/1‑80 (MQ=255)
                                                                                         aGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTatcccatcc                                                                                                                        <  5:199263/80‑1 (MQ=255)
                                                                                         aGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCACTCATCGCATCGCCTAAATGGAAAATTGGCAGATACATTGCCACCACCAGCGTACCAATAATtcctc                                                   >  2:190125/1‑149 (MQ=255)
                                                                                           gAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCACTCATCGCATCGCCTAAATGGAAAATTGGCAGATACATTGCCACCACCAGCGTACCAATAATTCCTccc                                                 <  8:23935/149‑1 (MQ=255)
                                                                                            aaCAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCACTCATCGCATCGCCTAAATGGAAAATTGGCAGATACAtt                                                                               >  7:199525/1‑118 (MQ=255)
                                                                                            aaCAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCACTCATCGCATCGCCTAAATGGAAAATTGGCAGATACAtt                                                                               <  8:199525/118‑1 (MQ=255)
                                                                                            aaCAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCACTCATCGCATCGCCTAAATGGAAAATTGGCAGATACATTGCCACCACCAGCGTACCAATAATTCCTCCCg                                                >  6:354476/1‑149 (MQ=255)
                                                                                             aCAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCACTCATCGCATCGCCTAAATGGAAAATTGGCAGATACATTGCCACCACCAGCGTACCAATAATTCCTCCCGt                                               <  3:347259/149‑1 (MQ=255)
                                                                                             aCAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCACTCATCGCATCGCCTAAATGGAAAATTGGCAGATACATTGCCACCACCAGCGTACCAATAATTCCTCCCGt                                               >  7:212172/1‑149 (MQ=255)
                                                                                             aCAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCACTCATCGCATCGCCTAAATGGAAAATTGGCAGATACATTGCCACCACCAGCGTACCAATAATTCCTCCCGt                                               >  8:250597/1‑149 (MQ=255)
                                                                                                aGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCACTCATCGCATCGCCTAAATGGAAAATTGGCAGATACATTGCCACCACCAGCGTACCAATAATTCCTCCCGTTAt                                            <  8:212172/149‑1 (MQ=255)
                                                                                                    aaCCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCACTCATCGCATCGCCTAAATGGAAAATTGGCAGATACATTGCCACCACCAGCGTACCaataa                                                         >  3:335997/1‑132 (MQ=255)
                                                                                                    aaCCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCACTCATCGCATCGCCTAAATGGAAAATTGGCAGATACATTGCCACCACCAGCGTACCAATAATtcc                                                     >  7:328926/1‑136 (MQ=255)
                                                                                                    aaCCGCATCTTCAACAACCTGTAATCTCACAACGCTGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCACTCATCGCATCGCCTAAATGGAAAATTGGCAGATACATTGCCACCACCAGCGTACCaataa                                                         <  4:335997/132‑1 (MQ=255)
                                                                                                     aCCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCACTCATCGCATCGCCTAAATGGAAAATTGGCAGATACATTGCCACCACCAGCGTACCAATAATTCCTCCCGTTATGAt                                         >  5:10965/1‑147 (MQ=255)
                                                                                                      ccGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCACTCATCGCATCGCCTAAATGGAAAATTGGCAGATACATTGCCACCACCAGCGTACCAATAATTCCTCCCGTTATGATCAg                                      >  7:328596/1‑149 (MQ=255)
                                                                                                         cATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCACTCATCGCATCGCCTAAATGGAAAATTGGCAGATACATTGCCACCACCAGCGTACCAATAATTCCTCCCGTTATGATCAGCaa                                   <  5:354476/149‑1 (MQ=255)
                                                                                                            cTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCACTCATCGCATCGCCTAAATGGAAAATTGGCAGATACATTGCCACCACCAGCGTACCAATAATTCCTCCCGTTATGATCAGCAACGc                                >  2:102252/1‑149 (MQ=255)
                                                                                                            cTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCACTCATCGCATCGCCTAAATGGAAAATTGGCAGATACATTGCCACCACCAGCGTACCAATAATTCCTCCCGTTATGATCAGCAACGc                                >  6:407085/1‑149 (MQ=255)
                                                                                                              tCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCACTCATCGCATCGCCTAAATGGAAAATTGGCAGATACATTGCCACCACCAGCGTACCAATAATTCCTCCCGTTATGATCAGCAACGCCg                              >  1:166133/1‑149 (MQ=255)
                                                                                                               caacaaCCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCACTCATCGCATCGCCTAAATGGAAAATTGGCAGATACATTGCCACCACCAGCGTACCAATAATTCCTCCCGTTATGATCAGCAACGc                                <  1:190125/146‑1 (MQ=255)
                                                                                                                   aacCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCACTCATCGCATCGCCTAAATGGAAAATTGGCAGATACATTGCCACCACCAGCGTACCAATAATTCCTCCCGTTATGATCAGCAACGCCGGTTc                          >  1:214578/1‑148 (MQ=255)
                                                                                                                    acCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCACTCATCGCATCGCCTAAATGGAAAATTGGCAGATACATTGCCACCACCAGCGTACCAATAATTCCTCCCGTTATGATCAGCAACGCCGGTTCCa                        <  5:251708/149‑1 (MQ=255)
                                                                                                                     ccTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCTTTATCCCATCCGACTCATCGCATCGCCTAAATGGAAAATTGGCAGATACATTGCCACCACCAGCGTACCAATAATTCCTCCCGTTATGATCAGCAACGCCGGTTc                          >  4:46006/1‑146 (MQ=255)
                                                                                                                     ccTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCACTCATCGCATCGCCTAAATGGAAAATTGGCAGATACATTGCCACCACCAGCGTACCAATAATTCCTCCCGTTATGATCAGCAACGCCGGTTCCAg                       <  1:308474/149‑1 (MQ=255)
                                                                                                                      cTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCACTCATCGCATCGCATAAATGGAAAATTGGCAGATACATTGCCACCACCAGCGTACCAATAATTCCTCCCGTTATGATCAGCAACGCCGGTTCCAg                       >  3:391563/1‑148 (MQ=255)
                                                                                                                       tGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCACTCATCGCATCGCCTAAATGGAAAATTGGCAGATACATTGCCACCACCAGCGTACCAATAATTCCTCCCGTTATGATCAGCAACGCCGGTTCCAGt                      >  3:350229/1‑148 (MQ=255)
                                                                                                                         tAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCACTCATCGCATCGCCTAAATGGAAAATTGGCAGATACATTGCCACCACCAGCGTACCAATAATTCCTCCCGTTATGATCAGCAACGCCGGTTCCAGTaa                    <  1:102252/148‑1 (MQ=255)
                                                                                                                         tAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCACTCATCGCATCGCCTAAATGGAAAATTGGCAGATACATTGCCACCACCAGCGTACCAATAATTCCTCCCGTTATGATCAGCAACGCCGGTTCCAGTaa                    <  5:242679/148‑1 (MQ=255)
                                                                                                                           aTCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCACTCATCGCATCGCCTAAATGGAAAATTGGCAGATACATTGCCACCACCAGCGTACCAATAATTCCTCCCGTTATGATCAGCAACGCCGGTTCCAGTAAGGc                 <  4:350229/149‑1 (MQ=255)
                                                                                                                              tcCCAACGCTGGCGTGGAGCAACCCGCCAGCGTTATCCCATCCCCCTCATCGCATCGCCTAAATGGAAAATTGGCAGATACATTGCCACCACCAGCGTACCAATAATTCCTCACGTTATGATCAGCAACGCCGGTTCCAGTAAGGCTGc              >  6:223435/1‑149 (MQ=255)
                                                                                                                                ccAACGCTGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCACTCATCGCATCGCCTAAATGGAAAATTGGCAGATACATTGCCACCACCAGCGTACCAATAATTCCTCCCGTTATGATCAGCAACGCCGGTTCCAGTAAGGCt                <  2:214578/145‑1 (MQ=255)
                                                                                                                                ccAACGCTGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCACTCATCGCATCGCCTAAATGGAAAATTGGCAGATACATTGCCACCACCAGCGTACCAATAATTCCTCCCGTTATGATCAGCAACGCCGGTTCCAGTAAGGCTg               <  8:328596/146‑1 (MQ=255)
                                                                                                                                ccAACGCTGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCACTCATCGCATCGCCTAAATGGAAAATTGGCAGATACATTGCCACCACCAGCGTACCAATAATTCCTCCCGTTATGATCAGCAACGCCGGTTCCAGTAAGGCTGc              >  6:133880/1‑147 (MQ=255)
                                                                                                                                  aaCGCTGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCACTCATCGCATCGCCTAAATGGAAAATTGGCAGATACATTGCCACCACCAGCGTACCAATAATTCCTCCCGTTATGATCAGCAACGCCGGTTCCAGTAAGGCTGCGAgg          <  6:10965/149‑1 (MQ=255)
                                                                                                                                    cGCTGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCACTCATCGCATCGCCTAAATGGAAAATTGGCAGATACATTGCCACCACCAGCGTACCAATAATTCCTCCCGTTATGATCAGCAACGCCGGTTCCAGTAAGGCTGCGAGGtt        <  2:166133/149‑1 (MQ=255)
                                                                                                                                     gCTGGCGTGTTGCTCCACGCCAGCGTTATCCCATCCCACTCATCGCATCGCCTAAATGGAAAATTGGCAGATACATTGCCACCACCAGCGTACCaat                                                           <  4:182440/97‑1 (MQ=39)
                                                                                                                                     gCTGGCGTGTTGCTCCACGCCAGCGTTATCCCATCCCACTCATCGCATCGCCTAAATGGAAAATTGGCAGATACATTGCCACCACCAGCGTACCaat                                                           >  3:182440/1‑97 (MQ=39)
                                                                                                                                       tGGCGTGTTGCTCCACGCCAGCGTTATCCCATCCCACTCATCGCATCGCCTAAATGGAAAATTGGCAGATACATTGCCACCACCAGCGTACCAATAATTCCTCCCGTTATGATCAGCAACGCCGGTTCCa                        <  2:158925/130‑1 (MQ=255)
                                                                                                                                       tGGCGTGTTGCTCCACGCCAGCGTTATCCCATCCCACTCATCGCATCGCCTAAATGGAAAATTGGCAGATACATTGCCACCACCAGCGTACCAATAATTCCTCCCGTTATGATCAGCAACGCCGGTTCCa                        >  1:158925/1‑130 (MQ=255)
                                                                                                                                       tGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCACTCATCGCATCGCCTAAATGGAAAATTGGCAGATACATTGCCACCACCAGCGTACCAATAATTCCTCCCGTTATGATCAGCAACGCCGGTTCCAGTAAGGCTGCGAGGTTATc     >  6:94106/1‑149 (MQ=255)
                                                                                                                                       tGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCACTCATCGCATCGCCTAAATGGAAAATTGGCAGATACATTGCCACCACCAGCGTACCAATAATTCCTCCCGTTATGATCAGCAACGCCGGTTCCAGTAAGGCTGCGAGGTTATc     >  5:163578/1‑149 (MQ=255)
                                                                                                                                       tGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCACTCATCGCATCGCCTAAATGGAAAATTGGCAGACACATTGCCACCACCAGCGTACCAATAATTCCTCCCGTTATGATCAGCAACGCCGGTTCCAGTAAGGCTGCGAGGTTATc     >  7:260012/1‑149 (MQ=255)
                                                                                                                                            tGGAGCAACACGCCAGCGTTATCCCATCCCACTCATCGCATCGCCTAAATGGAAAATTGGCAGATACATTGCCACCACCAGCGTACCAATAATTCCTCCCGTTATGATCAGCAACGCCGGTTCCAGTAAGGCTGCGAGGTTATCCGc  >  6:73004/1‑147 (MQ=255)
                                                                                                                                                   acacGCCAGCGTTATCCCATCCCACTCATCGCATCGCCTAAATGGAAAATTGGCAGATACATTGCCACCACCAGCGTACCAATAATTCCTCCCGTTATGATCAGc                                     >  7:171607/1‑105 (MQ=255)
                                                                                                                                                   acacGCCAGCGTTATCCCATCCCACTCATCGCATCGCCTAAATGGAAAATTGGCAGATACATTGCCACCACCAGCGTACCAATAATTCCTCCCGTTATGATCAGc                                     <  8:171607/105‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                           
TATTTTGGGCGGCCTGCGTTCAGCTTGTACATACGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCACTCATCGCATCGCCTAAATGGAAAATTGGCAGATACATTGCCACCACCAGCGTACCAATAATTCCTCCCGTTATGATCAGCAACGCCGGTTCCAGTAAGGCTGCGAGGTTATCCGC  >  NC_000913/114364‑114650

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