Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,181,279 G→T 8.2% R12L (CGT→CTT)  rpoB → RNA polymerase, beta subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,181,2790GT 8.2% 301.4 / 7.3 92R12L (CGT→CTT) rpoBRNA polymerase, beta subunit
Reads supporting (aligned to +/- strand):  ref base G (60/24);  new base T (4/4);  total (64/28)
Fisher's exact test for biased strand distribution p-value = 2.40e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGCGGAAAGTGTTCCATTTTCCGGTCAACAAAATAGTGTTGCACAAACTGTCCGCTCAATGGACAGATGGGTCGACTTGTCAGCGAGCTGAGGAACCCTATGGTTTACTCCTATACCGAGAAAAAACGTATTCGTAAGGATTTTGGTAAACGTCCACAAGTTCTGGATGTACCTTATCTCCTTTCTATCCAGCTTGACTCGTTTCAGAAATTTATCGAGCAAGATCCTGAAGGGCAGTATGGTCTGGAAGCTGCTTTCCGTTCCGT  >  NC_000913/4181146‑4181411
                                                                                                                                     |                                                                                                                                    
tGCGGAAAGTGTTCCATTTTCCGGTCAACAAAATAGTGTTGCACAAACTGTCCGCTCAATGGACAGATGGGTCGACTTGTCAGCGAGCTGAGGAACCCTATGGTTTACTCCTATACCGAGAAAAAACGTATTCGTAAGGa                                                                                                                                >  2:66588/1‑140 (MQ=255)
           ttCCATTTTCCGGTCAACAAAATAGTGTTGCACAAACTGTCCGCTCAATGGACAGATGGGTCGACTTGTCAGCGAGCTGAGGAACCCTATGGTTTACTCCTATACCGAGAAAAAACGTATTCGTAAGGATTTTGGTAAAc                                                                                                                     >  1:679539/1‑140 (MQ=255)
            tCCATTTTCCGGTCAACAAAATAGTGTTGCACAAACTGTCCGCTCAATGGACAGATGGGTCGACTTGTCCGCGAGCTGAGGAACCCTATGGTTTACTCCTATACCGAGAAAAAACGTATTCGTAAGGATTTTGGTAAAct                                                                                                                    >  2:1275549/1‑139 (MQ=255)
                 tttCCGGTCAACAAAATAGTGTTGCACAAACTGTCCGCTCAATGGACAGCTGGGTCGACTTGTCAGCGAGCTGAGGAACCCTATGGTTTACTCCTATACCGAGAAAAAACGTATTCGTAAGGATTTTGGTAAACGTCcac                                                                                                               <  2:1023747/140‑1 (MQ=255)
                 tttCCGGTCAACAAAATAGTGTTGCACAAACTGTCCGCTCAATGGACAGATGGGTCGACTTGTCAGCGAGCTGAGGAACCCTATGGTTTACTCCTATACCGAGAAAAAACGTATTCGTAAGGATTTTGGTAAACGTCcac                                                                                                               <  2:1023831/140‑1 (MQ=255)
                 tttCCGGTCAACAAAATAGTGTTGCACAAACTGTCCGCTCAATGGACAGATGGGTCGACTTGTCAGCGAGCTGAGGAACCCTATGGTTTACTCCTATACCGAGAAAAAACGTATTCGTAAGGATTTTGGTAAACGTCcac                                                                                                               <  1:1242392/140‑1 (MQ=255)
                  ttCCGGTCAACAAAATAGTGTTGCACAAACTGTCCGCTCAATGGACAGATGGGTCGACTTGTCAGCGAGCTGAGGACCCCTATGGTTTACTCCTATACCGAGAAAAAACGTATTCGTAAGGATTTTGGTAAACGTCcaca                                                                                                              >  2:1078321/1‑140 (MQ=255)
                  ttCCGGTCAACAAAATAGTGTTGCACAAACTGTCCGCTCAATGGACAGATGGGTCGACTTGTCAGCGAGCTGAGGAACCCTATGGTTTACTCCTATACCGAGAAAAAACGTATTCGTAAGGATTTTGGTAAACGTCcaca                                                                                                              >  2:1078357/1‑140 (MQ=255)
                   tCCGGTCAACAAAATAGTGTTGCACAAACTGTCCGCTCAATGGACAGATGGGTCGACTTGTCAGCGAGCTGAGGAACCCTATGGTTTACTCCTATACCGAGAAAAAACGTATTCGTAAGGATTTTGGTAAACGTCCACaa                                                                                                             >  1:660626/1‑140 (MQ=255)
                   tCCGGTCAACAAAATAGTGTTGCACAAACTGTCCGCTCAATGGACAGATGGGTCGACTTGTCAGCGAGCTGAGGAACCCTATGGTTTACTCCTATACCGAGAAAAAACGTATTCGTAAGGATTTTGGTAAACGTCCACaa                                                                                                             >  1:634890/1‑140 (MQ=255)
                   tCCGGTCAACAAAATAGTGTTGCACAAACTGTCCGCTCAATGGACAGATGGGTCGACTTGTCAGCGAGCTGAGGAACCCTATGGTTTACTCCTATACCGAGAAAAAACGTATTCGTAAGGATTTTGGTAAACGTCCACaa                                                                                                             >  2:522880/1‑140 (MQ=255)
                   tCCGGTCAACAAAATAGTGTTGCACAAACTGTCCGCTCAATGGACAGATGGGTCGACTTGTCAGCGAGCTGAGGAACCCTATGGTTTACTCCTATACCGAGAAAAAACCTATTCGTAAGGATTTTGGTAAACGTCCACaa                                                                                                             >  1:1150012/1‑140 (MQ=255)
                        tCAACAAAATAGTGTTGCACAAACTGTCCGCTCAAAGGACAGATGGGTCGACTTGTCAGGGAGCTGAGGAACCCTATGGTTTTCTCCTATACCGAGAAAAAACTTATACGTAAGGATTTTGGTAAACGTCGACAAGATCt                                                                                                        >  2:820544/1‑140 (MQ=255)
                             aaaaTAGTGTTGCACAAACTGTCCGCTCAATGGACAGATGGGTCGACTTGTCAGCGAGCTGAGGAACCCTATGGTTTACTCCTATACCGAGAAAAAACGTATTCGTAAGGATTTTGGTAAACGTCCACAAGTTCTTGATg                                                                                                   <  1:1436903/140‑1 (MQ=255)
                              aaaTAGTGTTGCACAAACTGTCCGCTCAATGGACAGATGGGTCGACTTGTCAGCGAGCTGAGGAACCCTATGGTTTACTCCTATACCGAGAAAAAACGTATTCGTAAGGATTTTGGTAAACGTCCACAAGTTCTGGATGt                                                                                                  <  2:205192/140‑1 (MQ=255)
                              aaaTAGTGTTGCACAAACTGTCCGCTCAATGGACAGATGGGTCGACTTGTCAGCGAGCTGAGGAACCCTATGGTTTACTCCTATACCGAGAAAAAACGTATTCGTAAGGATTTTGGTAAACGTCCACAAGTTCTGGATGt                                                                                                  <  2:552393/140‑1 (MQ=255)
                              aaaTAGTGTTGCACAAACTGTCCGCTCAATGGACAGATGGGTCGACTTGTCAGCGAGCTGAGGAACCCTATGGTTTACTCCTATACCGAGAAAAAACGTATTCGTAAGGATTTTGGTAAACGTCCACAAGTTCTGGATGt                                                                                                  <  2:552331/140‑1 (MQ=255)
                                    tgtTGCACAAACTGTCCGCTCAATGGACAGATGGGTCGACTTGTCAGCGAGCTGAGGAACCCTATGGTTTACTCCTATACCGAGAAAAAACGTATTCTTAAGGATTTTGGTAAACGTCCACAAGTTCTGGATGTACCTTa                                                                                            <  2:1097565/140‑1 (MQ=255)
                                    tgtTGCACAAACTGTCCGCTCAATGGACAGATGGGTCGACTTGTCAGCGAGCTGAGGAACCCTATGGTTTACTCCTATACCGAGAAAAAACGTATTCGTAAGGATTTTGGTAAACGTCCACAAGTTCTGGATGTACCTTa                                                                                            >  2:1048429/1‑140 (MQ=255)
                                    tgtTGCACAAACTGTCCGCTCAATGGACAGATGGGTCGACTTGTCAGCGAGCTGAGGAACCCTATGGTTTACTCCTATACCGAGAAAAAACGTATTCGTAAGGATTTTGGTAAACGTCCACAAGTTCTGGATGTACCTTa                                                                                            >  2:1048392/1‑140 (MQ=255)
                                    tgtTGCACAAACTGTCCGCTCAATGGACAGATGGGTCGACTTGTCAGCGAGCTGAGGAACCCTATGGTTTACTCCTATACCGAGAAAAAACGTATTCGTAAGGATTTTGGTAAACGTCCACAAGTTCTGGATGTACCTTa                                                                                            >  2:534105/1‑140 (MQ=255)
                                    tgtTGCACAAACTGTCCGCTCAATGGACAGATGGGTCGACTTGTCAGCGAGCTGAGGAACCCTATGGTTTACTCCTATACCGAGAAAAAACGTATTCGTAAGGATTTTGGTAAACGTCCACAAGTTCTGGATGTACCTTa                                                                                            >  2:534195/1‑140 (MQ=255)
                                    tgtTGCACAAACTGTACGCTAAATGGACAGATGGGTCGACTTGTCAGCGAGCTGAGGAACCCTATGGTTTACTCCTATACCGAGAAAAAACGTATTCGTAAGGATTTTGGTAAACGTCCACAAGTTCTGGATGTACCTTa                                                                                            >  1:1489576/1‑140 (MQ=255)
                                    tgtTGCACAAACAGTCCGCTCAATGGACAGATGGGTCGACTTGTCAGCGAGCTGATGAACCCTAAGGTTTACTCCTAAACCGAGAAAAAACGTATTCGTAAGGATTTTGTTAAACGTCCACAAGTTCAGGATTTAACTTa                                                                                            >  2:1481764/1‑140 (MQ=255)
                                       tGCACAAACTGTCCGCTCAATGGACAGATGGGTCGACTTGTCAGCGAGCTGAGGAACCCTATGGTTTACTCCTATACCGAGAAAAAACGTATTCTTAAGGATTTTGGTAAACGTCCACAAGTTCTGGATGTACCTTAtct                                                                                         <  1:727413/140‑1 (MQ=255)
                                       tGCACAAACTGTCCGCTCAATGGACAGATGGGTCGACTTGTCAGCGAGCTGAGGAACCCTATGGTTTACTCCTATACCGAGAAAAAACGTATTCGTAAGGATTTTGGTAAACGTCCACAAGTTCTGGATGTACCTTAtct                                                                                         <  1:1111090/140‑1 (MQ=255)
                                       tGCACAAACTGTCCGCTCAATGGACAGATGGGTCGACTTGTCAGCGAGATGAGGAACCCTATGGTTTACTCCTATACCGAGAAAAAACGTATTCGTAAGGATTTTGGTAAACGTCCACAAGTTCTGGATGTACCTTAtct                                                                                         <  1:936937/140‑1 (MQ=255)
                                              aCTGTCCGCTCAATGGACAGATGGGTCGACTTGTCAGCGAGCTGAGGAACCCTATGGTTTACTCCTATACCGAGTAAAAACGTATTCGTAAGGATTTTGGTAAACGTCCACAAGTTCTGGATGTACCTTATCTCCTTTCt                                                                                  >  1:1178697/1‑140 (MQ=255)
                                               cTGTCCGCTCAATGGACAGATGGGTCGACTTGTCAGCGAGCTGAGGAACCCTATGGTTTACTCCTATACCGAGAAAAAACGTATTCGTAAGGATTTTGGTAAACGTCCACAAGTTCTGGATGTACCTTATCTCCTTTCTa                                                                                 <  2:825134/140‑1 (MQ=255)
                                                tGTCCGCTCAATGGACAGATGGGTCGACTTGTCAGCGAGCTGAGGAACCCTATGGTTTACTCCTATACCGAGAAAAAACGTATTCGTAAGGATTTTGGTAAACGTCCACAAGTTCTGGATGTACCTTATCTCCTTTCTAt                                                                                <  2:526261/140‑1 (MQ=255)
                                                  tCCGCTCAATGGACAGATGGGTCGACTTGTCAGCGAGCTGGGGAACCCTATGGTTTACTCCTATACCGAGAAAAAACGTATTCTTAAGGATTTTGGTAAACGTCCACAAGTTCTGGATGTACCTTATCTCCTTTCTATcc                                                                              >  1:202790/1‑140 (MQ=255)
                                                  tCCGCTCAATGGACAGATGGGTCGACTTGTCAGCGAGCTGAGGAACCCTATGGTTTACTCCTATACCGAGAAAAAACGTATTCTTAAGGATTTTGGTAAACGTCCACAAGTTCTGGATGTACCTTATCTCCTTTCTATcc                                                                              >  1:201034/1‑140 (MQ=255)
                                                  tCCGCTCAATGGACAGATGGGTCGACTTGTCAGCGAGCTGAGGAACCCTATGGTTTACTCCTATACCGAGAAAAAACGTATTCGTAAGGATTTTGGTAAACGTCCACAAGTTCTGGATGTACAGTATCTCCTATCTATcc                                                                              >  1:1685963/1‑140 (MQ=255)
                                                   ccGCTCAATGGACAGATGGGTCGACTTGTCAGCGAGCTGAGGAACCCTATGGTTTACTCCTATACCGAGAAAAAACGTATTCGTAAGGATTTTGGTAAACGTCCACAAGTTCTGGATGTACCTTATCTCCTTTCTATCCa                                                                             >  2:596553/1‑140 (MQ=255)
                                                     gCTCAATGGACAGATGGGTCGACTTGTCAGCGAGCTGAGGAACCCTATGGTTTACTCCTATACCGAGAAAAAACGTATTCGTAAGGATTTTGGTAAACGTCCACACGTTCTGGATGTACCTTATCTCCTTTCTATCCAGc                                                                           >  1:1679724/1‑140 (MQ=255)
                                                     gCTCAATGGACAGATGGGTCGACTTGTCAGCGAGCTGAGGAACCCTATGGTTTACTCCTATACCGAGAAAAAACGTATTCGTAAGGATTTTGGTAAACGTCCACAAGTTCTGGATGTACCTTATCTCCTTTCTATCCAGc                                                                           >  1:290112/1‑140 (MQ=255)
                                                          aTGGACAGATGGGTCGACTTGTCAGCGAGCTGAGGAACCCTATGGTTTACTCCTATACCGAGAAAAAACGTATTCGTAAGGATTTTGGTAAACGTCCACAAGTTCTGGATGTACCTTATCTCCTTTCTATCCAGCTTGAc                                                                      >  1:1041944/1‑140 (MQ=255)
                                                            ggACAGATGGGTCGACTTGTCAGCGAGCTGAGGAACCCTATGGTTTACTCCTATACCGAGAAAAAACGTATTCGTAAGGATTTTGGTAAACGTCCACAAGTTCTGGATGTGCCTTATCTCCTTTCTATCCAGCTTGACTc                                                                    >  2:210105/1‑140 (MQ=255)
                                                            ggACAGATGGGTCGACTTGTCAGCGAGCTGAGGAACCCTATGGTTTACTCCTATACCGAGAAAAAACGTATTCGTAAGGATTTTGGTAAACGTCCACAAGTTCTGGATGTACCTTATCTCCTTTCTATCCAGCTTGACTc                                                                    >  2:1560934/1‑140 (MQ=255)
                                                            ggACAGATGGGTCGACTTGTCAGCGAGCTGAGGAACCCTATGGTTTACTCCTATACCGAGAAAAAACGTATTCGTAAGGATTTTGGTAAACGTCCACAAGTTCTGGATGTACCTTATCTCCTTTCTATCCAGCTTGACTc                                                                    >  2:1560985/1‑140 (MQ=255)
                                                                 gATGGGTCGACTTGACAGCGAGCTGAGGAACCCTATGGATTACTCCTAAAACGAGAAAAAACGTATTCGTAAGGATTATGGTAAACGTCCACAAGTTCTGGATGTACCTAATCTCCTTTCTATCCAGCATGACTCTTTTc                                                               >  2:538648/1‑140 (MQ=255)
                                                                          aCTTGTCAGCGAGCTGAGGAACCCTATGGTTTACTCCTATACCGAGAAAAAACGTATTCTTAAGGATTTTGGTAAACGTCCACAAGTTCTGGATGTACCTTATCTCCTTTCTATCCAGCTTGACTCGTTTCAGAAATTTa                                                      <  1:1474722/140‑1 (MQ=255)
                                                                          aCTTGTCAGCGAGCTGAGGAACCCTATGGTTTACTCCTATACCGAGAAAAAACGTATTCGTAAGGATTTTGGTAAACGTCCACAAGTTCTGGATGTACCTTATCTCCTTTCTATCCAGCTTGACTCGTTTCAGAAATTTa                                                      <  2:1041944/140‑1 (MQ=255)
                                                                            ttGTCAGCGAGCTGAGGAACCCTATGGTTTACTCCTATACCGAGAAAAAACGTATTCTTAAGGATTTTGGTAAACGTCCACAAGTTCTGGATGTACCTTATCTCCTTTCTATCCAGCTTGACTCGTTTCAGAAATTTATc                                                    >  2:1523059/1‑140 (MQ=255)
                                                                            ttGTCAGCGAGCTGAGGAACCCTATGGTTTACTCCTATACCGAGAAAAAACGTATTCGTAAGGATTTTGGTAAACGTCCACAAGTTCTGGATGTACCTTATCTCCTTTCTATCCAGCTTGACTCGTTTCAGAAATTTATc                                                    >  1:725339/1‑140 (MQ=255)
                                                                            ttGTCAGCGAGCTGAGGAACCCTATGGTTTACTCCTATACCGAGAAAAAACGTATTCGTAAGGATTTTGGTAAACGTCCACAAGTTCTGGATGTACCTTATCTCCTTTCTATCCAGCTTGACTCGTTTCAGAAATTTATc                                                    >  1:174249/1‑140 (MQ=255)
                                                                            ttGTCAGCGAGCTGAGGAACCCTATGGTTTACTCCTATACCGAGAAAAAACGTATTCGTAAGGATTTTGGTAAACGTCCACAAGTTCTGGATGTACCTTATCTCCTTTCTATCCAGCTTGACTCGTTTCAGAAATTTATc                                                    >  1:724888/1‑140 (MQ=255)
                                                                            ttGTCAGCGAGCTGAGGAACCCTATGGTTTACTCCTATACCGAGAAAAAACGTATTCGTAAGGATTTTGGTAAACGTCCACAAGTTCTGGATGTACCTTATCTCCTTTCTATCCAGCTTGACTCGTTTCAGAAATTTATc                                                    >  1:1013791/1‑140 (MQ=255)
                                                                                 aGCGAGCTGAGGAACCCTATGGTTTACTCCTATACCGAGAAAAAACGTATTCGTAAGGATTTTGGTAAACGTCCACAAGTTCTGGATGTACCTTATCTCCTTTCTATCCAGCTTGACTCGTTTCAGAAATTTATCGAGCa                                               <  1:617390/140‑1 (MQ=255)
                                                                                 aGCGAGCTGAGGAACCCTATGGTTTACTCCTATACCGAGAAAAAACGTATTCGTAAGGATTTTGGTAAACGTCCACAAGTTCTGGATGTACCTTATCTCCTTTCTATCCAGCTTGACTCGTTTCAGAAATTTATCGAGCa                                               >  1:1051899/1‑140 (MQ=255)
                                                                                           ggAACCCTATGGTTTACTCCTATACCGAGAAAAAACGTATTCTTAAGGATTTTGGTAAACGTCCACAAGTTCTGGATGTACCTTATCTCCTTTCTATCCAGCTTGACTCGTTTCAGAAATTTATCGAGCAAGATCCTGaa                                     <  2:1506255/140‑1 (MQ=255)
                                                                                           ggAACCCTATGGTTTACTCCTATACCGAGAAAAAACGTATTCGTAAGGATTTTGGTAAACGTCCACAAGTTCTGGATGTACCTTATCTCCTTTCTATCCAGCTTGACTCGTTTCAGAAATTTATCGAGCAAGATCCTGaa                                     <  2:1428289/140‑1 (MQ=255)
                                                                                              aCCCTATGGTTTACTCCTATACCGAGAAAAAACGTATTCGTAAGGATTTTGGTAAACGTCCACAAGTTCTGGATGTACCTTATCTCCTTTCTATCCAGCTTGACTCGTTTCAGAAATTTATCGAGCAAGATCCTGAAggg                                  <  1:1221592/140‑1 (MQ=255)
                                                                                              aCCCTATGGTTTACTCCTATACCGAGAAAAAACGTATTCGTAAGGATTTTGGTAAACGTCCACAAGTTCTGGATGTACCTTATCTCCTTTCTATCCAGCTTGACTCGTTTCAGAAATTTATCGAGCAAGATCCTGAAggg                                  <  2:263379/140‑1 (MQ=255)
                                                                                                  tATGGTTTACTCCTATACCGAGAAAAAACGTATTCGTAAGGATTTTGGTAAACGTCCACAAGTTCTGGATGTACCTTATCTCCTTTCTATCCAGCTTGACTCGTTTCAGAAATTTATCGAGCAAGATCCTGAAGGGCAGt                              >  1:376864/1‑140 (MQ=255)
                                                                                                    tGGTTTACTCCTATACCGAGAAAAAACGTATTCGTAAGGATTTTGGTAAACGTCCACAAGTTCTGGATGTACCTCATCTCCTTTCTATCCAGCTTGACTCGTTTCAGAAATTTATCGAGCAAGATCCTGAAGGGCAGTAt                            >  1:609456/1‑140 (MQ=255)
                                                                                                    tGGTTTACTCCTATACCGAGAAAAAACGTATTCGTAAGGAGTTTGGTAAACGTCCACAAGTTCTGGATGTACCTTATCTCCTTTCTATCCAGCTTGACTCGTTTCAGAAATTTATCGAGCAAGATCCTGAAGGGCAGTAt                            >  1:609415/1‑140 (MQ=255)
                                                                                                         tACTCCTATACCGAGAAAATACGTATTCGTAAGGATTTTGGTAAACGTCCACAAGTTCTGGATGTACGTTATCTCCTTTCTATCCAGCTTGTCTCGTTTCTGTAATTTATAGAGCTAGATCCTGTAGGGCAGTATGGTCt                       >  1:925334/1‑140 (MQ=255)
                                                                                                         tACTCCTATACCGAGAAAAAACGTATTCGTAAGGATTTTGGTAAACGTCCACAAGTTCTGGATGTACCTTATCTCCTTTCTATCCAGCTTGACTCGTTTCAGAAATTTATCGAGCAAGATCCTGAAGGGCAGTATGGTCt                       >  1:76930/1‑140 (MQ=255)
                                                                                                         tACTCCTATACCGAGAAAAAACGTATTCGTAAGGATTTTGGTAAACGTCCACAAGTTCTGGATGTACCTTATCTCCTTTCTATCCAGCTTGACTCGTTTCAGAAATTTATCGAGCAAGATCCTGAAGGGCAGTATGGTCt                       >  1:1077116/1‑140 (MQ=255)
                                                                                                         tACTCCTATACCGAGAAAAAACGTATTCGTAAGGATTTTGGTAAACGTCCACAAGTTCTGGATGTACCTTATCTCCTTTCTATCCAGCTTGACTCGTTTCAGAAATTTATCGAGCAAGATCCTGAAGGGCAGTATGGTCt                       >  1:89779/1‑140 (MQ=255)
                                                                                                         tACTCCTATACCGAGAAAAAACGTATTCGTAAGGATTTTGGTAAACGTCCACAAGTTCTGGATGTACCTTATCTCCTTTCTATCCAGCTTGACTCGTTTCAGAAATTTATCGAGCAAGATCCTGAAGGGCAGTATGGTCt                       <  2:1017439/140‑1 (MQ=255)
                                                                                                         tACTCCTATACCGAGAAAAAACGTAGTCGTAAGGATTTTGGTAAACGTCCACAAGTTCTGGATGTACCTTATCTCCTTTCTATCCAGCTTGACTCGTTTCAGAAATTTATCGAGCAAGATCCTGAAGGGCAGTATGGTCt                       <  2:1017433/140‑1 (MQ=255)
                                                                                                          aCTCCTATACCGAGAAAAAACGTATTCGTAAGGATTTTGGTAAACGTCCACAAGTTCTGGATGTACCTTATCTCCTTTCTATCCAGCTTGACTCGTTTCAGAAATTTATCGAGCAAGATCCTGAAGGGCAGTa                             >  1:1154857/1‑133 (MQ=255)
                                                                                                          aCTCCTATACCGAGAAAAAACGTATTCGTAAGGATTTTGGTAAACGTCCACAAGTTCTGGATGTACCTTATCTCCTTTCTATCCAGCTTGACTCGTTTCAGAAATTTATCGAGCAAGATCCTGAAGGGCAGTa                             <  2:1154857/133‑1 (MQ=255)
                                                                                                          aCTCCTATACCGAGAAAAAACGTATTCGTAAGGATTTTGGTAAACGTCCACAAGTTCTGGATGTACCTTATCTCCTTTCTATCCAGCTTGACTCGTTTCAGAAATTTATCGAGCAAGATCCTGAAGGGCAGTATGGTCTg                      >  2:749001/1‑140 (MQ=255)
                                                                                                          aCTCCTATACCGAGAAAAAACGTATTCGTAAGGATTTTGGTAAACGTCCACAAGTTCTGGATGTACCTTATCTCCTTTCTATCCAGCTTGACTCGTTTCAGAAATTTATCGAGCAAGATCCTGAAGGGCAGTATGGTCTg                      >  2:18781/1‑140 (MQ=255)
                                                                                                          aCTCCTATACCGAGAAAAAACGTATTCGTAAGGATTTTGGTAAACGTCCACAAGTTCTGGATGTACCTTATCTCCTTTCTATCCAGCTTGACTCGTTTCAGAAATTTATCGAGCAAGATCCTGAAGGGCAGTATGGTCTg                      >  1:1576965/1‑140 (MQ=255)
                                                                                                          aCTCCTATACCGAGAAAAAACGTATTCGTAAGGATTTTGGTAAACGTCCACAAGTTCTGGATGTACCTTATCTCCTTTCTATCCAGCTTGACTCGTGTCAGAAATTTATCGAGCAAGATCCTGAAGGGCAGTATGGTCTg                      >  1:40120/1‑140 (MQ=255)
                                                                                                                ataCCGAGAAAAAACGTATTCGTAAGGATTTTGGTAAACGTCCACAAGTTCTGGATGTACCTTGTCTCCTTTCTATCCAGCTTGACTCGTTTCAGAAATTTATCGAGCAAGATCCTGAAGGGCAGTATGGTCTGGAAgct                >  1:173905/1‑140 (MQ=255)
                                                                                                                ataCCGAGAAAAAACGTATTCGTAAGGATTTTGGTAAACGTCCACAAGTTCTGGATGTACCTTATCTCCTTTCTATCCAGCTTGACTCGTTTCAGAAATTTATCGAGCAAGATCCTGAAGGGCAGTATGGTCTGGAAgct                <  1:430186/140‑1 (MQ=255)
                                                                                                                ataCCGAGAAAAAACGTATTCGTAAGGATTTTGGTAAACGTCCACAAGTTCTGGATGTACCTTATCTCCTTTCTATCCAGCTTGACTCGTTTCAGAAATTTATCGAGCAAGATCCTGAAGGGCAGTATGGTCTGGAAgct                <  1:430185/140‑1 (MQ=255)
                                                                                                                ataCCGAGAAAAAACGTATTCGTAAGGATTTTGGTAAACGTCCACAAGTTCTGGATGTACCTTATCTCCTTTCTATCCAGCTTGACTCGTTTCAGAAATTTATCGAGCAAGATCCTGAAGGGCAGTATGGTCTGGAAgct                <  1:430213/140‑1 (MQ=255)
                                                                                                                 taCCGAGAAAAAACGTATTCGTAAGGATTTTGGTAAACGTCCACAAGTTCTGGATGTACCTTATCTCCTTTCTATCCAGCTTGACTCGTTTCAGAAATTTATCGAGCAAGATCCTGAAGGGCAGTATGGTCTGGAAgctg               >  2:1320235/1‑140 (MQ=255)
                                                                                                                      agaAAAAACGTATTCTTAAGGATTTTGGTAAACGTCCACAAGTTCTGGTTGTACGTTATGTCCTTTCTATCCAGCTTGAGTCGTTTCTGAAATTTATGGAGCAAGATCGTGTAGGGCAGTATGGTCTGGTAGATGCtttg          >  2:579899/1‑139 (MQ=255)
                                                                                                                      agaAAAAACGTATTCGTAAGGATTTTGGTAAGCGTCCACAAGTTCTGGATGTACCTTATCTCCTTTCTATCCAGCTTGACTCGTTTCAGAAATTTATCGAGCAAGATCCTGAAGGGCAGTATGGTCTGGAAGCTGCTttc          >  2:1383935/1‑140 (MQ=255)
                                                                                                                      agaAAAAACGTATTCGTAAGGATTTTGGTAAACGTCCACAAGTTCTGGATGTACCTTATCTCCTTTCTATCCAGCTTGACTCGTTTCAGAAATTTATGGAGCAAGTTCCTGAAGGGCAGTATGGTCTGGAAGCTGCTttc          >  2:579967/1‑140 (MQ=255)
                                                                                                                      agaAAAAACGTATTCGTAAGGATTTTGGTAAACGTCCACAAGTTCTGGATGTACCTTATCTCCTTTCTATCCAGCTTGACTCGTTTCAGAAATTTATCGAGCAAGATCCTGAAGGGCAGTATGGTCTGGAAGCTGCTttc          >  2:782636/1‑140 (MQ=255)
                                                                                                                      agaAAAAACGTATTCGTAAGGATTTTGGTAAACGTCCACAAGTTCTGGATGTACCTTATCTCCTTTCTATCCAGCTTGACTCGTTTCAGAAATTTATCGAGCAAGATCCTGAAGGGCAGTATGGTCTGGAAGCTGCTttc          >  2:578984/1‑140 (MQ=255)
                                                                                                                      agaAAAAACGTATTCGTAAGGATTTTGGTAAACGTCCACAAGTTCTGGATGTACCTTATCTCCTTTCTATCCAGCTTGACTCGTTTCAGAAATTTATCGAGCAAGATCCTGAAGGGCAGTATGGTCTGGAAGCTGCTttc          >  2:994757/1‑140 (MQ=255)
                                                                                                                      agaAAAAACGTATTCGTAAGGATTTTGGTAAACGTCCACAAGTTCTGGATGTACCTTATCTCCTTTCTATCCAGCTTGACTCGTTTCAGAAATTTATCGAGCAAGATCCTGAAGGGCAGTATGGTCTGGAAGCTGCTttc          >  2:1251274/1‑140 (MQ=255)
                                                                                                                      agaAAAAACGTATTCGTAAGGATTTTGGTAAACGTCCACAAGTTCTGGATGTACCTTATCTCCTTTCTATCCAGCTTGACTCGTTTCAGAAATTTATCGAGCAAGATCCTGAAGGGCAGTATGGTCTGGAAGCTGCTttc          >  2:1245922/1‑140 (MQ=255)
                                                                                                                      agaAAAAACGTATTCGTAAGGATTTTGGTAAACGTCCACAAGTTCTGGATGTACCTTATCTCCTTTCTATCCAGCTTGACTCGTTTCAGAAATTTATCGAGCAAGATCCTGAAGGGCAGTATGGTCTGGAAGCTGCTttc          >  1:96743/1‑140 (MQ=255)
                                                                                                                      agaAAAAACGTATTCGTAAGGATTTTGGTAAACGTCCACAAGTTCTGGATGTACCTTATCTCCTTTCTATCCAGCTTGACTCGTTTCAGAAATTTATCGAGCAAGATCCTGAAGGGCAGTATGGTCTGGAAGCTGCTttc          >  1:96710/1‑140 (MQ=255)
                                                                                                                      agaAAAAACGTATTCGTAAGGATTTTGGTAAACGTCCACAAGTTCTGGATGTACCTTATCTCCTTTCTATCCAGCTTGACTCGTTTCAGAAATTTATCGAGCAAGATCCTGAAGGGCAGTATGGTCTGGAAGCTGCTttc          >  1:449479/1‑140 (MQ=255)
                                                                                                                      agaAAAAACGTATTCGTAAGGATTTTGGTAAACGTCCACAAGTTCTGGATGTACCTTATCTCCTTTCTATCCAGCTTGACTCGTTTCAGAAATTTATCGAGCAAGATCCTGAAGGGCAGTATGGTCTGGAAGCTGCTttc          >  1:33083/1‑140 (MQ=255)
                                                                                                                      agaAAAAACGTATTCGTAAGGATTTTGGTAAACGTCCACAAGTTCTGGATGTACCTTATCTCCTTTCTATCCAGCTTGACTCGTTTCAGAAATTTATCGAGCAAGATCCTGAAGGGCAGTATGGTCTGGAAGCTGCTttc          >  1:33020/1‑140 (MQ=255)
                                                                                                                      agaAAAAACGTATTCGTAAGGATTTTGGTAAACGTCCACAAGTTCTGGATGTACCTTATCTCCTTTCTATCCAGCTTGACTCGTTTCAGAAATTTATCGAGCAAGATCCTGAAGGGCAGTATGGTCTGGAAGCTGCTttc          >  1:1298343/1‑140 (MQ=255)
                                                                                                                      agaAAAAACGTATTCGTAAGGATTTTGGTAAACGTCCACAAGTTCTGGATGTACCTTATCTCCTTTCTATCCAGCTTGACTCGTTTCAGAAATTTATCGAGCAAGATCCTGAAGGGCAGTATGGTCTGGAAGCTGCTttc          >  1:1282556/1‑140 (MQ=255)
                                                                                                                      agaAAAAACGTATTCGTAAGGATTTTGGTAAACGTCCACAAGTTCTGGATGTACCTTATCTCCTTTCTATCCAGCTTGACTCGTTTCAGAAATTTATCGAGCAAGATCCTGAAGGGCAGTATGGTCTGGAAGCTGCTttc          >  1:1215627/1‑140 (MQ=255)
                                                                                                                        aaaaaaCGTATTCGTAAGGATTTTGGTAAACGTCCACAAGTTCTGGATGTACCTTATCTCCTTTCTATCCAGCTTGACTCGTTTCAGAAATTTATCGAGCAAGATCCTGAAGGGCAGTATGGTCTGGAAGCTGCTttccg        <  1:565487/140‑1 (MQ=255)
                                                                                                                              cGTATTCGTAAGGATTTTGGTAAACGTCCACAAGTTCTGGATGTACCTTATCTCCTTTCTATCCAGCTTGACTCGTTTCAGAAATTTATCGAGCAAGATCCTGAAGGGCAGTATGGTCTGGAAGCTGCTTtccgttccgt  <  2:1432778/140‑1 (MQ=255)
                                                                                                                                     |                                                                                                                                    
TGCGGAAAGTGTTCCATTTTCCGGTCAACAAAATAGTGTTGCACAAACTGTCCGCTCAATGGACAGATGGGTCGACTTGTCAGCGAGCTGAGGAACCCTATGGTTTACTCCTATACCGAGAAAAAACGTATTCGTAAGGATTTTGGTAAACGTCCACAAGTTCTGGATGTACCTTATCTCCTTTCTATCCAGCTTGACTCGTTTCAGAAATTTATCGAGCAAGATCCTGAAGGGCAGTATGGTCTGGAAGCTGCTTTCCGTTCCGT  >  NC_000913/4181146‑4181411

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: