Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 642,862 G→C 5.5% D259H (GAT→CAT)  ybdR → putative Zn‑dependent NAD(P)‑binding oxidoreductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913642,8620GC 5.5% 261.2 / 6.0 73D259H (GAT→CAT) ybdRputative Zn‑dependent NAD(P)‑binding oxidoreductase
Reads supporting (aligned to +/- strand):  ref base G (34/35);  new base C (2/2);  total (36/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.21e-01

GGTGCGGAACAGATTTTTGTTGTTGATCATCATCCCTACCGCTTGCATTTCGCCGCCGACCGCTACGGCGCGATCCCGATTAATTTTGATGAAGACAGCGATCCGGCACAGTCAATTATTGAACAAACGGCAGGTCACCGGGGCGTGGATGCAGTAATAGACGCCGTCGGTTTTGAAGCGAAAGGCAGCACCACGGAAACGGTGCTGACTAACCTGAAACTGGAGGGCAGCAGCGGTAAAGCGTTGCGTCAGTGTATTGCGGCGGTCAGGCGTGG  >  NC_000913/642715‑642989
                                                                                                                                                   |                                                                                                                               
ggTGCGGAACAGATTTTTGTTGTTGATCATCATCCCTACCGCTTGCATTTCGCCGCCGACCGCTACGGCGCGATCCCGATTAATTTTGATGAAGACAGCGATCCGGCACAGTCAATTATTGAACAAACGGCAGGTCACCGGGGCGTgg                                                                                                                                 <  4:924498/148‑1 (MQ=255)
     ggAACAGATTTTTGTTGTTGATCATCATCCCTACCGCTTGCATTTCGCCGCCGACCGCTACGGCGCGATCCCGATTAATTTTGATGAAGACAGCGATCCGGCACAGTCAATTATTGAACAAACGGCAGGTCACCGGGGCGTGGATGCa                                                                                                                            <  4:68344/148‑1 (MQ=255)
            aTTTTTGTTGTTGATCATCATCCCTACCGCTTGCATTTCGCCGCCGACCGCTACGGCGCGATCCCGATTAATTTTGATGAAGACAGCGATCCGGCACAGTCAATTATTGAACAAACGGCAGGTCACCGGGGCGTGGATGCAGTAATAg                                                                                                                     >  4:152641/1‑148 (MQ=255)
                   ttgttgATCATCATCCCTACCGCTTGCATTTCGCCGCCGACCGCTACGGCGCGATCCCGATTAATTTTGATGAAGACAGCGATCCGGCACAGTCAATTATTGAACAAACGGCAGGTCACCGGGGCGTGGATGCAGTAATAGACGCCGt                                                                                                              <  3:152641/148‑1 (MQ=255)
                      ttgATCATCATCCCTACCGCTTGCATTTCGCCGCCGACCGCTACGGCGCGATCCCGATTAATTTTGATGAAGACAGCGATCCGGCACAGTCAATTATTGAACAAACGGCAGGTCACCGGGGCGTGGATGCAGTAATAGACGCCGTCgg                                                                                                           >  4:776070/1‑148 (MQ=255)
                         atcatcatcCCTACCGCTTGCATTTCGCCGCCGACCGCTACGGCGCGATCCCGATTAATTTTGATGAAGACAGCGATCCGGCACAGTCAATTATTGAACAAACGGCAGGTCAACGGGGCGTGGATGCAGTAAt                                                                                                                       <  1:691223/133‑1 (MQ=255)
                         aacatcatcCCTACCGCTTGCATTTCGCCGCCGACCGCTACGGCGCGATCCCGATTAATTTTGATGAAGACAGCGATCCGGCACAGTCAATTATTGAACAAACGGCAGGTCACCGGGGCGTGGATGCAGTAAt                                                                                                                       >  2:691223/3‑133 (MQ=255)
                            atcatcCCTACCGCTTGCATTTCGCCGCCGACCGCTACGGCGCGATCCCGATTAATTTTGATGAAGACAGCGATCCGGCACAGTCAATTATTGAACAAACGGCAGGTCACCGGGGCGTGGATGCAGt                                                                                                                          >  1:1030026/1‑127 (MQ=255)
                            atcatcCCTACCGCTTGCATTTCGCCGCCGACCGCTACGGCGCGATCCCGATTAATTTTGATGAAGACAGCGATCCGGCACAGTCAATTATTGAACAAACGGCAGGTCACCGGGGCGTGGATGCAGt                                                                                                                          <  2:1030026/127‑1 (MQ=255)
                              catcCCTACCGCTTGCATTTCGCCGCCGACCGCTACGGCGCGATCCCGATTAATTTTGATGAAGACAGCGATCCGGCACAGTCAATTATTGAACAAACGGCAGGTCACCGGGGCGTGGATGCAGTATTAGACGCCGTCGGTTTTGa                                                                                                     <  1:895095/146‑1 (MQ=255)
                              catcCCTACCGCTTGCATTTCGCCGCCGACCGCTACGGCGCGATCCCGATTAATTTTGATGAAGACAGCGATCCGGCACAGTCAATTATTGAACAAACGGCAGGTCACCGGGGCGTGGATGCAGTAATAGACGCCGTCGGTTTTGa                                                                                                     >  2:895095/1‑146 (MQ=255)
                                  ccTACCGCTTGCATTTCGCCGCCGACCGCTACGGCGCGATCCCGATTAATTTTGATGAAGACAGCGATCCGGCACAGTCAATTATTGAACAAACGGCAGGTCACCGGGGCGTGGATGCAGTAATAGa                                                                                                                    <  1:884585/127‑1 (MQ=255)
                                  ccTACCGCTTGCATTTCGCCGCCGACCGCTACGGCGCGATCCCGATTAATTTTGATGAAGACAGCGATCCGGCACAGTCAATTATTGAACAAACGGCAGGTCACCGGGGCGTGGATGCAGTAATAGa                                                                                                                    >  2:884585/1‑127 (MQ=255)
                                       cGCTTGCATTTCGCCGCCGACCGCTACGGCGCGATCCCGATTAATTTTGATGAAGACAGCGATCCGGCACAGTCAATTATTGAACAAACGGCAGGTCACCGGGGCGTGGAt                                                                                                                               >  4:66053/1‑111 (MQ=255)
                                       cGCTTGCATTTCGCCGCCGACCGCTACGGCGCGATCCCGATTAATTTTGATGAAGACAGCGATCCGGCACAGTCAATTATTGAACAAACGGCAGGTCACCGGGGCGTGGAt                                                                                                                               <  3:66053/111‑1 (MQ=255)
                                           tGCATTTCGCCGCCGACCGCTACGGCGCGATCCCGATTAATTTTGATGAAGACAGCGATCCGGCACAGTCAATTATTGAACAAACGGCAGGTCACCGGGGCGTGGATGCAGTAATAg                                                                                                                     <  4:1022318/117‑1 (MQ=255)
                                           tGCATTTCGCCGCCGACCGCTACGGCGCGATCCCGATTAATTTTGATGAAGACAGCGATCCGGCACAGTCAATTATTGAACAAACGGCAGGTCACCGGGGCGTGGATGCAGTAATAg                                                                                                                     >  3:1022318/1‑117 (MQ=255)
                                                     cgccgACCGCTATGGCGCGATCCCGATTAATTTTGATGAAGACAGCGCTCCGGCACAGTCAATTATTGAACAAACGGCAGGTCACCGGGGCGTGGATGCAGTAATAGACGCCGTCGGTTTTGAAGCGAAAGGCAGCACCACGGAAACg                                                                            >  2:792479/1‑148 (MQ=255)
                                                      gccgACCGCTATGGCGCGATCCCGATTAATTTTGATGAAGACAGCGATCCGGCACAGTCAATTATTGAACAAACGGCAGGTCACCGGGGCGTGGATGCAGTAATAGACGCCGTCGGTTTTGAAGCGAAAGGCAGCACCACGGAAACgg                                                                           <  1:792479/148‑1 (MQ=255)
                                                        cgACCGCTACGGCGCGATCCCGATTAATTTTGATGAAGACAGCGATCCGGCACAGTCAATTATTGAACAAACGGCAGGTCACCGGGGCGTGGATGCAGTAATAGACGCCGTCGGTTTTGAAGCGAAAGGCAGCACCACGGAAACGGTg                                                                         <  3:776070/148‑1 (MQ=255)
                                                         gACCGCTACGGCGCGATCCCGATTAATTTTGATGAAGACAGCGATCCGGCACAGTCAATTATTGAACAAACGGCAGGTCACCGGGGCGTGGATGCAGTAATAGACGCCGTCGGTTTTGAAGCGAAAGGCAGCACCACGGAAACGGTGc                                                                        >  1:946800/1‑148 (MQ=255)
                                                           ccGCTACGGCGCGATCCCGATTAATTTTGATGAAGACAGCGATCCGGCACAGTCAATTATTGAACAAACGGCAGGTCACCGGGGCGTGGATGCAGTAATAGACGCCGTCGGTTTTGAAGCGAAAGGCAGCACCACGGAAACGGTGCTg                                                                      <  2:946800/148‑1 (MQ=255)
                                                              cTACGGCGCGATCCCGATTAATTTTGATGAAGACAGCGATCCGGCACAGTCAATTATTGAACAAACGGCAGGTCACCGGGGCGTGGATGCAGTAATAGACGCCGTCGGt                                                                                                          <  3:639599/109‑1 (MQ=255)
                                                              cTACGGCGCGATCCCGATTAATTTTGATGAAGACAGCGATCCGGCACAGTCAATTATTGAACAAACGGCAGGTCACCGGGGCGTGGATGCAGTAATAGACGCCGTCGGt                                                                                                          >  4:639599/1‑109 (MQ=255)
                                                                              aTTAATTTTGATGAAGACAGCGATCCGGCACAGTCAATTATTGAACAAACGGCAGGTCACCGGGGCGTGGATGCAGTAATAGACGCCGTCGGTTTTGAAGCGAAAGGCAGCACCACGGAAACGGTGCTGACTaa                                                                 <  3:547271/134‑1 (MQ=255)
                                                                              aTTAATTTTGATGAAGACAGCGATCCGGCACAGTCAATTATTGAACAAACGGCAGGTCACCGGGGCGTGGATGCAGTAATAGACGCCGTCGGTTTTGAAGCGAAAGGCAGCACCACGGAAACGGTGCTGACTaa                                                                 >  4:547271/1‑134 (MQ=255)
                                                                               ttAATTTTGATGAAGACAGCGATCCGGCACAGTCAATTATTGAACAAACGGCAGGTCACCGGGGCGTGGATGCAGTAATAGACGCCGTCGGTTTTGAAGCGAAAGGCAGCACCACGGAAACGGTGCTGACTAACCTGAAACTGGAggg                                                  >  3:957012/1‑148 (MQ=255)
                                                                                tAATTTTGATGAAGACAGCGATCCGGCACAGTCAATTATTGAACAAACGGCAGGTCACCGGGGCGTGGATGCAGTAATAGACGCCGTCGGTTTTGAAGCGAAAGGCAGCACCACGGAAACGGTGc                                                                        >  3:1015084/1‑125 (MQ=255)
                                                                                tAATTTTGATGAAGACAGCGATCCGGCACAGTCAATTATTGAACAAACGGCAGGTCACCGGGGCGTGGATGCAGTAATAGACGCCGTCGGTTTTGAAGCGAAAGGCAGCACCACGGAAACGGTGc                                                                        <  4:1015084/125‑1 (MQ=255)
                                                                                    tttGATGAAGACAGCGATCCGGCACAGTCAATTATTGAACAAACGGCAGGTCACCGGGGCGTGGATGCAGTAATAGACGCCGTCGGTTTTGAAGCGAAAGGCAGCACCACGGAAACGGTGCTGACTaa                                                                 >  4:167469/1‑128 (MQ=255)
                                                                                    tttGATGAAGACAGCGATCCGGCACAGTCAATTATTGAACAAACGGCAGGTCACCGGGGCGTGGATGCAGTAATAGACGCCGTCGGTTTTGAAGCGAAAGGCAGCACCACGGAAACGGTGCTGACTaa                                                                 <  3:167469/128‑1 (MQ=255)
                                                                                      tgatgaAGACAGCGATCCGGCACAGTCAATTATTGAACAAACGGCAGGTCACCGGGGCGTGGATGCAGTAATAGACGCCGTCGGTTTTGAAGCGAAAGGCAGCACCACGGAAACGGTGCTg                                                                      <  3:142107/121‑1 (MQ=255)
                                                                                      tgatgaAGACAGCGATCCGGCACAGTCAATTATTGAACAAACGGCAGGTCACCGGGGCGTGGATGCAGTAATAGACGCCGTCGGTTTTGAAGCGAAAGGCAGCACCACGGAAACGGTGCTg                                                                      >  4:142107/1‑121 (MQ=255)
                                                                                          gaAGACAGCGATCCGGCACAGTCAATTATTGAACAAACGGCAGGTCACCGGGGCGTGGATGCAGTAATAGACGCCGTCGGTTTTGAAGCGAAAGGCAGCACCACGGAAACGGTg                                                                         >  3:774031/1‑114 (MQ=255)
                                                                                          gaAGACAGCGATCCGGCACAGTCAATTATTGAACAAACGGCAGGTCACCGGGGCGTGGATGCAGTAATAGACGCCGTCGGTTTTGAAGCGAAAGGCAGCACCACGGAAACGGTg                                                                         <  4:774031/114‑1 (MQ=255)
                                                                                          gaAGACAGCGATCCGGCACAGTCAATTATTGAACAAACGGCAGGTCACCGGGGCGTGGATGCAGTAATAGACGCCGTCGGTTTTGAAGCGAAAGGCAGCACCACGGAAACGGTGCTGACTAACCTGAAACTGGAGGGCAGCAGCGGTa                                       >  1:168486/1‑148 (MQ=255)
                                                                                          gaAGACAGCGATCCGGCACAGTCAATTATTGAACAAACGGCAGGTCACCGGGGCGTGGATGCAGTAATAGACGCCGTCGGTTTTGAAGCGAAAGGCAGCACCACGGAAACGGTGCTGACTAACCTGAAACTGGAGGGCAGCAGCGGTa                                       >  4:956751/1‑148 (MQ=255)
                                                                                           aaGACAGCGATCCGGCACAGTCAATTATTGAACAAACGGCAGGTCACCGGGGCGTGGATGCAGTAATAGACGCCg                                                                                                               >  1:557914/1‑75 (MQ=255)
                                                                                           aaGACAGCGATCCGGCACAGTCAATTATTGAACAAACGGCAGGTCACCGGGGCGTGGATGCAGTAATAGACGCCg                                                                                                               <  2:557914/75‑1 (MQ=255)
                                                                                            aGACAGCGATCCGGCACAGTCAATTATTGAACAAACGGCAGGTCACCGGGGCGTGGATGCAGTAATAGACGCCGTCGGTTTTGAAGCGAAAGGCAGCACCACGGAAACGGTg                                                                         >  3:324813/1‑112 (MQ=255)
                                                                                            aGACAGCGATCCGGCACAGTCAATTATTGAACAAACGGCAGGTCACCGGGGCGTGGATGCAGTAATAGACGCCGTCGGTTTTGAAGCGAAAGGCAGCACCACGGAAACGGTg                                                                         <  4:324813/112‑1 (MQ=255)
                                                                                                aGCGATCCGGCACAGTCAATTATTGAACAAACGGCAGGTCACCGGGGCGTGGATGCAGTAATAGACGCCGTCGGTTTTGAAGCGAAAGGCAGCACCACGGAAACGGTGCTGACTAACCTGAAACTGGAGGGCAGCAGCGGTAAAGCGt                                 <  2:168486/148‑1 (MQ=255)
                                                                                                aGCGATCCGGCACAGTCAATTATTGAACAAACGGCAGGTCACCGGGGCGTGGATGCAGTAATAGACGCCGTCGGTTTTGAAGCGAAAGGCAGCACCACGGAAACGGTGCTGACTAACCTGAAACTGGAGGGCAGCAGCGGTAAAGCGt                                 <  3:956751/148‑1 (MQ=255)
                                                                                                 gCGATCCGGCACAGTCAATTATTGAACAAACGGCAGGTCACCGGGGCGTGGATGCAGTAATAGACGCCGTCGGTTTTGAAGCGAAAGGCAGCACCACGGAAACGGTGCTGACTAACCTGAAACTg                                                       >  2:996524/1‑125 (MQ=255)
                                                                                                 gCGATCCGGCACAGTCAATTATTGAACAAACGGCAGGTCACCGGGGCGTGGATGCAGTAATAGACGCCGTCGGTTTTGAAGCGAAAGGCAGCACCACGGAAACGGTGCTGACTAACCTGAAACTg                                                       <  1:996524/125‑1 (MQ=255)
                                                                                                  cGATCCGGCACAGTCAATTATTGAACAAACGGCAGGTCACCGGGGCGTGGATGCAGTAATAGACGCCGTCGGTTTTGAAGCGAAAGGCAGCACCACGGAAACGGTGCTGACTAACCTGAAACTGGAGGGCAGCAGCGGTAAAGCGTTg                               <  4:957012/148‑1 (MQ=255)
                                                                                                    aTCCGGCACAGTCAATTATTGAACAAACGGCAGGTCACCGGGGCGTGGATGCAGTAATAGACGCCGTCGGt                                                                                                          <  2:175773/71‑1 (MQ=255)
                                                                                                    aTCCGGCACAGTCAATTATTGAACAAACGGCAGGTCACCGGGGCGTGGATGCAGTAATAGACGCCGTCGGt                                                                                                          >  1:175773/1‑71 (MQ=255)
                                                                                                          cacaGTCAATTATTGAACAAACGGCAGGTCACCGGGGCGTGGATGCAGTAATAGACGCCGTCGGTTTTGAAGCGAAAGGCAGCACCACGGAAACGGt                                                                          <  4:249613/97‑1 (MQ=255)
                                                                                                          cacaGTCAATTATTGAACAAACGGCAGGTCACCGGGGCGTGGATGCAGTAATAGACGCCGTCGGTTTTGAAGCGAAAGGCAGCACCACGGAAACGGt                                                                          >  3:249613/1‑97 (MQ=255)
                                                                                                          cacaGTCAATTATTGAACAAACGGCAGGTCACCGGGGCGTGGATGCAGTAATAGACGCCGTCGGTTTTGAAGCGAAAGGCAGCACCACGGAAACGGTGCTGACTAACCTGAAACTGGAGGGCAGCAGCGGTAAAGCGTTGCGTCAgtg                       >  3:636032/1‑148 (MQ=255)
                                                                                                              gTCAATTATTGAACAAACGGCAGGTCACCGGGGCGTGGATGCAGTAATAGACGCCGTCGGTTTTGAAGCGAAAGGCAGCACCACGGAAACGGTGCTGACTa                                                                  >  3:869770/1‑101 (MQ=255)
                                                                                                              gTCAATTATTGAACAAACGGCAGGTCACCGGGGCGTGGATGCAGTAATAGACGCCGTCGGTTTTGAAGCGAAAGGCAGCACCACGGAAACGGTGCTGACTa                                                                  <  4:869770/101‑1 (MQ=255)
                                                                                                                cAATTATTGAACAAACGGCAGGTCACCGGGGCGTGGATGCAGTAATAGACGCCGTCGGTTTTGAAGCGAAAGGCAGCACCACGGAAACGGt                                                                          <  2:368284/91‑1 (MQ=255)
                                                                                                                cAATTATTGAACAAACGGCAGGTCACCGGGGCGTGGATGCAGTAATAGACGCCGTCGGTTTTGAAGCGAAAGGCAGCACCACGGAAACGGt                                                                          >  1:368284/1‑91 (MQ=255)
                                                                                                                 aaTTATTGAACAAACGGCAGGTCACCGGGGCGTGGATGCAGTAATAGACGCCGTCGGTTTTGAAGCGAAAGGCAGCACCACGGAAACGGTGCTGACTAACCTGAAACTGGAGGGCAGCAGCGGTAAAGCGTTGCGTCAGTGTATTgcg                <  4:636032/148‑1 (MQ=255)
                                                                                                                  attattGAACAAACGGCAGGTCACCGGGGCGTGGATGCAGTAATAGACGCCGTCGGTTTTGAAGCGaaa                                                                                              <  4:144360/69‑1 (MQ=255)
                                                                                                                  attattGAACAAACGGCAGGTCACCGGGGCGTGGATGCAGTAATAGACGCCGTCGGTTTTGAAGCGaaa                                                                                              >  3:144360/1‑69 (MQ=255)
                                                                                                                  attattGAACAAACGGCAGGTCACCGGGGCGTGGATGCAGTAATAGACGCCGTCGGTTTTGAAGCGAAAGGCAGCACCACGGAAACGGTGCTGACTAACCTGAAACTGGAggg                                                  >  2:564626/1‑113 (MQ=255)
                                                                                                                  attattGAACAAACGGCAGGTCACCGGGGCGTGGATGCAGTAATAGACGCCGTCGGTTTTGAAGCGAAAGGCAGCACCACGGAAACGGTGCTGACTAACCTGAAACTGGAggg                                                  <  1:564626/113‑1 (MQ=255)
                                                                                                                  attattGAACAAACGGCAGGTCACCGGGGCGTGGATGCAGTAATAGACGCCGTCGGTTTTGAAGCGAAAGGCAGCACCACGGAAACGGTGCTGACTAACCTGAAACTGGAGGGCAGCAGCGGTAAAGCGTTGCGTCAGTGTATTgcgg               >  3:815010/1‑148 (MQ=255)
                                                                                                                    tattGAACAAACGGCAGGTCACCGGGGCGTGGATGCAGTAATAGACGCCGTCGGTTTTGAAGCGAAAGGCAGCACCACGGAAACGGTGCTGACTaa                                                                 <  2:910372/96‑1 (MQ=255)
                                                                                                                    tattGAACAAACGGCAGGTCACCGGGGCGTGGATGCAGTAATAGACGCCGTCGGTTTTGAAGCGAAAGGCAGCACCACGGAAACGGTGCTGACTaa                                                                 >  1:910372/1‑96 (MQ=255)
                                                                                                                       tGAACAAACGGCAGGTCACCGGGGCGTGGATGCAGTAATAGACGCCGTCGGTTTTGAAGCGAAAGGCAGCACCACGGAAACGGTGCTGACTAACCTGAAACt                                                        <  1:105467/102‑1 (MQ=255)
                                                                                                                       tGAACAAACGGCAGGTCACCGGGGCGTGGATGCAGTAATAGACGCCGTCGGTTTTGAAGCGAAAGGCAGCACCACGGAAACGGTGCTGACTAACCTGAAACt                                                        >  2:105467/1‑102 (MQ=255)
                                                                                                                         aaCAAACGGCAGGTCACCGGGGCGTGGATGCAGTAATAGACGCCGTCGGTTTTGAAGCGAAAGGCAGCACCACGGAAAAGGTGCTGACTAACCTGAAACTGGAGGGCAGCAGCGGTAAAGCGt                                 <  4:393849/123‑1 (MQ=255)
                                                                                                                         aaCAAACGGCAGGTCACCGGGGCGTGGATGCAGTAATAGACGCCGTCGGTTTTGAAGCGAAAGGCAGCACCACGGAAAAGGTGCTGACTAACCTGAAACTGGAGGGCAGCAGCGGTAAAGCGt                                 >  3:393849/1‑123 (MQ=255)
                                                                                                                               cGGCAGGTCACCGGGGCGTGGATGCAGTAATAGACGCCGTCGGTTTTGAAGCGAAAGGCAGCACCACGGAAACGGTGCTGACTAACCTGAAACTGGAGGGCAGCAGCGGTAAAGCGTTGCGTCAGTGTATTGCGGCGGTCAGGCGTgg  <  4:815010/148‑1 (MQ=255)
                                                                                                                               cGGCAGGTCACCGGGGCGTGCATGCAGTAATAGACGCCGTCGGTTTTGAAGCGAAAGGCAGCACCACGGAAACGGTGCTGACTAACCTGAAACTGGAGGgcagca                                             <  4:248063/105‑1 (MQ=255)
                                                                                                                               cGGCAGGTCACCGGGGCGTGCATGCAGTAATAGACGCCGTCGGTTTTGAAGCGAAAGGCAGCACCACGGAAACGGTGCTGACTAACCTGAAACTGGAGGgcagca                                             <  4:216599/105‑1 (MQ=255)
                                                                                                                               cGGCAGGTCACCGGGGCGTGCATGCAGTAATAGACGCCGTCGGTTTTGAAGCGAAAGGCAGCACCACGGAAACGGTGCTGACTAACCTGAAACTGGAGGgcagca                                             >  3:248063/1‑105 (MQ=255)
                                                                                                                               cGGCAGGTCACCGGGGCGTGCATGCAGTAATAGACGCCGTCGGTTTTGAAGCGAAAGGCAGCACCACGGAAACGGTGCTGACTAACCTGAAACTGGAGGgcagca                                             >  3:216599/1‑105 (MQ=255)
                                                                                                                                 gCAGGTCACCGGGGCGTGGATGCAGTAATAGACGCCGTCGGTTTTGAAGCGAAAGGCAGCACCACg                                                                                  <  4:554095/66‑1 (MQ=255)
                                                                                                                                 gCAGGTCACCGGGGCGTGGATGCAGTAATAGACGCCGTCGGTTTTGAAGCGAAAGGCAGCACCACg                                                                                  >  3:554095/1‑66 (MQ=255)
                                                                                                                                                   gATGCAGTAATAGACGCCGTCGGTTTTGAAGCGAAAGGCAGCACCACGGAAACGGTGCTGACTAACCTGAAACTGGAGGGCAGCAGCGGTaa                                      >  4:526982/1‑92 (MQ=255)
                                                                                                                                                   gATGCAGTAATAGACGCCGTCGGTTTTGAAGCGAAAGGCAGCACCACGGAAACGGTGCTGACTAACCTGAAACTGGAGGGCAGCAGCGGTaa                                      <  3:526982/92‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                               
GGTGCGGAACAGATTTTTGTTGTTGATCATCATCCCTACCGCTTGCATTTCGCCGCCGACCGCTACGGCGCGATCCCGATTAATTTTGATGAAGACAGCGATCCGGCACAGTCAATTATTGAACAAACGGCAGGTCACCGGGGCGTGGATGCAGTAATAGACGCCGTCGGTTTTGAAGCGAAAGGCAGCACCACGGAAACGGTGCTGACTAACCTGAAACTGGAGGGCAGCAGCGGTAAAGCGTTGCGTCAGTGTATTGCGGCGGTCAGGCGTGG  >  NC_000913/642715‑642989

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