Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,336,763 A→C 6.4% intergenic (‑119/+30) yfaA ← / ← gyrA DUF2138 family protein, putative host defense protein/DNA gyrase (type II topoisomerase), subunit A

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,336,7630AC 6.4% 455.3 / 20.2 125intergenic (‑119/+30)yfaA/gyrADUF2138 family protein, putative host defense protein/DNA gyrase (type II topoisomerase), subunit A
Reads supporting (aligned to +/- strand):  ref base A (60/57);  new base C (4/4);  total (64/61)
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.31e-01

Only 100 of 129 total aligned reads displayed.

AGCCTTTCGCCTTTTTTTCACCACTCATTGCCACATTCCTTGTGTATAGCCAGCCATTTTTTACGGGCACAGCCAAACTTTACCGTGCCCTAATACGACAAAAGCCCAGACTTTGCAGCCTGGACTTTTCAATTCAAACAAGGGAGATAGCTCCCTTTTGGCATGAAGAAGTAAAATTATTCTTCTTCTGGCTCGTCGTCAACGTCCACTTCCGGAGCGATTTCATCGTCCCCTTCCGCGGCACTGCCGTCGATGGTATCCAGATCTTCCTCGTCAACCGGTTCAGCAACACGTT  >  NC_000913/2336617‑2336911
                                                                                                                                                  |                                                                                                                                                    
aGCCTTTCGCCATTTTTTCACCACTCATTGCCACATTCCTTGTGTATAGCCAGCCATTTTTTACGGGCACAGCCAAACTTTACCGTGCCCTAATACGACAAAAGCCCAGACTTTGCAGCCTGGACTTTTCAATTCAAACAAGGGAGATa                                                                                                                                                    >  3:1635358/1‑149 (MQ=255)
    tttCGCCTTTTTTTCACCACTCATTGCCACATTCCTTGTGTATAGCCAGCCATTTTTTACGGGCACAGCCAAACTTTACCGTGCCCTAATACGACAAAAGCCCAGACTTTGCAGCCTGGACTTTTCAATTCAAACAAGGGAGATAGCTc                                                                                                                                                >  1:1700840/1‑149 (MQ=255)
     ttCGCCTTTTTTTCACCACTCATTGCCACATTCCTTGTGTATAGCCAGCCATTTTTTACGGGCACAGCCAAACTTTACCGTGCCCTAATACGACAAAAGCCCAGACTTTGCAGCCTGGACTTTTCAATTCAAACAAGGGAGATAGCTcc                                                                                                                                               <  2:1524672/149‑1 (MQ=255)
        gCCTTTTTTTCACCACTCATTGCCACATTCCTTGTGTATAGCCAGCCATTTTTTACGGGCACAGCCAAACTTTACCGTGCCCTAATACGACAAAAGCCCAGACTTTGCAGCCTGGACTTTTCAATTCAAACAAGGGAGCTATCTCCCtt                                                                                                                                            >  4:1391537/1‑149 (MQ=255)
         ccTTTTTTTCACCACTCATTGCCACATTCCTTGTGTATAGCCAGCCATTTTTTACGGGCACAGCCAAACTTTACCGTGCCCTAATACGACAAAAGCCCAGACTTTGCAGCCTGGACTTTTCAATTCAAACAAGGGAGATAGCTCCCttt                                                                                                                                           >  1:803710/1‑149 (MQ=255)
           tttttttCACCACTCATTGCCACATTCCTTGTGTATAGCCAGCCATTTTTTACGGGCACAGCCAAACTTTACCGTGCCCTAATACGACAAAAGCCCAGACTTTGCAGCCTGGACTTTTCAATTCAAACAAGGGAGATAGCTCCCTTTTg                                                                                                                                         <  4:781379/149‑1 (MQ=255)
             tttttCACCACTCATTGCCACATTCCTTGTGTATAGCCAGCCATTTTTTACGGGCACAGCCAAACTTTACCGTGCCCTAATACGACAAAAGCCCAGACTTTGCAGCCTGGACTTTTCAATTCAAACAAGGGAGATAGCTCCCTTTTGGc                                                                                                                                       <  2:803710/149‑1 (MQ=255)
               tttCACCACTCATTGCCACATTCCTTGTGTATAGCCAGCCATTTTTTACGGGCACAGCCAAACTTTACCGTGCCCTAATACGACAAAAGCCCAGACTTTGCAGCCTGGACTTTTCAATTCAAACAAGGGAGCTATCTCCCttgtttgaa                                                                                                                                     <  3:1391537/149‑8 (MQ=255)
               tttCACCACTCATTGCCACATTCCTTGTGTATAGCCAGCCATTTTTTACGGGCACAGCCAAACTTTACCGTGCCCTAATACGACAAAAGCCCAGACTTTGCAGCCTGGACTTTTCAATTCAAACAAGGGAGATAGCTCCCTTTTGGCAt                                                                                                                                     <  3:1734879/149‑1 (MQ=255)
                ttCACCACTCATTGCCACATTCCTTGTGTATAGCCAGCCATTTTTTACGGGCACAGCCAAACTTTACCGTGCCCTAATACGACAAAAGCCCAGACTTTGCAGCCTGGACTTTTCAATTCAAACAAGGGAGATAGCTCCCTTTTGGCATg                                                                                                                                    >  2:663655/1‑149 (MQ=255)
                 tCACCACTCATTGCCACATTCCTTGTGTATAGCCAGCCATTTTTTACGGGTACAGCCAAACTTTACCGTGCCCTAATACGACAAAAGCCCAGACTTTGCAGCCTGGACTTTTCAATTCAAACAAGGGAGCTATCTCCCtt                                                                                                                                            <  1:1283632/140‑1 (MQ=255)
                 tCACCACTCATTGCCACATTCCTTGTGTATAGCCAGCCATTTTTTACGGGTACAGCCAAACTTTACCGTGCCCTAATACGACAAAAGCCCAGACTTTGCAGCCTGGACTTTTCAATTCAAACAAGGGAGCTATCTCCCtt                                                                                                                                            >  2:1283632/1‑140 (MQ=255)
                  caccacTCATTGCCACATTCCTTGTGTATAGCCAGCCATTTTTTACGGGCACAGCCAAACTTTACCGTGCCCTAATACGACAAAAGCCCAGACTTTGCAGCCTGGACTTTTCAATTCAAACAAGGGAGATAGCTCCCtttt                                                                                                                                          >  3:875899/1‑141 (MQ=255)
                  caccacTCATTGCCACATTCCTTGTGTATAGCCAGCCATTTTTTACGGGCACAGCCAAACTTTACCGTGCCCTAATACGACAAAAGCCCAGACTTTGCAGCCTGGACTTTTCAATTCAAACAAGGGAGATAGCTCCCtttt                                                                                                                                          <  4:875899/141‑1 (MQ=255)
                  caccacTCATTGCCACATTCCTTGTGTATAGCCAGCCATTTTTTACGGGCACAGCCAAACTTTACCGTGCCCTAATACGACAAAAGCCCAGACTTTGCAGCCTGGACTTTTCAATTCAAACAAGGGAGATAGCTCCCTTTTGGCATgaa                                                                                                                                  <  3:1133482/149‑1 (MQ=255)
                    ccacTCATTGCCACATTCCTTGTGTATAGCCAGCCATTTTTTACGGGCACAGCCAAACTTTACCGTGCCCTAATACGACAAAAGCCCAGACTTTGCAGCCTGGACTTTTCAATTCAAACAAGGGAGATa                                                                                                                                                    <  4:1427981/129‑1 (MQ=255)
                    ccacTCATTGCCACATTCCTTGTGTATAGCCAGCCATTTTTTACGGGCACAGCCAAACTTTACCGTGCCCTAATACGACAAAAGCCCAGACTTTGCAGCCTGGACTTTTCAATTCAAACAAGGGAGATAGCTCCCTTTTGGCATgaa                                                                                                                                  >  3:560855/1‑147 (MQ=255)
                    ccacTCATTGCCACATTCCTTGTGTATAGCCAGCCATTTTTTACGGGCACAGCCAAACTTTACCGTGCCCTAATACGACAAAAGCCCAGACTTTGCAGCCTGGACTTTTCAATTCAAACAAGGGAGATAGCTCCCTTTTGGCATgaa                                                                                                                                  <  4:560855/147‑1 (MQ=255)
                    ccacTCATTGCCACATTCCTTGTGTATAGCCAGCCATTTTTTACGGGCACAGCCAAACTTTACCGTGCCCTAATACGACAAAAGCCCAGACTTTGCAGCCTGGACTTTTCAATTCAAACAAGGGAGATAGCTCCCTTTTGGCATgaa                                                                                                                                  <  2:1884274/147‑1 (MQ=255)
                     cacTCATTGCCACATTCCTTGTGTATAGCCAGCCATTTTTTACGGGCACAGCCAAACTTTACCGTGCCCTAATACGACAAAAGCCCAGACTTTGCAGCCTGGACTTTTCAATTCAAACAAGGGAGATAGCTCCCTTTTGGCATgaagaa                                                                                                                               <  4:375049/149‑1 (MQ=255)
                       cTCATTGCCACATTCCTTGTGTATAGCCAGCCATTTTTTACGGGCACAGCCAAACTTTACCGTGCCCTAATACGACAAAAGCCCAGACTTTGCAGCCTGGACTTTTCAATTCAAACAAGGGAGATAGCTc                                                                                                                                                >  3:810277/1‑130 (MQ=255)
                           ttGCCACATTCCTTGTGTATAGCCAGCCATTTTTTACGGGCACAGCCAAACTTTACCGTGCCCTAATACGACAAAAGCCCAGACTTTGCAGCCTGGACTTTTCAATTCAAACAAGGGAGATAGCTCCCTTTTGGCATGAAGAAGTaaaa                                                                                                                         <  4:612502/149‑1 (MQ=255)
                                  aTTCCTTGTGTATAGCCAGCCATTTTTTACGGGCACAGCCAAACTTTACCGTGCCCTAATACGACAAAAGCCCAGACTTTGCAGCCTGGACTTTTCAATTCAAACAAGGGAGATAGCTCCCTTTTGGCATGAAGAAGTAAAattatt                                                                                                                    <  3:34019/147‑1 (MQ=255)
                                  aTTCCTTGTGTATAGCCAGCCATTTTTTACGGGCACAGCCAAACTTTACCGTGCCCTAATACGACAAAAGCCCAGACTTTGCAGCCTGGACTTTTCAATTCAAACAAGGGAGATAGCTCCCTTTTGGCATGAAGAAGTAAAattatt                                                                                                                    >  4:34019/1‑147 (MQ=255)
                                   ttCCTTGTGTATAGCCAGCCATTTTTTACGGGCACAGCCAAACTTTACCGTGCCCTAATACGACAAAAGCCCAGACTTTGCAGCCTGGACTTTTCAATTCAAACAAGGGAGATAGCTCCCTTTTGGCATGAAGAAGTAAAATTAttctt                                                                                                                 >  4:1518543/1‑149 (MQ=255)
                                       ttGTGTATAGCCAGCCATTTTTTACGGGCACAGCCAAACTTTACCGTGCCCTAATACGACAAAAGCCCGGACTTTGCAGCCTGGACTTTTCAATTCAAACAAGGGAGATAGCTCCCtttt                                                                                                                                          <  2:797487/120‑1 (MQ=255)
                                       ttGTGTATAGCCAGCCATTTTTTACGGGCACAGCCAAACTTTACCGTGCCCTAATACGACAAAAGCCCAGACTTTGCAGCCTGGACTTTTCAATTCAAACAAGGGAGATAGCTCCCtttt                                                                                                                                          >  1:797487/1‑120 (MQ=255)
                                       ttGTGTATAGCCAGCCATTTTTTACGGGCACAGCCAAACTTTACCGTGCCCTAATACGACAAAAGCCCAGACTTTGCAGCCTGGACTTTTCAATTCAAACAAGGGAGATAGCTCCCTTTTGGCATGAAGAAGTAAAATTAttcttcttc                                                                                                             >  4:798883/1‑149 (MQ=255)
                                           gtATAGCCAGCCATTTTTTACGGGCACAGCCAAACTTTACCGTGCCCTAATACGACAAAAGCCCAGACTTTGCAGCCTGGACTTTTCAATTCAAACAAGGGAGATAGCTCCCTTTTGGCAt                                                                                                                                     <  1:96155/121‑1 (MQ=255)
                                               agccagccATTTTTTACGGGCACAGCCAAACTTTACCGTGCCCTAATACGACAAAAGCCCAGACTTTGCAGCCTGGACTTTTCAATTCAAACAAGGGAGATAGCTCCCTTTTg                                                                                                                                         >  4:1530483/1‑113 (MQ=255)
                                               agccagccATTTTTTACGGGCACAGCCAAACTTTACCGTGCCCTAATACGACAAAAGCCCAGACTTTGCAGCCTGGACTTTTCAATTCAAACAAGGGAGATAGCTCCCTTTTg                                                                                                                                         <  3:1530483/113‑1 (MQ=255)
                                                   agccaTTTTTTACGGGCACAGCCAAACTTTACCGTGCCCTAATACGACAAAAGCCCAGACTTTGCAGCCTGGACTTTTCAATTCAAACAAGGGAGATAGCTCCCTTTTGGCATgaagaa                                                                                                                               <  1:314536/119‑1 (MQ=255)
                                                      caTTTTTTACGGGCACAGCCAAACTTTACCGTGCCCTAATACGACAAAAGCCCAGACTTTGCAGCCTGGACTTTTCAATTCAAACAAGGGAGATAGCTCCCTTTTGGCATGAAGAAGTAAAattatt                                                                                                                    >  4:1423052/1‑127 (MQ=255)
                                                      caTTTTTTACGGGCACAGCCAAACTTTACCGTGCCCTAATACGACAAAAGCCCAGACTTTGCAGCCTGGACTTTTCAATTCAAACAAGGGAGATAGCTCCCTTTTGGCATGAAGAAGTAAAattatt                                                                                                                    <  3:1423052/127‑1 (MQ=255)
                                                       aTTTTTTACGGGCACAGCCAAACTTTACCGTGCCCTAATACGACAAAAGCCCAGACTTTGCAGCCTGGACTTTTCAATTCAAACAAGGGAGATAGCTCCCTTTTGGCATGAAGAAGTAAAattat                                                                                                                     <  3:525719/125‑1 (MQ=255)
                                                       aTTTTTTACGGGCACAGCCAAACTTTACCGTGCCCTAATACGACAAAAGCCCAGACTTTGCAGCCTGGACTTTTCAATTCAAACAAGGGAGATAGCTCCCTTTTGGCATGAAGAAGTAAAattat                                                                                                                     >  4:525719/1‑125 (MQ=255)
                                                       aTTTTTTACGGGCACAGCCAAACTTTACCGTGCCCTAATACGACAAAAGCCCAGACTTTGCAGCCTGGACTTTTCAATTCAAACAAGGGAGATAGCTCCCTTTTGGCATGAAGAAGTAAAATTATTCTTCTTCTGGCTCGTCGTCAACg                                                                                             <  1:1856947/149‑1 (MQ=255)
                                                        ttttttACGGGCACAGCCAAACTTTACCGTGCCCTAATACGACAAAAGCCCAGACTTTGCAGCCTGGACTTTTCAATTCAAACAAGGGAGATAGCTCCCTTTTGGCATGAAGAAGTAAAATTATTCTTCTTCTGGCtcgtc                                                                                                    <  3:1507309/141‑1 (MQ=255)
                                                                  gCACAGCCAAACTTTACCGTGCCCTAATACGACAAAAGCCCAGACTTTGCAGCCTGGACTTTTCAATTCAAACAAGGGAGATAGCTCCCTTTTGGCATgaag                                                                                                                                 <  3:1020558/102‑1 (MQ=255)
                                                                  gCACAGCCAAACTTTACCGTGCCCTAATACGACAAAAGCCCAGACTTTGCAGCCTGGACTTTTCAATTCAAACAAGGGAGATAGCTCCCTTTTGGCATgaag                                                                                                                                 >  4:1020558/1‑102 (MQ=255)
                                                                  gCACAGCCAAACTTTACCGTGCCCTAATACGACAAAAGCCCAGACTTTGCAGCCTGGACTTTTCAATTCAAACAAGGGAGATAGCTCCCTTTTGGCATgaa                                                                                                                                  <  2:1535325/101‑1 (MQ=255)
                                                                        ccAAACTTTACCGTGCCCTAATACGACAAAAGCCCAGACTTTGCAGCCTGGACTTTTCAATTCAAACAAGGGAGATAGCTCCCTTTTGGCATGAAGAAGTAAAATTATTCTTCTTCTGGCTCGTCGTCAACGTCCACTTCCGGAGCGAt                                                                            <  3:1518543/149‑1 (MQ=255)
                                                                         cAAACTTTACCGTGCCCTAATACGACAAAAGCCCAGACTTTGCAGCCTGGACTTTTCAATTCAAACAAGGGAGATAGCTCCCTTTTGGCATGAAGAAGTAAAATTAttctt                                                                                                                 <  3:1866633/111‑1 (MQ=255)
                                                                         cAAACTTTACCGTGCCCTAATACGACAAAAGCCCAGACTTTGCAGCCTGGACTTTTCAATTCAAACAAGGGAGATAGCTCCCTTTTGGCATGAAGAAGTAAAATTAttctt                                                                                                                 >  4:1866633/1‑111 (MQ=255)
                                                                          aaaCTTTACCGTGCCCTAATACGACAAAAGCCCAGACTTTGCAGCCTGGACTTTTCAATTCAAACAAGGGAGATAGCTCCCTTTTGGCATGAAGAAGTAAAATTATTCTTCTTCTGGCTCGTCGTCAACGTCCa                                                                                         <  2:1294533/134‑1 (MQ=255)
                                                                          aaaCTTTACCGTGCCCTAATACGACAAAAGCCCAGACTTTGCAGCCTGGACTTTTCAATTCAAACAAGGGAGATAGCTCCCTTTTGGCATGAAGAAGTAAAATTATTCTTCTTCTGGCTCGTCGTCAACGTCCa                                                                                         >  1:1294533/1‑134 (MQ=255)
                                                                             cTTTACCGTGCCCTAATACGACAAAAGCCCAGACTTTGCAGCCTGGACTTTTCAATTCAAACAAGGGAGCTATCTCCCttgt                                                                                                                                          >  4:1546912/1‑80 (MQ=255)
                                                                              tttACCGTGCCCTAATACGACAAAAGCCCAGACTTTGCAGCCTGGACTTTTCAATTCAAACAAGGGAGATAGCTCCCTTTTGGCATGAAGAAGTAAAATTATTCTTCTTCTGGCTCGTCGTCAACGTCCACTTCCGGAGCGATTTCATc                                                                      <  1:663655/149‑1 (MQ=255)
                                                                                   cGTGCCCTAATACGACAAAAGCCCAGACTTTGCAGCCTGGACTTTTCAATTCAAACAAGGGAGATAGCTCCCTTTTGGCATgaagaa                                                                                                                               <  1:1723789/87‑1 (MQ=255)
                                                                                   cGTGCCCTAATACGACAAAAGCCCAGACTTTGCAGCCTGGACTTTTCAATTCAAACAAGGGAGATAGCTCCCTTTTGGCATgaagaa                                                                                                                               >  2:1723789/1‑87 (MQ=255)
                                                                                        ccTAATACGACAAAAGCCCAGACTTTGCAGCCTGGACTTTTCAATTCAAACAAGGGAGATAGCTCCCtttt                                                                                                                                          <  4:1199452/71‑1 (MQ=255)
                                                                                        ccTAATACGACAAAAGCCCAGACTTTGCAGCCTGGACTTTTCAATTCAAACAAGGGAGATAGCTCCCTTTTGGCATGAAGAAGTAAAATTATTCTTCTTCTGGCTCGTCGTCAACGTCCACTTCCGGAGCGATTTCATCGTc                                                                   <  1:266994/142‑1 (MQ=255)
                                                                                        ccTAATACGACAAAAGCCCAGACTTTGCAGCCTGGACTTTTCAATTCAAACAAGGGAGATAGCTCCCTTTTGGCATGAAGAAGTAAAATTATTCTTCTTCTGGCTCGTCGTCAACGTCCACTTCCGGAGCGATTTCATCGTc                                                                   >  2:266994/1‑142 (MQ=255)
                                                                                          tAATACGACAAAAGCCCAGACTTTGCAGCCTGGACTTTTCAATTCAAACAAGGGAGATAGCTCCCtt                                                                                                                                            >  4:385951/1‑67 (MQ=255)
                                                                                          tAATACGACAAAAGCCCAGACTTTGCAGCCTGGACTTTTCAATTCAAACAAGGGAGATAGCTCCCtt                                                                                                                                            <  3:385951/67‑1 (MQ=255)
                                                                                                         ccAGACTTTGCAGCCTGGACTTTTCAATTCAAACAAGGGAGATAGCTCCCTTTTGGCATGAAGAAGTaaaa                                                                                                                         >  3:416289/1‑71 (MQ=255)
                                                                                                         ccAGACTTTGCAGCCTGGACTTTTCAATTCAAACAAGGGAGATAGCTCCCTTTTGGCATGAAGAAGTAAAATTATTCTTCTTCTGGCTCGTCGTCAACGTCCACTTCCGGAGCGATTTCATCGTCCCCTTCCGCGGCACTGCCGTCGAt                                           >  1:1065397/1‑149 (MQ=255)
                                                                                                         ccAGACTTTGCAGCCTGGACTTTTCAATTCAAACAAGGGAGATAGCTCCCTTTTGGCATGAAGAAGTAAAATTATTCTTCTTCTGGCTCGTCGTCAACGTCCACTTCCGGAGCGATTTCATCGTCCCCTTCCGCGGCACTGCCGTCGAt                                           >  2:1129248/1‑149 (MQ=255)
                                                                                                         ccAGACTTTGCAGCCTGGACTTTTCAATTCAAACAAGGGAGATAGCTCCCTTTTGGCATGAAGAAGTAAAATTATTCTTCTTCTGGCTCGTCGTCAACGTCCACTTCCGGAGCGATTTCATCGTCCCCTTCCGCGGCACTGCCGTCGAt                                           >  1:643796/1‑149 (MQ=255)
                                                                                                          cAGACTTTGCAGCCTGGACTTTTCAATTCAAACAAGGGAGATAGCTCCCTTTTGGCATGAAGAAGTAAAATTATTCTTCTTCTGGCTCGTCGTCAACGTCCACTTCCGGAGCGATTTCATCGTCCCCTTCCGCGGCACTGCCGTc                                              <  3:900233/145‑1 (MQ=255)
                                                                                                          cAGACTTTGCAGCCTGGACTTTTCAATTCAAACAAGGGAGATAGCTCCCTTTTGGCATGAAGAAGTAAAATTATTCTTCTTCTGGCTCGTCGTCAACGTCCACTTCCGGAGCGATTTCATCGTCCCCTTCCGCGGCACTGCCGTCGATg                                          >  3:195405/1‑149 (MQ=255)
                                                                                                           aGACTTTGCAGCCTGGACTTTTCAATTCAAACAAGGGAGATAGCTCCCTTTTGGCATGAAGAAGTAAAATTATTCTTCTTCTGGCTCGTCGTCAACGTCCACTTCCGGAGCGATTTCATCGTCCCCTTCCGCGGCACTGCCGTCGa                                            <  3:697070/146‑1 (MQ=255)
                                                                                                            gACTTTGCAGCCTGGACTTTTCAATTCAAACAAGGGAGATAGCTCCCTTTTgg                                                                                                                                        >  3:1001467/1‑53 (MQ=255)
                                                                                                                   cAGCCTGGACTTTTCAATTCAAACAAGGGAGATAGCTCCCTTTTGGCATGAAGAAGTAAAATTATTCTTCTTCTGGCTCGTCGTCAACGTCCACTTCCGGAGCGATTTCATCGTCCCCTTCCGCGGCACTGCCGTCGATGGTATCCAGa                                 <  2:643796/149‑1 (MQ=255)
                                                                                                                     gtcTGGGCTTTTCAATTCAAACAAGGGAGATAGCTCCCTTTTGGCATGAAGAAGTAAAATTATTCTTCTTCTGGCTCGTCGTCAACGTCCACTTCCGGAGCGATTTCATCGTCCCCTTCCGCGGCACTGCCGTCGATGGTATCCAGATc                               <  4:304751/147‑1 (MQ=255)
                                                                                                                       cTGGACTTTTCAATTCAAACAAGGGAGATAGCTCCCTTTTGGCATGAAGAAGTAAAATTATTCTTCTTCTGGCTCGTCGTCAACGTCCACTTCCGGAGCGATTTCATCGTcc                                                                  <  4:1476957/112‑1 (MQ=255)
                                                                                                                       cTGGACTTTTCAATTCAAACAAGGGAGATAGCTCCCTTTTGGCATGAAGAAGTAAAATTATTCTTCTTCTGGCTCGTCGTCAACGTCCACTTCCGGAGCGATTTCATCGTcc                                                                  >  3:1476957/1‑112 (MQ=255)
                                                                                                                        tGGACTTTTCAATTCAAACAAGGGAGATAGCTCCCTTTTGGCATGAAGAAGTAAAATTAttcttct                                                                                                               >  1:1285261/1‑66 (MQ=255)
                                                                                                                        tGGACTTTTCAATTCAAACAAGGGAGATAGCTCCCTTTTGGCATGAAGAAGTAAAATTATTCTTCTTCTGGCTTGTCGTCAACGTCCACTTCCGGAGCGATTTCATCGTCCCCTTCCGCGGCACTGCCGTCGATGGTATCCAGATCTTc                            <  1:1146729/149‑1 (MQ=255)
                                                                                                                        tGGACTTTTCAATTCAAACAAGGGAGATAGCTCCCTTTTGGCATGAAGAAGTAAAATTATTCTTCTTCTGGCTCGTCGTCAACGTCCACTTCCGGAGCGATTTCATCGTCCCCTTCCGCGGCACTGCCGTCGATGGt                                        <  3:698040/137‑1 (MQ=255)
                                                                                                                        tGGACTTTTCAATTCAAACAAGGGAGATAGCTCCCTTTTGGCATGAAGAAGTAAAATTATTCTTCTTCTGGCTCGTCGTCAACGTCCACTTCCGGAGCGATTTCATCGTCCCCTTCCGCGGCACTGCCGTCGATGGt                                        >  4:698040/1‑137 (MQ=255)
                                                                                                                        tGGACTTTTCAATTCAAACAAGGGAGATAGCTCCCTTTTGGCATGAAGAAGTAAAATTATTCTTCTTCTGGCTCGTCGTCAACGTCCACTTCCGGAGCGATTTCATCGTCCCCTTCCGCGGCACTGCCGTCGATGGTATCCAGAt                                >  1:1769986/1‑145 (MQ=255)
                                                                                                                        tGGACTTTTCAATTCAAACAAGGGAGATAGCTCCCTTTTGGCATGAAGAAGTAAAATTATTCTTCTTCTGGCTCGTCGTCAACGTCCACTTCCGGAGCGATTTCATCGTCCCCTTCCGCGGCACTGCCGTCGATGGTATCCAGAt                                <  2:1769986/145‑1 (MQ=255)
                                                                                                                        tGGACTTTTCAATTCAAACAAGGGAGATAGCCCCCTTTTGGCATGAAGAAGTAAAATTATTCTTCTTCTGGCTCGTCGTCAACGTCCACTTCCGGAGCGATTTCATCGTcc                                                                  >  4:1641759/1‑111 (MQ=255)
                                                                                                                        tGGACTTTTCAATTCAAACAAGGGAGATAGCCCCCTTTTGGCATGAAGAAGTAAAATTATTCTTCTTCTGGCTCGTCGTCAACGTCCACTTCCGGAGCGATTTCATCGTcc                                                                  <  3:1641759/111‑1 (MQ=255)
                                                                                                                          gACTTTTCAATTCAAACAAGGGAGATAGCTCCCTTTTGGCATGAAGAAGTAAAATTATTCTTCTTCTGGCTCGTCGTCAACGTCCACTTCCGGAGCGATTTCATCGTCCCCTTCCGCGGCACTGCCGTCGATGGTATCCAGATCTTCCt                          <  2:1065397/149‑1 (MQ=255)
                                                                                                                             ttttCAATTCAAACAAGGGAGATAGCTCCCTTTTGGCATGAAGAAGTAAAATTATTCTTCTTCTGGCTCGTCGTCAACGTCCACTTCCGGAGCGATTTCATCGTc                                                                   <  4:54813/105‑1 (MQ=255)
                                                                                                                               ttcaattcaaACAAGGGAGATAGCTCCCTTTTGGCATGAAGAAGTAAAATTATTCTTCTTCTGGCTCGTCGTCAACGTCCACTTCCGGAGCGATTTCATCGTCCCCTTCCGCGGCACTGCCGTCGATGGTATCCAGATCTTCCTCGTCa                     >  1:462344/1‑149 (MQ=255)
                                                                                                                                   attcaaACAAGGGAGATAGCTCCCTTTTGGCATGAAGAAGTAAAATTATTCTTCTTCTGGCTCGTCGTCAACGTCCACTTCCGGAGCGATTTCATCGTCCCCTTCCGCGGCACTGCCGTCGATGGTATCCAGATCTTCCTCGTCAACCg                 >  2:1408487/1‑149 (MQ=255)
                                                                                                                                    ttcaaACAAGGGAGATAGCTCCCTTTTGGCATGAAGAAGTAAAATTATTCTTCTTCTGGCTCGTCGTCAACGTCCACTTCCGGAGCGATTTCATCGTCCCCTTCCGCGGCACTGCCGTCGATg                                          <  2:1272117/123‑1 (MQ=255)
                                                                                                                                    ttcaaACAAGGGAGATAGCTCCCTTTTGGCATGAAGAAGTAAAATTATTCTTCTTCTGGCTCGTCGTCAACGTCCACTTCCGGAGCGATTTCATCGTCCCCTTCCGCGGCACTGCCGTCGATGGTATCCAGATCTTCCTCGTCAACCgg                >  4:887666/1‑149 (MQ=255)
                                                                                                                                    ttcaaACAAGGGAGATAGCTCCCTTTTGGCATGAAGAAGTAAAATTATTCTTCTTCTGGCTCGTCGTCAACGTCCACTTCCGGAGCGATTTCATCGTCCCCTTCCGCGGCACTGCCGTCGATGGTATCCAGATCTTCCTCGTCAACCgg                >  4:1843812/1‑149 (MQ=255)
                                                                                                                                    ttcaaACAAGGGAGATAGCTCCCTTTTGGCATGAAGAAGTAAAATTATTCTTCTTCTGGCTCGTCGTCAACGTCCACTTCCGGAGCGATTTCATCGTCCCCTTCCGCGGCACTGCCGTCGATGGTATCCAGATCTTCCTCGTCAACCgg                <  1:1129248/149‑1 (MQ=255)
                                                                                                                                    ttcaaACAAGGGAGATAGCTCCCTTTTGGCATGAAGAAGTAAAATTATTCTTCTTCTGGCTCGTCGTCAACGTCCACTTCCAGAGCGATTTCATCGTCCCCTTCCGCGGCACTGCCGTCGATGGTATCCAGATCTTc                            >  3:891665/1‑137 (MQ=255)
                                                                                                                                    ttcaaACAAGGGAGATAGCTCCCTTTTGGCATGAAGAAGTAAAATTATTCTTCTTCTGGCTCGTCGTCAACGTCCACTTCCAGAGCGATTTCATCGTCCCCTTCCGCGGCACTGCCGTCGATGGTATCCAGATCTTc                            <  4:891665/137‑1 (MQ=255)
                                                                                                                                       aaaCAAGGGAGATAGCTCCCTTTTGGCATGAAGAAGTAAAATTATTCTTCTTCTGGCTCGTCGTCAACGTCCACTTCCGGAGCGATTTCATCGTCCCCTTCCGCGGCACTGCCGTCGATGGTATCCAGATCTTCCTCGTCAACCGGTTc             >  2:1768970/1‑149 (MQ=255)
                                                                                                                                       aaaCAAGGGAGATAGCTCCCTATTGGCATGAAGAAGTAAAATTATTCTTCTTCTGGCTCGTCGTCAACGTCCACTTCCGGAGCGATTCCATCGTCCCCTTCCGCGGCACTGCCGTCGATGGTATCCAGATCTTcc                           <  3:1291494/135‑1 (MQ=255)
                                                                                                                                         aaaaGGGAGCTATCTCCCTTTTGGCATGAAGAAGTAAAATTATTCTTCTTCTGGCTCGTCGTCAACGTCCACTTCCGGAGCGATTTCATCGTCCCCTTCCGCGGCACTGCCGTCGATGGTATCCAGATCTTCCTCGTCAACCGGTTCa            >  3:305018/3‑148 (MQ=255)
                                                                                                                                           aaGGGAGATAGCTCCCTTTTGGCATGAAGAAGTAAAATTATTCTTCTTCTGGCTCGTCGTCAACGTcc                                                                                          <  2:433735/68‑1 (MQ=255)
                                                                                                                                           aaGGGAGATAGCTCCCTTTTGGCATGAAGAAGTAAAATTATTCTTCTTCTGGCTCGTCGTCAACGTcc                                                                                          <  2:1244793/68‑1 (MQ=255)
                                                                                                                                           aaGGGAGATAGCTCCCTTTTGGCATGAAGAAGTAAAATTATTCTTCTTCTGGCTCGTCGTCAACGTcc                                                                                          >  1:1244793/1‑68 (MQ=255)
                                                                                                                                           aaGGGAGATAGCTCCCTTTTGGCATGAAGAAGTAAAATTATTCTTCTTCTGGCTCGTCGTCAACGTCCACTTCCGGAGCGATTTCATCGTc                                                                   >  2:647128/1‑91 (MQ=255)
                                                                                                                                           aaGGGAGATAGCTCCCTTTTGGCATGAAGAAGTAAAATTATTCTTCTTCTGGCTCGTCGTCAACGTCCACTTCCGGAGCGATTTCATCGTc                                                                   <  1:647128/91‑1 (MQ=255)
                                                                                                                                           aaGGGAGATAGCTCCCTTTTGGCATGAAGAAGTAAAATTATTCTTCTTCTGGCTCGTCGTCAACGTCCACTTCCGGAGCGATTTCATCGTCCCCTTCCGCGGCACTGCCGTCGATGGTATCCAGATCTTCCTCGTCAACCGGTTCAGCa         >  2:187723/1‑149 (MQ=255)
                                                                                                                                             gggAGATAGCTCCCTTTTGGCATGAAGAAGTAAAATTATTCTTCTTCTGGCTCGTCGTCAACGTCCACTTCCGGAGCGATTTCATCGTCCCCTTCCGCGGCACTGCCGTCGATGGTATCCAGATCTTCCTCg                        >  3:1032617/1‑132 (MQ=255)
                                                                                                                                             gggAGATAGCTCCCTTTTGGCATGAAGAAGTAAAATTATTCTTCTTCTGGCTCGTCGTCAACGTCCACTTCCGGAGCGATTTCATCGTCCCCTTCCGCGGCACTGCCGTCGATGGTATCCAGATCTTCCTCg                        <  4:1032617/132‑1 (MQ=255)
                                                                                                                                             gggAGATAGCTCCCTTTTGGCATGAAGAAGTAAAATTATTCTTCTTCTGGCTCGTCGTCAACGTCCACTTCCGGAGCGATTTCATCGTCCCCTTCCGCGGCACTGCCGTCGATGGTATCCAGATCTTCCTCGTCAACCGGTTCAGCAac       >  1:391392/1‑149 (MQ=255)
                                                                                                                                               gagaTAGCTCCCTTTTGGCATGAAGAAGTAAAATTATTCTTCTTCTGGCTCGTCGTCAACGTCCACTTCCGGAGCGATTTCATCGTCCCCTTCCGCGGCACTGCCGTCGATGGTATCCAGATCTTCCTCGTCAACCGGTTCAGCAacac     >  1:904668/1‑149 (MQ=255)
                                                                                                                                                 gaTAGCTCCCTTTTGGCATGAAGAAGTAAAATTATTCTTCTTCTGGCTCGTCGTCAACGTCCACTTCCGGAGCGATTTCATCGTCCCCTTCCGCGGCACTGCCGTCGATGGTATCCAGATCTTc                            <  2:850812/124‑1 (MQ=255)
                                                                                                                                                 gaTAGCTCCCTTTTGGCATGAAGAAGTAAAATTATTCTTCTTCTGGCTCGTCGTCAACGTCCACTTCCGGAGCGATTTCATCGTCCCCTTCCGCGGCACTGCCGTCGATGGTATCCAGATCTTc                            >  1:850812/1‑124 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
AGCCTTTCGCCTTTTTTTCACCACTCATTGCCACATTCCTTGTGTATAGCCAGCCATTTTTTACGGGCACAGCCAAACTTTACCGTGCCCTAATACGACAAAAGCCCAGACTTTGCAGCCTGGACTTTTCAATTCAAACAAGGGAGATAGCTCCCTTTTGGCATGAAGAAGTAAAATTATTCTTCTTCTGGCTCGTCGTCAACGTCCACTTCCGGAGCGATTTCATCGTCCCCTTCCGCGGCACTGCCGTCGATGGTATCCAGATCTTCCTCGTCAACCGGTTCAGCAACACGTT  >  NC_000913/2336617‑2336911

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