Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,521,344 C→T 7.3% intergenic (+16/+249) lysV → / ← xapR tRNA‑Lys/DNA‑binding transcriptional activator XapR

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,521,3440CT 7.3% 122.3 / 3.3 56intergenic (+16/+249)lysV/xapRtRNA‑Lys/DNA‑binding transcriptional activator XapR
Reads supporting (aligned to +/- strand):  ref base C (23/28);  new base T (2/2);  total (26/30)
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.57e-01

CCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAAT  >  NC_000913/2521202‑2521486
                                                                                                                                              |                                                                                                                                               
cccTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGTTCGCATGTTCGAATCCTGCACGACCCACCAATGTCAACACGCGccc                                                                                                                                                 >  2:175466/1‑143 (MQ=14)
        agggggGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGgcg                                                                                                                                         <  2:675995/143‑1 (MQ=255)
              gggTCGGCGGTTCGATCCCGTCATCACCCACCACTTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCTTTAGGGCGCttttt                                                                                                                                   <  1:9339/143‑1 (MQ=255)
              gggTCGGCGGTTCGATCCCGTCATCACCCACCACTGCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTGTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCttttt                                                                                                                                   <  2:599336/143‑1 (MQ=255)
                gTCGGCGGTTCGATCCCGTCATCACCCACCACTTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTa                                                                                                                                 <  1:423827/143‑1 (MQ=255)
                  cggcggTTCGATCCCGTCATCACCCACCACTTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTAGCAGGGTCGAGTCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACt                                                                                                                               >  2:791665/1‑143 (MQ=255)
                    gcggTTCGATCCCGTCATCACCCACCACTTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCTTTAGGGCGCTTTTTTAcatt                                                                                                                             <  2:53115/143‑4 (MQ=255)
                    gcggTTCGATCCCGTCATCACCCACCACTTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAACGGCGCTTTTTTACTAt                                                                                                                             >  2:379867/1‑143 (MQ=255)
                        ttCGATCCCGTCATCACCCACCACTTCGGGTAGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTga                                                                                                                         >  1:832045/1‑142 (MQ=255)
                          cGATCCCGTCATCACCCACCACTTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGa                                                                                                                       <  2:386059/143‑1 (MQ=255)
                           gATCCCGTCATCACCAACAACTTCGGGACGTTAGCTCAGTTGGTAGAACAGTTGACTTTTAATCAATTGGTAGCAGGTTCGAATCCTAAACGACCCACCAATGTAAAAAAGCGCCATAAAGGAGCTTTTTTACTATCTGcact                                                                                                                      >  2:235073/1‑140 (MQ=14)
                                  tcatcaCCCACCACTTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCaaa                                                                                                               <  1:428311/143‑1 (MQ=255)
                                     tcaCCCACCACTTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGAGCTTTTTTACTATCTGCGATACTCAAAGAt                                                                                                            <  2:555065/143‑1 (MQ=255)
                                       acccaccACTTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTc                                                                                                          <  2:549192/143‑1 (MQ=255)
                                         ccacccCTTCGGGTCGGTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGa                                                                                                        <  2:160937/143‑1 (MQ=255)
                                           accacTTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTCCTATCTGCGATACTCAACGATTCGCAc                                                                                                      >  2:879067/1‑143 (MQ=255)
                                              acTTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTg                                                                                                   >  2:898024/1‑143 (MQ=35)
                                                 tCGGGTCGTTAGCTGAGTTGGCTGAGCGGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTgcag                                                                                                <  2:590742/143‑1 (MQ=255)
                                                 tCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTgcag                                                                                                <  1:726013/143‑1 (MQ=33)
                                                     gTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTTCAGCAgg                                                                                            >  2:35285/1‑143 (MQ=17)
                                                        gTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGttt                                                                                         >  2:661109/1‑143 (MQ=32)
                                                          tAGCTCTGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTCTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCTGCAGGTtcga                                                                                       <  1:602377/143‑1 (MQ=255)
                                                          tAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTttga                                                                                       <  1:680432/143‑1 (MQ=32)
                                                          tAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTttga                                                                                       <  1:325546/143‑1 (MQ=32)
                                                               cAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTttgagttga                                                                                  <  1:311583/143‑1 (MQ=32)
                                                                aggtGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTttgagttgag                                                                                 <  2:700988/140‑1 (MQ=17)
                                                                  ttGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAgcg                                                                               >  2:156006/1‑143 (MQ=32)
                                                                    ggTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCa                                                                             <  1:545484/143‑1 (MQ=32)
                                                                       agagCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCg                                                                          <  1:789259/143‑1 (MQ=32)
                                                                         agCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGaa                                                                        <  2:897430/143‑1 (MQ=32)
                                                                           cAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGATACAGGTGAGCGCAGCGAaac                                                                      >  2:710872/1‑143 (MQ=255)
                                                                             gTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAGAAGAGCGCACGAAAGGAGCTCTTTTACTCTCTGCGATACTCAACGGTTCGAACCTGCAGCAGGTTCGGGTTGAGCGGCGCGaaccga                                                                    >  2:315879/1‑139 (MQ=11)
                                                                             gTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGcaacaa                                                                    >  1:263473/1‑143 (MQ=17)
                                                                              ttGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAaacaac                                                                   >  2:47401/1‑143 (MQ=32)
                                                                                        aTCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGCTTTGAGTTGAGCGCAGCGCACCCACGGAGCCGCTc                                                         >  2:422231/1‑143 (MQ=255)
                                                                                          cAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTcgc                                                       <  1:379867/143‑1 (MQ=32)
                                                                                           aaTTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTcgcg                                                      <  1:47401/143‑1 (MQ=32)
                                                                                            aTTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTCAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGGGCAGCGAAACAACGGAGCCGCTCgcgg                                                     >  1:176420/1‑143 (MQ=17)
                                                                                            aTTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGCAACAACGGAGCCGCTCgcgg                                                     >  2:501880/1‑143 (MQ=17)
                                                                                                gTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAGAGAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGCGCTGCGCGCAACGACACAACTGTGACCCTCGCGGCGAc                                                 >  2:63040/1‑143 (MQ=11)
                                                                                                 tCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGATCGCGGCGACg                                                >  1:581512/1‑143 (MQ=17)
                                                                                                   gCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGc                                              >  2:127428/1‑143 (MQ=32)
                                                                                                         tCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCACAGATTCGAATGTGTCGCGAGTTTGAGTTGAGCTCAGCGAACCAACAGAGCCGCTCGCGGCGACGGGCCGacg                                        >  2:792598/1‑141 (MQ=11)
                                                                                                              tCCTGTGCGACCCACCATGGTCACGGAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGgcga                                  <  2:836174/143‑1 (MQ=11)
                                                                                                              tCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGgcga                                  <  2:135602/143‑1 (MQ=32)
                                                                                                                   cACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCGAAGATTCGAGCGTGGAGCCGGTTTGAGGAGAGCGCAGCGAGACAACGGAGCCGCTCGCGACGACGGCGCGAAGGgcgagcgca                             >  2:305937/1‑141 (MQ=11)
                                                                                                                   cACGACCCACCAACGTAAGAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGAGTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGagcgaa                             <  2:773957/143‑1 (MQ=12)
                                                                                                                       ccccacccaTGTAACACAGCGCCCTAAGGGCGCTTTTTTACTATCTGCGATACGCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGagcgaagcga                         <  2:760412/142‑1 (MQ=11)
                                                                                                                       aCCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTGTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGagcgaagcga                         <  1:170340/143‑1 (MQ=17)
                                                                                                                                  tAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAAAAACGGAGCCGCTCGCGGCGACAGCCCGAAGGGCGAGAGAAGCGAGTCATCCTGCa              >  2:542825/1‑143 (MQ=12)
                                                                                                                                     aaaaGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCTATGATTCAAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGa           <  1:661109/143‑1 (MQ=12)
                                                                                                                                      acaGCGCCCTAAACGCGCGTTTTTACAAGCTGCGATACTCACAGATTCGAACCTGCAGCAGGTTTGAGTAGAGCTCAGCGAAACAACGTAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGAc          <  1:169923/141‑1 (MQ=11)
                                                                                                                                      aaaGCGCCTTTAGGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGAc          >  1:146977/1‑143 (MQ=255)
                                                                                                                                      aaaGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCCCGAc          >  2:277975/1‑143 (MQ=17)
                                                                                                                                       aaGCGCCTTTAGGGCGCTTTTTTACTATCTGCGATACCCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGCGCGCAGCGCAACAACGGAGCCGCTCGCGGCGACGGCCCGCAGGGCGAGCGAAGCGAGTCATCCTGCACGAcc         >  1:366393/1‑143 (MQ=11)
                                                                                                                                       aaGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACGCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGAcc         <  2:515637/143‑1 (MQ=17)
                                                                                                                                           gcCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACccacc     >  2:471397/1‑143 (MQ=32)
                                                                                                                                            cccTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACccacca    >  2:766926/1‑143 (MQ=32)
                                                                                                                                             ccTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCaa   <  1:898024/143‑1 (MQ=32)
                                                                                                                                             ccTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCAGCTCGCGGCGACGGCCCGACGGACGCCGCGAGCGAGTCATACGGCACGACccaccac   >  2:576566/1‑142 (MQ=11)
                                                                                                                                              cTAAGGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGCGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGTCGACGGCCCGGAGGGAGATCGAAGCGAGTCATCCTGCACGACccaccgag  >  2:884871/1‑140 (MQ=11)
                                                                                                                                              |                                                                                                                                               
CCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAAT  >  NC_000913/2521202‑2521486

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: