Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 457,664 G→T 5.6% G80V (GGC→GTC)  clpX → ATPase and specificity subunit of ClpX‑ClpP ATP‑dependent serine protease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913457,6640GT 5.6% 254.3 / 4.9 73G80V (GGC→GTC) clpXATPase and specificity subunit of ClpX‑ClpP ATP‑dependent serine protease
Reads supporting (aligned to +/- strand):  ref base G (32/36);  new base T (2/2);  total (35/38)
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.23e-01

ATCCGTGTATATCTGCGACGAATGTGTTGATTTATGTAACGACATCATTCGCGAAGAGATTAAAGAAGTTGCACCGCATCGTGAACGCAGTGCGCTACCGACGCCGCATGAAATTCGCAACCACCTGGACGATTACGTTATCGGCCAGGAACAGGCGAAAAAAGTGCTGGCGGTCGCGGTATACAACCATTACAAACGTCTGCGCAACGGCGATACCAGCAATGGCGTCGAGTTGGGCAAAAGTAACATTCTGCTGATCGGTCCGACCGGTTCCGGTAAAACGC  >  NC_000913/457521‑457804
                                                                                                                                               |                                                                                                                                            
aTCCGTGTATATCTGCGACGAATGTGTTGATTTATGTAACGACATCATTCGCGAAGAGATTAAAGAAGTTGCACCGCATCGTGAACGCAGTGCGCTACCGACGCCGCATGAAATTCGCAACCACCTGGCCGATAACGTAATCGTcc                                                                                                                                            <  1:158558/146‑1 (MQ=255)
aTCCGTGTATATCTGCGACGAATGTGTTGATTTATGTAACGACATCATTCGCGAAGAGATTAAAGAAGTTGCACCGCATCGTGAACGCAGTGCGCTACCGACGCCGCATGAAATTCGCAACCACCTGGCCGATAACGTAATCGTcc                                                                                                                                            >  2:158558/1‑146 (MQ=255)
     tgtATATCTGCGACGAATGTGTTGATTTATGTAACGACATCATTCGCGAAGAGATTAAAGAAGTTGCACCGCATCGTGAACGCAGTGCGCTACCGACGCCGCATGAAATTCGCAACCACCTGGACGATTACGTTATCGGCCAGGAACa                                                                                                                                     <  2:369718/148‑1 (MQ=255)
      gtATATCTGCGACGAATGTGTTGATTTATGTAACGACATCATTCGCGAAGAGATTAAAGAAGTTGCACCGCATCGTGAACGCAGTGCGCTACCGACGCCGCATGAAATTCGCAACCACCTGGACGATTACGTTATCgg                                                                                                                                              <  1:1124758/138‑1 (MQ=255)
      gtATATCTGCGACGAATGTGTTGATTTATGTAACGACATCATTCGCGAAGAGATTAAAGAAGTTGCACCGCATCGTGAACGCAGTGCGCTACCGACGCCGCATGAAATTCGCAACCACCTGGACGATTACGTTATCgg                                                                                                                                              >  2:1124758/1‑138 (MQ=255)
            cTGCGACGAATGTGTTGATTTATGTAACGACATCATTCGCGAAGAGATTAAAGAAGTTGCACCGCATCGTGAACGCAGTGCGCTACCGACGCCGCATGAAATTCGCAACCACCTGGACGATTACGTTATCGGCCAGGAACAGGCGaaa                                                                                                                              <  3:1132759/148‑1 (MQ=255)
                 acgaATGTGTTGATTTATGTAACGACATCATTCGCGAAGAGATTAAAGAAGTTGCACCGCATCGTGAACGCAGTGCGCTACCGACGCCGCATGAAATTCGCAACCACCTGGACGATTACGTTATCGGCCAGGAACAGGCGAAAAAAGt                                                                                                                         <  3:626162/148‑1 (MQ=255)
                     aTGTGTTGATTTATGTAACGACATCATTCGCGAAGAGATTAAAGAAGTTGCACCGCATCGTGAACGCAGTGCGCTACCGACGCCGCATGAAATTCGCAACCACCTGGACGATTACGTTATCGGcc                                                                                                                                            >  4:1021941/1‑125 (MQ=255)
                     aTGTGTTGATTTATGTAACGACATCATTCGCGAAGAGATTAAAGAAGTTGCACCGCATCGTGAACGCAGTGCGCTACCGACGCCGCATGAAATTCGCAACCACCTGGACGATTACGTTATCGGcc                                                                                                                                            <  3:1021941/125‑1 (MQ=255)
                          ttGATTTATGTAACGACATCATTCGCGAAGAGATTAAAGAAGTTGCACCGCATCGTGAACGCAGTGCGCTACCGACGCCGCATGAAATTCGCAACCACCTGGACGATTACGTTATCGGCCAGGAACAGGCGAAAAAAGTGCTGGCGGt                                                                                                                >  1:652704/1‑148 (MQ=255)
                           tGATTTATGTAACGACATCATTCGCGAAGAGATTAAAGAAGTTGCACCGCATCGTGAACGCAGTGCGCTACCGACGCCGCATGAAATTCGCAACCACCTGGACGATTACGTTATCGGCCAGGAACAGGCGaaaaa                                                                                                                            >  2:933043/1‑135 (MQ=255)
                           tGATTTATGTAACGACATCATTCGCGAAGAGATTAAAGAAGTTGCACCGCATCGTGAACGCAGTGCGCTACCGACGCCGCATGAAATTCGCAACCACCTGGACGATTACGTTATCGGCCAGGAACAGGCGaaaaa                                                                                                                            <  1:933043/135‑1 (MQ=255)
                             aTTTATGTAACGACATCATTCGCGAAGAGATTAAAGAAGTTGCACCGCATCGTGAACGCAGTGCGCTACCGACGCCGCATGAAATTCGCAACCACCTGGACGATTACGTTATCGGCCAGGAACAGGCGAAAAAAGTGCTGGCGGTcgc                                                                                                             <  2:652704/148‑1 (MQ=255)
                              tttATGTAACGACATCATTCGCGAAGAGATTAAAGAAGTTGCACCGCATCGTGAACGCAGTGCGCTACCGACGCCGCATGAAATTCGCAACCACCTGGACGATTACGTTATCGGCCAGGAACAg                                                                                                                                    >  1:881717/1‑124 (MQ=255)
                              tttATGTAACGACATCATTCGCGAAGAGATTAAAGAAGTTGCACCGCATCGTGAACGCAGTGCGCTACCGACGCCGCATGAAATTCGCAACCACCTGGACGATTACGTTATCGGCCAGGAACAg                                                                                                                                    <  2:881717/124‑1 (MQ=255)
                                                tCGCGAAGAGATTAAAGAAGTTGCACCGCATCGTGAACGCAGTGCGCTACCGACGCCGCATGAAATTCGCAACCACCTGGACGATTACGTTATCGGcc                                                                                                                                            <  1:708475/98‑1 (MQ=255)
                                                tCGCGAAGAGATTAAAGAAGTTGCACCGCATCGTGAACGCAGTGCGCTACCGACGCCGCATGAAATTCGCAACCACCTGGACGATTACGTTATCGGcc                                                                                                                                            >  2:708475/1‑98 (MQ=255)
                                                        agaTTAAAGAAGTTGCACCGCATCGTGAACGCAGTGCGCTACCGACGCCGCATGAAATTCGCAACCACCTGGACGATTACGTTATCGGCCAGGAACAGGCGAAAAAAGTGCt                                                                                                                      >  2:1019656/1‑112 (MQ=255)
                                                        agaTTAAAGAAGTTGCACCGCATCGTGAACGCAGTGCGCTACCGACGCCGCATGAAATTCGCAACCACCTGGACGATTACGTTATCGGCCAGGAACAGGCGAAAAAAGTGCt                                                                                                                      <  1:1019656/112‑1 (MQ=255)
                                                           ttAAAGAAGTTGCACCGCATCGTGAACGCAGTGCGCTACCGACGCCGCATGAAATTCGCAACCACCTGGACGATTACGTTATCGGCCAGGAACAGGCGAAAAAAGTGCTGGCGGTCGCGGTATACAACCATTACAAACGTCTGCGCaa                                                                               >  1:991418/1‑148 (MQ=255)
                                                              aagaagTTGCACCGCATCGTGAACGCAGTGCGCTACCGACGCCGCATGAAATTCGCAACCACCTGGACGATTACGTTATCGGCCAGGAACAGGCGAAAAAAGTGCTGGCGGTCGCGGTATACa                                                                                                     >  3:601705/1‑123 (MQ=255)
                                                              aagaagTTGCACCGCATCGTGAACGCAGTGCGCTACCGACGCCGCATGAAATTCGCAACCACCTGGACGATTACGTTATCGGCCAGGAACAGGCGAAAAAAGTGCTGGCGGTCGCGGTATACa                                                                                                     <  4:601705/123‑1 (MQ=255)
                                                               agaagTTGCACCGCATCGTGAACGCAGTGCGCTACCGACGCCGCATGAAATTCGCAACCACCTGGACGATTACGTTATCGGCCAGGAACAGGCGAAAAAAGTGCTGGCGGTCGCGGTATACAACCATTACAAACGTCt                                                                                     <  3:571818/138‑1 (MQ=255)
                                                               agaagTTGCACCGCATCGTGAACGCAGTGCGCTACCGACGCCGCATGAAATTCGCAACCACCTGGACGATTACGTTATCGGCCAGGAACAGGCGAAAAAAGTGCTGGCGGTCGCGGTATACAACCATTACAAACGTCt                                                                                     >  4:571818/1‑138 (MQ=255)
                                                                gaagTTGCACCGCATCGTGAACGCAGTGCGCTACCGACGCCGCATGAAATTCGCAACCACCTGGACGATTACGTTATCGGCCAGGAACAGGCGAAAAAAGTGCTGGCGGTcgc                                                                                                             <  2:559113/113‑1 (MQ=255)
                                                                gaagTTGCACCGCATCGTGAACGCAGTGCGCTACCGACGCCGCATGAAATTCGCAACCACCTGGACGATTACGTTATCGGCCAGGAACAGGCGAAAAAAGTGCTGGCGGTcgc                                                                                                             >  1:559113/1‑113 (MQ=255)
                                                                gaagTTGCACCGCATCGTGAACGCAGTGCGCTACCGACGCCGCATGAAATTCGCAACCACCTGGACGATTACGTTATCGGCCAAGAACAGGCGAAAAAAGTGCTGGCGGTcgc                                                                                                             <  4:340003/113‑1 (MQ=255)
                                                                gaagTTGCACCGCATCGTGAACGCAGTGCGCTACCGACGCCGCATGAAATTCGCAACCACCTGGACGATTACGTTATCGGCCAAGAACAGGCGAAAAAAGTGCTGGCGGTcgc                                                                                                             >  3:340003/1‑113 (MQ=255)
                                                                 aagTTGCACCGCATCGTGAACGCAGTGCGCTACCGACGCCGCATGAAATTCGCAACCACCTGGACGATTACGTTATCGGCCAGGAACAGGCGAAAAAAGTGCTGGCGGTCGCGGTATACAACCATTACaa                                                                                           <  3:20691/130‑1 (MQ=255)
                                                                 aagTTGCACCGCATCGTGAACGCAGTGCGCTACCGACGCCGCATGAAATTCGCAACCACCTGGACGATTACGTTATCGGCCAGGAACAGGCGAAAAAAGTGCTGGCGGTCGCGGTATACAACCATTACaa                                                                                           >  4:20691/1‑130 (MQ=255)
                                                                    ttGCACCGCATCGTGAACGCAGTGCGCTACCGACGCCGCATGAAATTCGCAACCACCTGGACGATTACGTTATCGGCCAGGAACAgg                                                                                                                                   <  3:30473/87‑1 (MQ=255)
                                                                    ttGCACCGCATCGTGAACGCAGTGCGCTACCGACGCCGCATGAAATTCGCAACCACCTGGACGATTACGTTATCGGCCAGGAACAgg                                                                                                                                   >  4:30473/1‑87 (MQ=255)
                                                                    ttGCACCGCATCGTGAACGCAGTGCGCTACCGACGCCGCATGAAATTCGCAACCACCTGGACGATTACGTTATCGGCCAGGAACAGGCGAAAAAAGTGc                                                                                                                       <  3:761701/99‑1 (MQ=255)
                                                                    tggCACCGCATCGTGAAGGCAGTGCGCTACCGACGCCGCATGAAATTCGCAACCACCTGGACGATTACGTTATCGGCCAGGAACAGGCGAAAAAAGTGc                                                                                                                       >  4:761701/3‑99 (MQ=255)
                                                                        aCCGCATCGTGAACGCAGTGCGCTACCGACGCCGCATGAAATTCGCAACCACCTGGCCGATAACGTAATCGTCCAgg                                                                                                                                         >  2:916964/1‑77 (MQ=37)
                                                                        aCCGCATCGTGAACGCAGGGCGCTACCGACGCCGCATGAAATTCGCAACCACCTGGCCGATAACGTAATCGTCCAgg                                                                                                                                         <  1:916964/77‑1 (MQ=37)
                                                                                gTGAACGCAGTGCGCTACCGACGCCGCATGAAATTCGCAACCACCTGGACGATTACGTTATCGGCCAGGAACAGGCGAAAAAAGTGCTGGCGGTCGCGGTATACAACCATTACAAACGTCTGCGCAACGGCGATACCAGCAATGGCGt                                                          <  4:125640/148‑1 (MQ=255)
                                                                                    aCGCAGTGCGCTACCGACGCCGCATGAAATTCGCAACCACCTGGACGATTACGTTATCGGCCAGGAACAGGCGAAAAAAGTGCTGGCGGTCGCGGTATACAACCATTACAAACGTCTGCGCAACGGCGATACCAGCAATGGCGTCGAg                                                      <  2:991418/148‑1 (MQ=255)
                                                                                     cGCAGTGCGCTACCGACGCCGCATGAAATTCGCAACCACCTGGACGATTACGTTATCGGCCAGGAACAGGCGAAAAAAGTGCTGGCGGTCGCGGTATACAACCATTACAAACGTCTGCGCAACGGCGATACCAGCAATGGCGTCGAGt                                                     >  3:1047284/1‑148 (MQ=255)
                                                                                                aCCGACGCCGCATGAAATTCGCAACCACCTGGACGATTACGTTATCGGCCAGGAACAGGCGAAAAAAGTGCTGGCGGTCGCGGTATACAACCATTACAAACGTCTGCGCAACGGCGATACCAGCAATGGCGTCGAGTTGGGCAAAAGt                                          >  3:181338/1‑148 (MQ=255)
                                                                                                    aCGCCGCATGGAATTCGCAACCACCTGGACGATTACGTTATCGGCCAGGAACAGGCGAAAAAAGTGCTGGCGGTCGCGGTATACAACCATTACAAACGTCTGCGCAACGGCGATACCAGCAATGGCGTCGAGTTGGGCaaa                                             >  1:594858/1‑141 (MQ=255)
                                                                                                    aCGCCGCATGAAATTCGCAACCACCTGGACGATTACGTTATCGGCCAGGAACAGGCGAAAAAAGTGCTGGCGGTCGCGGTATACAACCATTACAAACGTCTGCGCAACGGCGATACCAGCAATGGCGTCGAGTTGGGCaaa                                             <  2:594858/141‑1 (MQ=255)
                                                                                                            tGAAATTCGCAACCACCTGGACGATTACGTTCTCGCCCAGGAACAGGCGAAAAAAGTGCTGGCGGTCGCGGTATACAACCATTACAAACGTCTGCGCAACGGCGATACCAGCAATGGCGTCGAGTTGGGCAAAAGTAACATTctgctg                              >  1:726815/1‑148 (MQ=255)
                                                                                                             gAAATTCGCAACCACCTGGACGATTACGTTATCGGCCAGGAACAGGCGAAAAAAGTGCTGGCGGTCGCGGTATACAACCATTACAAACGTCTGCGCAACGGCGATACCAGCAATGGCGTCGAGTTGGGCAAAAGTAACATTCTGCTGa                             >  1:563140/1‑148 (MQ=255)
                                                                                                             gAAATTCGCAACCACCTGGACGATTACGTTATCGGCCAGGAACAGGCGAAAAAAGTGCTGGCGGTCGCGGTATACAACCATTACAAACGTCTGCGCAACGGCGATACCAGCAATGGCGTCGAGTTGGGCAAAAGTAACATTCTGCTGa                             <  4:1047284/148‑1 (MQ=255)
                                                                                                              aaaTTCGCAACCACCTGGACGATTACGTTATCGGCCAGGAACAGGCGAAAAAAGTGCt                                                                                                                      <  4:174724/58‑1 (MQ=255)
                                                                                                              aaaTTCGCAACCACCTGGACGATTACGTTATCGGCCAGGAACAGGCGAAAAAAGTGCt                                                                                                                      >  3:174724/1‑58 (MQ=255)
                                                                                                              aaaTTCGCAACCACCTGGACGATTACGTTATCGGCCAGGAACAGGCGAAAAAAGTGCTGGCGGTCGCGGTATACAACCATTACAAACGTCTGCGCAACGGCGATACCa                                                                    >  1:1036291/1‑108 (MQ=255)
                                                                                                              aaaTTCGCAACCACCTGGACGATTACGTTATCGGCCAGGAACAGGCGAAAAAAGTGCTGGCGGTCGCGGTATACAACCATTACAAACGTCTGCGCAACGGCGATACCa                                                                    <  2:1036291/108‑1 (MQ=255)
                                                                                                              aaaTTCGCAACCACCTGGACGATTACGTTATCGGCCAGGAACAGGCGAAAAAAGTGCTGGCGGTCGCGGTATACAACCATTACAAACGTCTGCGCAACGGCGATACCAGCaa                                                                <  1:1094888/112‑1 (MQ=255)
                                                                                                              aaaTTCGCAACCACCTGGACGATTACGTTATCGGCCAGGAACAGGCGAAAAAAGTGCTGGCGGTCGCGGTATACAACCATTACAAACGTCTGCGCAACGGCGATACCAGCaa                                                                >  2:1094888/1‑112 (MQ=255)
                                                                                                                aTTCGCAACCACCTGGACGATTACGTTATCGGCCAGGAACAGGCGAAAAAAGTGCTGGCGGTCGCGGTATACAACCATTACAAACGTCTGCGCAACGGCGATACCAGCa                                                                 <  3:1134311/109‑1 (MQ=255)
                                                                                                                aTTCGCAACCACCTGGACGATTACGTTATCGGCCAGGAACAGGCGAAAAAAGTGCTGGCGGTCGCGGTATACAACCATTACAAACGTCTGCGCAACGGCGATACCAGCa                                                                 >  4:1134311/1‑109 (MQ=255)
                                                                                                                aTTCGCAACCACCTGGACGATTACGTTATCGGCCAGGAACAGGCGAAAAAAGTGCTGGCGGTCGCGGTATACAACCATTACAAACGTCTGCGCAACGGCGATACCAGCAATGGCGTCGAGt                                                     <  2:270835/121‑1 (MQ=255)
                                                                                                                aTTCGCAACCACCTGGACGATTACGTTATCGGCCAGGAACAGGCGAAAAAAGTGCTGGCGGTCGCGGTATACAACCATTACAAACGTCTGCGCAACGGCGATACCAGCAATGGCGTCGAGt                                                     >  1:270835/1‑121 (MQ=255)
                                                                                                                         caccTGGACGATTACGTTATCGGCCAGGAACAGGCGAAAAAAGTGCTGGCGGTCGCGGTATACAACCATTACAAACGTCTGCGCAACGGCGATACCAGCAATGGCGTCGAGTTGGGCAAAAGTAACATTCTGCTGATCGGTCCGACCg                 <  4:181338/148‑1 (MQ=255)
                                                                                                                           ccTGGACGATTACGTTATCGGCCAGGAACAGGCGAAAAAAGTGCTGGCGGTCGCGGTATACAACCATTACAAACGTCTGCGCAACgg                                                                            >  2:549501/1‑87 (MQ=255)
                                                                                                                           ccTGGACGATTACGTTATCGGCCAGGAACAGGCGAAAAAAGTGCTGGCGGTCGCGGTATACAACCATTACAAACGTCTGCGCAACgg                                                                            <  1:549501/87‑1 (MQ=255)
                                                                                                                            cTGGACGATTACGTTATCGGCCAGGAACAGGCGAAAAAAGTGCTGGCGGTCGCGGTATACAACCATTACAAACGTCTGCGCAACGGCGATACCAGCAATGGCGTCGAGTTGGGCAAAAGTAACATTCTGCTGATCGGTCCGACCGGtt              <  2:563140/148‑1 (MQ=255)
                                                                                                                              ggACGATTACGTTATCGGCCAGGAACAGGCGAAAAAAGTGCTGGCGGTCGCGGTATACAACCATTACAAACGTCTGCGCAACGGCGAt                                                                        >  4:1022401/1‑88 (MQ=255)
                                                                                                                              ggACGATTACGTTATCGGCCAGGAACAGGCGAAAAAAGTGCTGGCGGTCGCGGTATACAACCATTACAAACGTCTGCGCAACGGCGAt                                                                        <  3:1022401/88‑1 (MQ=255)
                                                                                                                              ggACGATTACGTTATCGGCCAGGAACAGGCGAAAAAAGTGCTGGCGGTCGCGGTATACAACCATTACAAACGTCTGCGCAACGGCGATa                                                                       >  2:999878/1‑89 (MQ=255)
                                                                                                                              ggACGATTACGTTATCGGCCAGGAACAGGCGAAAAAAGTGCTGGCGGTCGCGGTATACAACCATTACAAACGTCTGCGCAACGGCGATa                                                                       <  1:999878/89‑1 (MQ=255)
                                                                                                                              ggACGATTACGTTATCGGCCAGGAACAGGCGAAAAAAGTGCTGGCGGTCGCGGTATACAACCATTACAAACGTCTGCGCAACGGCGATACCAGCAATGGCGTCGAGTTGGGCAAAAGTAACATTCTGCTGATCGGTCCGACCGGTTcc            >  4:934236/1‑148 (MQ=255)
                                                                                                                                aCGATTACGTTATCGGCCAGGAACAGGCGAAAAAAGTGCTGGCGGTCGCGGTATACAACCATTACAAACGTCTGCGCAACGGCGATACCAGCAATGGCGTGGAGTTGGGCaaa                                             <  4:735065/113‑1 (MQ=255)
                                                                                                                                aCGATTACGTAATCGGCCAGGAACAGGCGAAAAAAGTGCTGGCGGTCGCGGTATACAACCATTACAAACGTCTGCGCAACGGCGATACCAGCAATGGCGTCGAGTTGGGCaaa                                             >  3:735065/1‑113 (MQ=255)
                                                                                                                                   aTTACGTTATCGGCCAGGAACAGGCGAAAAAAGTGCTGGCGGTCGCGGTATACAACCATTACAAACGTCTGCGCAACGGCGATAcc                                                                     >  3:926214/1‑86 (MQ=255)
                                                                                                                                   aTTACGTTATCGGCCAGGAACAGGCGAAAAAAGTGCTGGCGGTCGCGGTATACAACCATTACAAACGTCTGCGCAACGGCGATAcc                                                                     <  4:926214/86‑1 (MQ=255)
                                                                                                                                    ttACGTTATCGGCCAGGAACAGGCGAAAAAAGTGCTGGCGGTCGCGGTATACAACCATTACAAACGTCTGCGCAACGGCGATACCAGCAATGGCGTCGAGTTGGGCAAAAGTAACATTCTGCTGATc                           >  1:279250/1‑127 (MQ=255)
                                                                                                                                    ttACGTTATCGGCCAGGAACAGGCGAAAAAAGTGCTGGCGGTCGCGGTATACAACCATTACAAACGTCTGCGCAACGGCGATACCAGCAATGGCGTCGAGTTGGGCAAAAGTAACATTCTGCTGATc                           <  2:279250/127‑1 (MQ=255)
                                                                                                                                     tACGTTATCGGCCAGGAACAGGCGAAAAAAGTGCTGGCGGTCGCGGTATACAACCATTACAAACGTCTGCGCAACGGCGATACCAGCAATGGCGTCGAGTTGGGCaaa                                             <  4:487684/108‑1 (MQ=255)
                                                                                                                                     tACGTTATCGGCCAGGAACAGGCGAAAAAAGTGCTGGCGGTCGCGGTATACAACCATTACAAACGTCTGCGCAACGGCGATACCAGCAATGGCGTCGAGTTGGGCaaa                                             >  3:487684/1‑108 (MQ=255)
                                                                                                                                        gTTATCGGCCAGGAACAGGCGAAAAAAGTGCTGGCGGTCGCGGTATACAACCATTACAAACGTCTGCGCAACGGCGATACCAGCAATGGCGTCGAGTTGGGCAAAAGTAACATTCTGCTGATCGGTCCGACCGGTTCCGGTAAAACgc  <  2:726815/148‑1 (MQ=255)
                                                                                                                                           aTCGGCCAGGAACAGGCGAAAAAAGTGCTGGCGGTCGCGGTATACAACCATTACAAACGTCTGCGCAACGGCGATACCAGCAATGGCGTCGAGTTGGGCaaa                                             <  4:296343/102‑1 (MQ=255)
                                                                                                                                           aTCGGCCAGGAACAGGCGAAAAAAGTGCTGGCGGTCGCGGTATACAACCATTACAAACGTCTGCGCAACGGCGATACCAGCAATGGCGTCGAGTTGGGCaaa                                             >  3:296343/1‑102 (MQ=255)
                                                                                                                                               |                                                                                                                                            
ATCCGTGTATATCTGCGACGAATGTGTTGATTTATGTAACGACATCATTCGCGAAGAGATTAAAGAAGTTGCACCGCATCGTGAACGCAGTGCGCTACCGACGCCGCATGAAATTCGCAACCACCTGGACGATTACGTTATCGGCCAGGAACAGGCGAAAAAAGTGCTGGCGGTCGCGGTATACAACCATTACAAACGTCTGCGCAACGGCGATACCAGCAATGGCGTCGAGTTGGGCAAAAGTAACATTCTGCTGATCGGTCCGACCGGTTCCGGTAAAACGC  >  NC_000913/457521‑457804

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: