Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,336,763 A→C 10.0% 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,7630AC10.0% 222.7 / 18.1 70intergenic (‑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 (34/29);  new base C (3/4);  total (37/33)
Fisher's exact test for biased strand distribution p-value = 6.99e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.58e-01

CAGCCTTTCGCCTTTTTTTCACCACTCATTGCCACATTCCTTGTGTATAGCCAGCCATTTTTTACGGGCACAGCCAAACTTTACCGTGCCCTAATACGACAAAAGCCCAGACTTTGCAGCCTGGACTTTTCAATTCAAACAAGGGAGATAGCTCCCTTTTGGCATGAAGAAGTAAAATTATTCTTCTTCTGGCTCGTCGTCAACGTCCACTTCCGGAGCGATTTCATCGTCCCCTTCCGCGGCACTGCCGTCGATGGTATCCAGATCTTCCTCGTCAACCGGTTCAGCAACACGTT  >  NC_000913/2336616‑2336911
                                                                                                                                                    |                                                                                                                                                    
cAGCCTTTCGCCTTTTTTTCACCACTCATTGCCACATTCCTTGTGTATAGCCAGCCATTTTTTACGGGCACAGCCAAACTTTACCGTGCCCTAATACGACAAAAGCCCAGACTTTGCAGCCTGGACTTTTCAATTCAAACAAGGGAGAt                                                                                                                                                     <  1:47899/149‑1 (MQ=255)
          ccTTTTTTTCACCACTCATTGCCACATTCCTTGTGTATAGCCAGCCATTTTTTACGGGCACAGCCAAACTTTACCGTGCCCTAATACGACAAAAGCCCAGACTTTGCAGCCTGGACTTTTCAATTCAAACAAGGGAGATAGCTCCCttt                                                                                                                                           <  2:260147/149‑1 (MQ=255)
           cTTTTTTTCACCACTCATTGCCACATTCCTTGTGTATAGCCAGCCATTTTTTACGGGCACAGCCAAACTTTACCGTGCCCTAATACGACAAAAGCCCAGACTTTGCAGCCTGGACTTTTCAATTCAAACAAGGGAGATAGCTCCCtttt                                                                                                                                          <  4:656996/149‑1 (MQ=255)
            tttttttCACCACTCATTGCCACATTCCTTGTGTATAGCCAGCCATTTTTTACGGGCACAGCCAAACTTTACCGTGCCCTAATACGACAAAAGCCCAGACTTTGCAGCCTGGACTTTTCAATTCAAACAAGGGAGCTATCTCCCttgtt                                                                                                                                         <  3:289701/149‑4 (MQ=255)
            tttttttCACCACTCATTGCCACATTCCTTGTGTATAGCCAGCCATTTTTTACGGGCACAGCCAAACTTTACCGTGCCCTAATACGACAAAAGCCCAGACTTTGCAGCCTGGACTTTTCAATTCAAACAAGGGAGCTATCTCCCttgtt                                                                                                                                         >  3:30931/1‑146 (MQ=255)
            tttttttCACCACTCATTGCCACATTCCTTGTGTATAGCCAGCCATTTTTTACGGGCACAGCCAAACTTTACCGTGCCCTAATACGACAAAAGCCCAGACTTTGCAGCCTGGACTTTTCAATTCAAACAAGGGAGCTATCTCCCttgtt                                                                                                                                         >  3:184550/1‑146 (MQ=255)
            tttttttCACCACTCATTGCCACATTCCTTGTGTATAGCCAGCCATTTTTTACGGGCACAGCCAAACTTTACCGTGCCCTAATACGACAAAAGCCCAGACTTTGCAGCCTGGACTTTTCAATTCAAACAAGGGAGATAGCTCCCTTTTg                                                                                                                                         >  3:300956/1‑149 (MQ=255)
            tttttttCACCACTCATTGCCACATTCCTTGTGTATAGCCAGCCATTTTTTACGGGCACAGCCAAACTTTACCGTGCCCTAATACGACAAAAGCCCAGACTTTGCAGCCTGGACTTTTCAATTCAAACAAGGGAGATAGCTCCCTTTTg                                                                                                                                         >  1:359728/1‑149 (MQ=255)
            tttttttCACCACTCATTGCCACATTCCTTGTGTATAGCCAGCCATTTTTTACGGGCACAGCCAAACTTTACCGTGCCCTAATACGACAAAAGCCCAGACTTTGCAGCCTGGACTTTTCAATTCAAACAAGGGAGATAGCTCCCTTTTg                                                                                                                                         >  2:482152/1‑149 (MQ=255)
             ttttttCACCACTCATTGCCACATTCCTTGTGTATAGCCAGCCATTTTTTTACGGGCACAGCCAAACTTTACCGTGCCCTAATACGACAAAAGCCCAGACTTTGCAGCCTGGACTTTTCAATTCAAACAAGGGAGCTATCTCCCttgtt                                                                                                                                         <  1:390180/149‑4 (MQ=255)
                tttCACCACTCATTGCCACATTCCTTGTGTATAGCCAGCCATTTTTTACGGGCACAGCCAAACTTTACCGTGCCCTAATACGACAAAAGCCCAGACTTTGCAGCCTGGACTTTTCAATTCAAACAAGGGAGCTATCTCCCttgtttgaa                                                                                                                                     <  1:3662/149‑8 (MQ=255)
                tttCACCACTCATTGCCACATTCCTTGTGTATAGCCAGCCATTTTTTACGGGCACAGCCAAACTTTACCGTGCCCTAATACGACAAAAGCCCAGACTTTGCAGCCTGGACTTTTCAATTCAAACAAGGGAGCTATCTCCCttgtttgaa                                                                                                                                     <  4:30931/149‑8 (MQ=255)
                      cacTCATTGCCACATTCCTTGTGTATAGCCAGCCATTTTTTACGGGCACAGCCAAACTTTACCGTGCCCTAATACGACAAAAGCCCAGACTTTGCAGCCTGGACTTTTCAATTCAAACAAGGGAGATAGCTCCCTTTTGGCATgaagaa                                                                                                                               <  2:671291/149‑1 (MQ=255)
                        cTCATTGCCACATTCCTTGTGTATAGCCAGCCATTTTTTACGGGCACAGCCAAACTTTACCGTGCCCTAATACGACAAAAGCCCAGACTTTGCAGCCTGGACTTTTCAATTCAAACAAGGGAGATAGCTCCCTTTTGGCATgaagaa                                                                                                                               <  2:629116/147‑1 (MQ=255)
                        cTCATTGCCACATTCCTTGTGTATAGCCAGCCATTTTTTACGGGCACAGCCAAACTTTACCGTGCCCTAATACGACAAAAGCCCAGACTTTGCAGCCTGGACTTTTCAATTCAAACAAGGGAGATAGCTCCCTTTTGGCATgaagaa                                                                                                                               >  1:629116/1‑147 (MQ=255)
                           aTTGCCACATTCCTTGTGTATAGCCAGCCATTTTTTACGGGCACAGCCAAACTTTACCGTGCCCTAATACGACAAAAGCCCAGACTTTGCAGCCTGGACTTTTCAATTCAAACAAGGGAGATAGCTCCCTTTTGGCATGAAGAAGTaaa                                                                                                                          >  4:491477/1‑149 (MQ=255)
                             tGCCACATTCCTTGTGTATAGCCAGCCATTTTTTACGGGCACAGCCAAACTTTACCGTGCCCTAATACGACAAAAGCCCAGACTTTGCAGCCTGGACTTTTCAATTCAAACAAGGGAGATAGCTCCCTTTTGGCATGAAGAAGTAAAat                                                                                                                        <  1:482152/149‑1 (MQ=255)
                              gCCACATTCCTTGTGTATAGCCAGCCATTTTTTACGGGCACAGCCAAACTTTACCGTGCCCTAATACGACAAAAGCCCAGACTTTGCAGCCTGGACTTTTCAATTCAAACAAGGGAGATAGCTCCCTTTTGGCATGAAGAAGTAAAatt                                                                                                                       >  2:416454/1‑149 (MQ=255)
                               ccACATTCCTTGTGTATAGCCAGCCATTTTTTACGGGCACAGCCAAACTTTACCGTGCCCTAATACGACAAAAGCCCAGACTTTGCAGCCTGGACTTTTCAATTCAAACAAGGGAGATAGCTCCCTTTTGGCATgaa                                                                                                                                  <  3:572847/137‑1 (MQ=255)
                               ccACATTCCTTGTGTATAGCCAGCCATTTTTTACGGGCACAGCCAAACTTTACCGTGCCCTAATACGACAAAAGCCCAGACTTTGCAGCCTGGACTTTTCAATTCAAACAAGGGAGATAGCTCCCTTTTGGCATgaa                                                                                                                                  >  4:572847/1‑137 (MQ=255)
                                cacaTTCCTTGTGTATAGCCAGCCATTTTTTACGGGCACAGCCAAACTTTACCGTGCCCTAATACGACAAAAGCCCAGACTTTGCAGCCTGGACTTTTCAATTCAAACAAGGGAGATAGCTCCCTTTTGGCATGAAGAAGTAAAattat                                                                                                                     <  3:491477/149‑1 (MQ=255)
                                   aTTCCTTGTGTATAGCCAGCCATTTTTTACGGGCACAGCCAAACTTTACCGTGCCCTAATACGACAAAAGCCCAGACTTTGCAGCCTGGACTTTTCAATTCAAACAAGGGAGATAGCTCCCTTTTGGCATGAAGAAGTAAAATTAttct                                                                                                                  >  2:318251/1‑149 (MQ=255)
                                            gtATAGCCAGCCATTTTTTACGGGCACAGCCAAACTTTACCGTGCCCTAATACGACAAAAGCCCAGACTTTGCAGCCTGGACTTTTCAATTCAAACAAGGGAGATAGCTCCCTTTTGGCATGAAGAAGTAAAATTATTCTTCTTCTGGc                                                                                                         >  2:66474/1‑149 (MQ=255)
                                                  ccagccaTTTTTTACGGGCACAGCCAAACTTTACCGTGCCCTAATACGACAAAAGCCCAGACTTTGCAGCCTGGACTTTTCAATTCAAACAAGGGAGATAGCTCCCTTTTGGCATGAAGAAGTAAAATTATTCTTCTTCTGGCtcgtcg                                                                                                   <  4:300956/149‑1 (MQ=255)
                                                  ccagccaTTTTTTACGGGCACAGCCAAACTTTACCGTGCCCTAATACGACAAAAGCCCAGACTTTGCAGCCTGGACTTTTCAATTCAAACAAGGGAGATAGCTCCCTTTTGGCATGAAGAAGTAAAATTATTCTTCTTCTGGCtcgtcg                                                                                                   >  2:275683/1‑149 (MQ=255)
                                                    agccaTTTTTTACGGGCACAGCCAAACTTTACCGTGCCCTAATACGACAAAAGCCCAGACTTTGCAGCCTGGACTTTTCAATTCAAACAAGGGAGATAGCTCCCTTTTGGCATGAAGAAGTAAAATTATTCTTCTTCTGGCTCgtcgtc                                                                                                 <  3:443340/149‑1 (MQ=255)
                                                       caTTTTTTACGGGCACAGCCAAACTTTACCGTGCCCTAATACGACAAAAGCCCAGACTTTGCAGCCTGGACTTTTCAATTCAAACAAGGGAGATAGCTCCCTTTTGGCATGAAGAAGTAAAATTAttcttcttc                                                                                                             <  1:289822/134‑1 (MQ=255)
                                                       caTTTTTTACGGGCACAGCCAAACTTTACCGTGCCCTAATACGACAAAAGCCCAGACTTTGCAGCCTGGACTTTTCAATTCAAACAAGGGAGATAGCTCCCTTTTGGCATGAAGAAGTAAAATTAttcttcttc                                                                                                             >  2:289822/1‑134 (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                                                                   <  3:374082/149‑1 (MQ=255)
                                                                                   aCCGTGCCCTAATACGACAAAAGCCCAGACTTTGCAGCCTGGACTTTTCAATTCAAACAAGGGAGATAGCTCCCTTTTGGCATGAAGAAGTAAAATTATTCTTCTTCTGGCTCGTCGTCAACGTCCACTTCCGGAGCGATTTCATCGTc                                                                   <  1:339853/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)
                                                                                                           ccAGACTTTGCAGCCTGGACTTTTCAATTCAAACAAGGGAGAt                                                                                                                                                     >  2:334173/1‑43 (MQ=255)
                                                                                                           ccAGACTTTGCAGCCTGGACTTTTCAATTCAAACAAGGGAGAt                                                                                                                                                     <  1:334173/43‑1 (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)
                                                                                                                                                    |                                                                                                                                                    
CAGCCTTTCGCCTTTTTTTCACCACTCATTGCCACATTCCTTGTGTATAGCCAGCCATTTTTTACGGGCACAGCCAAACTTTACCGTGCCCTAATACGACAAAAGCCCAGACTTTGCAGCCTGGACTTTTCAATTCAAACAAGGGAGATAGCTCCCTTTTGGCATGAAGAAGTAAAATTATTCTTCTTCTGGCTCGTCGTCAACGTCCACTTCCGGAGCGATTTCATCGTCCCCTTCCGCGGCACTGCCGTCGATGGTATCCAGATCTTCCTCGTCAACCGGTTCAGCAACACGTT  >  NC_000913/2336616‑2336911

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: