Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,521,346 A→T 8.7% intergenic (+18/+247) 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,3460AT 8.7% 132.3 / 7.1 60intergenic (+18/+247)lysV/xapRtRNA‑Lys/DNA‑binding transcriptional activator XapR
Reads supporting (aligned to +/- strand):  ref base A (23/30);  new base T (2/3);  total (27/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.58e-01

TTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGT  >  NC_000913/2521205‑2521488
                                                                                                                                             |                                                                                                                                              
ttACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACGCACCAACGTAAAAAAGCGCCCTaa                                                                                                                                               >  2:122373/1‑143 (MQ=255)
  accaggaggGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAg                                                                                                                                             <  1:406229/140‑1 (MQ=255)
         gggggTCGGCGGTTCGATCCCGTCATCACCCACCACTTCGGGTCGTTAGCACAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCTTTAGGGCGCttt                                                                                                                                      <  1:36910/143‑1 (MQ=17)
           gggTCGGCGGTTCGATCCCGTCATCACCCACCACTTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCttttt                                                                                                                                    <  1:59705/143‑1 (MQ=255)
           gggCCGGCGGTTCGAGCCCGTCATCACCCACCACTTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCttttt                                                                                                                                    <  1:590310/143‑1 (MQ=255)
              tCGGCGGTTCGATCCCGTCATCACCCACCACTTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCTTTAGGGCGCtttt                                                                                                                                     <  1:465796/139‑1 (MQ=255)
              tCGGCGGTTCGATCCCGTCATCACCCACCACTTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCTTTAGGGCGCtttt                                                                                                                                     >  2:465796/1‑139 (MQ=255)
              tCGGCGGTTCGATCCCGTCATCACCCACCACTTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAACGGCGCTTTTTTag                                                                                                                                 >  1:539540/1‑142 (MQ=255)
              tCGGCGGTTCGATCCCGTCATCACCCACCACTTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTAc                                                                                                                                 <  1:364751/143‑1 (MQ=255)
              tCGGCGGTTCGATCCCGTCATCACCCACCACTTCGGGTCGTTAGCTCAGGTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCGTTTGTAc                                                                                                                                 >  2:711452/1‑143 (MQ=255)
                      tCGATCCCGTCATCACCCACCACTTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCg                                                                                                                         <  1:338734/143‑1 (MQ=255)
                         aTCCCGTCATCACCCACCACTTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATa                                                                                                                      >  2:903377/1‑143 (MQ=255)
                         aTCCCGTCATCACCCACCACTTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATa                                                                                                                      >  2:757944/1‑143 (MQ=255)
                           cccGTCATCACCCACCACTTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTCACTTTTAATCCATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACt                                                                                                                    <  2:393833/143‑1 (MQ=255)
                             cGTCATCACCCACCACTTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCa                                                                                                                  >  1:1089194/1‑143 (MQ=255)
                              gTCATCACCCACCACTTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGCCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCac                                                                                                                 >  2:1129859/1‑142 (MQ=255)
                               tcatcaCCCACCACTTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGCATCCTGCACGACCTACTAATGAAAAAAAACGCCCTAAAGGCGCTTTTTTACTATCTACGATACTCaaa                                                                                                                >  1:882796/1‑143 (MQ=14)
                                catcaCCCACCACTTCGGGTCGTTAGCTCAGTGGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAGTCGTGCACGGCCCACTAATGTGAAAAAGCGCCCTGAAGGCGCTTTTGTACTATCTGCGAGACTCGAAg                                                                                                               >  2:190852/1‑143 (MQ=14)
                                  tcaCCCACCACTTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTGGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGACACTCaacgct                                                                                                             >  2:186739/1‑139 (MQ=255)
                                   cacccagcAATTCGGGTCTGTAGCTCAGGTGGTAGAGCAGTTGACTTTCAATCAATTGGTCGCAGGTTCGAATCCTGCAAGACACAACAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGAtt                                                                                                            <  1:1035884/143‑1 (MQ=14)
                                   cacccaccACTTCGGGTCGTTAGCTCAGTTTGTAGAGCAGTTGATTTTTTATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTACAAAAGCGCCCTAATGGCGCTTTTTTACTATCTGCGATACTCAAAGAtt                                                                                                            <  2:738479/143‑1 (MQ=14)
                                   cacccaccACTTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGAtt                                                                                                            <  2:1093188/143‑1 (MQ=255)
                                       caccgcTGCGGGTCGTTAGCTCAGTTGGTAGAGCCGTTGACTTTTAATCAAGTGGTGTTAAGGTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGaa                                                                                                        <  1:1125375/143‑1 (MQ=11)
                                       caccacTTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGaa                                                                                                        <  1:243217/143‑1 (MQ=40)
                                        accacTTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAAc                                                                                                       >  2:857444/1‑143 (MQ=39)
                                                 ggTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGcagcag                                                                                              <  1:797606/143‑1 (MQ=32)
                                                  gttgTTAGCTCATTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAgg                                                                                             <  1:1106004/140‑1 (MQ=17)
                                                   tCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTAGAACCTGCAGCAggg                                                                                            >  1:78313/1‑142 (MQ=17)
                                                       tAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTttga                                                                                        <  1:124804/143‑1 (MQ=32)
                                                       tAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTtcga                                                                                        >  2:1088013/1‑143 (MQ=17)
                                                         gCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTttgagt                                                                                      <  1:373219/143‑1 (MQ=32)
                                                          cTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTttgagtt                                                                                     <  2:383533/143‑1 (MQ=32)
                                                            cAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGATCCTGCAGCAGGTttgagttga                                                                                   >  1:947896/1‑143 (MQ=17)
                                                            cAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTttgagttga                                                                                   >  2:485226/1‑143 (MQ=32)
                                                                 ggTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCa                                                                              <  1:546305/143‑1 (MQ=32)
                                                                      agCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAACCAGGCACTAGAGTGGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGaa                                                                         <  1:261796/143‑1 (MQ=11)
                                                                        cagtgttcTTTTAATCAATTGGTCGCAGGTTCGAATCGTGCACTACCCACCAATGTAAAAAAGGGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAaac                                                                       <  2:704742/136‑1 (MQ=11)
                                                                          gTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCATCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGACCGCAGCGAaacaa                                                                     >  1:1081614/1‑143 (MQ=17)
                                                                            tGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACg                                                                   <  2:32653/143‑1 (MQ=32)
                                                                             gACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACgg                                                                  <  2:1089194/143‑1 (MQ=32)
                                                                              aCTTTCATTCAATTGGTAGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGa                                                                 <  1:11177/143‑1 (MQ=12)
                                                                               cTTTTAATCAATGGGACGCAGGGTCGAATCCTTCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGTAGGTTTGAGTTGAGCGCAGCGAAACAACGGAg                                                                <  1:331983/143‑1 (MQ=11)
                                                                                ttttAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGCTTTGAGCTGAGCGCAGCGACACAACGGAgg                                                               >  2:768873/1‑142 (MQ=12)
                                                                                   tAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTCTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGCTTGAGTTGAGCGCAGCGAAACAACGGAGCCGc                                                            >  2:75178/1‑143 (MQ=17)
                                                                                     aTCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTc                                                          >  2:895625/1‑143 (MQ=32)
                                                                                      tCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTcg                                                         >  2:254770/1‑143 (MQ=32)
                                                                                          ttGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCgcggc                                                     >  1:1103089/1‑143 (MQ=32)
                                                                                          ggtagagCAGTTTCGAACCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCgcggc                                                     <  1:531487/137‑1 (MQ=255)
                                                                                                gCAGGTTCGACTCCTGCACGACCCACCAATGTAAAGAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATCCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGc                                               <  2:374281/143‑1 (MQ=12)
                                                                                                    gTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCCGGTTTGCGTTCAGCTCAGCGAAACAACGGAGCCCCCCGCGGCGACGGCCCgc                                           >  2:159046/1‑142 (MQ=11)
                                                                                                              tGCACGACCCACCAATGTAAAAAAGCGCCCTGCGCTGCCGCATTTACTATCTGCGATACTCAAAAATGCGAGGCGCCCGCCACTTGGAAGCCAGCGCAGCGGAACAACAGATCGGCTCGCTGCGACGGCCCGAAGGacgaggc                                 >  2:859843/1‑141 (MQ=2)
                                                                                                                caacactcaccacTCTAAAAAAGCGCCCTAAAGCCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGAAGCCGCTCGCGGCGACGGCCCGAAGGGCGagcgaa                               <  2:380445/136‑1 (MQ=11)
                                                                                                                    aCCCACCAATGTAAATAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATGCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGagcgaagcga                           <  1:240639/143‑1 (MQ=17)
                                                                                                                    aCCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGagcgaagcga                           <  1:895625/143‑1 (MQ=32)
                                                                                                                       caccaATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCAGATCGCGGGGACGGCCCGAAGGGCGAGCGAAGCGAGTc                        >  1:571281/1‑143 (MQ=12)
                                                                                                                        accaATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCa                       >  1:977371/1‑143 (MQ=32)
                                                                                                                                 aaaaaGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAGGCGAGTCATCCTGCACg              >  2:141214/1‑143 (MQ=17)
                                                                                                                                  aaaaGCGCCTTTAGGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTttgagttga                                                                                   <  2:248536/73‑1 (MQ=11)
                                                                                                                                  aaaaGCGCCTTTAGGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTttgagttga                                                                                   >  1:248536/1‑73 (MQ=11)
                                                                                                                                         cccGAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACccacca      <  1:803849/143‑1 (MQ=17)
                                                                                                                                           cTAACGGCGCTCTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAAt    <  1:916587/143‑1 (MQ=17)
                                                                                                                                             aaaGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGCCATCCTGCACGACCCACCAATGt  >  1:1061876/1‑143 (MQ=17)
                                                                                                                                             |                                                                                                                                              
TTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGT  >  NC_000913/2521205‑2521488

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