Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 246,747 G→T S12S (TCG→TCT yafL → NlpC/P60 family protein YafL

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913246,7470GT97.6% 146.7 / NA 42S12S (TCG→TCTyafLNlpC/P60 family protein YafL
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base T (22/19);  minor base C (1/0);  total (23/19)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

ATTCAGCTTGAATTGTAGCTCAAATGAGCAACAGCGACAATCCATGAGTAATTTAATAAACGGATGTTTTTACATAATTAACTGAAATATATATGTTAATTTTATAATAATGGTTTATTAAGTTTGTTCAGAAACGATCCGGGATACACTGCAACATTAAAGCATGACCAGCCATTATTGAGTATGTCCTTGCCGTCGATTCCATCCTTTGTATTGTCGGGATTACTGTTGATTTGTTTGCCGTTTTCTTCATTTGCCAGCGCCACCACATCACATATCTCTTTCAGCTACGCCGCCCGCCAGCGGATGCAAAACCGTGCGCGTTTATTAAAACAGTACCAAACTCATCTGAAAAAGCAGGCCAGCTATATTGTGGAAGGCAATGCCGAAAGCAAAAGGGCGCTACGCCAGCACAACCGGGAGCAGATAAAACAGCATCCA  >  NC_000913/246529‑246969
                                                                                                                                                                                                                            |                                                                                                                                                                                                                              
attCAGCTTGAATTGTAGCTCAAATGAGCAACAGCGACAATCCATGAGTAATTTAATAAACGGATGTTTTTACATAATTAACTGAAATATATATGTTAATTTTATAATAATGGTTTATTAAGTTTGTTCAGAAACGATCCGGGATACACTGCAACATTAAAGCATGACCAGCCATTATTGAGTATGTCCTTGCCGTCGATTCCATCCTTTGTATTGTCTGGATTACTGTTGATTTGTTTGCCGTTttcttc                                                                                                                                                                                                <  1:305935/251‑1 (MQ=255)
  tCAGCTTGAATTGTAGCTCAAATGAGCAACAGCGACAATCCATGAGTAATTTAATAAACGGATGTTTTTACATAATTAACTGAAATATATATGTTAATTTTATAATAATGGTTTATTAAGTTTGTTCAGAAACGATCCGGGATACACTGCAACATCAAAGCATGACCAGCCATTATTGAGTATGTCCTTGCCGTCGATTCCATCCTTTGTATTGTCTGGATTACTGTTGATTTGTTTGCCGTTTTCTTCa                                                                                                                                                                                               <  2:170601/250‑1 (MQ=255)
         gAATTGTAGCTCAAATGAGCAACAGCGACAATCCATGAGTAATTTAATAAACGGATGTTTTTACATAATTAACTGAAATATATATGTTAATTTTATAATAATGGTTTATTAAGTTTGTTCAGAAACGATCCGGGATACACTGCAACATTAAAGCATGACCAGCCATTATTGAGTATGTCCTTGCCGTCGATTCCATCCTTTGTATTGTCTGGATTACTGTTGATTTGTTTGCCGTTTTCTTCATTTGCCAg                                                                                                                                                                                       <  1:180655/251‑1 (MQ=255)
         gAATTGTAGCTCAAATGAGCAACAGCGACAATCCATGAGTAATTTAATAAACGGATGTTTTTACATAATTAACTGAAATATATATGTTAATTTTATAATAATGGTTTATTAAGTTTGTTCAGAAACGAGCCGGGATACACTGCAACATTAAAGCATGACCAGCCATTATTGAGTATGTCCTTGCCGTCGATTCCATCCTTTGTATTGTCTGGATTACTGTTGATTTGTTTGGCGGTTTCTTCACTTGCCa                                                                                                                                                                                        <  2:88811/250‑1 (MQ=255)
           aTTGTAGCTCAAATGAGCAACAGCGACAATCCATGAGTAATTTAATAAACGGATGTTTTTACATAATTAACTGAAATATATATGTTAATTTTATAATAATGGTTTATTAAGTTTGTTCAGAAACGATCCGGGATACACTGCAACATTAAAGCATGACCAGCCATTATTGAGTATGTCCTTGCCGGCGATTCCATCCTTTGTATTGTCTGGATTACTGTTGATTTGTTTGCCGTTTTCTTCATTTGCCAgcg                                                                                                                                                                                     <  1:70889/251‑1 (MQ=255)
                aGCTCAAATGAGCAACAGCGACAATCCATGAGTAATTTAATAAACGGATGTTTTTACATAATTAACTGAAATATATATGTTAATTTTATAATAATGGTTTATTAAGTTTGTTCAGAAACGATCCGGGATACACTGCAACATTAAAGCATGACCAGCCATTATTGAGTATGTCCTTGCCGTCGATTCCATCCTTTGTATTGTCTGGATTACTGTTGATTTGTTTGCCGTTTTCTTCAGTTGCCAGCGccac                                                                                                                                                                                 <  2:193782/250‑1 (MQ=255)
                           gCAACAGCGACAATCCATGAGTAATTTAATAAACGGATGTTTTTACATAATTAACTGAAATATATATGTTAATTTTATAATAATGGTTTATTAAGTTTGTTCAGAAACGATCCGGGATACACTGCAACATTAAAGCATGACCAGCCATTATTGAGTATGTCCTTGCCGTCGATTCCATCCTTTGTATTGTCTGGATTACTGTTGATTTGTTTGCCGTTTTCTTCATTTGCCAGCGCCACCACATCACatat                                                                                                                                                                     <  1:167791/251‑1 (MQ=255)
                                   gACAATCCATGAGTAATTTAATAAACGGATGTTTTTACATAATTAACTGAAATATATATGTTAATTTTATAATAATGGTTTATTAAGTTTGTTCAGAAACGATCCGGGATACACTGCAACATTAAAGCATGACCAGCCATTATTGAGTATGTCCTTGCCGTCGATTCCATCCTTTGTATTGTCTGGATTACTGTTGATTTGTTTGCCGTTTTCTTCATTTGCCAGCGCCACCACATCACATATCTCTTTCa                                                                                                                                                             <  1:425624/251‑1 (MQ=255)
                                     cAATCCATGAGTAATTTAATAAACGGATGTTTTTACATAATTAACTGAAATATATATGTTAATTTTATAATAATGGTTTATTAAGTTTGTTCAGAAACGATCCGGGATACACTGCAACATTAAAGCATGACCAGCCATTATTGAGTATGTCCTTGCCGTCGATTCCATCCTTTGTATTGTCTGGATTACTGTTGATTTGTTTGCCGTTTTCTTCATTTGCCAGCGCCACCACATCACATATCTCTTTCAGc                                                                                                                                                           <  2:157878/251‑1 (MQ=255)
                                        accATGAGTAATTTAATAAACGTATGTTTTTACATAATTAACTGAAATATATATGTTAATTTTATAATAATGGTTTATTAAGTTTGTTCAGAAACGATCCGGGATACACTGCAACATTAAAGCATGACCAGCCATTATTGAGTATGTCCTTGCCGTCGATTCCATCCTTTGTATTGTCTGGATTACTTTTGATTTGTTTGCCGTTTTCTTCATTTGCCAGCGCCACCACATCACATACCTCTTTCAGCTAc                                                                                                                                                        <  2:192356/250‑1 (MQ=255)
                                            tGAGTAATTTAATAAACGGATGTTTTTACATAATTAACTGAAATATATATGTTAATTTTATAATAATGGTTTATTAAGTTTGTTCAGAAACGATCCGGTATACACTGCAACATTAAAGCATGACCAGCCATTATTGAGTATGTCCTTGCCGTCGATTCCATCCTTTGTATTGTCTGGATTTCTGTTGCTTTGTTTGCCGTTTTCTTCATTTGCCAGCGCCACCACTTCACATATCTCTTTCTGCtcctccg                                                                                                                                                    >  2:315223/1‑245 (MQ=255)
                                                tactttAATACACGGATGTTTTTACATAATTAACTGAAATATATATGTTAATTTTATCATAATGGTTTATTAAGTTTGATCAGAAACGATCCGGGATAAACTGCAACATTAAAGCATGACAAGCCATTATTGGGTATGTCCTTGCCGTCGATTCCATCCTTTGTATTGTCTGGATTATTGTTGCTTTGTTTGCCGTTTTCTTCCTTTGACAGCGACACCAAATCACATATTTCTTTTAGCtcctctgcccg                                                                                                                                                >  2:383665/4‑241 (MQ=255)
                                                          aaCGGATGTTTTTACATACTTAACTGAAATATATATGTTAATTTTATAATAATGGTTTATTAAGTTTGTTAAGAAACGATCTGTTATACACTTCAACTTTAAAGCATGACCATCCATTATAGATTATATCCTTGCCGTCGATTCCATCCTTTGTATTGTCTGGGTTACTGTTGATTTGCTTGCCGTTTTCTTCATTTGCTAGCGCCACCACATCACTTATTTCTTTCATCTACGCTGTCCGCCAGCgcctc                                                                                                                                      >  2:3091/1‑247 (MQ=255)
                                                           aCGGATGTTTTTACATAATTAACTGAAATATATATGTTAATTTTATAATAATGGTTTATTAAGTTTGTTCAGAAACGATCCGGGATACACTGCAACATTAAAGCATGACCAGCCATTATTGAGTATGTCCTTGCCGTCGATTCCATCCTTTGTATTGTCTGGATTACTGTTGATTTGTTTGCCGTTTTCTTCATTTGCCAGCGCCACCACATCACATATCTCTTTCAGCTACGCCGCCCGCCAGCGGATGc                                                                                                                                     >  1:312048/1‑251 (MQ=255)
                                                                        cATAATTAACTGAAATATATATGTTAATTTTATAATAATGGTTTATTAAGTTTGTTCAGAAACGATCCGGGATACACTGCAACATTAAAGCATGACCAGCCATTATTGAGTATGTCCTTGCCGTCGATTCCATCCTTTGTATTGTCTGGATTACTGTTGATTTGTTTGCCGTTTTCTTCATTTGCCAGCGCCACCACATCACATATCTCTTTCAGCTACGCCGCCCGCCAGCGGATGCAAAACCGTGcgcg                                                                                                                        >  1:454211/1‑251 (MQ=255)
                                                                        cATAATTAACTGAAAGATATATGTTAATTTTATAATAATGGTTTATTAAGTTTGTTCAGAAACGATCCGGGATACACTGCAACATTAAAGCATGACCAGCCATTATTGAGTATGTCCTTGCCGTCGATTCCATCCTTTGTATTGTCTGGATTACTGTTGATTTGTTTGCCGTTTTCTTCATTTGCCAGCGCCACCACATCACATATCTCTTTCAGCTACGCCGCCCGCCAGCGGATGCAAAACCGTGcgcg                                                                                                                        <  1:498601/251‑1 (MQ=255)
                                                                               aaCTGAAATATATATGTTAATTTTATAATAATGGTTTATTAAGTTTGTTCAGAAACGATCCGGGATACACTGCAACATTAAAGCATGACCAGCCATTATTGCGTATGTCCTTGCCGTCGATTCCATCCTTTGTATTGTCTGGATTACTGTTGATTTGTTTGCCGTTTTCTTCATTTGCCAGCGCCACCACATCACATATCTCTTTCAGCTACGCCGCCCGCCAGCGGATGCACAACCGTCCGCGTttatt                                                                                                                  >  2:124427/1‑250 (MQ=255)
                                                                                          atatGTTAATTTTATAATAATGGTTTATTAAGTTTGTTCAGAAACGATCCGGGATACACTGCAACATTAAAGCATGACCAGCCATTATTGAGTATGTCCTTGCCGTCGATTCCATCCTTTGTATTGTCTGGAGTACGGTTGATTTGTGTGCCGTTTTCTTCATTTGCCAGCGCCACCACAGCACATATCTCTTTCAGCTACGCCGGCCGCCAGCGGATGCAAAACCGTGCGCGTTTATTAAAACAGTACCa                                                                                                      <  1:376782/251‑1 (MQ=255)
                                                                                                     ttATAATAATGGTTTATTAAGTTTGTTCAGAAACGATCCGGGATACACTGCAACATTAAAGCATGACCAGCCATTATTGAGTATGTCCTTGCCGTCGATTCCATCCTTTGTATTGTCTGGATTACTGTTGATTTGTTTGCCGTTTTCTTCATTTGCCAGCGCCACCACATCACATATCTCTTTCAGCTACGCCGCCCGCCAGCGGATGCAAAACCGTGCGCGTTTATTAAAACAGTACCAAACTCATCTg                                                                                            <  1:214817/250‑1 (MQ=255)
                                                                                                         aataatGGTTTATTAAGTTTGTTCCGAAAAGATCCGGGATCCACTGCAACATTAAAGCATGACCAGCCATTATTGAGTATGTCCTTGCCGTCGATTCCATCCTTTGTATTGTCTGGATTACTGTTGATTTGTTTGCCGTTTTCTTCATTTGCCAGCGCCACCACATCACATATCTCTTTCAGCTACGCCGCCCGCCAGCGGATGCAAAACCGTGCGCGTTTATTAAAACAGTACCAAACTCATCTGcaaa                                                                                        >  2:396350/1‑250 (MQ=255)
                                                                                                            aatGGTTTATTAAGTTTGTTCAGAAACGATCCGGGATACACTGCAACATTAAAGCATGACCAGCCATTATTGAGTATGTCCTTGCCGTCGATTCCATCCTTTGTATTGTCTGGATTACTGTTGATTTGTTTGCCGTTTTCTTCATTTGCCAGCGCCACCACATCACATATCTCTTTCAGCTACGCCGCCCGCCAGCGGATGCAAAACCGTGCGCGTTTATTAAAACAGTACCAAACTCATCTGAAAAAGc                                                                                     >  1:430508/1‑250 (MQ=255)
                                                                                                             atGGTTTATTAAGTTTGTTCAGAAACGATCCGGGCTACACTGCAACATTAAAGCATGACCAGCCATTATTGAGTATGTCCTTGCCGTCGATTCCATCCTTTGTATTGTCTGGATTACTGTTGATTTGTTTGCCGTTTTCTTCATTTGCCAGCGCCACCACATCACATATCTCTTTCAGCTACGCCGCCCGCCAGCGGATGCAAAACCGTGCGCGTTTATTAAAACAGTACCAAACTCATCTGAAAAAGc                                                                                     >  2:325412/1‑249 (MQ=255)
                                                                                                             atGGTTTATTAAGTTTGTTCAGAAACGATCCGGGATACACTGCAACATTAAAGCATGACCAGCCATTATTGAGTATGTCCTTGCCGTCGATTCCATCCTTTGTATTGTCTGGATTACTGTTGATTTGTTTGCCGTTTTCTTCATTTGCCAGCGCCACCACATCACATATCTCTTTCAGCTACGCCGCCCGCCAGCGGATGCAAAACCGTGCGCGTTTATTAAAACAGTACCAAACTCATCTGAAAAAGCAg                                                                                   <  1:214831/251‑1 (MQ=255)
                                                                                                                 tttATTAAGTTTGTTCAGAAACGATCCGGGATACACTGCAACATTAAAGCATGACCAGCCATTATTGAGTATGTCCTTGCCGTCGATTCCATCCTTTGTATTGTCTGGATTACTGTTGATTTGTTTGCCGTTTTCTTCATTTGCCAGCGCCACCACATCACATATCTCTTTCAGCTACGCCGCCCGCCAGCGGATGCAAAACAGTGCGCGTTTATTAAAACAGTACCAAACTCATCTGAAAAAGCAGGcc                                                                                >  2:21722/1‑250 (MQ=255)
                                                                                                                   tattaAGTTTGTTCAGAAACGATCCGGGATACACTGCAACATTAAAGCATGACCAGCCATTATTGAGTATGTCATTGCCGTCGATTCCATCCTTTGTATTGTCTGGATTACTGTTGATTTGTTTGCCGTTTTCTTCATTTGCCAGCGCCACCACATCACATATCTCTTTCAGCTACGCCGCCCGCCAGCGGATGCAAAACCGTGCGCGTTTATTAAAACAGTACCAAACTCATCTGAAAAAGCAGGCAAGc                                                                             >  1:424168/1‑251 (MQ=255)
                                                                                                                          tttGTTCAGAAACGATCCGGGATACACTGCAACATTAAAGCATGACCAGCCATTATTGAGTATGTCCTTGCCGTCGATTCCATCCTTTGTATTGTCTGGATTACTGTTGATTTGTTTGCCGTTTTCTTCATTTGCCAGCGCCACCACATCACATATCTCTTTCAGCTACGCCGCCCGCCAGCGGATGCAAAACCGTGCGCGTTTATTAAAACAGTACCAAACTCATCTGAAAAAGCAGGCCAGCTATATtg                                                                      >  2:314446/1‑251 (MQ=255)
                                                                                                                             gTTCAGAAACGATCCGGGAAACACTGCACCATTAACTCATGACCCGCTATTCTTGAGTATATCCTTGCCGTCGAATCCATACTATGTCTTGTCCGGATTACTGGTGATTTGTTTGCCGTTATCTTCATTTGCTAGCGCCACCACATCACATATCTCTTTCACCTCCGCCGCCCGCCAACGGATGCAAAACCGCGCGCGTTTATTAAAACAGAAACAAACTCATCTGAAAAAGCAGGCCAGCTATATTGTgg                                                                   >  2:563903/1‑251 (MQ=255)
                                                                                                                              ttCAGAAACGATCCGGGATACAATGACAAATTACAGCATGACCAGCCATTATTGAGTATGTCCTAGCCGTCGATTCCATCCTTTGTATTGTCTGGATTACTGTTGATTTGTTTGCCGTTTTCTTCATTTGCCAGCGCCACCACATCACATATCTCTTTTAGCTACGCCGCCCGCCAGCGGATGCAAAACCGTGCGCGTTTATTAAAACAGTACCAAACTCCTACGAAAAAGCAGGCAAGCTAAATtgagg                                                                   >  2:270096/1‑247 (MQ=255)
                                                                                                                                      cGATCCGGGATACTCTGCAACATAAAAGCATGCCCAGCCATTATTGATTATGTCCTTGCCGTCGATTCCATCCTTTTGATTGTCTGGATTACTGTTGATTTGTTTGCCGTTTTCTTCATTTGCCGGCGCCACCACATCACATATCTCTGTCAGCTACGCCTCCCGCCAGCGGAGGCAAAACCGTGCGCGTTTATTAAAACAGTACCAAACTCATCTGAAAAAGCAGGCCAGCGATATTGTGGAAGGCAATg                                                          <  2:84247/251‑1 (MQ=255)
                                                                                                                                        aTCCGGGATACACTGCAACATTAAAGCATGACCAGCCATTATTGAGTATGTCCTTGCCGTCGATTCAATACTTTGTATTGTCTGGATTACTGTTGATTTGATTGCCGTTTTCTTCATTTGCCAGCGCCCCCACATCACATATCTCTTTCAGCTACGCCGCCCGCCAGCGGATGCAAAACCGTGCGCGTTTATTAACACAGTACCAAACTCATCTGAAAAAGCAGGCCAGCTATATTGTGGAAGGCAATGcc                                                        >  1:364213/1‑251 (MQ=255)
                                                                                                                                                  accatgcaaccatAAAGCATTGCCAGCCATTATTGAGTATTCCCTTGCGGGCGATTCCATCCTTTGTAGTGTGTGGATTACTGTTGATTTGTTTGCCGTTTTCTTCATTTGCCAGCGCCCCCACATCACATATCTCTTTCAGCTACGCCGCCCGCCAGCGGATGCAAAACCGTGCGCGTTTATTAAAACAGTACCAAACTCATCTGAAAAAGCAGGCCAGCTATATTGTGGAAGGCTATTCCGAAAGCaa                                                <  2:242036/238‑1 (MQ=255)
                                                                                                                                                    cTGCAACATTAAAGCATGACCAGCCATTATTGAGTATGTCCTTGCCGTCGATTCCATCCTTTGTATTGTCTGGATTACTGTTGATTTGTTTGCCGTTTTCTTCCTTTGCCAGCGCCACCACATCACATATCTCTTTCAGCTACGCCGCCCGCCAGCGGATGCAAAACCGTGCGCGTTTATTAAACCAGTACCAAACTCATCTGAAAAAGCAGGCCCGCTATATTGTGGAAGGCAATGCAGAAAGCAAAAgg                                            >  1:526954/1‑251 (MQ=255)
                                                                                                                                                            aTTAAAGCATGACCCGCGCATATTGAGTATGTCCTTGCCGTCGATTCCATCCTTTGTATTGTCTGGATTACTGTTGATTTGTTTGCCGTTTTCTTCATTTGCCAGCGCCACCACATCACATATCTCTTTCAGCTACGCCGCCAGCCACCGGATGCAACAACGTGCGCGTTTATTAAAACAGTACCAAACTCATCTGAAAAAGCAGGCCAGCTATATTGTGGAAGGCAATGCCGAAAGCAAAAGGGCGCTAc                                     >  2:379596/1‑251 (MQ=255)
                                                                                                                                                                   cATGACCAGCCATTATTGAGTATGTCCTTGCCGTCGATTCCATCCTTTGTATTGTCTGGATTACTGTTGATTTGTTTGCCGTTTTCTTCATTTGCCAGCGCCACCACATCACATATCTCTTTCAGCTACGCCGCCCGCCAGCGGATGCAAAACCGTGCGCGTTTATTAAAACAGTACCAAACTCATCTGAAAAAGCAGGCCAGCTATATTGTGGAAGGCAATGCCGAAAGCAAAAGGGCGCTACGCCAGc                               >  1:446812/1‑250 (MQ=255)
                                                                                                                                                                                   ttGAGTATGTCCTTGCCGTCGATTCCATCCTTTGTATTGTCTGGATTACTGTTGATTTGTTTGCCGTTTTCTTCATTTGCCAGCGCCACCACATCACATATCTCTTTCAGCTACGCCGCCCGCCAGCGGATGCAAAACCGTGCGCGTTTATTAAAACAGTACCAAACTCATCTGAAAAAGCAGGCCAGCTATATTGTGGAAGGCAATGCCGAAAGCAAAAGGGCGCTACGCCAGCACAACCGGGAGc                   >  1:319572/1‑247 (MQ=255)
                                                                                                                                                                                      aGTATGTCCTTGCCGTAGATTCCATCCTTTGTATTGTCTGGATTACTGTTGATTTGTTTGCCGTTTTCTTCATTTGCCAGCGCCAACACATCACATATCTCTTTCAGCTACGACGCCCGCCAGCGGATGCAAAACCGTGCGCGTTTATTACAACAGGACCAAACTCATCTGAAAAAGCAGGCCAGCTATATTGTGGAAGGCAATGCCGAAAGCAAAAGGGCGCTACGCCAGCACCACCGGGAGCAGATaa             >  1:540119/1‑250 (MQ=255)
                                                                                                                                                                                        tATGTCCTTGCCGTCGATTCCATCCTTTGTATTGTCTGGATTACTGTTGATTTGTTTGCCGTTTTCTTCATTTGCCAGCGCCACCACATCACATATCTCTTTCAGCTACGCCGCCCGCCAGCGGATGCAAAACCGTGCGCGTTTATTAAAACAGTACCAAACTCATCTGAAAAAGCAGGCCAGCTATATTGTGGAAGGCAATGCCGAAAGCAAAAGGGCGCTACGCCAGCACAACCGGGAGCAGATAAAAc          >  2:363695/1‑251 (MQ=255)
                                                                                                                                                                                            tccttgctggtgtgtCCATCATTTTTATTGTCTGGATTCCTGTTGATTTGTTGTCCGTTTTTTTCATTTGCCTACGCCGCCCCAGCACATATCTCTTGCCTCTACGCCGTCCGCCAGCGGATGCAAAACCGTGCGCGTTTATTAAAACAGTACCAAACTCATCTGAAAAAGCAGTGCAGCTATATTGTGGAAGGCAATGCTGAAATCAAAAGGGCGCCACGGCAGCAAAACCGTGTGCAGATtaaagacc       <  2:413529/236‑5 (MQ=255)
                                                                                                                                                                                            tCCTTGCCGTCGATTCCATCCTTTGTATTGTCTGGATTACTGTTGATTTGTTTGCCGTTTTCTTCATTTGCCAGCGCCACCACATCACATATCTCTTTCAGCTACGCCGCCCGCCAGCGGATGCAAAACCGTGCGCGTTTATTAAAACAGTACCAAACTCATCTGAAAAAGCAGGCCAGCTATATTGTGGAAGGCaa                                                            <  2:144362/197‑1 (MQ=255)
                                                                                                                                                                                            tCCTTGCCGTCGATTCCATCCTTTGTATTGTCTGGATTACTGTTGATTTGTTTGCCGTTTTCTTCATTTGCCAGCGCCACCACATCACATATCTCTTTCAGCTACGCCGCCCGCCAGCGGATGCAAAACCGTGCGCGTTTATTAAAACAGTACCAAACTCATCTGAAAAAGCAGGCCAGCTATATTGTGGAAGGCaa                                                            >  1:144362/1‑197 (MQ=255)
                                                                                                                                                                                                 gCCGTCGATTCCATCCTTTGTATTGTCTGGATTACTGTTGATTTTTTTGCCGTTTCCTTCATTTGCCAGCGCCACCACATCACATATCTCTTTCAGCTACGCCGCCCGCCAGCGGATGCAAAACCGTGCGCGTTTATTAAAACAGTACCAAACTCATCTGAAAAAGCAGGCCAGCTATATTGTGGAAGGCAATGCCGAAAGCAAAAGGGCGCTACGCCAGCACAACCGGGAGCAGATAAAACAGCATCCa  <  1:262653/250‑1 (MQ=255)
                                                                                                                                                                                                  ccGTCGATTCCATCCTTTGTATTGTCTGGATTACTGTTGATTTGTTTGCCGTTTTCTTCATTTGCCAGCGCCACCACATCACATATCTCTTTCAGCTACGCCGCCCGCCAGCGGATGCAAAACCGTGCGCGTTTATTAAAACAGTACCAAACTCATCTGAAAAAGCAGGCCAGCTATATTGTGGAAGGCAATGCCGAAAGCAAAAGGGCGCTACGCCAGCACAACCGGGAGCAGATAAAACAGCATcc   >  1:140405/1‑248 (MQ=255)
                                                                                                                                                                                                                            |                                                                                                                                                                                                                              
ATTCAGCTTGAATTGTAGCTCAAATGAGCAACAGCGACAATCCATGAGTAATTTAATAAACGGATGTTTTTACATAATTAACTGAAATATATATGTTAATTTTATAATAATGGTTTATTAAGTTTGTTCAGAAACGATCCGGGATACACTGCAACATTAAAGCATGACCAGCCATTATTGAGTATGTCCTTGCCGTCGATTCCATCCTTTGTATTGTCGGGATTACTGTTGATTTGTTTGCCGTTTTCTTCATTTGCCAGCGCCACCACATCACATATCTCTTTCAGCTACGCCGCCCGCCAGCGGATGCAAAACCGTGCGCGTTTATTAAAACAGTACCAAACTCATCTGAAAAAGCAGGCCAGCTATATTGTGGAAGGCAATGCCGAAAGCAAAAGGGCGCTACGCCAGCACAACCGGGAGCAGATAAAACAGCATCCA  >  NC_000913/246529‑246969

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