Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,336,766 G→T 10.1% intergenic (‑122/+27) 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,7660GT10.1% 236.2 / 14.8 69intergenic (‑122/+27)yfaA/gyrADUF2138 family protein, putative host defense protein/DNA gyrase (type II topoisomerase), subunit A
Reads supporting (aligned to +/- strand):  ref base G (34/28);  new base T (3/4);  total (37/32)
Fisher's exact test for biased strand distribution p-value = 6.96e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.27e-01

CCTTTTTTTCACCACTCATTGCCACATTCCTTGTGTATAGCCAGCCA‑TTTTTTACGGGCACAGCCAAACTTTACCGTGCCCTAATACGACAAAAGCCCAGACTTTGCAGCCTGGACTTTTCAATTCAAACAAGGGAGATAGCTCCCTTTTGGCATGAAGAAGTAAAATTATTCTTCTTCTGGCTCGTCGTCAACGTCCACTTCCGGAGCGATTTCATCGTCCCCTTCCGCGGCACTGCCGTCGATGGTATCCAGATCTTCCTCGTCAACCGGTTCAGCAACACGTTGCA  >  NC_000913/2336626‑2336914
                                                                                                                                             |                                                                                                                                                    
ccTTTTTTTCACCACTCATTGCCACATTCCTTGTGTATAGCCAGCCA‑TTTTTTACGGGCACAGCCAAACTTTACCGTGCCCTAATACGACAAAAGCCCAGACTTTGCAGCCTGGACTTTTCAATTCAAACAAGGGAGATAGCTCCCttt                                                                                                                                              <  2:260147/149‑1 (MQ=255)
 cTTTTTTTCACCACTCATTGCCACATTCCTTGTGTATAGCCAGCCA‑TTTTTTACGGGCACAGCCAAACTTTACCGTGCCCTAATACGACAAAAGCCCAGACTTTGCAGCCTGGACTTTTCAATTCAAACAAGGGAGATAGCTCCCtttt                                                                                                                                             <  4:656996/149‑1 (MQ=255)
  tttttttCACCACTCATTGCCACATTCCTTGTGTATAGCCAGCCA‑TTTTTTACGGGCACAGCCAAACTTTACCGTGCCCTAATACGACAAAAGCCCAGACTTTGCAGCCTGGACTTTTCAATTCAAACAAGGGAGCTATCTCCCttgtt                                                                                                                                            <  3:289701/149‑4 (MQ=255)
  tttttttCACCACTCATTGCCACATTCCTTGTGTATAGCCAGCCA‑TTTTTTACGGGCACAGCCAAACTTTACCGTGCCCTAATACGACAAAAGCCCAGACTTTGCAGCCTGGACTTTTCAATTCAAACAAGGGAGCTATCTCCCttgtt                                                                                                                                            >  3:30931/1‑146 (MQ=255)
  tttttttCACCACTCATTGCCACATTCCTTGTGTATAGCCAGCCA‑TTTTTTACGGGCACAGCCAAACTTTACCGTGCCCTAATACGACAAAAGCCCAGACTTTGCAGCCTGGACTTTTCAATTCAAACAAGGGAGCTATCTCCCttgtt                                                                                                                                            >  3:184550/1‑146 (MQ=255)
  tttttttCACCACTCATTGCCACATTCCTTGTGTATAGCCAGCCA‑TTTTTTACGGGCACAGCCAAACTTTACCGTGCCCTAATACGACAAAAGCCCAGACTTTGCAGCCTGGACTTTTCAATTCAAACAAGGGAGATAGCTCCCTTTTg                                                                                                                                            >  1:359728/1‑149 (MQ=255)
  tttttttCACCACTCATTGCCACATTCCTTGTGTATAGCCAGCCA‑TTTTTTACGGGCACAGCCAAACTTTACCGTGCCCTAATACGACAAAAGCCCAGACTTTGCAGCCTGGACTTTTCAATTCAAACAAGGGAGATAGCTCCCTTTTg                                                                                                                                            >  3:300956/1‑149 (MQ=255)
  tttttttCACCACTCATTGCCACATTCCTTGTGTATAGCCAGCCA‑TTTTTTACGGGCACAGCCAAACTTTACCGTGCCCTAATACGACAAAAGCCCAGACTTTGCAGCCTGGACTTTTCAATTCAAACAAGGGAGATAGCTCCCTTTTg                                                                                                                                            >  2:482152/1‑149 (MQ=255)
   ttttttCACCACTCATTGCCACATTCCTTGTGTATAGCCAGCCATTTTTTTACGGGCACAGCCAAACTTTACCGTGCCCTAATACGACAAAAGCCCAGACTTTGCAGCCTGGACTTTTCAATTCAAACAAGGGAGCTATCTCCCttgtt                                                                                                                                            <  1:390180/149‑4 (MQ=255)
      tttCACCACTCATTGCCACATTCCTTGTGTATAGCCAGCCA‑TTTTTTACGGGCACAGCCAAACTTTACCGTGCCCTAATACGACAAAAGCCCAGACTTTGCAGCCTGGACTTTTCAATTCAAACAAGGGAGCTATCTCCCttgtttgaa                                                                                                                                        <  1:3662/149‑8 (MQ=255)
      tttCACCACTCATTGCCACATTCCTTGTGTATAGCCAGCCA‑TTTTTTACGGGCACAGCCAAACTTTACCGTGCCCTAATACGACAAAAGCCCAGACTTTGCAGCCTGGACTTTTCAATTCAAACAAGGGAGCTATCTCCCttgtttgaa                                                                                                                                        <  4:30931/149‑8 (MQ=255)
            cacTCATTGCCACATTCCTTGTGTATAGCCAGCCA‑TTTTTTACGGGCACAGCCAAACTTTACCGTGCCCTAATACGACAAAAGCCCAGACTTTGCAGCCTGGACTTTTCAATTCAAACAAGGGAGATAGCTCCCTTTTGGCATgaagaa                                                                                                                                  <  2:671291/149‑1 (MQ=255)
              cTCATTGCCACATTCCTTGTGTATAGCCAGCCA‑TTTTTTACGGGCACAGCCAAACTTTACCGTGCCCTAATACGACAAAAGCCCAGACTTTGCAGCCTGGACTTTTCAATTCAAACAAGGGAGATAGCTCCCTTTTGGCATgaagaa                                                                                                                                  >  1:629116/1‑147 (MQ=255)
              cTCATTGCCACATTCCTTGTGTATAGCCAGCCA‑TTTTTTACGGGCACAGCCAAACTTTACCGTGCCCTAATACGACAAAAGCCCAGACTTTGCAGCCTGGACTTTTCAATTCAAACAAGGGAGATAGCTCCCTTTTGGCATgaagaa                                                                                                                                  <  2:629116/147‑1 (MQ=255)
                 aTTGCCACATTCCTTGTGTATAGCCAGCCA‑TTTTTTACGGGCACAGCCAAACTTTACCGTGCCCTAATACGACAAAAGCCCAGACTTTGCAGCCTGGACTTTTCAATTCAAACAAGGGAGATAGCTCCCTTTTGGCATGAAGAAGTaaa                                                                                                                             >  4:491477/1‑149 (MQ=255)
                   tGCCACATTCCTTGTGTATAGCCAGCCA‑TTTTTTACGGGCACAGCCAAACTTTACCGTGCCCTAATACGACAAAAGCCCAGACTTTGCAGCCTGGACTTTTCAATTCAAACAAGGGAGATAGCTCCCTTTTGGCATGAAGAAGTAAAat                                                                                                                           <  1:482152/149‑1 (MQ=255)
                    gCCACATTCCTTGTGTATAGCCAGCCA‑TTTTTTACGGGCACAGCCAAACTTTACCGTGCCCTAATACGACAAAAGCCCAGACTTTGCAGCCTGGACTTTTCAATTCAAACAAGGGAGATAGCTCCCTTTTGGCATGAAGAAGTAAAatt                                                                                                                          >  2:416454/1‑149 (MQ=255)
                     ccACATTCCTTGTGTATAGCCAGCCA‑TTTTTTACGGGCACAGCCAAACTTTACCGTGCCCTAATACGACAAAAGCCCAGACTTTGCAGCCTGGACTTTTCAATTCAAACAAGGGAGATAGCTCCCTTTTGGCATgaa                                                                                                                                     >  4:572847/1‑137 (MQ=255)
                     ccACATTCCTTGTGTATAGCCAGCCA‑TTTTTTACGGGCACAGCCAAACTTTACCGTGCCCTAATACGACAAAAGCCCAGACTTTGCAGCCTGGACTTTTCAATTCAAACAAGGGAGATAGCTCCCTTTTGGCATgaa                                                                                                                                     <  3:572847/137‑1 (MQ=255)
                      cacaTTCCTTGTGTATAGCCAGCCA‑TTTTTTACGGGCACAGCCAAACTTTACCGTGCCCTAATACGACAAAAGCCCAGACTTTGCAGCCTGGACTTTTCAATTCAAACAAGGGAGATAGCTCCCTTTTGGCATGAAGAAGTAAAattat                                                                                                                        <  3:491477/149‑1 (MQ=255)
                         aTTCCTTGTGTATAGCCAGCCA‑TTTTTTACGGGCACAGCCAAACTTTACCGTGCCCTAATACGACAAAAGCCCAGACTTTGCAGCCTGGACTTTTCAATTCAAACAAGGGAGATAGCTCCCTTTTGGCATGAAGAAGTAAAATTAttct                                                                                                                     >  2:318251/1‑149 (MQ=255)
                                  gtATAGCCAGCCA‑TTTTTTACGGGCACAGCCAAACTTTACCGTGCCCTAATACGACAAAAGCCCAGACTTTGCAGCCTGGACTTTTCAATTCAAACAAGGGAGATAGCTCCCTTTTGGCATGAAGAAGTAAAATTATTCTTCTTCTGGc                                                                                                            >  2:66474/1‑149 (MQ=255)
                                        ccagcca‑TTTTTTACGGGCACAGCCAAACTTTACCGTGCCCTAATACGACAAAAGCCCAGACTTTGCAGCCTGGACTTTTCAATTCAAACAAGGGAGATAGCTCCCTTTTGGCATGAAGAAGTAAAATTATTCTTCTTCTGGCtcgtcg                                                                                                      >  2:275683/1‑149 (MQ=255)
                                        ccagcca‑TTTTTTACGGGCACAGCCAAACTTTACCGTGCCCTAATACGACAAAAGCCCAGACTTTGCAGCCTGGACTTTTCAATTCAAACAAGGGAGATAGCTCCCTTTTGGCATGAAGAAGTAAAATTATTCTTCTTCTGGCtcgtcg                                                                                                      <  4:300956/149‑1 (MQ=255)
                                          agcca‑TTTTTTACGGGCACAGCCAAACTTTACCGTGCCCTAATACGACAAAAGCCCAGACTTTGCAGCCTGGACTTTTCAATTCAAACAAGGGAGATAGCTCCCTTTTGGCATGAAGAAGTAAAATTATTCTTCTTCTGGCTCgtcgtc                                                                                                    <  3:443340/149‑1 (MQ=255)
                                             ca‑TTTTTTACGGGCACAGCCAAACTTTACCGTGCCCTAATACGACAAAAGCCCAGACTTTGCAGCCTGGACTTTTCAATTCAAACAAGGGAGATAGCTCCCTTTTGGCATGAAGAAGTAAAATTAttcttcttc                                                                                                                >  2:289822/1‑134 (MQ=255)
                                             ca‑TTTTTTACGGGCACAGCCAAACTTTACCGTGCCCTAATACGACAAAAGCCCAGACTTTGCAGCCTGGACTTTTCAATTCAAACAAGGGAGATAGCTCCCTTTTGGCATGAAGAAGTAAAATTAttcttcttc                                                                                                                <  1:289822/134‑1 (MQ=255)
                                                  ttttACGGGCACAGCCAAACTTTACCGTGCCCTAATACGACAAAAGCCCAGACTTTGCAGCCTGGACTTTTCAATTCAAACAAGGGAGATAGCTCCCTTTTGGCATGAAGAAGTAAAATTATTCTTCTTCTGGCTCGTCGTCAACGTcc                                                                                             <  2:359728/149‑1 (MQ=255)
                                                        gggCACAGCCAAACTTTACCGTGCCCTAATACGACAAAAGCCCAGACTTTGCAGCCTGGACTTTTCAATTCAAACAAGGGAGATAGCTCCCTTTTGGCATGAAGAAGTAAAATTATTCTTCTTCTGGCTCGTCGTCAACGTCCACTTcc                                                                                       <  1:653383/149‑1 (MQ=255)
                                                          gCACAGCCAAACTTTACCGTGCCCTAATACGACAAAAGCCCAGACTTTGCAGCCTGGACTTTTCAATTCAAACAAGGGAGATAGCTCCCTTTTg                                                                                                                                            <  3:397428/94‑1 (MQ=255)
                                                          gCACAGCCAAACTTTACCGTGCCCTAATACGACAAAAGCCCAGACTTTGCAGCCTGGACTTTTCAATTCAAACAAGGGAGATAGCTCCCTTTTg                                                                                                                                            >  4:397428/1‑94 (MQ=255)
                                                            acaGCCAAACTTTACCGTGCCCTAATACGACAAAAGCCCAGACTTTGCAGCCTGGACTTTTCAATTCAAACAAGGGAGATAGCTCCCTTTTGGCATGAAGAAGTaaaa                                                                                                                            <  2:668984/108‑1 (MQ=255)
                                                            acaGCCAAACTTTACCGTGCCCTAATACGACAAAAGCCCAGACTTTGCAGCCTGGACTTTTCAATTCAAACAAGGGAGATAGCTCCCTTTTGGCATGAAGAAGTaaaa                                                                                                                            >  1:668984/1‑108 (MQ=255)
                                                                ccAAACTTTACCGTGCCCTAATACGACAAAAGCCCAGACTTTGCAGCCTGGACTTTTCAATTCAAACAAGGGAGATAGCTCCCTTTTGGCATGAAGAAGTAAAATTATTCTTCTTCTGGCTCGTCGTCAACGTCCACTTCCGGAGCGAt                                                                               >  2:336167/1‑149 (MQ=255)
                                                                 cAAACTTTACCGTGCCCTAATACGACAAAAGCCCAGACTTTGCAGCCTGGACTTTTCAATTCAAACAAGGGAGATAGCTCCCTTTTGGCATGAAGAAGTAAAATTATTCTTCTTCTGGCTCGTCGTCAACGTCCACTTCCGGAGCGAtt                                                                              >  4:182205/1‑149 (MQ=255)
                                                                   aaCTTTACCGTGCCCTAATACGACAAAAGCCCAGACTTTGCAGCCTGGACTTTTCAATTCAAACAAGGGAGATAGCTCCCTTTTGGCATGAAGAAGTAAAATTATTCTTCTTCTGGCTCGTCGTCAACGTCCACTTCCGGAGCGATTTc                                                                            >  4:45639/1‑149 (MQ=255)
                                                                         aCCGTGCCCTAATACGACAAAAGCCCAGACTTTGCAGCCTGGACTTTTCAATTCAAACAAGGGAGATAGCTCCCTTTTGGCATGAAGAAGTAAAATTATTCTTCTTCTGGCTCGTCGTCAACGTCCACTTCCGGAGCGATTTCATCGTc                                                                      <  1:339853/149‑1 (MQ=255)
                                                                         aCCGTGCCCTAATACGACAAAAGCCCAGACTTTGCAGCCTGGACTTTTCAATTCAAACAAGGGAGATAGCTCCCTTTTGGCATGAAGAAGTAAAATTATTCTTCTTCTGGCTCGTCGTCAACGTCCACTTCCGGAGCGATTTCATCGTc                                                                      <  3:374082/149‑1 (MQ=255)
                                                                             tGCCCTAATACGACAAAAGCCCAGACTTTGCAGCCTGGACTTTTCAATTCAAACAAGGGAGATAGCTCCCTTTTGGCATGAAGAAGTAAAATTATTCTTCTTCTGGCTCGTCGTCAACGTCCACTTCCGGAGCGATTTCATCGTCCCCt                                                                  >  2:48433/1‑149 (MQ=255)
                                                                              gCCCTAATACGACAAAAGCCCAGACTTTGCAGCCTGGACTTTTCAATTCAAACAAGGGAGATAGCTCCCTTTTGGCATGAAGAAGTAAAATTATTCTTCTTCTGGCTCGTCGTCAACGTCCACTTCCGGAGCGATTTCATCGTCCCCtt                                                                 <  1:230744/149‑1 (MQ=255)
                                                                                 cTAATACGACAAAAGCCCAGACTTTGCAGCCTGGACTTTTCAATTCAAACAAGGGAGATAGCTCCCTTTTGGCATGAAGAAGTAAAATTATTCTTCTTCTGGCTCGTCGTCAACGTCCACTTCCGGAGCGATTTCATCGTCCCCTTCcg                                                              >  4:466354/1‑149 (MQ=255)
                                                                                    aTACGACAAAAGCCCAGACTTTGCAGCCTGGACTTTTCAATTCAAACAAGGGAGATAGCTCCCTTTTGGCATGAAGAAGTAAAATTATTCTTCTTCTGGCTCGTCGTCAACGTCCACTTCCGGAGCGATTTCATCGTCCCCTTCCGCgg                                                           >  1:286719/1‑149 (MQ=255)
                                                                                        gACAAAAGCCCAGACTTTGCAGCCTGGACTTTTCAATTCAAACAAGGGAGATAGCTCCCTTTTGGCATGAAGAAGTAAAATTATTCTTCTTCTGGCTCGTCGTCAACGTCCACTTCCTGAGCGATTTCATCGTCCCCTTCCGCGGCACt                                                       >  4:473280/1‑149 (MQ=255)
                                                                                          cAAAAGCCCAGACTTTGCAGCCTGGACTTTTCAATTGAAACAAGAGAGATAGCTCCCTTTTGAAATGAAGAAGTAAAATTATTCTTCTTCTGGCTCGTCGTCAACGTCCACTTCCGAAGCGATTTCATCGTCCCCTTCCGCGGCACTGc                                                     >  3:431601/1‑149 (MQ=255)
                                                                                           aaaaGCCCAGACTTTGCAGCCTGGACTTTTCAATTCAAACAAGGGAGATAGCTCCCTTTTGGCATGAAGAAGTAAAATTATTCTTCTTCTGGCTCGTCGTCAACGTCCACTTCCGGAGCGATTTCATCGTCCCCTTCCGCGGCACTGcc                                                    >  4:425185/1‑149 (MQ=255)
                                                                                                  cAGACTTTGCAGCCTGGACTTTTCAATTCAAACAAGGGAGATAGCTCCCTTTTGGCATGAAGAAGTAAAATTATTCTTCTTCTGGCTCGTCGTCAACGTCCACTTCCGGAGCGATTTCATCGTCCCCTTCCGCGGCACTGCCGTCGATg                                             >  3:244650/1‑149 (MQ=255)
                                                                                                  cAGACTTTGCAGCCTGGACTTTTCAATTCAAACAAGGGAGATAGCTCCCTTTTGGCATGAAGAAGTAAAATTATTCTTCTTCTGGCTCGTCGTCAACGTCCACTTCCGGAGCGATTTCATAGTCCCCTTCCGCGGCACTGCCGTCGATg                                             <  3:45639/149‑1 (MQ=255)
                                                                                                   aGACTTTGCAGCCTGGACTTTTCAATTCAAACAAGGGAGATAGCTCCCTTTTGGCATGAAGAAGTAAAATTATTCTTCTTCTGGCTCGTCGTCAACGTCCACTTCCGGAGCGATTTCATCGTCCCCTTCCGCGGCACTGCCGTCGATgg                                            >  3:47319/1‑149 (MQ=255)
                                                                                                        ttGCAGCCTGGACTTTTCAATTCAAACAAGGGAGATAGCTCCCTTTTGGCATGAAGAAGTAAAATTATTCTTCTTCTGGCTCGTCGTCAACGTCCACTTCCGGAGCGATTTCATCGTCCCCTTCCGCGGCACTGCCGTCGATGGTATcc                                       <  1:307759/149‑1 (MQ=255)
                                                                                                             gCCTGGACTTTTCAATTCAAACAAGGGAGATAGCTCCCTTTTGGCATGAAGAAGTAAAATTATTCTTCTTCTGGCTCGTCGTCAACGTCCACTTCCGGAGCGATTTCATCGTCCCCTTCCGCGGCACTGCCGTCGATGGTATCCAGATc                                  <  3:182205/149‑1 (MQ=255)
                                                                                                               cTGGACTTTTCAATTCAAACAAGGGAGATAGCTCCCTTTTGGCATGAAGAAGTAAAATTATTCTTCTTCTGGCTCGTCGTCAACGTCCACTTCCGGAGCGATTTCATCGTCCCCTTCCGCGGCACTGCCGTCGATGGTATCCAGATCtt                                >  3:194250/1‑149 (MQ=255)
                                                                                                               cTGGACTTTTCAATTCAAACAAGGGAGATAGCTCCCTTTTGGCATGAAGAAGTAAAATTATTCTTCTTCTGGCTCGTCGTCAACGTCCACTTCCGGAGCGATTTCATCGTCCCCTTCCGCGGCACTGCCGTCGATGGTATCCAGATCtt                                >  4:648398/1‑149 (MQ=255)
                                                                                                                     ttttCAATTCAAACAAGGGAGATAGCTCCCTTTTGGCATGAAGAAGTAAAATTATTCTTCTTCTGGCTCGTCGTCAACGTCCACTTCCGGAGCGATTTCATCGTCCCCTTCCGCGGCACTGCCGTCGATg                                             >  1:167753/1‑130 (MQ=255)
                                                                                                                     ttttCAATTCAAACAAGGGAGATAGCTCCCTTTTGGCATGAAGAAGTAAAATTATTCTTCTTCTGGCTCGTCGTCAACGTCCACTTCCGGAGCGATTTCATCGTCCCCTTCCGCGGCACTGCCGTCGATg                                             <  2:167753/130‑1 (MQ=255)
                                                                                                                       ttcaattcaaACAAGGGAGATAGCTCCCTTTTGGCATGAAGAAGTAAAATTATTCTTCTTCTGGCTCGTCGTCAACGTCCACTTCCGGAGCGATTTCATCGTCCCCTTCCGCGGCACTGCCGTCGATGGTATCCAGATCTTCCTCGTCa                        >  2:217780/1‑149 (MQ=255)
                                                                                                                        tcaattcaaACAAGGGAGATAGCTCCCTTTTGGCATGAAGAAGTAAAATTATTCTTCTTCTGGCTCGTCGTCAACGTCCACTTCCGGAGCGATTTCATCGTCCCCTTCCGCGGCACTGCCGTCGATGGTATCCAGATCTTCCTCGTCaa                       >  4:622933/1‑149 (MQ=255)
                                                                                                                         caattcaaACAAGGGAGATAGCTCCCTTTTGGCATGAAGAAGTAAAATTATTCTTCTTCTGGCTCGTCGTCAACGTCCACTTCCGGAGCGATTTCATCGTCCCCTTCCGCGGCACTGCCGTCGATGGTATCCAGATCTTCCTCGTCAAc                      <  3:425185/149‑1 (MQ=255)
                                                                                                                             tcaaACAAGGGAGATAGCTCCCTTTTGGCATGAAGAAGTAAAATTATTCTTCTTCTGGCTCGTCGTCAACGTCCACTTCCGGAGCGATTTCATCGTCCCCCTCCGCGGCACTGCCGTCGATGGTATCCAGATCTTCCTCGTCAACCGGt                  <  4:503524/149‑1 (MQ=255)
                                                                                                                              caaACAAGGGAGATAGCTCCCTTTTGGCATGAAGAAGTAAAATTATTCTTCTTCTGGCTCGTCGTCAACGTCCACTTCCGGAGCGATTTCATCGTCCCCTTCCGCGGCACTGCCGTCGATGGTATCCAGATCTTCCTCGTCAACCGGtt                 <  1:66474/149‑1 (MQ=255)
                                                                                                                               aaaCAAGGGAGATAGCTCCCTTTTGGCATGAAGAAGTAAAATTATTCTTCTTCTGGCTCGTCGTCAACGTCCACTTCCGGAGCGATTTCATCGTCCCCTTCCGCGGCACTGCCGTCGATGGTATCCAGATCTTCCTCGTCAACCGGTTc                >  3:177569/1‑149 (MQ=255)
                                                                                                                                aaCAAGGGAGATAGCTCCCTTTTGGCATGAAGAAGTAAAATTATTCTTCTTCTGGCTCGTCGTCAACGTcc                                                                                             >  3:596663/1‑71 (MQ=255)
                                                                                                                                aaCAAGGGAGATAGCTCCCTTTTGGCATGAAGAAGTAAAATTATTCTTCTTCTGGCTCGTCGTCAACGTcc                                                                                             <  4:596663/71‑1 (MQ=255)
                                                                                                                                 aaaaGGGAGCTATCTCCCTTTTGGCATGAAGAAGTAAAATTATTCTTCTTCTGGCTCGTCGTCAACGTCCACTTCCGGAGCGATTTCATCGTCCCCTTCCGCGGCACTGCCGTCGATGGTATCCAGATCTTCCTCGTCAACCGGTTCAg              >  4:473274/3‑149 (MQ=255)
                                                                                                                                 aCAAGGGAGATAGCTCCCTTTTGGCATGAAGAAGTAAAATTATTCTTCTTCTGGCTCGTCGTCAACGTCCACTTCCGGAGCGATTTCATCGTCCCCTTCCGCGGCACTGCCGTCGAt                                              <  1:350189/117‑1 (MQ=255)
                                                                                                                                 aCAAGGGAGATAGCTCCCTTTTGGCATGAAGAAGTAAAATTATTCTTCTTCTGGCTCGTCGTCAACGTCCACTTCCGGAGCGATTTCATCGTCCCCTTCCGCGGCACTGCCGTCGAt                                              >  2:350189/1‑117 (MQ=255)
                                                                                                                                     gggAGATAGCTCCCTTTTGGCATGAAGAAGTAAAATTATTCTTCTTCTGGCTCGTCGTCAACGTCCACTTCCGGAGCGATTTCATCGTCCCCTTCCGCGGCACTGCCGTCGATGGTATCCAGATCTTCCTCGTCAACCGGTTCAGCAac          >  2:80518/1‑149 (MQ=255)
                                                                                                                                      ggAGATAGCTCCCTTTTGGCATGAAGAAGTAAAATTATTCTTCTTCTGGCTCGTCGTCAACGTCCACTTCCGGAGCGATTTCATCGTCCCCTTCCGCGGCACTGCCGTCGATGGTATCCAGATCTTCCTCGTCAACCGGTTCAGCAaca         <  4:244650/149‑1 (MQ=255)
                                                                                                                                          aTAGCTCCCTTTTGGCATGAAGAAGTAAAATTATTCTTCTTCTGGCTCGTCGTCAACGTCCACTTCCGGAGCGATTTCATCGTCCCCTTCCGCGGCACTGCCGTCGATGGTATCCAGATCTTCCTCGTCAACCGGTTCAGCAACACGtt     >  4:290864/1‑149 (MQ=255)
                                                                                                                                            aGCTCCCTTTTGGCATGAAGAAGTAAAATTATTCTTCTTCTGGCTCGTCGTCAACGTCCACTTCCGGAGCGATTTCATCGTCCCCTTCCGCGGCACTGCCGTCGATGGTATCCAGATCTTCCTCGTCAACCGGTTCAGCAACACGTTGc   <  1:275683/149‑1 (MQ=255)
                                                                                                                                             gCTCCCTTTTGGCATGAAGAAGTAAAATTATTCTTCTTCTGGCTCGTCGTCAACGTCCACTTCCGGAGCGATTTCATCGTCCCCTTCCGCGGCACTGCCGTCGATGGTATCCAGATCTTCCTCGTCAACCGGTTCAGCAACACGTTGCa  >  3:180095/1‑149 (MQ=255)
                                                                                                                                             |                                                                                                                                                    
CCTTTTTTTCACCACTCATTGCCACATTCCTTGTGTATAGCCAGCCA‑TTTTTTACGGGCACAGCCAAACTTTACCGTGCCCTAATACGACAAAAGCCCAGACTTTGCAGCCTGGACTTTTCAATTCAAACAAGGGAGATAGCTCCCTTTTGGCATGAAGAAGTAAAATTATTCTTCTTCTGGCTCGTCGTCAACGTCCACTTCCGGAGCGATTTCATCGTCCCCTTCCGCGGCACTGCCGTCGATGGTATCCAGATCTTCCTCGTCAACCGGTTCAGCAACACGTTGCA  >  NC_000913/2336626‑2336914

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: ‑