Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 826,359 G→A 6.3% N298N (AAC→AAT)  ybhS ← putative ABC exporter membrane subunit YbhS

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913826,3590GA 6.3% 196.2 / 2.8 64N298N (AAC→AAT) ybhSputative ABC exporter membrane subunit YbhS
Reads supporting (aligned to +/- strand):  ref base G (26/33);  new base A (2/2);  total (29/35)
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.48e-01

CTGGAATATTCCCGGCGAGGAACAGGCTTTGCAGGGTGCTGACGAAATAACGAGCGGGAATAATGTACGTCACCGCGCGGATCACCGCGGGCATACTGTCGATCTGAAAAATAAAGCCGGAAAGCATAATCGACGGCAGAAAAGCGGCGTTCAGGGCGACCTGAGCGGCATTGAACTGGTTGCGGGTAATCGTGGAAATCAGCAGCCCCATCCCCAGGGTACTGAGTAAAAACAGGCTGGAGATAAAAAACAGAATCAGCAGCGACCCGCGATACGGCACGCCGAGA  >  NC_000913/826211‑826497
                                                                                                                                                    |                                                                                                                                          
cTGGAATATTACCTGAGAGGAACAGGCTTTGCAGGTTTAAGAATAAATAACGAGCCGTAATAATGTACGTCAACGAGCTGATAACAGCGGGCATAATGTCGATCTGAAAAATAAATCCTTAAATCATAATCTAAGGCAGAAAAGCGGAGt                                                                                                                                           >  2:1308330/1‑150 (MQ=255)
  ggAATATTCCCGGCGAGGAACAGGCTTTGCAGGGTGCTGACGAAATAACGAGCGGGAATAATGTACGTCACCGCGCGGATCACCGCGGGCATACTGTCGATCTGAAAAATAAAGCCGGAAAGCATAATCGACGGCAGAAAAGCGGCGTTc                                                                                                                                         >  1:628209/1‑150 (MQ=255)
     atatTCCCGGCGAGGAACAGGCTTTGCAGGGTGCTGACGAAATAACGAGCGGGAATAATGTACGTCACCGCGCGGATCACCGCGGGCATACTGTCGATCTGAAAAATAAAGCCGGAAAGCATAATCGACGGCAGAAAAGCGGCGTTCAgg                                                                                                                                      <  2:1736981/150‑1 (MQ=255)
          cccGGCGAGGAACAGGCTTTGCAGGGTGCTGACGAAATAACGAGCGGGAATAATGTACGTCACCGCGCGGATCACCGCGGGCATACTGTCGATCTGAAAAATAAAGCCGGAAAGCATAATCGACGGTAGAAAAGCGGCGTTCAGGGCGAc                                                                                                                                 >  2:1102556/1‑150 (MQ=255)
           ccGGCGAGGAACAGGCTTTGCAGGGTGCTGACGAAATAACGAGCGGGAATAATGTACGTCACCGCGCGGATCACCGCGGGCATACTGTCGATCTGAAAAATAAAGCCGGAAAGCATAATCGACGGCAGAAAAGCGGCGt                                                                                                                                           >  1:112005/1‑139 (MQ=255)
           ccGGCGAGGAACAGGCTTTGCAGGGTGCTGACGAAATAACGAGCGGGAATAATGTACGTCACCGCGCGGATCACCGCGGGCATACTGTCGATCTGAAAAATAAAGCCGGAAAGCATAATCGACGGCAGAAAAGCGGCGt                                                                                                                                           <  2:112005/139‑1 (MQ=255)
           ccGGCGAGGAACAGGCTTTGCAGGGTGCTGACGAAATAACGAGCGGGAATAATGTACGTCACCGCGCGGATCACCGCGGGCATACTGTCGATCTGAAAAATAAAGCCGGAAAGCATAATCGACGGCAGAAAAGCGGCGTTCAGGGCGAcc                                                                                                                                >  1:154653/1‑150 (MQ=255)
              gCGAGGAACAGGCTTTGCAGGGTGCTGACGAAATAACGAGCGGGAATAATGTACGTCACCGCGCGGATCACCGCGGGCATACTGTCGATCTGAAAAATAAAGCCGGAAAGCATAATCGACGGCAGAAAAGCGGCGTTCAGGGCGACCTGa                                                                                                                             <  2:1933542/150‑1 (MQ=255)
                    aaCAGGCTTTGCAGGGTGCTGACGAAATAACGAGCGGGAATAATGTACGTCACCGCGCGGATCACCGCGGGCATACTGTCGATCTGAAAAATAAAGCCGGAAAGCATAATCGACGCCAGAAAAGCGGCCTTCAGGGCGACCTGAGCCGCa                                                                                                                       >  2:962058/1‑150 (MQ=255)
                     aCAGGCTTTGCAGGGTGCTGACGAAATAACGAGCGGGAATAATGTACGTCACCGCGCGGATCACCGCGGGCATACTGTCGATCTGAAAAATAAAGCCGGAAAGCATAATCGACGGCAGAAAAGCGGCGTTCAGGGCGACCTGAGCGGCAt                                                                                                                      >  2:1741770/1‑150 (MQ=255)
                          cTTTGCAGGGTGCTGACGAAATAACGAGCGGGAATAATGTACGTCACCGCGCGGATCACCGCGGGCATACTGTCGATCTGAAAAATAAAGCCGGAAAGCATAATCGACGGCAGAAAAGCGGCGTTCAGGGCGACCTGAGCGGCATTGAAc                                                                                                                 <  2:1423263/150‑1 (MQ=255)
                          cTTTGCAGGGTGCTGACGAAATAACGAGCGGGAATAATGTACGTCACCGCGCGGATCACCGCGGGCATACTGTCGATCTGAAAAATAAAGCCGGAAAGCATAATCGACGGCAGAAAAGCGGCGTTCAGGGCGACCTGAGCGGCATTGAAc                                                                                                                 <  1:1639636/150‑1 (MQ=255)
                              gCAGGGTGCTGACGAAATAACGAGCGGGAATAATGTACGTCACCGCGCGGATCACTGCGGGCATACTGTCGATCTGAAAAATAAAGCCGGAAAGCATAATCGACGGCAGAAAAGCGGCGTTCAGGGCGACCTGAGCGGCATTGAACTGGt                                                                                                             <  2:1184272/150‑1 (MQ=255)
                              gCAGGGTGCTGACGAAATAACGAGCGGGAATAATGTACGTCACCGCGCGGATCACCGCGGGCATACTGTCGATCTGAAAAATAAAGCCGGAAAGCATAATCGACGGCAGAAAAGCGGCGTTCAGGGCGACCTGAGCGGCATTGAACTGGt                                                                                                             >  2:1500575/1‑150 (MQ=255)
                              ccAGGATGCTAACAAAAAACAGAGAGGTAATAATGTACGTCAACGCGCGGATAACAGCGGGCATACTGACGATCTGAAAAATAAAGCCGGAAAGCATAATCGACGGCAGAAAAGCGGCGTTCAGGGCGACCTGAGCGGCATTGAACTGGt                                                                                                             >  2:1821712/2‑150 (MQ=255)
                                 gggTGCTGACGAAATAACGAGCGGGAATAATGTACGTCACCGCGCGGATCACCGCGGGCATACTGTCGATCTGAAAAATAAAGCCGGAAAGCATAATCGACGGCAGAAAAGCTGCGTTCAGGGCGACCTGAGCGGCATTGAACTGGTTGc                                                                                                          <  2:1703046/150‑1 (MQ=255)
                                  ggTGCTGACGAAATAACGAGCGTGAATAATGTACGTCACCGCGCGGATCACCGCGGGCATACTGTCGATCTGAAAAATAAAGCCGGAAAGCATAATCGACGGCAGAAAAGCGGCGTTCAGGGCGACCTGAGCGGCATTGAACTGGTTGCg                                                                                                         <  1:1308330/150‑1 (MQ=255)
                                     gCTGACGAAATAACGAGCGGGAATAATGTACGTCACCTCGCGGATCACCGCGGGCATACTGTCGATCTGAATAATAAAGCCGGAAAGCATAATCGACGGCATAAAAGCGGCGTTCAGGGCGACCTGAGCGGCATTGAACTGGTTGCGGGt                                                                                                      <  2:1453985/150‑1 (MQ=255)
                                     gCTGACGAAATAACGAGCGGGAATAATGTACGTCACCGCGCGGATCACCGCGGGCATACTGTCGATCTGAAAAATAAAGCCGGAAAGCATAATCGACGGCAGAAAAGCGGCGTTCAGGGCGACCTGAGCGGCATTGAACTGGTTGCGGGt                                                                                                      >  1:341837/1‑150 (MQ=255)
                                             aaTAACGAGCGGGAATAATGTACGTCACCGCGCGGATCACCGCGGGCATACTGTCGATCTGAAAAATAAAGCCGGAAAGCATAATCGACGGCAGAAAAGCGGCGTTCAGGGCGACCTGAGCGGCATTGAACTGGTTGCGGGTAATCGTgg                                                                                              >  1:900356/1‑150 (MQ=255)
                                                  cGAGCGGGAATAATGTACGTCACCGCGCGGATCACCGCGGGCATACTGTCGATCTGAAAAATAAAGCCGGAAAGCATAATCGACGGTAGAAAAGCGGCGTTCAGGGCGACCTGAGCGGCATTGAACTGGTTGCGGGTAATCGTGGAAATc                                                                                         <  1:1102556/150‑1 (MQ=255)
                                                  cGAGCGGGAATAATGTACGTCACCGCGCGGATCACCGCGGGCATACTGTCGATCTGAAAAATAAAGCCGGAAAGCATAATCGACGGCAGAAAAGCGGCGTTCAGGGCGACCTGAGCGGCATTGAACTGGTTGCGGGTAATCGTGGAAATc                                                                                         >  1:1213982/1‑150 (MQ=255)
                                                   gAGCGGGAATAATGTACGTCACCGCGCGGATCACCGCGGGCATACTGTCGATCTGAAAAATAAAGCCGGAAAGCATAATCGACGGCAGAAAAGCGGCGTTCAGGGCGACCTGAGCGGCATTGAACTGGTTGCGGGTAATCGTGGAAATca                                                                                        >  1:1909640/1‑150 (MQ=255)
                                                              atGTACGTCACCGCGCGGATCACCGCGGGCATACTGTCGATCTGAAAAATAAAGCCGGAAAGCATAATCGACGGCAGAAAAGCGGCGTTCAGGGCGACCTGAGCGGCATTGAACTGGTTGCGGGTAATCg                                                                                                 >  1:1114941/1‑130 (MQ=255)
                                                              atGTACGTCACCGCGCGGATCACCGCGGGCATACTGTCGATCTGAAAAATAAAGCCGGAAAGCATAATCGACGGCAGAAAAGCGGCGTTCAGGGCGACCTGAGCGGCATTGAACTGGTTGCGGGTAATCg                                                                                                 <  2:1114941/130‑1 (MQ=255)
                                                              atGTACGTCACCGCGCGGATCACCGCGGGCATACTGTCGATCTGAAAAATAAAGCCGGAAAGCATAATCGACGGCAGAAAAGCGGCGTTCAGGGCGACCTGAGCGGCATTGAACTGGTTGCGGGTAATCGTGGAAATCAGCAGCCCCATc                                                                             <  2:1426567/150‑1 (MQ=255)
                                                               tGTACGTCACCGCGCGGATCACCGCGGGCATACTGTCGATCTGAAAAATAAAGCCGGAAAGCATAATCGACGGCAGAAAAGCGGCGTTCAGGGCGACCTGAGCGGCATTGAACTGGTTGCGGGTAATCGTGGAAATCAGCAGCCCCATcc                                                                            <  2:1762290/150‑1 (MQ=255)
                                                                gTACGTCACCGCGCGGATCACCGCGGGCATACTGTCGATCTGAAAAATAAAGCCGGAAAGCATAATCGACGGCAGAAAAGCGGCGTTCAGGGCGACCTGAGCGGCATTGAACTGGTTGCGGGTAATCGTGGAAATCAGCAGCCCCATccc                                                                           <  2:1410925/150‑1 (MQ=255)
                                                                  aCGTCACCGCGCGGATCACCGCGGGCATACTGTCGATCTGAAAAATAAAGCCGGAAAGCATAATCGACGGCAGAAAAGCGGCGTTCAGGGCGACCTGAGCGGCATTGAACTGGTTGCGGGTAATCGTGGAAATCAGCAGCCCCATCCCCa                                                                         >  2:616356/1‑150 (MQ=255)
                                                                   cGTCACCGCGCGGATCACCGCGGGCATACTGTCGATCTGAAAAATAAAGCCGGAAAGCATAATCGACGGCAGAAAAGCGGCGTTCAGGGCGACCTGAGCGGCATTGAACTGGTTGCGGGTAATCGTGGAAATCAGCAGCCCCATCCCCAg                                                                        >  1:29796/1‑150 (MQ=255)
                                                                        ccGCGCGGATCACCGCGGGCATACTGTCGATCTGAAAAATAAAGCCGGAAAGCATAATCGACGGCAGAAAAGCGGCGTTCAGGGCGACCTGAGCGGCATTGAACTGGTTGCGGGTAATCGTGGAAATCAGCAGCCCCATCCCCAGGGTAc                                                                   <  1:1741770/150‑1 (MQ=255)
                                                                        ccGCGCGGATCACCGCGGGCATACTGTCGATCTGAAAAATAAAGCCGGAAAGCATAATCGACGGCAGAAAAGCGGCGTTCAGGGCGACCTGAGCGGCATTGAACTGGTTGCGGGTAATCGTGGAAATCAGCAGCCCCATCCCCAGGGTAc                                                                   >  1:1407710/1‑150 (MQ=255)
                                                                             cgGATCACCGCGGGCATACTGTCGATCTGAAAAATAAAGCCGGAAAGCATAATCGACGGCAGAAAAGCGGCGTTCAGGGCGACCTGAGCGGCATTGAACTGGTTGCGGGTAATCGTGGAAATCAGCAGCCCCATCCCCAGGGTACTGAGt                                                              <  2:1654383/150‑1 (MQ=255)
                                                                                aTCACCGCGGGCATACTGTCGATCTGAAAAATAAAGCCGGAAAGCATAATCGACGGCAGAAAAGCGGCGTTCAGGGCGACCTGAGCGGCATTGAACTGGTTGCGGGTAATCGTGGAAATCAGCAGCCCCATCCCCAGGGTACTGAGTaaa                                                           >  2:291075/1‑150 (MQ=255)
                                                                                 tCACCGCGGGCATACTGTCGATCTGAAAAATAAAGCCGGAAAGCATAATCGACGGCAGAAAAGCGGCGTTCAGGGCGACCTGAGCGGCATTGAACTGGTTGCGGGTAATCGTGGAAATCAGCAGCCCCATCCCCAGGGTACTGAGTaaaa                                                          <  2:154653/150‑1 (MQ=255)
                                                                                       cgGGCATACTGTCGATCTGAAAAATAAAGCCGGAAAGCATAATCGACGGCAGAAAAGCGTCGTTCAGGGCGACCTGAGCGGCATTGAACTGGTTGCGGGTAATCGTGGAAATCAGCAGCCCCATCCCCAGGGTACTGAGTAAAAACAGGc                                                    <  1:1932262/150‑1 (MQ=255)
                                                                                          gCATACTGTCGATCTGAAAAATAAAGCCGGAAAGCATAATCGACGGCAGAAAAGCGGCGTTCAGGGCGACCTGAGCGGCATTGAACTGGTTGCGGGTAATCGTGGAAATCAGCAGCCCCATCCCCAGGGTACTGAGTAAAAACAGGCTgg                                                 >  1:1429993/1‑150 (MQ=255)
                                                                                          gCATACTGTCGATCTGAAAAATAAAGCCGGAAAGCATAATCGACGGCAGAAAAGCGGCGTTCAGGGCGACCTGAGCGGCATTGAACTGGTTGCGGGTAATCGTGGAAATCAGCAGCCCCATCCCCAGGGTACTGAGTAAAAACAGGCTgg                                                 <  1:1008196/150‑1 (MQ=255)
                                                                                            aTACTGTCGATCTGAAAAATAAAGCCGGAAAGCATAATCGACGGCAGAAAAGCGGCGTTCAGGGCGACCTGAGCGGCATTGAACTGGTTGCGGGTAATCGTGGAAATCAGCAGTCCCATCCCCTGGTTACTGAGTAAATACAGGCTGgag                                               <  2:341837/150‑1 (MQ=255)
                                                                                            aTACTGTCGATCTGAAAAATAAAGCCGGAAAGCATAATCGACGGCAGAAAAGCGGCGTTCAGGGCGACCTGAGCGGCATTGAACTGGTTGCGGGTAATCGTGGAAATCAGCAGCCCCATCCCCAGGGTACTGAGTAAAAACAGGGTGgag                                               <  2:628209/150‑1 (MQ=255)
                                                                                            aTACTGTCGATCTGAAAAATAAAGCCGGAAAGCATAATCGACGGCAGAAAAGCGGCGTTCAGGGCGACCTGAGCGGCATTGAACTGGTTGCGGGTAATCGTGGAAATCAGCAGCCCCATCCCCAGGGTACTGAGTAAAAACAGGCTGgag                                               <  1:1628101/150‑1 (MQ=255)
                                                                                              aCTGTCGATCTGAAAAATAAAGCCGGAAAGCATAATCGACGGCAGAAAAGCGGCGTTCAGGGCGACCTGAGCGGCATTGAACTGGTTGCGGGTAATCGTGGAAATCAGCAGCCCCATCCCCAGGGTACTGAGTAAAAACAGGCTGGAGAt                                             >  1:541023/1‑150 (MQ=255)
                                                                                               cTGGCGATCTGAAAAATAAAGCCGGAAAGCATAATCGACGGCAGAAAAGCGGCGTTCAGGGCGACCTGAGCGGCATTGAACTGGTTGCGGGTAATCGTGGAAATCAGCAGCCCCATCCCCAGGGTACTGAGTAAAAACAGGCTGGAGATa                                            <  2:1325036/150‑1 (MQ=255)
                                                                                                  tCGATCTGAAAAATAAAGCCGGAAAGCATAATCGACGGCAGAAAAGCGGCGTTCAGGGCGACCTGAGCGGCATTGAACTGGTTGCGGGTAATCGTGGAAATCAGCAGCCCCATCCCCAGGGTACTGAGTAAAAACAGGCTGGAGATaaaa                                         <  2:29796/150‑1 (MQ=255)
                                                                                                   cGATCTGAAAAATAAAGCCGGAAAGCATAATCGACGGCAGAAAAGCGGCGTTCAGGGCGACCTGAGCGGCATTGAACTGGTTGCGGGTAATCGTGGAAATCAGCAGCCCCATCCCCAGGGTACTGAGTAAAAACAGGCTGGAGATaaaaa                                        >  2:677742/1‑150 (MQ=255)
                                                                                                     aTCTGAAAAATAAAGCCGGAAAGCATAATCGACGGCAGAAAAGCGGCGTTCAGGGCGACCTGAGCGGCATTGAACTGGTTGCGGGTAATCGTGGAAATCAGCAGCCCCATCCCCAGGGTACTGAGTAAAAACAGGCTGGAGATAAAAAAc                                      <  1:1158061/150‑1 (MQ=255)
                                                                                                      tCTGAAAAATAAAGCCGGAAAGCATAATCGACGGCAGAAAAGCGGCGTTCAGGGCGACCTGAGCGGCATTGAACTGGTTGCGGGTAATCGTGGAAATCAGCAGCCCCATCCCCAGGGTACTGAGTAAAAACAGGCTGGAGATAAAAAACa                                     <  1:616356/150‑1 (MQ=255)
                                                                                                        tGAAAAATAAAGCCGGAAAGCATAATCGACGGCAGAAAAGCGGCGTTCAGGGCGACCTGAGCGGCATTGAACTGGTTGCGGGTAATCGTGGAAATCAGCAGCCCCATCCCCAGGGTACTGAGTAAAAACAGGCTGGAGATAAAAAACAGa                                   <  1:962058/150‑1 (MQ=255)
                                                                                                         gAAAAATAAAGCCGGAAAGCATAATCGACGGCAGAAAAGCGGCGTTCAGGGCGACCTGAGCGGCATTGAACTGGTTGCGGGTAATCGTGGAAATCAGCAGCCCCATCCCCAGGGTACTGAGTAAAAACAGGCTGGAGATAAAAAACAGaa                                  >  2:18114/1‑150 (MQ=255)
                                                                                                            aaaTAAAGCCGGAAAGCATAATCGACGGCAGAAAAGCGGCGTTCAGGGCGACCTGAGCGGCATTGAACTGGTTGCGGGTAATCGTGGAAATCAGCAGCCCCATCCCCAGGGTACTGAGTAAAAACAGGCTGGAGATAAAAAACAGAATca                               <  1:1435416/150‑1 (MQ=255)
                                                                                                            aaaTAAAGCCGGAAAGCATAATCGACGGCAGAAAAGCGGCGTTCAGGGCGACCTGAGCGGCATTGAACTGGTTGCGGGTAATCGTGGAAATCAGCAGCCCCATCCCCAGGGTACTGAGTAAAAACAGGCTGGAGATAAAAAACAGAATca                               <  1:18114/150‑1 (MQ=255)
                                                                                                             aaTAAAGCCGGAAAGCATAATCGACGGCAGAAAAGCGGCGTTCAGGGCGACCTGAGCGGCATTGAACTGGTTGCGGGTAATCGTGGAAATCAGCAGCCCCATCCCCAGGGTACTGAGTAAAAACAGGCTGGAGATAAAAAACAGAATcag                              <  1:1500575/150‑1 (MQ=255)
                                                                                                                 aaGCCGGAAAGCATAATCGACGGCAGAAAAGCGGCGTTCAGGGCGACCTGAGCGGCATTGAACTGGTTGCGGGTAATCGTGGAAATCAGCAGCCCCATCCCCAGGGTACTGAGTAAAAACAGGCTGGAGATAAAAAACAGAATCagcagc                          <  2:1909640/150‑1 (MQ=255)
                                                                                                                   gCCGGAAAGCATAATAGACGGCAGAAAAGCGGCGATCAGGGCGACCTGAGCGGCATTGAACTGGTTGCGGGTAATCGTGGAAATCAGCAGCCCCATCCCCAGGGTACTGAGTAAAAACAGGCTGGAGATAAAAAACAGAATCAGCAGCGa                        >  2:1660692/1‑150 (MQ=255)
                                                                                                                        aaaGCATAATCGACGGCAGAAAAGCGGCGTTCAGGGCGACCTGAGCGGCATTGAACTGGTTGCGGGTAATCGTGGAAATCAGCAGCCCCATCCCCAGGGTACTGAGTAAAAACAGGCTGGAGATAAAAAACAGAATCAGCAGCGACCcgc                   >  1:1915298/1‑150 (MQ=255)
                                                                                                                         aaGCATAATCGACGGCAGAAAAGCGGCGTTCAGGGCGACCTGAGCGGCATTGAACTGGTTGCGGGTAATCGTGGAAATCAGCAGCCCCATCCCCAGGGTACTGAGTAAAAACAGGCTGGAGATAAAAAACAGAATCAGCAGCGACCcgcg                  >  2:370642/1‑150 (MQ=255)
                                                                                                                                aTCGACGGCAGAAAAGCGGCGTTCAGGGCGACCTGAGCGGCATTGAACTGGTTGCGGGTAATCGTGGAAATCAGCAGCCCCATCCCCAGGGTACTGAGTAAAAACAGGCTGGAGATAAAAAACAGAATCAGCAGCGACCCGCGATACGGc           >  1:422449/1‑150 (MQ=255)
                                                                                                                                        cAGAAAAGCGGCATTCAGGGCGACCTGAGCGGCATTGAACTGGTTGCGGGTAATCGTGGAAATCAGCAGCCCCATCCCCAGGGTACTGAGTAAAAACAGGCTGGAGATAAAAAACAGAATCAGCAGCGAcc                      >  1:721677/1‑131 (MQ=255)
                                                                                                                                        cAGAAAAGCGGCATTCAGGGCGACCTGAGCGGCATTGAACTGGTTGCGGGTAATCGTGGAAATCAGCAGCCCCATCCCCAGGGTACTGAGTAAAAACAGGCTGGAGATAAAAAACAGAATCAGCAGCGAcc                      >  1:728530/1‑131 (MQ=255)
                                                                                                                                        cAGAAAAGCGGCATTCAGGGCGACCTGAGCGGCATTGAACTGGTTGCGGGTAATCGTGGAAATCAGCAGCCCCATCCCCAGGGTACTGAGTAAAAACAGGCTGGAGATAAAAAACAGAATCAGCAGCGAcc                      <  2:721677/131‑1 (MQ=255)
                                                                                                                                        cAGAAAAGCGGCATTCAGGGCGACCTGAGCGGCATTGAACTGGTTGCGGGTAATCGTGGAAATCAGCAGCCCCATCCCCAGGGTACTGAGTAAAAACAGGCTGGAGATAAAAAACAGAATCAGCAGCGAcc                      <  2:728530/131‑1 (MQ=255)
                                                                                                                                         aGAAAAGCGGCGTTCAGGGCGACCTGAGCGGCATTGAACTGGTTGCGGGTAATCGTGGAAATCAGCAGCCCCATCCCCAGGGTACTGAGTAAAAACAGGCTGGAGATAAAAAACAGAATCAGCAGCGACCCGCGATACGGCACGCCgaga  <  1:1981595/150‑1 (MQ=255)
                                                                                                                                              aGCGGCGTTCAGGGCGACCTGAGCGGCATTGAACTGGTTGCGGGTAATCGTGGAAATCAGCAGCCCCATCCCCAGGGTACTGAGTAAAAACAGGCTGGAGATaaa                                          >  1:1012344/1‑105 (MQ=255)
                                                                                                                                              aGCGGCGTTCAGGGCGACCTGAGCGGCATTGAACTGGTTGCGGGTAATCGTGGAAATCAGCAGCCCCATCCCCAGGGTACTGAGTAAAAACAGGCTGGAGATaaa                                          <  2:1012344/105‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                          
CTGGAATATTCCCGGCGAGGAACAGGCTTTGCAGGGTGCTGACGAAATAACGAGCGGGAATAATGTACGTCACCGCGCGGATCACCGCGGGCATACTGTCGATCTGAAAAATAAAGCCGGAAAGCATAATCGACGGCAGAAAAGCGGCGTTCAGGGCGACCTGAGCGGCATTGAACTGGTTGCGGGTAATCGTGGAAATCAGCAGCCCCATCCCCAGGGTACTGAGTAAAAACAGGCTGGAGATAAAAAACAGAATCAGCAGCGACCCGCGATACGGCACGCCGAGA  >  NC_000913/826211‑826497

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