Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,657,599 T→A 20.2% intergenic (+32/+387) hdeD → / ← arrS acid‑resistance membrane protein/antisense sRNA ArrS, function unknown

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,657,5990TA20.2% 227.1 / 55.6 94intergenic (+32/+387)hdeD/arrSacid‑resistance membrane protein/antisense sRNA ArrS, function unknown
Reads supporting (aligned to +/- strand):  ref base T (35/40);  new base A (7/12);  total (42/52)
Fisher's exact test for biased strand distribution p-value = 6.06e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.96e-01

TTGGTGCAACGCCGATGGTGTCTGTAACGCTGGTTTCCACACTGGTGGGAATTGAACTGATATTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTACCTGCTAGTGATTATTTCAACCTACTGAATTTCATCTAATTTTTTTCACTCTATGGCAAATTAGCCA  >  NC_000913/3657452‑3657746
                                                                                                                                                   |                                                                                                                                                   
ttGGTGCAACGCCGATGGTGTCTGTAACGCTGGTTTCCACACTGGTGGGAATTGAACTGATATTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACtt                                                                                                                                                    <  1:1069841/149‑1 (MQ=255)
ttGGTGCAACGCCGATGGTGTCTGTAACGCTGGTTTCCACACTGGTGGGAATTGAACTGATATTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACtt                                                                                                                                                    <  1:487257/149‑1 (MQ=255)
  ggTGCAACGCCGATGGTGTCTGTAACGCTGGTTTCCACACTGGTGGGAATTGAACTGATATTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGa                                                                                                                                                  <  4:927048/149‑1 (MQ=255)
   gTGCAACGCCGATGGTGTCTGTAACGCTGGTTTCCACACTGGTGGGAATTGAACTGATATTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGAt                                                                                                                                                 <  1:777245/149‑1 (MQ=255)
   gTGCAACGCCGATGGTGTCTGTAACGCTGGTTTCCACACTGGTGGGAATTGAACTGATATTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGAt                                                                                                                                                 <  2:633430/149‑1 (MQ=255)
     gCAACGCCGATGGTGTCTGTAACGCTGGTTTCCACACTGGTGGGAATTGAACTGATATTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGAtgt                                                                                                                                               <  4:1448982/149‑1 (MQ=255)
       aaCGCCGATGGTGTCTGTAACGCTGGTTTCCTCACTGGTGGGAATTGAACTGATATTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGa                                                                                                                                             <  2:747403/149‑1 (MQ=255)
       aaCGCCGATGGTGTCTGTAACGCTGGTTTCCACACTGGTGGGAATTGAACTGATATTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACATCAAGTGa                                                                                                                                             <  3:513345/149‑1 (MQ=255)
          gCCGATGGTGTCTGTAACGCTGGTTTCCACACTGGTGGGAATTGAACTGATATTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTa                                                                                                                                          <  1:612633/149‑1 (MQ=255)
          gCCGATGGTGTCTGTAACGCTGGTTTCCACACTGGTGGGAATTGAACTGATATTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTa                                                                                                                                          <  2:573904/149‑1 (MQ=255)
          gCCGATGGTGTCTGTAACGCTGGTTTCCACACTGGTGGGAATTGAACTGATATTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACATCAAGTGACTa                                                                                                                                          <  1:1245785/149‑1 (MQ=255)
          gCCGATGGTGTCTGTAACGCTGGTTTCCACACTGGTGGGAATTGAACTGATATTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACATCAAGTGACTa                                                                                                                                          <  3:140558/149‑1 (MQ=255)
           ccGATGGTGTCTGTAACGCTGGTTTCCACACTGGTGGGAATTGAACTGATATTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACATCAAGTGACTAt                                                                                                                                         <  1:1106115/149‑1 (MQ=255)
            cGATGGTGTCTGTAACGCTGGTTTCCACACTGGTGGGAATTGAACTGATATTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACATCAAGTGACTATg                                                                                                                                        <  3:1213835/149‑1 (MQ=255)
                      tGTAACGCTGGTTTCCACACTGGTGGGAATTGAACTGATATTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGa                                                                                                                                             <  3:1338870/134‑1 (MQ=255)
                      tGTAACGCTGGTTTCCACACTGGTGGGAATTGAACTGATATTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGa                                                                                                                                             >  4:1338870/1‑134 (MQ=255)
                      tGTAACGCTGGTTTCCACACTGGTGGGAATTGAACTGATATTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACATCAAGTGACTATGa                                                                                                                                       <  2:446130/140‑1 (MQ=255)
                      tGTAACGCTGGTTTCCACACTGGTGGGAATTGAACTGATATTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACATCAAGTGACTATGa                                                                                                                                       >  1:446130/1‑140 (MQ=255)
                          aCGCTGGTTTCCACACTGGTGGGAATTGAACTGATATTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGAt                                                                                                                                 <  2:395659/142‑1 (MQ=255)
                          aCGCTGGTTTCCACACTGGTGGGAATTGAACTGATATTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGAt                                                                                                                                 >  1:395659/1‑142 (MQ=255)
                            gCTGGTTTCCACACTGGTGGGAATTGAACTGATATTTAGCGCCGCCAGCCTGTTCGGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGa                                                                                                                                             <  2:987265/128‑1 (MQ=255)
                            gCTGGTTTCCACACTGGTGGGAATTGAACTGATATTTAGCGCCGCCAGCCTGTTCGGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGa                                                                                                                                             >  1:987265/1‑128 (MQ=255)
                                gTTTCCACACTGGTGGGAATTGAACTGATATTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGAtgt                                                                                                                                               >  3:953624/1‑122 (MQ=255)
                                gTTTCCACACTGGTGGGAATTGAACTGATATTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGAtgt                                                                                                                                               <  4:953624/122‑1 (MQ=255)
                                gTTTCCACACTGGTGGGAATTGAACTGATATTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTcc                                                                                                                    >  3:771311/1‑149 (MQ=255)
                                gTTTCCACACTGGTGGGAATTGAACTGATATTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACATCAAGTGACTAt                                                                                                                                         >  1:772837/1‑128 (MQ=255)
                                gTTTCCACACTGGTGGGAATTGAACTGATATTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACATCAAGTGACTAt                                                                                                                                         <  2:772837/128‑1 (MQ=255)
                                    ccACACTGGTGGGAATTGAACTGATATTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACATCAAGTGACTAt                                                                                                                                         <  1:788799/124‑1 (MQ=255)
                                    ccACACTGGTGGGAATTGAACTGATATTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACATCAAGTGACTAt                                                                                                                                         >  2:788799/1‑124 (MQ=255)
                                                 aaTTGAACTGATATTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGa                                                                                                                                                  <  4:1342690/102‑1 (MQ=255)
                                                 aaTTGAACTGATATTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGa                                                                                                                                                  >  3:1342690/1‑102 (MQ=255)
                                                 aaTTGAACTGATATTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGa                                                                                                                                                  >  1:1447399/1‑102 (MQ=255)
                                                 aaTTGAACTGATATTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGa                                                                                                                                                  <  2:1447399/102‑1 (MQ=255)
                                                 aaTTGAACTGATATTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGa                                                                                                                                  >  1:1018610/1‑118 (MQ=255)
                                                 aaTTGAACTGATATTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGa                                                                                                                                  <  2:1018610/118‑1 (MQ=255)
                                                    tGAACTGATATTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATc                                                                                                <  4:771311/149‑1 (MQ=255)
                                                       aCTGATATTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACATCAAGTGACTAt                                                                                                                                         <  1:1408732/105‑1 (MQ=39)
                                                       aCTGATATTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACATCAAGTGACTAt                                                                                                                                         >  2:1408732/1‑105 (MQ=39)
                                                         tGATATTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACATCAAGTGACTa                                                                                                                                          <  4:816484/102‑1 (MQ=255)
                                                         tGATATTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACATCAAGTGACTATGa                                                                                                                                       >  4:536657/1‑105 (MQ=39)
                                                         tGATATTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACATCAAGTGACTATGa                                                                                                                                       <  3:536657/105‑1 (MQ=255)
                                                         tGATATTTAGCGCCGCCAGCCTGTACAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACATCAAGTGACTa                                                                                                                                          >  3:816484/1‑102 (MQ=39)
                                                          gATATTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTgc                                                                                          >  1:821395/1‑149 (MQ=255)
                                                            tatTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGAt                                                                                                                                      <  3:70245/103‑1 (MQ=255)
                                                            tatTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGAt                                                                                                                                      >  4:70245/1‑103 (MQ=255)
                                                             atTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCa                                                                                       <  2:821395/149‑1 (MQ=255)
                                                               ttAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACtt                                                                                                                                                    <  4:1450553/86‑1 (MQ=255)
                                                               ttAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACtt                                                                                                                                                    >  3:1450553/1‑86 (MQ=255)
                                                                  gcgcCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATc                                                                                  >  3:1148081/1‑149 (MQ=255)
                                                                    gccgccAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAgtgt                                                                                                        >  2:325040/1‑125 (MQ=255)
                                                                    gccgccAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAgtat                                                                                                        <  1:325040/125‑3 (MQ=255)
                                                                    gccgccAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCAt                                                                                <  4:1148081/149‑1 (MQ=255)
                                                                            ccTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAaataa                                                                        >  1:136529/1‑149 (MQ=255)
                                                                            ccTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAaataa                                                                        <  2:136529/149‑1 (MQ=255)
                                                                                ttCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATc                                                                                                                                     <  1:1071499/84‑1 (MQ=255)
                                                                                ttCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATc                                                                                                                                     >  2:1071499/1‑84 (MQ=255)
                                                                                          ccAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTACCTGCTAGTGa                                                          <  2:346668/149‑1 (MQ=255)
                                                                                          ccAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTACCTGCTAGTGa                                                          >  1:346668/1‑149 (MQ=255)
                                                                                          ccAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTACCTGCTAGTGa                                                          >  1:929917/1‑149 (MQ=255)
                                                                                           cAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCaa                                                                                              >  1:1214515/1‑112 (MQ=255)
                                                                                           cAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCaa                                                                                              <  2:1214515/112‑1 (MQ=255)
                                                                                            aGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGCTTAATACTCTCTCCGCTACGCAGtgttg                                                                                                      <  4:17912/103‑1 (MQ=255)
                                                                                            aGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGtgttg                                                                                                      >  3:17912/1‑103 (MQ=255)
                                                                                            aGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGa                                                                                                  <  4:1377835/107‑1 (MQ=255)
                                                                                            aGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGa                                                                                                  >  3:1377835/1‑107 (MQ=255)
                                                                                                     gTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGAt                                                                                                                                                 >  4:566744/1‑51 (MQ=255)
                                                                                                     gTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGAt                                                                                                                                                 <  3:566744/51‑1 (MQ=255)
                                                                                                        aaGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCAtt                                                                               >  2:202169/1‑114 (MQ=255)
                                                                                                        aaGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCAtt                                                                               <  1:202169/114‑1 (MQ=255)
                                                                                                          gcagcaATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTa                                                                                                    <  2:136553/91‑1 (MQ=255)
                                                                                                          gcagcaATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTa                                                                                                    >  1:136553/1‑91 (MQ=255)
                                                                                                            agcaATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTACCTGCTAGTGATTATTTCAACCTACTGaa                                        <  2:929917/149‑1 (MQ=255)
                                                                                                                  aatTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATc                                                                                                                                     >  2:755242/1‑50 (MQ=255)
                                                                                                                  aatTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATc                                                                                                                                     <  1:755242/50‑1 (MQ=255)
                                                                                                                        cccGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAgtgt                                                                                                        >  1:26461/1‑73 (MQ=255)
                                                                                                                        cccGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAgtgt                                                                                                        <  2:26461/73‑1 (MQ=255)
                                                                                                                          cGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTACCTGCTAGTGa                                                          <  1:485641/117‑1 (MQ=255)
                                                                                                                          cGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTACCTGCTAGTGa                                                          >  2:485641/1‑117 (MQ=255)
                                                                                                                                  aCCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAAt                                                                                                                             >  1:163814/1‑42 (MQ=255)
                                                                                                                                  aCCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAAt                                                                                                                             <  2:163814/42‑1 (MQ=255)
                                                                                                                                  aCCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTa                                                                     >  3:1009818/1‑98 (MQ=255)
                                                                                                                                  aCCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTa                                                                     <  4:1009818/98‑1 (MQ=255)
                                                                                                                                  aCCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTACCTGCTa                                                              <  3:1041310/105‑1 (MQ=255)
                                                                                                                                  aCCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTACCTGCTa                                                              >  4:1041310/1‑105 (MQ=255)
                                                                                                                                   cccGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCa                                                                                 >  2:403171/1‑85 (MQ=255)
                                                                                                                                   cccGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCa                                                                                 <  1:403171/85‑1 (MQ=255)
                                                                                                                                    ccGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCAtt                                                                               >  1:1339297/1‑86 (MQ=255)
                                                                                                                                    ccGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCAtt                                                                               <  2:1339297/86‑1 (MQ=255)
                                                                                                                                            tAGTCACATCAAGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTACCTGCTAGTGATTATTTCAACCTACTGAATTTCAt                                  >  2:1172936/1‑123 (MQ=255)
                                                                                                                                            tAGTCACATCAAGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTACCTGCTAGTGATTATTTCAACCTACTGAATTTCAt                                  <  1:1172936/123‑1 (MQ=255)
                                                                                                                                              gTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTACCTGCTa                                                              >  2:545383/1‑93 (MQ=255)
                                                                                                                                              gTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTACCTGCTa                                                              <  1:545383/93‑1 (MQ=255)
                                                                                                                                               tCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTACCTGCTAGTGATTATTTCAACCTACTGAATTTc                                    <  1:315437/118‑1 (MQ=255)
                                                                                                                                               tCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTACCTGCTAGTGATTATTTCAACCTACTGAATTTc                                    >  2:315437/1‑118 (MQ=255)
                                                                                                                                                 aCTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTACCTGCTAGTGATTATTTCAACCTACTGAATTTCATCTAATTTTTTTCACTCTATGGCAAATTAGcc   >  2:164801/1‑149 (MQ=255)
                                                                                                                                                  cTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTACCTGCTAGTGATTATTTCAACCTACTGAATTTCATCTAATTTTTTTCACTCTATGGCaaa         >  3:670079/1‑142 (MQ=255)
                                                                                                                                                  cTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTACCTGCTAGTGATTATTTCAACCTACTGAATTTCATCTAATTTTTTTCACTCTATGGCaaa         <  4:670079/142‑1 (MQ=255)
                                                                                                                                                  cTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTACCTGCTAGTGATTATTTCAACCTACTGAATTTCATCTAATTTTTTTCACTCTATGGCAAATTAGCCa  >  1:1305718/1‑149 (MQ=255)
                                                                                                                                                   ttGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCAt                                                                                >  4:1105009/1‑70 (MQ=255)
                                                                                                                                                   ttGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCAt                                                                                <  3:1105009/70‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
TTGGTGCAACGCCGATGGTGTCTGTAACGCTGGTTTCCACACTGGTGGGAATTGAACTGATATTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTACCTGCTAGTGATTATTTCAACCTACTGAATTTCATCTAATTTTTTTCACTCTATGGCAAATTAGCCA  >  NC_000913/3657452‑3657746

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: