Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 826,369 A→C 9.5% V295G (GTC→GGC)  ybhS ← putative ABC transporter permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913826,3690AC 9.5% 149.4 / 5.8 60V295G (GTC→GGC) ybhSputative ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base A (25/29);  new base C (2/4);  total (27/33)
Fisher's exact test for biased strand distribution p-value = 6.81e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.12e-01

CCGGCGAGGAACAGGCTTTGCAGGGTGCTGACGAAATAACGAGCGGGAATAATGTACGTCACCGCGCGGATCACCGCGGGCATACTGTCGATCTG‑AAAAATAAAGCCGGAAAGCATAATCGACGGCAGAAAAGCGGCGTTCAGGGCGACCTGAGCGGCATTGAACTGGTTGCGGGTAATCGTGGAAATCAGCAGCCCCATCCCCAGGGTACTGAGTAAAAACAGGCTGGAGATAAAAAACAGAATCAGCAGCGACCCGCGATACGGCACGCCGAGAATAAACACTGACACCAGCAT  >  NC_000913/826222‑826517
                                                                                                                                                    |                                                                                                                                                    
ccGGCGAGGAACAGGCTTTGCAGGGTGCTGACGAAATAACGAGCGGGAATAATGTACGTCACCGCGCGGATCACCGCGGGCATACTGTCGATCTG‑AAAAATAAAGCCGGAAAGCATAATCGACGGCAGAAAAGCGGCGTTCAGGGCGAcc                                                                                                                                                    <  2:845632/150‑1 (MQ=255)
  ggCGATTAACAGTCTTTGCAGGTTGCTGACGAAATAACGAGCGGGAATAATGTACGTCACCGCGCGGATCACCTCGGGCATACTGTCGATCTG‑AAAAATAAAGCCGGAAAGCATAATCGACGTCAGAAAAGCGGCGTTCAGGGCGACCTg                                                                                                                                                  <  1:106630/150‑1 (MQ=255)
       ggAACAGGCTTTGCAGGGTGCTGACGAACTAACGCGCGGGAATAATGTACGGCACCGCGCGGATCACCGCGGGCATAGTGTCGATCTG‑AAAAATAAAGCCGTAAAGCATAATCGACGGCAGAAAAGCGGCGTTCAGGGCGACCTGAGCgg                                                                                                                                             <  1:1108672/150‑1 (MQ=255)
        gAACAGGCTTTGCAGGCTGCTGACGAAATAACGAGCGGGAATAATGTACGTCACCGCGCGGATCACCGCGGGCATACTGTCGATCTG‑AAAAATAAAGCCGGAAAGCATAATCGACGGCAGAAAAGCGGCGCTCAGGTCGCCCTGAACgcc                                                                                                                                            <  2:1231759/150‑3 (MQ=255)
        gAACAGGCGTTGCAGGGTGCGGACGAAATAACGAGCGGGAATAATGTACGTCACCGCGCGGATCACCGCGGGCATACTGTCGATCTG‑AAAAATAAAGCCGGAAAGCATAATCGACGGCAGAAAAGCGGCGCTCAGGTCGCCCTGAACgcc                                                                                                                                            <  2:1226842/150‑3 (MQ=255)
         aaCAGGCTTTGCAGGGTGCTGACGAAATAACGAGCGGGAATAATGTACGTCACCGCGCGGATCACCGCGGGCATACTGTCGATCTG‑AAAAATAAAGCCGGAAAGCATAATCGACGGCAGAAAAGCGGCGCTCAGGTCGCCCTGAACgccg                                                                                                                                           <  2:771907/150‑4 (MQ=255)
            aGGCTTTGCAGGGTGCTGACGAAATAACGAGCGGGAATAATGTACGTCACCGCGCGGATCACCGCGGGCATACTGTCGATCTG‑AAAAATAAAGCCGGAAAGCATAATCGACGGCAGAAAAGCGGCGTTCAGGGCGACCTGAGCGGCATTg                                                                                                                                        <  2:853025/150‑1 (MQ=255)
                     aGGGTGCTGACGAAATAACGAGCGGGAATAATGTACGTCACCGCGCGGATCACCGCGGGCATACTGTCGATCTG‑AAAAATAAAGCCGGAAAGCATAATCGACGGCAGAAAAGCGGCGTTCAGGGCGACCTGAGCGGCATTGAACTGGTTg                                                                                                                               <  2:1136254/150‑1 (MQ=255)
                       ggggCTGACGAAATAACGCGCGGGAATAATGTACGTCACCGCGCGGATCACCGCGGGCATACTGTCGATCGG‑AAAATAAAAGCCGGAAAGCATACTCGACGGCAGAATAGCGGCGTTCAGGGCGACCTGAGCGGCATTGAACTGGTTGCg                                                                                                                             <  1:1228006/147‑1 (MQ=255)
                       ggTGCTGACGAAATAACGAGCGGGAATAATGTACGTCACCGCGCGGATCACCGCGGGCATACTGTCGATCTG‑AAAAATAAAGCCGGAAAGCATAATCGACGGCAGAAAAGCGGCGTTCAGGGCGACCTGAGCGGCATTGAACTGGTTGCg                                                                                                                             <  1:1137559/150‑1 (MQ=255)
                         tGCTGACGAAATAACGAGCGGGAATAATGTACGTCACCGCGCGGATCACCGCGGGCATACTGTCGATCTG‑AAAAATAAAGCCGGAAAGCATAATCGACGGCAGAAAAGCGGCGTTCAGGGCGACCTGAGCGGCATTGAACTGGTTGCggg                                                                                                                           <  1:238674/150‑1 (MQ=255)
                         ggCTGACCAAATAACGAGCGTGAATAATGTACGTCACCGCGCGGATCTCCGCGGGCATACTGTCGATCTG‑AAAAATAAAGCCGGAAAGCATAATCGACGGCAGAAAAGCGGCGTTCAGGGCGACCTGAGCGGCATTGAACTGGTTGCggg                                                                                                                           <  2:106763/149‑1 (MQ=255)
                           cTGACGAAATAACGAGCGGGAATAATGTACGTCACCGCGCGGATCACCGCGGGCATACTGTCGATCTG‑AAAAATAAAGCCGGAAAGCATAATCGACGGCAGAAAAGCGGCGTTCAGGGCGACCTGAGCGGCATTGAACTGGTTGCGGGTa                                                                                                                         >  2:586645/1‑150 (MQ=255)
                            tGACGAAATAACGAGCGGGAATAATGTACGTCACCGCGCGGATCACCGCGGGCATACTGACGATCTG‑AAAAATAAAGCCGGAAATCATAATCGACGGCAGAAACGCGGCGTTCAGGGCGACCTGAGCGGCATTGAACTGGTTGCGGGTaa                                                                                                                        >  1:21683/1‑150 (MQ=255)
                               cGAAATAACGAGCGGGAATAATGTACGTCACCGCGCGGATCACCGCGGGCATACTGTCGATCTG‑AAAAATAAAGCCGGAAAGCATAATCGACGGCAGAAAAGCGGCGTTCAGGGCGACCTGCGCGGCATTGAACTGGTTGCGGGTAATCg                                                                                                                     >  1:118871/1‑150 (MQ=255)
                                 aaaTAACGAGCGGGAATAATGTACGTCACCGCGCGGATCACCGCGGGCATACTGTCGATCTG‑AAAAATAAAGCCGGAAAGCATAATCGACGGCAGAAAAGCGGCGCTCAGGTCGCCCTg                                                                                                                                                  >  1:395443/1‑119 (MQ=255)
                                 aaaTAACGAGCGGGAATAATGTACGTCACCGCGCGGATCACCGCGGGCATACTGTCGATCTG‑AAAAATAAAGCCGGAAAGCATAATCGACGGCAGAAAAGCGGCGCTCAGGTCGCCCTg                                                                                                                                                  <  2:395443/119‑1 (MQ=255)
                                     aaCGAGCGGGAATAATGTACGTCACCGCGCGGATCACCGCGGGCATACTGTCGATCTG‑AAAAATAAAGCCGGAAAGCATAATCGACGGCAGAAAAGCGGCGTTCAGGGCGACCTGAGCGGCATTGAACTGGTTGCGGGTAATCGTGGaaa                                                                                                               >  2:753465/1‑150 (MQ=255)
                                         aGCGGGAATAATGTACGTCACCGCGCGGATCACCGCGGGCATACTGTCGATCTG‑AAAAATAAAGCCGGAAAGCATAATCGACGGCAGAAAAGCGGCGTTCAGGGCGACCTGAGCGGCATTGAACTGGTTGCGGGTAATCGTGGAAATcag                                                                                                           >  1:1230481/1‑150 (MQ=255)
                                           cGGGAATAATGTACGTCACCGCGCGGATCACCGCGGGCATACTGTCGATCTG‑AAAAATAAAGCCGGAAAGCATAATCGCCGGCAGAAAAGCGGCGTTCAGGGCGACCTGAGCGGCATTGAACTGGTTGC‑GGTAATCGTGGAAATcagcag                                                                                                        >  1:524752/1‑150 (MQ=255)
                                             ggAATAATGTACGTCACCGCGCGGATCACCGCGGGCATACTGTCGATCTG‑AAAAATAAAGCCGGAAAGCATAATCGACGGCAGAAAAGCGGCGTTCAGGGCGCCCTGAGCGGCATTGAACTGGTTGCGGGTTATCGTGGAAATCagcagc                                                                                                       >  2:583892/1‑150 (MQ=255)
                                              gAATAATGTACGTCACCGCGCGGATCACCGCGGGCATACTGTCGATCTG‑AAAAATAAAGCCGGAAAGCATAATCGACGGCAGAAAAGCGGCGTTCAGGGCGACCTGAGCGGCATTGAACTGGTTGCGGGTAATCGTGGAAATCAGCAGcc                                                                                                      <  2:782874/150‑1 (MQ=255)
                                                ataatGTACGTCACCGCGCGGATCACCGCGGGCATACTGTCGATCTG‑AAAAATAAAGCCGGAAAGCATAATCGACGGCAGCAAAGCGGCGTTCAGGGCGACCTGAGCGGCATTGAACTGGTTGCGGGTAATCGTGGAAATCAGCAGcccc                                                                                                    >  2:973989/1‑150 (MQ=255)
                                                 taatGTACGTCACCGCGCGGATCACCGCGGGCATACTGTCGATCTG‑AAAAATAAAGCCGGAAAGCATAATCGACGGCAGAAAAGCGGCGTTCAGGGCGACCTGAGCGGCATTGAACTGGTTGCGGGTAATCGTGGAAATCAGCAGCCCCa                                                                                                   >  1:1131103/1‑150 (MQ=255)
                                                     gTACGTCACCGCGCGGATCACCGCGGGCATACTGTCGATCTG‑AAAAATAAAGCCGGAAAGCATAATCGACGGCAGAAAAGCGGCGTTCAGGGCGACCTGAGCGGCATTGAACTGGTTGCGGGTAATCGTGGAAATCAGCAGCCCCATccc                                                                                               <  1:524615/150‑1 (MQ=255)
                                                       aCGTCACCGCGCGGATCACCGCGGGCATACTGGCGATCTG‑AAAAATAAAGCCGGAAAGCATAATCGACGGCAGAAAAGCGGCGTTCAGGGCGACCTGAGCGGCATTGAACTGGTTGCGGGTAATCGTGGAAATCAGCAGCCCCATCCCCa                                                                                             <  1:248501/150‑1 (MQ=255)
                                                        cGTCACCGCGCGGATCACCGCGGGCATACTGTCGATCTG‑AAAAATAAAGCCGGAAAGCATAATCGACGGCAGAAAAGCGGCGTTCAGGGCGACCTGAGCGGCATTGAACTGGTTGCGGGTAATCGTGGAAATCAGCAGCCCCATCCCCAg                                                                                            >  2:874943/1‑150 (MQ=255)
                                                          tCACCGCTCGGATCACCGCGGGCATACTGTCGATCTG‑AAAAATAAACCAGGAAAGCATAATCGACGGCAGAAAAGCGGCGTTCAGGGCGACCTGAGCGGCATTGAACTGGTTGCGGGTAATCGTGGAAATCAGCAGCCCCATCCCCAggg                                                                                          <  2:167852/150‑1 (MQ=255)
                                                          tCACCGCGCGGATCACCGCGGGCATACTGTCGATCTG‑AAAAATAAAGCCGGAAAGCATAATCGACGGCAGAAAAGCGGCGTTCAGGGCGACCTGAGCGGCATTGAACTGGTTGCGGGTAATCGTGGAAATCAGCAGCCCCATCCCCAggg                                                                                          <  1:689288/150‑1 (MQ=255)
                                                          tCACCGCGCGGATCACCGCGGGCATACTGTCGATCTG‑AAAAATAAAGCCGGAAAGCATAATCGACGGCAGAAAAGCGGCGTTCAGGGCGACCTGAGCGGCATTGAACTGGTTGCGGGTAATCGTGGAAATCAGCAGCCCCATCCCCAggg                                                                                          >  2:417400/1‑150 (MQ=255)
                                                               gcgcgGATCACCGCGGGCATACTGTCGATCTG‑AAAAATAAAGCCGGAAAGCATAATCGACGGCAGAAAAGCGGCGTTCAGGGCGACCTGAGCGGCATTGAACTGGTTGCGGGTAATCGTGGAAATCAGCAGCCCCATCCCCAGGGTACTg                                                                                     >  1:1109385/1‑150 (MQ=255)
                                                                cgcgGATCACCGCGGGCATACTGTCGATCTG‑AAAAATAAAGCCGGAACGCATAATCGACGGCAGAAAAGCGGCGTTCAGGGCGACCTGAGCGGCATTGATCTGGTTGCGGGTAATCGTGGAAATCAGCAGCCCCAACCCCAGGGTACTGa                                                                                    >  1:10434/1‑150 (MQ=255)
                                                                  cgGATCACCGCGGGCATACTGTCGATCTG‑AAAAATAAAGCCGGAAAGCATAATCGACGGCAGAAAAGCGGCGTTCAGGGCGACCTGAGCGGCATTGAACTGGTTGCGGGTAATCGTGGAAATCAGCAGCCCCATCCCCAGGGTACTGAGt                                                                                  >  2:47687/1‑150 (MQ=255)
                                                                   ggATCACCGCGGGCATACTGTCGATCTGAAAAAATACAGCCGGAAAGCATAATCGACGGCAGAAAAGCGGCGTTCAGGGCGACCTGAGCGGCATTGAACTGGGTGCGGGTAATCGTGGAAATCAGCAGCCCCATCCCCAGGGTACTGAGt                                                                                  <  1:998653/150‑1 (MQ=255)
                                                                   ggATCACCGCGGGCATACTGTCGATCTG‑AAAAATACAGCCGGAAAGCATAATCGACGGCAGAAACGCCGCGGTCAGGGCGACCTGAGCGGCATTGAACTGGTTGCGGGTAATCGTGGAAATCAGCAGCCCCATCCCCAGGGTACTGAGTa                                                                                 <  2:36856/150‑1 (MQ=255)
                                                                       cACCGCGGGCATACTGTCGATCTG‑AAAAATAAAGCCGGAAAGCATAATCGACGGCAGAAAAGCGGCGTTCAGGGCGACCTGAGCGGCATTGAACTGGTTGCGGGTAATCGTGGAAATCAGCAGCCCCATCCCCAGGGAACTGAGTaaaaa                                                                             >  2:1197032/1‑150 (MQ=255)
                                                                                        cGATCTG‑AAAAATAAATCCGTAAATCATAATCGATGGCAGAAAAGCGGCGTTCAGGGCGACCTGAGCGGCATTGAACTGGTTGCGGGTAATCGTGGAAATCAGCAGCCCCATCCCCAGGGTACTGAGTAAAAACAGGCTGGAGATaaaaa                                                            <  2:408848/150‑1 (MQ=255)
                                                                                        cGATCTG‑AAAAATAAAGCCGGAAAGCATAATCGACGGCAGAAAAGCGGCGTTCAGGGCGACCTGAGCGGCATTGAACTGGTTGCGGGTAATCGTGGAAATCAGCAGCCCCATCCCCAGGGTACTGAGTAAAAACAGGCTGGAGATaaaaa                                                            <  2:1094314/150‑1 (MQ=255)
                                                                                          aTCTG‑AAAAATAAAGCCGGAAAGCATCAGCGAGTGCAGAATAGCGGCGTTCAGGGCGACCTGAGCGGCATTGAACTGGTTGCGGGTAATCGTGGAAATCAGCAGCCCCATCCCCAGGGTACTGAGTAAAAACAGGCTGGAGATAAAAAAc                                                          <  1:902822/150‑1 (MQ=255)
                                                                                          aTCTG‑AAAAATAAAGCCGGAAAGCATAATCGACGGCAGAAAAGCGGCGTTCAGGGCGACCTGAGCGGCATTGAACTGGTTGCGGGTAATCGTGGAAATCAGCAGCCCCATCCCCAGGGTACTGAGTAAAAACAGGCTGGAGATAAAAAAc                                                          <  1:1127641/150‑1 (MQ=255)
                                                                                           tCTG‑AAAAATAAAGCCGGAAATCATAATCGACGGCAGAAAAGCGGCGTTCAGGGCGACATGAGCGTCATTGAACTGGTTGCGGGTAATCGTGGAAATCAGCAGCCCCATCCCCAGCGTACTGAGTAAAAACAGGCTGGAGATAAAAAACa                                                         <  2:639978/150‑1 (MQ=255)
                                                                                            cTG‑AAAAATAAAGCCGGAAAGCATAATCGACGGCAGAAAAGCGGCGTTCAGGGCGACCTGAGCGGCATTGAACTGGTTGCGGGTAATCGTGGAAATCAGCAGCCCCATCCCCAGGGTACTGAGTAAAAACAGGCTGGAGATAAAAAACAg                                                        >  1:912847/1‑150 (MQ=255)
                                                                                                aaaaaTAAAGCCGGAAAGCATAATCGACGGCAGAAAAGCGGCGTTCAGGGCGACCTGAGCGGCATTGAACTGGTTGCGGGTAATCGTGGAAATCAGCAGCCCCATCCCCAGGGTACTGAGGAAAAACAGGCTGGAGATaaaaaaaaaaaa                                                     >  2:130120/1‑144 (MQ=255)
                                                                                                   aaTAAAGCCGGAAAGCATAATCGACGGCAGAAAAGCGGCGTTCAGGGCGACCTGAGCGGCATTGAACTGGTTGCGGGTAATCGTGGAAATCAGCAGCCCCATCCCCAGGGTACTGAGTAAAAACAGGCTGGAGATAAAAACCAGAATcag                                                  >  2:878709/1‑150 (MQ=255)
                                                                                                     tAAAGCCGGAAAGCATAATCGACGGCAGAAAAGCGGCGTTCAGGGCGACCTGAGCGGCATTGAACTGGTTGCGGGTAATCGTGTAAATCAGCAGCCCCATCCCCAGGGTACTGAGTAAAAACAGGCGGGAGATAAAAAACAGAATcagcc                                                >  2:1207295/1‑149 (MQ=255)
                                                                                                          ccGGAAAGCATAATCGACGGCAGAAAAGCGGCGTTCAGGGCGACCTGAGCGGCATTGAACTGGTTGCGGGTAATCGTGGAAATCAGCAGCCCCATCCCCAGGGTACTGAGTAAAAACAGGCTGGAGATAAAAAACAGAATCAGCAGCGAc                                           >  2:1286730/1‑150 (MQ=255)
                                                                                                            ggAAAGCATAATCGACGGCAGAAAAGCGGCGTTCAGGGCGACCTGAGCGGCATTGAACTGGTTGCGGGTAATCGTGGAAATCAGCAGCCCCATCCCCAGGGTACTGAGTAAAAACAGGATGGAGATAAAAAACAGAATCAGCAGCGcccc                                         >  1:131654/1‑150 (MQ=255)
                                                                                                                  cATAATCGACGGCAGAAAAGCGGCGTTCAGGGCGACCTGAGCGGCATTGAACTGGTTGCGGGTAATCGTGGAAATCAGCAGCCCCATCCCCAGGGTACTGAGTAAAAACAGGCTGGAGATAAAAAACAGAATCAGCAGCGACCCGCGATa                                   >  1:1106318/1‑150 (MQ=255)
                                                                                                                    tAATCGACGGCAGAAAAGCGGCGTTCAGGGCGACCTGAGCGGCATTGAACTGGTTGCGGGTAATCGTGGAAATCAGCAGCCCCATCCCCAGGGTACTGAGTAAAAACAGGCTGGAGATAAAACACAGAATCAGCAGCGCCCCGCGATACg                                 >  1:1151465/1‑150 (MQ=255)
                                                                                                                      aTCGACGGCAGAAAAGCGGCGTTCAGGGCGACCTGAGCGGCATTGAACTGGTTGCGGGTAATCGTGGAAATCAGCAGCCCCATCCCCAGGGTACTGAGTAAAAACAGGCGGGAGATAAAAAACAGAATCAGCAGCGACCCGCGATACGGc                               <  1:126141/150‑1 (MQ=255)
                                                                                                                      aTCGACGGCAGAAAAGCGGCGTTCAGGGCGACCTGAGCAGCATTGAACTGGTTGCGGGTAATCGTGGAAATCAGCAGCCCCATCCCCAGGGTACTGAGTAAAAACAGGCTGGAGATAAAAAACAGAATCAGCAGCGACCCGCGATACGGc                               <  2:141314/150‑1 (MQ=255)
                                                                                                                       tCGACGGCAGAAAAGCGGCGTTCAGGGCGACCTGAGCGGCATTGAACTGGTTGCGGGTAATCGTGGAAATCAGCAGCCCCATCCCCAGGGTACTGAGTAAAAACAGGCTGGAGATAAAAAACAGAATCAGCAGCGACCCGCGATACGGCa                              <  1:716744/150‑1 (MQ=255)
                                                                                                                          aCGGCAGAAAAGCGGCGTTCAGGGCGACCTGAGCGGCATTGAACTGGTTGCGGGTAATCGTGGAAATCAGCAGCCCCATCCCCAGGGTACTGAGTAAAAACAGGCTGGAGATATAAAACAGAATCAGCAGCGCCGCGCGATACGGCCCGc                           >  1:230970/1‑150 (MQ=255)
                                                                                                                              cAGAAAAGCGGCGTTCAGGGCGACCTGAGCGGCATTGAACTGGTTGCGGGTAATCGTGGAAATCAGCAGCCCCATCCCCAGGGTACTGAGTAAAAACAGGCTGGAGATAAAAAACAGAATCAGCAGCGACCCGCGATACGGCACGCCgag                       <  1:756499/150‑1 (MQ=255)
                                                                                                                                      cggcgTTCAGGGCGACCTGAGCGGCATTGAACTGGTTGCGGGTAATCGTGGAAATCAGCAGCCCCATCCCCAGGGTACTGAGTAAAAACAGGCTAGAGATAAAAAACAGAATCAGCAGCGACCCGCGATACGGCACGCCGAGAATaatcc               >  2:708196/1‑147 (MQ=255)
                                                                                                                                          gTTCAGGGCGACCTGAGCGGCATTGAACTGGTTGCGGGTAATCGTGGAAATCAGCAGCCCCATCCCCAGGGTACTGAGTAAAAACAGGCTGGAGATAAAAAACAGAATCAGCAGCGACCCGCGATACGGCACGCCGAGAATAAACACTGa           <  1:417400/150‑1 (MQ=255)
                                                                                                                                                ggCGACCTGAGCGTCATTGAACTGGTTGCGGGTAATCGGGTAAATCAGCAGCCCCATCCCCAGGGTACTGAGTAAAAACAGGCTGGAGATAAAAAACAGAATCAGCAGCGACCCGCGATACGGCACGCCGAGAATAAACACTGACACCAg     <  2:711019/150‑1 (MQ=255)
                                                                                                                                                  cGACCTGAGCGGCATTGAACTGGTTGCGGGTAATCGTGGAAATCAGCAGCCCCATCCCCAGGGTACTGAGTAAAAACAGGCTGGAGATAAAAAACAGAATCAGCAGCGACCCGCGATACGGCACGCCGAGAATAAACACTGACACCAGCa   <  2:912847/150‑1 (MQ=255)
                                                                                                                                                  cGACCTGAGCGGCATTGAACTGGTTGCGGGTAATCGTGGAAATCAGCAGCCCCATCCCCAGGGTACTGAGTAAAAACAGGCTGGAGATAAAAAACAGAATCAGCAGCGACCCGCGATACGGCACGCCGAGAATAAACACTGACACCAGCa   <  1:855157/150‑1 (MQ=255)
                                                                                                                                                   gACCTGAGCGGCATTGAACTGGTTGCGGGTAATCGTGGAAATCAGCAGCCCCATCCCCAGGGTACTGAGTAAAAACAGGCTGGAGATAAAAAACAGAATCAGCAGCGACCCGCGATACGGCACGCCGAGAATAAACAATGACACCAGCAt  >  1:650660/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
CCGGCGAGGAACAGGCTTTGCAGGGTGCTGACGAAATAACGAGCGGGAATAATGTACGTCACCGCGCGGATCACCGCGGGCATACTGTCGATCTG‑AAAAATAAAGCCGGAAAGCATAATCGACGGCAGAAAAGCGGCGTTCAGGGCGACCTGAGCGGCATTGAACTGGTTGCGGGTAATCGTGGAAATCAGCAGCCCCATCCCCAGGGTACTGAGTAAAAACAGGCTGGAGATAAAAAACAGAATCAGCAGCGACCCGCGATACGGCACGCCGAGAATAAACACTGACACCAGCAT  >  NC_000913/826222‑826517

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