Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 457,659 T→A 9.0% V78V (GTT→GTA 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,6590TA 9.0% 263.1 / 12.2 77V78V (GTT→GTAclpXATPase and specificity subunit of ClpX‑ClpP ATP‑dependent serine protease
Reads supporting (aligned to +/- strand):  ref base T (33/37);  new base A (4/3);  total (37/40)
Fisher's exact test for biased strand distribution p-value = 7.05e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

CCGGTCCATCCGTGTATATCTGCGACGAATGTGTTGATTTATGTAACGACATCATTCGCGAAGAGATTAAAGAAGTTGCACCGCATCGTGAACGCAGTGCGCTACCGACGCCGCATGAAATTCGCAACCACCTGGACGATTACGTTATCGGCCAGGAACAGGCGAAAAAAGTGCTGGCGGTCGCGGTATACAACCATTACAAACGTCTGCGCAACGGCGATACCAGCAATGGCGTCGAGTTGGGCAAAAGTAACATTCTGCTGATCGGTCCGACCGGTTCCGGTAAAACGC  >  NC_000913/457514‑457804
                                                                                                                                                 |                                                                                                                                                 
ccGGTCCATCCGTGTATATCTGCGACGAATGTGTTGATTTATGTAACGACATCATTCGCGAAGAGATTAAAGAAGTTGCACCGCATCGTGAACGCAGTGCGCTACCGACGCCGCATGAAATTCGCAACCACCTGGACGATTACGTTAt                                                                                                                                                 <  4:66490/148‑1 (MQ=255)
  ggTCCATCCGTGTATATCTGCGACGAATGTGTTGATTTATGTAACGACATCATTCGCGAAGAGATTAAAGAAGTTGCACCGCATCGTGAACGCAGTGCGCTACCGACGCCGCATGAAATTCGCAACCACCTGGACGATTACGTTATCg                                                                                                                                               >  4:626162/1‑148 (MQ=255)
       aTCCGTGTATATCTGCGACGAATGTGTTGATTTATGTAACGACATCATTCGCGAAGAGATTAAAGAAGTTGCACCGCATCGTGAACGCAGTGCGCTACCGACGCCGCATGAAATTCGCAACCACCTGGCCGATAACGTAATCGTcc                                                                                                                                            >  2:158558/1‑146 (MQ=255)
       aTCCGTGTATATCTGCGACGAATGTGTTGATTTATGTAACGACATCATTCGCGAAGAGATTAAAGAAGTTGCACCGCATCGTGAACGCAGTGCGCTACCGACGCCGCATGAAATTCGCAACCACCTGGCCGATAACGTAATCGTcc                                                                                                                                            <  1:158558/146‑1 (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                                                                                                                                    <  2:881717/124‑1 (MQ=255)
                                     tttATGTAACGACATCATTCGCGAAGAGATTAAAGAAGTTGCACCGCATCGTGAACGCAGTGCGCTACCGACGCCGCATGAAATTCGCAACCACCTGGACGATTACGTTATCGGCCAGGAACAg                                                                                                                                    >  1:881717/1‑124 (MQ=255)
                                            aaCGACATCATTCGCGAAGAGATTAAAGAAGTTGCACCGCATCGTGAACGCAGTGCGCTACCGACGCCGCATGAAATTCGCAACCACCTGGCCGATAACGTAATCgtc                                                                                                                                             >  3:1033026/1‑106 (MQ=255)
                                            aaCGACATCATTCGCGAAGAGATTAAAGAAGTTGCACCGCATCGTGAACGCAGTGCGCTACCGACGCCGCATGAAATTCGCAACCACCTGGCCGATAACGTAATCgtc                                                                                                                                             <  4:1033026/108‑3 (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                                                                                     >  4:571818/1‑138 (MQ=255)
                                                                      agaagTTGCACCGCATCGTGAACGCAGTGCGCTACCGACGCCGCATGAAATTCGCAACCACCTGGACGATTACGTTATCGGCCAGGAACAGGCGAAAAAAGTGCTGGCGGTCGCGGTATACAACCATTACAAACGTCt                                                                                     <  3:571818/138‑1 (MQ=255)
                                                                       gaagTTGCACCGCATCGTGAACGCAGTGCGCTACCGACGCCGCATGAAATTCGCAACCACCTGGACGATTACGTTATCGGCCAGGAACAGGCGAAAAAAGTGCTGGCGGTcgc                                                                                                             <  2:559113/113‑1 (MQ=255)
                                                                       gaagTTGCACCGCATCGTGAACGCAGTGCGCTACCGACGCCGCATGAAATTCGCAACCACCTGGACGATTACGTTATCGGCCAGGAACAGGCGAAAAAAGTGCTGGCGGTcgc                                                                                                             >  1:559113/1‑113 (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                                                     <  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                           <  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)
                                                                                                                                               gTTATCGGCCAGGAACAGGCGAAAAAAGTGCTGGCGGTCGCGGTATACAACCATTACAAACGTCTGCGCAACGGCGATACCAGCAATGGCGTCGAGTTGGGCAAAAGTAACATTCTGCTGATCGGTCCGACCGGTTCCGGTAAAACgc  <  2:726815/148‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                 
CCGGTCCATCCGTGTATATCTGCGACGAATGTGTTGATTTATGTAACGACATCATTCGCGAAGAGATTAAAGAAGTTGCACCGCATCGTGAACGCAGTGCGCTACCGACGCCGCATGAAATTCGCAACCACCTGGACGATTACGTTATCGGCCAGGAACAGGCGAAAAAAGTGCTGGCGGTCGCGGTATACAACCATTACAAACGTCTGCGCAACGGCGATACCAGCAATGGCGTCGAGTTGGGCAAAAGTAACATTCTGCTGATCGGTCCGACCGGTTCCGGTAAAACGC  >  NC_000913/457514‑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: