Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 457,654 T→A 8.1% Y77N (TAC→AAC)  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,6540TA 8.1% 257.1 / 11.7 74Y77N (TAC→AAC) clpXATPase and specificity subunit of ClpX‑ClpP ATP‑dependent serine protease
Reads supporting (aligned to +/- strand):  ref base T (33/35);  new base A (3/3);  total (36/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 = 7.09e-01

CTGATTGCCGGTCCATCCGTGTATATCTGCGACGAATGTGTTGATTTATGTAACGACATCATTCGCGAAGAGATTAAAGAAGTTGCACCGCATCGTGAACGCAGTGCGCTACCGACGCCGCATGAAATTCGCAACCACCTGGACGATTACGTTATCGGCCAGGAACAGGCGAAAAAAGTGCTGGCGGTCGCGGTATACAACCATTACAAACGTCTGCGCAACGGCGATACCAGCAATGGCGTCGAGTTGGGCAAAAGTAACATTCTGCTGATCGGTCCGACCGGTTCC  >  NC_000913/457507‑457794
                                                                                                                                                   |                                                                                                                                            
cTGATTGCCGGTCCATCCGTGTATATCTGCGACGAATGTGTTGATTTATGTAACGACATCATTCGCGAAGAGATTAAAGAAGTTGCACCGCATCGTGAACGCAGTGCGCTACCGACGCCGCATGAAATTCGCAACCACCTGGACGAtt                                                                                                                                              <  4:413710/148‑1 (MQ=255)
       ccGGTCCATCCGTGTATATCTGCGACGAATGTGTTGATTTATGTAACGACATCATTCGCGAAGAGATTAAAGAAGTTGCACCGCATCGTGAACGCAGTGCGCTACCGACGCCGCATGAAATTCGCAACCACCTGGACGATTACGTTAt                                                                                                                                       <  4:66490/148‑1 (MQ=255)
         ggTCCATCCGTGTATATCTGCGACGAATGTGTTGATTTATGTAACGACATCATTCGCGAAGAGATTAAAGAAGTTGCACCGCATCGTGAACGCAGTGCGCTACCGACGCCGCATGAAATTCGCAACCACCTGGACGATTACGTTATCg                                                                                                                                     >  4:626162/1‑148 (MQ=255)
              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                                                                                                                  <  1:933043/135‑1 (MQ=255)
                                         tGATTTATGTAACGACATCATTCGCGAAGAGATTAAAGAAGTTGCACCGCATCGTGAACGCAGTGCGCTACCGACGCCGCATGAAATTCGCAACCACCTGGACGATTACGTTATCGGCCAGGAACAGGCGaaaaa                                                                                                                  >  2:933043/1‑135 (MQ=255)
                                           aTTTATGTAACGACATCATTCGCGAAGAGATTAAAGAAGTTGCACCGCATCGTGAACGCAGTGCGCTACCGACGCCGCATGAAATTCGCAACCACCTGGACGATTACGTTATCGGCCAGGAACAGGCGAAAAAAGTGCTGGCGGTcgc                                                                                                   <  2:652704/148‑1 (MQ=255)
                                            tttATGTAACGACATCATTCGCGAAGAGATTAAAGAAGTTGCACCGCATCGTGAACGCAGTGCGCTACCGACGCCGCATGAAATTCGCAACCACCTGGACGATTACGTTATCGGCCAGGAACAg                                                                                                                          <  2:881717/124‑1 (MQ=255)
                                            tttATGTAACGACATCATTCGCGAAGAGATTAAAGAAGTTGCACCGCATCGTGAACGCAGTGCGCTACCGACGCCGCATGAAATTCGCAACCACCTGGACGATTACGTTATCGGCCAGGAACAg                                                                                                                          >  1:881717/1‑124 (MQ=255)
                                                 gTAACGACATCATTCGCGAAGAGATTAAAGAAGTTGCACCGCATCGTGAACGCAGTGCGCTACCGACGCCGCATGAAATTCGCAACCACCTGGACGATTACGt                                                                                                                                          >  2:929323/1‑103 (MQ=255)
                                                 gTAACGACATCATTCGCGAAGAGATTAAAGAAGTTGCACCGCATCGTGAACGCAGTGCGCTACCGACGCCGCATGAAATTCGCAACCACCTGGACGATTACGt                                                                                                                                          <  1:929323/103‑1 (MQ=255)
                                                   aaCGACATCATTCGCGAAGAGATTAAAGAAGTTGCACCGCATCGTGAACGCAGTGCGCTACCGACGCCGCATGAAATTCGCAACCACCTGGCCGATAACGTAATCgtc                                                                                                                                   >  3:1033026/1‑106 (MQ=255)
                                                   aaCGACATCATTCGCGAAGAGATTAAAGAAGTTGCACCGCATCGTGAACGCAGTGCGCTACCGACGCCGCATGAAATTCGCAACCACCTGGCCGATAACGTAATCgtc                                                                                                                                   <  4:1033026/108‑3 (MQ=255)
                                                              tCGCGAAGAGATTAAAGAAGTTGCACCGCATCGTGAACGCAGTGCGCTACCGACGCCGCATGAAATTCGCAACCACCTGGACGATTACGTTATCGGcc                                                                                                                                  >  2:708475/1‑98 (MQ=255)
                                                              tCGCGAAGAGATTAAAGAAGTTGCACCGCATCGTGAACGCAGTGCGCTACCGACGCCGCATGAAATTCGCAACCACCTGGACGATTACGTTATCGGcc                                                                                                                                  <  1:708475/98‑1 (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                                                                           >  4:571818/1‑138 (MQ=255)
                                                                             agaagTTGCACCGCATCGTGAACGCAGTGCGCTACCGACGCCGCATGAAATTCGCAACCACCTGGACGATTACGTTATCGGCCAGGAACAGGCGAAAAAAGTGCTGGCGGTCGCGGTATACAACCATTACAAACGTCt                                                                           <  3:571818/138‑1 (MQ=255)
                                                                              gaagTTGCACCGCATCGTGAACGCAGTGCGCTACCGACGCCGCATGAAATTCGCAACCACCTGGACGATTACGTTATCGGCCAGGAACAGGCGAAAAAAGTGCTGGCGGTcgc                                                                                                   >  1:559113/1‑113 (MQ=255)
                                                                              gaagTTGCACCGCATCGTGAACGCAGTGCGCTACCGACGCCGCATGAAATTCGCAACCACCTGGACGATTACGTTATCGGCCAGGAACAGGCGAAAAAAGTGCTGGCGGTcgc                                                                                                   <  2:559113/113‑1 (MQ=255)
                                                                              gaagTTGCACCGCATCGTGAACGCAGTGCGCTACCGACGCCGCATGAAATTCGCAACCACCTGGACGATTACGTTATCGGCCAAGAACAGGCGAAAAAAGTGCTGGCGGTcgc                                                                                                   >  3:340003/1‑113 (MQ=255)
                                                                              gaagTTGCACCGCATCGTGAACGCAGTGCGCTACCGACGCCGCATGAAATTCGCAACCACCTGGACGATTACGTTATCGGCCAAGAACAGGCGAAAAAAGTGCTGGCGGTcgc                                                                                                   <  4:340003/113‑1 (MQ=255)
                                                                               aagTTGCACCGCATCGTGAACGCAGTGCGCTACCGACGCCGCATGAAATTCGCAACCACCTGGACGATTACGTTATCGGCCAGGAACAGGCGAAAAAAGTGCTGGCGGTCGCGGTATACAACCATTACaa                                                                                 >  4:20691/1‑130 (MQ=255)
                                                                               aagTTGCACCGCATCGTGAACGCAGTGCGCTACCGACGCCGCATGAAATTCGCAACCACCTGGACGATTACGTTATCGGCCAGGAACAGGCGAAAAAAGTGCTGGCGGTCGCGGTATACAACCATTACaa                                                                                 <  3:20691/130‑1 (MQ=255)
                                                                                  ttGCACCGCATCGTGAACGCAGTGCGCTACCGACGCCGCATGAAATTCGCAACCACCTGGACGATTACGTTATCGGCCAGGAACAgg                                                                                                                         >  4:30473/1‑87 (MQ=255)
                                                                                  ttGCACCGCATCGTGAACGCAGTGCGCTACCGACGCCGCATGAAATTCGCAACCACCTGGACGATTACGTTATCGGCCAGGAACAgg                                                                                                                         <  3:30473/87‑1 (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                                                          <  2:1036291/108‑1 (MQ=255)
                                                                                                                            aaaTTCGCAACCACCTGGACGATTACGTTATCGGCCAGGAACAGGCGAAAAAAGTGCTGGCGGTCGCGGTATACAACCATTACAAACGTCTGCGCAACGGCGATACCa                                                          >  1:1036291/1‑108 (MQ=255)
                                                                                                                            aaaTTCGCAACCACCTGGACGATTACGTTATCGGCCAGGAACAGGCGAAAAAAGTGCTGGCGGTCGCGGTATACAACCATTACAAACGTCTGCGCAACGGCGATACCAGCaa                                                      >  2:1094888/1‑112 (MQ=255)
                                                                                                                            aaaTTCGCAACCACCTGGACGATTACGTTATCGGCCAGGAACAGGCGAAAAAAGTGCTGGCGGTCGCGGTATACAACCATTACAAACGTCTGCGCAACGGCGATACCAGCaa                                                      <  1:1094888/112‑1 (MQ=255)
                                                                                                                              aTTCGCAACCACCTGGACGATTACGTTATCGGCCAGGAACAGGCGAAAAAAGTGCTGGCGGTCGCGGTATACAACCATTACAAACGTCTGCGCAACGGCGATACCAGCa                                                       <  3:1134311/109‑1 (MQ=255)
                                                                                                                              aTTCGCAACCACCTGGACGATTACGTTATCGGCCAGGAACAGGCGAAAAAAGTGCTGGCGGTCGCGGTATACAACCATTACAAACGTCTGCGCAACGGCGATACCAGCa                                                       >  4:1134311/1‑109 (MQ=255)
                                                                                                                              aTTCGCAACCACCTGGACGATTACGTTATCGGCCAGGAACAGGCGAAAAAAGTGCTGGCGGTCGCGGTATACAACCATTACAAACGTCTGCGCAACGGCGATACCAGCAATGGCGTCGAGt                                           >  1:270835/1‑121 (MQ=255)
                                                                                                                              aTTCGCAACCACCTGGACGATTACGTTATCGGCCAGGAACAGGCGAAAAAAGTGCTGGCGGTCGCGGTATACAACCATTACAAACGTCTGCGCAACGGCGATACCAGCAATGGCGTCGAGt                                           <  2:270835/121‑1 (MQ=255)
                                                                                                                                       caccTGGACGATTACGTTATCGGCCAGGAACAGGCGAAAAAAGTGCTGGCGGTCGCGGTATACAACCATTACAAACGTCTGCGCAACGGCGATACCAGCAATGGCGTCGAGTTGGGCAAAAGTAACATTCTGCTGATCGGTCCGACCg       <  4:181338/148‑1 (MQ=255)
                                                                                                                                         ccTGGACGATTACGTTATCGGCCAGGAACAGGCGAAAAAAGTGCTGGCGGTCGCGGTATACAACCATTACAAACGTCTGCGCAACgg                                                                  <  1:549501/87‑1 (MQ=255)
                                                                                                                                         ccTGGACGATTACGTTATCGGCCAGGAACAGGCGAAAAAAGTGCTGGCGGTCGCGGTATACAACCATTACAAACGTCTGCGCAACgg                                                                  >  2:549501/1‑87 (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                 <  2:279250/127‑1 (MQ=255)
                                                                                                                                                  ttACGTTATCGGCCAGGAACAGGCGAAAAAAGTGCTGGCGGTCGCGGTATACAACCATTACAAACGTCTGCGCAACGGCGATACCAGCAATGGCGTCGAGTTGGGCAAAAGTAACATTCTGCTGATc                 >  1:279250/1‑127 (MQ=255)
                                                                                                                                                   tACGTTATCGGCCAGGAACAGGCGAAAAAAGTGCTGGCGGTCGCGGTATACAACCATTACAAACGTCTGCGCAACGGCGATACCAGCAATGGCGTCGAGTTGGGCaaa                                   <  4:487684/108‑1 (MQ=255)
                                                                                                                                                   tACGTTATCGGCCAGGAACAGGCGAAAAAAGTGCTGGCGGTCGCGGTATACAACCATTACAAACGTCTGCGCAACGGCGATACCAGCAATGGCGTCGAGTTGGGCaaa                                   >  3:487684/1‑108 (MQ=255)
                                                                                                                                                   |                                                                                                                                            
CTGATTGCCGGTCCATCCGTGTATATCTGCGACGAATGTGTTGATTTATGTAACGACATCATTCGCGAAGAGATTAAAGAAGTTGCACCGCATCGTGAACGCAGTGCGCTACCGACGCCGCATGAAATTCGCAACCACCTGGACGATTACGTTATCGGCCAGGAACAGGCGAAAAAAGTGCTGGCGGTCGCGGTATACAACCATTACAAACGTCTGCGCAACGGCGATACCAGCAATGGCGTCGAGTTGGGCAAAAGTAACATTCTGCTGATCGGTCCGACCGGTTCC  >  NC_000913/457507‑457794

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: