Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,182,389 G→A P16S (CCA→TCA)  ygiD ← 4,5‑DOPA dioxygenase extradiol

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,182,3890GA98.1% 137.7 / NA 52P16S (CCA→TCA) ygiD4,5‑DOPA dioxygenase extradiol
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (26/25);  minor base T (0/1);  total (26/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.73e-01

CCGGAGCAGGATAATGCGTATCGTACAGCGCCTGCGGGAAGCCAC‑CAAAGTCATGAATCGTGGGCGGCGTCTCCATCGCGGTCACTCCTGTTCCACGGGTAAACCAGTGAGCCGAAACCACCACAATCGCTTGCGGGCGTGGCAATGTCATCCCCAACTTCTGCCAGCTGCGG‑GTATACAAATTATCTTCCAGCACGTTCATCGGACTACCGTGACCTAAAAACAATGCTGGCATACGTGTTGAAGACATGATGATATCCTTAACTAAAGGTGTCATTTTGATATCCTCACAATACGCTTGTTCGGCGGAGTAAGAACCCGGATAACAATGATGATGATCATCAGTTATTTTGACGATCTGCCTGAAGGTGAAGATTTATAAGGAGTTGTCGATGTCAGTACCTCTCATTCTGACCATACTGGCGGGGGCAGCCACGTTTATTGGCGCGTTTCTCGGCGTTCTCGGGCAAAAACCCT  >  NC_000913/3182158‑3182634
                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                     
ccGGAGCAGGATAATGAGTATCGTACAGCGCCTGCGGGAAGCCAC‑CAAAGTCATGAATCGTGGGCGGCGTCTCCATCGCGGTCACTCCTGTTCCACGGGTAAACCAGTGAGCCGAAACCACCACAATCGCTTGCGGGCGTGGCAATGTCATCCCCACCTTCTGCCAGCTGCGG‑GTATACAAATTATCTTCCAGCACGTTCATCGGACTACCGTGACCTAAACACAATGCTGACATACGTTTTGCAGACatg                                                                                                                                                                                                                                    >  2:377419/1‑251 (MQ=255)
      cAGGATAAGGCGTATGGTACAGCTCCTGCGGGAAGCCAC‑CAAAGTCATGAGTCGTTGGCGGCGGCACCAACGTGGTCACGCCTGGTCCTCGGGTAAACCAGTGAGCCGAAACCACCACAATCGCTTGCGGTATTGGCAATGTCATCCCCAACTTCTGCCAGCTGCGG‑GTATACAACTTAGCTTCCAGCACGTTCATCGGTCTACCGGGACCTAAAAACAATGCTTACATACTGGGTGAAGGCATGATGata                                                                                                                                                                                                                              <  2:391612/251‑1 (MQ=255)
         gATAATGCGTATCGTACAGCGCCTGCGGGAAGCCAC‑CAAAGTCATGAATCGTGGGCGGCGTCTCCATCGCGGTCACTCCTGTTCCACGGGTAAACCAGTGAGCCGAAACCACCACAATCGCTTGCGGGCGTGGCAATGTCATCCCCAACTTCTGCCAGCTGCGG‑GTATACAAATTATCTTCCAGCACGTTCATCGGACTACCGTGACCTAAAAACAATGCTGACATACGTGTTGAAGACATGATGATATcc                                                                                                                                                                                                                           >  1:93683/1‑251 (MQ=255)
            aatggtactagTACAGCGCCATCGGGAAGC‑ACACAAAGCCATTAATCGGTGGCGGTGTCTCCATCTCGGTCACTCCTTTTCCACGGGTAAACCAGTGAGCCGAAACCACCACAATCGCTTGCGGGCGTGTCAAGTTCATCCCCAACTTCTGCCAGCTGCGG‑GGATACAAATGATCTTACAGCACGTTCATCGGACTACCGTGACCTAAAAACAAACCTGACAGACGTGTTGAAGACATGATGATATCCTTa                                                                                                                                                                                                                        <  2:231916/241‑1 (MQ=255)
                   aTCGTACAGCGCCTGCGGGAAGCCAC‑CAAAGTCATGAATCGTGGGCGGCGTCTCCATCGCGGTCACTCCTGTTCCACGGGTAAACCAGTGAGCCGAAACCACCACAATCGCTTGCGGGCGTGGCAATGTCATCCCCAACTTCTGCCAGCTGCGG‑GTATACAAATTATCTTCCAGCACGTTCATCGGACTACCGTGACCTAAAAACAATGCTGACCTACGTGTTGAAGACATGATGATATCCTTCACTaa                                                                                                                                                                                                                   >  2:124683/1‑249 (MQ=255)
                     cGTACAGCGCCTGCGGGAAGCCAC‑CAAAGTCATGAATCGTGGGCGGCGTCTCCATCGCGGTCACTCCTGTTCCACGGGTAAACCAGTGAGCCGAAACCACCACAATCGCTTGCGGGCGTGGCAATGTCATCCCCAACTTCTGCCAGCTGCGG‑GTATACAAATTATCTTCCAGCACGTTCATCGGACTACCGTGACCTAAAAACAATGCTGACATACGTGTTGAAGACATGATGATATCCTTAACTAAAGgt                                                                                                                                                                                                               <  1:437553/251‑1 (MQ=255)
                      gTACAGCGCCTGCGGGAAGCCAC‑CAAAGTCATGAATCGTGGGCGGCGTCTCCATCGCGGTCACTCCTGTTCCACGGGTAAACCAGTGAGCCGAAACCACCACAATCGCTTGCGGGCGTGGCAATGTCATCCCCAACTTCTGCCAGCTGCGG‑GTATACAAATTATCTTCCAGCACGTTCATCGGACTACCGTGACCTAAAAACAATGCTGACATACGTGTTGAAGACATGATGATATCCTTAACTAAAGgtg                                                                                                                                                                                                              <  1:70261/251‑1 (MQ=255)
                       tACAGCGCCTGCGGGAAGCCAC‑CAAAGTCATGAATCGTGGGCGGCGTCTCCATCGCGGTCACTCCTGTTCCACGGGTAAACCAGTGAGCCGAAACCACCACAATCGCTTGCGGGCGTGGCAATGTCATCCCCAACTTCTGCCAGCTGCGG‑GTATACAAATTATCTTCCAGCACGTTCATCGGACTACCGTGACCTAAAAACAATGCTGACATACGTGTTGAAGACATGATGATATCCTTAACTAAAGgtgt                                                                                                                                                                                                             >  1:413403/1‑251 (MQ=255)
                        aCAGCGCCTGCGGGAAGCCAC‑CAAAGTCATGAATCGTGGGCGGCGTCTCCATCGCGGTCACTCCTGTTCCACGGGTAAACCAGTGAGCCGAAACCACCACAATCGCTTGCGGGCGTGGCAATGTCATCCCCAACTTCTGCCAGCTGCGG‑GTATACAAATTATCTTCCAGCACGTTCATCGGACTACCGTGACCTAAAAACAATGCTGACATACGTGTTGAAGACATGATGATATCCTTAACTAAAGgtgt                                                                                                                                                                                                             <  1:199955/250‑1 (MQ=255)
                         cAGCGCCTGCGGGAAGACAC‑AAAAGTCAGGAAGCGTGGGCGGCGTCTCCATCGCGGGCACTCCTGTTACACGGGTAAACCAGTGAGCCGAAACCACCACAATCGCGTGCGGGCGTGGCAATGTCATCCCCAACGTCTGCCAGCTGCGG‑GTATACAAATTAGCTTCCAGCACGTTCATCGGACTACCGTGATCTAAAAACAATGCTGACAGACGTGTTGAAGACATGATGATATCCTGAACTAAAGGTGTCa                                                                                                                                                                                                           <  2:84315/251‑1 (MQ=255)
                            cgcCTGCGGGAAGCCAC‑CAAAGTCATGAATCGTGGGCGGCGTCTCCACCGCGGTCACTCCTGATCCACGGGTACACCAGTGAGCCGAAACCACCACAATCGCTTGCGGGCGTGGCAATGTCATCCCCAACTACTGCCAGCTGCGG‑GTATACAACTTATCTTCCCGCACGTTCATCGGACTACCGTGACCTAAAAACAATGCTGACATACGTGTTGAAGACATGATGATATCCTTCACTAAAGGTGACAttt                                                                                                                                                                                                        >  2:287962/1‑251 (MQ=255)
                            cgcCTGCGGGAAGCAAC‑CAAAGGCAGGAATCGTGGGCGGCGCCTCCATCGCGGGCAGTCCTGTTCCACGGGTAAACCAGTGAGCCGAAACCACCACAATCGCTGGCGGGCGTGGCAATGTCATCCCCAACTTCTGCCAGCTGCGG‑GTATACAAATTAGCTTCCAGCACGTTCATCGGACTACCGTGACCTAAAAACAATGCTGACATACGTGTTGAAGACATGATGATATCCTTAACTAAATGTGTCAttt                                                                                                                                                                                                        <  2:191453/251‑1 (MQ=255)
                                   gggAAGCCAC‑CAAAGTCATGAATCGTGGGCGGCGTCTCCATCGCGGTCACTCCTGTTCCACGGGTAAACCAGTGAGCCGAAACCACCACAATCGCTTGCGGGCGTGGCAATGTCATCCCCAACTTCTGCCAGCTGCGG‑GTATACAAATTATCTTCCAGCACGTTCATCGGACTACCGTGACCTACAAACAATGCTGACATACGTGTTGAAGACATGATGATATCCTTAACTACAGGTGTCATTTTGATATc                                                                                                                                                                                                 >  2:91684/1‑251 (MQ=255)
                                                 aGTCATGAATCGTGGGCGGCGTCTCCATCGCGGTCACTCCTGTTCCACGGGTAAACCAGTGAGCCGAAACCACCACAATCGCTTGCGGGCGTGGCAAGGTCATCCCCAACTTCTGCCAGCTGCGG‑GTATACAAATTATCTTCCAGCACGTTCATCGGACTACCGTGACCTAAAAACAATGCTGACATACGTGTTGAAGACATGATGATATCCTTAACTAAAGGTGTCATTTTGATATCCTCACAATACGCt                                                                                                                                                                                    <  1:21285/251‑1 (MQ=255)
                                                  gTCATGAATCGTGGGCGGCGTCGCCATCGCGGTCACTCCTGTTCCACGGGTAAACCAGTGAGCCGAAGCCACCACAAGCGCTTGCGGGCGTGGCAATGTCATCCCCAACTTCTGCCAGCTGCTG‑GTATACAAATTATCTTCCAGCACGTGCATCTGACTACCGTGACGTAAAAACAACGGTGACATATGTGTTGAAGACATGATGATATGCTTAAGTAAAGGTGTCATTTTGagaccctgccagtacgctt                                                                                                                                                                                   <  2:21350/251‑19 (MQ=255)
                                                    cATGAAGCGTGGGCGGCGTCTCCATCGCGGTCACTCCTGTTCCACGGGTAAACCAGTGAGCCGAAACCACCACAATCGCTTGCGGGCGTGGCAATGTCATCCCCAACTTCTGCCAGCTGCGG‑GTATACAAATTAGCTTCCAGCACGTTCATCGGACTACCGTGACCTAAAAACAATGCTGACATACGTGTTTAAGACATGATGATATCCTTAACTAAAGTTGTCATTTTGATATGCTCACAATACGCTTg                                                                                                                                                                                  <  2:36811/250‑1 (MQ=255)
                                                       gAATCGTGGGCGGCGTCTCCATCGCTGTCACTCCTGTTCCACGGGTAATCCAGTGAGCCGAAACCACCACAATCGCTTGCGGGCGTGGCAATGTCATCCCCAACTTCTGCCAGCTGCGG‑GTATACAAATTATCTTCCAGCACGTGCATCGGACTACCGTGACCTAAAATCAAGGCTGACATTCGTGTTGAAGACATGATGATATCCTTAACTAAAGGTGTCATTTTGATATCCTCACAATACGCTGGTTcg                                                                                                                                                                              <  1:135498/251‑1 (MQ=255)
                                                        aaTCGAGGGCGGCGTCTCCATAGCGGTCACTCCTGTTCCACGGGTAAACCAGTGAGCCGAAACCACCACAATCGCTTGCGGGCGTGGCAATGTCATCCCCAACTTCTGCCAGCTGCGG‑GTATACAAATTATCTTCCAGCACGTTCATCGGACTACCGCGACCTAAAAACAATGCTGACATACGTGTTGAAGACATGATGATATCCTTAACTAAAGGTGTCATTTTGATATCCTCACAATACGCTTGTTcgg                                                                                                                                                                             >  2:60181/1‑251 (MQ=255)
                                                            gggggCGGCGTCTCCATCGCGGTCACTCCTGTTCCACGGGTAAACCAGTGAGCCGAAACCACCACAATCGCTTGCGGGCGTGGCAATGTCATCCCCAACTGCTGCCAGCTGCGG‑GTATACAAATTATCTTCCAGCACGTTCATCGGACTACCGTGACCTAAAAACAATGCTGTCATACGTGTTGAAGACATGATGATATCCTTAACTAAAGGTGTCATTTTGATATCCTCACAATACGCTTGTTCGGCGGa                                                                                                                                                                         <  2:190853/249‑1 (MQ=255)
                                                            gTGGGCGGCGTCTCAATCGCGGTCACTCCTGTTCCACGGGTAAACCAGTGAGCCGAAACCACCACAATCGCTTGCGGGCGTGGCAATGTCATCCCCAACTTCTGCCAGCTGCGG‑GTATACAAATTATCTTCCAGCACGTTCATCGGACTACCGTGACCTAAAAACAATGCTGACATACGTGTTGAAGACATGATGATATCCTTAACTAAAGGTGTCATTTTGATATCCTCACAATACGCTTGTTcggcgg                                                                                                                                                                          >  2:358883/1‑250 (MQ=255)
                                                               ggcggcGTCTCCATCGCGGTCACTCCTGTTCCAAGGGTCAACAAGCCAGCCGAACCCACCACAATCGCTTGCGGGCGTGGCCATGTAATCACCAACTTCTGCCAGCTGCGG‑GTATACAAATTATCTTCCAGCACGTTCATCGGACTACCGTGACCTAAAAACAATGCTGACATACGTGTTGAAGACATGATGATATCCTTAACTAAAGGAGTCATTTTGATATCCTCACAATACGCTTGTTCGGCGGAGt                                                                                                                                                                       >  1:412953/1‑250 (MQ=255)
                                                                    cgTCTCCATCGCGGTCACTCCTGTTCCACGGGTAAACCAGTGAGCCGAAACCACCACAATCGCTTGCGGGCGTGGCAATGTCATCCCCAACTTCTGCCAGCTGCGG‑GTATACAAATTATCTTCCAGCACGTTCATCGGACTACCGTGACCTAAAAACAATGCTGACATACGTGTTGAAGACATGATGATATCCTTAACTAAAGGTGTCATTTTGATATCCTCACAATACGCTTGTTCGGCGGAGTAAGAAc                                                                                                                                                                 <  1:150062/251‑1 (MQ=255)
                                                                         ccATCGCGGTCACTCCTGTTCCACTGGTAAACCAGTGATCCGAAACCACCACAAGCGCTTGCGGGCGTGGCAATGTCATCCCCAACTTCTGCCAGCTGCTG‑GTATACAAATTATCTTCCAGCACGTTCATCGGACTACCGTGACCTAAAAACAATGCTGACATACGTTTTGAAGGCATTATGATATCCTTAACTAAAGGTGTCATTTTGAGAGCCTCATAATACGCTTGTTCGGCTGAGTAATAACCCGGa                                                                                                                                                            <  2:49574/251‑1 (MQ=255)
                                                                             cgcgGTCACTCCTGTTCCACGGGTAAACCAGTGAGCCGAAACCACCACAATCGCTTGCGGGCGTGGCAATGTCATCCCCAACTTCTGCCAGCTGCGG‑GTATACAAATTATCTTCCAGCACGTTCATCGGACTACCGTGACCTAAAAACAATGCTGACATACGTGTTGAAGACATGATGATATCCTTAACTAAAGGTGTCATTTTGATATCCTCACAATACGCTTGTTCGGCGGAGTAAGAACCCGGATAAc                                                                                                                                                        <  1:72580/251‑1 (MQ=255)
                                                                                        ctgatCCACGGGTAAACCAGTGAGCCGAAACCACCACAAACGCTTGCGGGCGTGGCAATGTCATCCCCAACTTCTGCCAGCTGCGG‑GTATACAAATTATCTTCCAGCACGTTCATCGGACTACCGTGACCTAAAAACAATGCTGACATACGTGTTGAAGACATGATGATATCCTTAACTAAAGGTGTCATTTTGATATCCTCACAATACGCTTGTTCGGCGGAGTAAGAACACGGATAACAattctgatg                                                                                                                                              >  2:259007/5‑243 (MQ=255)
                                                                                            tCCACGGATAAACCAGTGAGCCGAAACCACCAAAATCGCTTGCGGGCCTGGCAATGTCATCCCCACCTTATGCCAGCTGCGG‑GTATACAAATTATCTCCAAGCCCGTTAATCGCACTACCGTGACATAAAACCAATGCTGACATAAGTGTTAAAGAAATGCTGATATCATTAACTAAAGGTGTCATTTTGATATCCTCACAATACGCTTGTTCGGCGGAGTAAGAACCCGGATAACAATGATGATGttcat                                                                                                                                         >  1:379810/1‑251 (MQ=255)
                                                                                                    tAAACCCGTGAGCCGAAACCACCACAATCGCTTGCGGGCGTGGCAATGTCATCCCCAACTTCTGCCAGCTGCGG‑GTATACAAATTATCTTCCAGCACGTTCATCGGACTACCGTGACCTAAAAACAATGCTGACATACGTGTTGAAGACATGATGATATCCTTAACTAAAGGTGTCATTTTGATATCCTCACAATACGCTTGTTCGGCGGAGTAAGAACCCGGATAACAATGATGATGATCATCAGTTAtt                                                                                                                                 >  1:299981/1‑251 (MQ=255)
                                                                                                        ccattgtgCCGAAACCACCACAAGCGCTTGCGGGCGTTGCAATTTCAGCCCCAACTGCATCCAGCTGTGG‑GTATACAAATTATCTTCCAGCACGTTCATCGGACTACCGTGACCTAAAAACAATGCGGACATACGGGTTGAAGACATGATGGTGTTCTGTAGTAAAGGTGTCATGGTGAGATCCTCACAATACGCTTGTGCGTCGGAGTAAGAACCCGGTTAAGAATGATGATGATCATCAGTTATTTTGa                                                                                                                             <  2:378307/244‑1 (MQ=255)
                                                                                                           gTGAGCCGAAACCACCACAATCGCTTGCGGGCGTGGCAATGTCATCCCCAACTTCTGCCAGCTGCGG‑GTATACAAATTATCTTCCAGCACGTTCATCGGACTACCGTGACCTAAAAACAATGCTGACATACGTGTTGAAGACATGATGATATCCTTAACTAAAGGTGTCATTTTGATATCCTCACAATACGCTTGTTCGGCGGAGTAAGAACCCGGATAACAATGATGATGATCATCAGTTATTTTGACGa                                                                                                                          <  1:446437/251‑1 (MQ=255)
                                                                                                             gAGCCGAAACCACCACAATCGCTTGCGGGCGTGGCAATGTCATCCCCAACGTCTGCCAGCTGCGG‑GTATACAAATTATCTTCCAGCACGTTCATCGGACTACCGTGACCTAAAAACAATGCTGACATACGTGTTGAAGACATGATGATATCCTTAACTAAAGGTGTCATTTTGATATCCTCACAATACGCTTGTTCGGCGGAGTAAGAACCCGGATAACAATGATGATGATCATCAGTTATTTTGACGAt                                                                                                                         <  2:24281/250‑1 (MQ=255)
                                                                                                             gAGCCGAAACCACCACAATCGCGTGCGGGCGTGTCAATGTCATCCCCAACTTCTGCCAGCTGCGG‑GTATACAAATTATCTTCCAGCTCGTTCATCGGACTACCGTGACCTAAAAACAATGCTGACATACGTGTTGAAGACATGATGATATCCTTAACTAAAGGTGGCATTTTGATATCCTCACAATACGCTTGTTCGGCGTAGTAAGAACCCGGATAACAATGATGATGATCAGCAGTTATTTTGACGAt                                                                                                                         <  2:149778/250‑1 (MQ=255)
                                                                                                                   aaaCCACCACAATCGCTTGCGGGCGTGGCAATGTCATCCCCAACTTCTGCCAGCTGCGG‑GTATACAAATTATCTTCCAGCACGTTCATCGGACTACCGTGACCTAAAAACAATGCTGACATACGTGTTGAAGACATGATGATATCCTTAACTAAAGGTGTCATTTTGATATCCTCACAATACGCTTGTTCGGCGGAGTAAGAACCCGGATAACAATGATGATGATCATCAGTTATTTTGACGATCTGCCTg                                                                                                                  >  2:424295/1‑251 (MQ=255)
                                                                                                                                             ggCAATGTCATCCCCAACTTCTGCCAGCTGCGG‑GTATACAAATTATCTTCCAGCACGTTCATCGGACTACCGTGACCTAAAAACAATGCTGACATACGTGTTGAAGACATGATGATATCCTTACCTAAAGGTGTCATTTTGATATCCTCACAATACGCTTGTTCGGCGGAGTAAGAACCCGGATAACAATGATGATGATCATCAGTTATTTTGACGATCTGCCTGAAGGTGAAGATTTATAAGGAGTTGTc                                                                                        >  2:170380/1‑251 (MQ=255)
                                                                                                                                               cAATGTCATCCCCAACTTCTGCCAGCTGCGG‑GTATACAAATTATCTTCCAGCACGTTCATCGGACTACCGTGACCTAAAAACAATGCTGACATACGTGTTGAAGACATGATGATATCCTTAACTAAAGGTGTCATTTTGATATCCTCACAATACGCTTGTTCGGCGGAGTAAGAACCCGGATAACAATGATGATGATCATCAGTTATTTTGACGATCTGCCTGAAGGTGAAGATTTATAAGGATATGTCGa                                                                                      <  1:124683/251‑1 (MQ=255)
                                                                                                                                                aaTGTCATCCCCAACTTCTGCCAGCTGCGG‑GTATACAAATTATCTTCCAGCACGTTCATCGGACTACCGTGACCTAAAAACAATGCTGACATACGTGTTGAAGACATGATGATATCCTTAACTAAAGGTGTCATTTTGATATCCTCTCAATACGCTTGTTCGGCGGAGGAAGAACCTGGATAACAAGGATTATTATCAGCAGTTATTTTTACGATCTGCCTGAAGGTGAAGATTTATAAGGAGTTGTCGAt                                                                                     <  2:207366/251‑1 (MQ=255)
                                                                                                                                                        ccccAACTTCTGCCAGCT‑CGGAGTATACAAATTATCTTCAAGCACGATCATCGGACTAACGTGACCTAAAAACAATGCTGACATACGTGATGAAGACATGATGATATCCTTAACTAAAGGAGTCATTTTGATATCCTCACAATACGCTTGTTCGGCGGCGTAAGAACCCGGATAACAATGATGATGATCATCAGATATTTTGACGATCTGCCTGAAGGAGAAGATTAATAAGGAGTTGTCGATGTCAGTAc                                                                             >  2:388408/1‑251 (MQ=255)
                                                                                                                                                            aaCTTCTGCCAGCTGCGG‑GTATACAAATTATCTTCCAGCACGTTCATCGGACTACCGTGACCTAAAAACAATGCTGACATACGTGTTGAAGACATGATGATATCCTTAACTAAAGGTGTCATTTTGATATCCTCACAATACGCTTGTTCGGCGGAGTAAGAACCCGGATAACAATGATGATGATCATCAGTTATTTTGACGATCTGCCTGAAGGTGAAGATTTATAAGGAGTTGTCGATGTCAGTACctct                                                                         >  2:89084/1‑251 (MQ=255)
                                                                                                                                                                    ccAGCTGCGG‑GAATAAAAATTATCTTCAAGCCCGCTCAACGGCCTACCGTGACCTCAAAACAATGCTGACATACGTGTTGAAGACATGGTGATATCCTTAACTAAAGGTGTCATTTTGATATCCTCACAATACGCTTGTTCGGCGAAGTAAGAAAACGGCTAACAATGAAGATGATCATCAGTTATTTTGCAGATCTGCCTGAAGGTGAAGATTTATAAGGAGTTGTCTATGTAAGTACCTCTCATTCTg                                                                  >  2:419232/1‑250 (MQ=255)
                                                                                                                                                                        cTGCGG‑GCATACAACTTATCTTCCCGCACGCTCATCGGACTACCGTGACCTAAAAACAATGCTGACATACGCGTTGAAGACATGATGATATCCTTAACTACAGGTGTCATTTTGATATCCTCCCAATACGCTTGTTCGGCGGAGTAAGAACCCGGATAACAATGATGATGATCATCAGTTATTTTGACGATCTGCCTGAAGGTGAAGATTTATAAGGAGTTGTCGATGTCAGTACCTCTCATTCTGACCAt                                                             >  1:164774/1‑251 (MQ=255)
                                                                                                                                                                          gCGG‑GTATACAAATTATCTTCCAGCACGTTCATCGGACTACCGTGACCTAAAAACAATGCTGACATACGTGTTGAAGACATGATGATATCCTTAACTAAAGGTGTCATTTTGATATCCTCACAATACGCTTGTTCGGCGGAGTAAGAACCCGGATAACAATGATGATGATCATCAGTTATTTTGACGATCTGCCTGAAGGTGAAGATTTATAAGGAGTTGTCGATGTCAGTACCTCTCATTCTGACCATa                                                            <  2:59015/250‑1 (MQ=255)
                                                                                                                                                                                  taCAACTTATCTTCCAGCACGCTAATCGGACTACCGTGACCTAAAAACAATGCTGACATACGTGTTGAAGACATGATGATATCCTTAACTAAAGGTGTCATTTTGATATCCTCACAATACGCTTGTTCGGCGGAGTAAGAACCCGGATAACAATGATGACGATCATCAGTTATTTTGACGATCTGCCTGAAGGTGAAGATTTATAAGGAGTTGTCGATGTCCGTACCTCTCATTCTGACCATACTGGCggg                                                    >  2:22500/1‑251 (MQ=255)
                                                                                                                                                                                       aTTATCTTCCAGCACGTTCATCGGACTACCGTGACCTAAAAACAATGCTGACATACGTGTTGAAGACATGATGATATCCTTAACTAAAGGAGTCATTTTGATATCCTCACAATACGCTTGTTCGGCGGAGTAAGAACCCGGATAACAATGATGATGATCATCAGTTATTTAGACGATCTGCCTGAAGGTGAAGATTTATAAGGAGTTGTCGATGTCAGAACCTCTCATTCTGACCATACTGGCgggggg                                                 >  2:31725/1‑248 (MQ=255)
                                                                                                                                                                                           tCTTCAAGCACGTTCATCGGACTCCCGTGACCTAAAAAAAATGCTGAAATCCGAGTTGAAGACCTCATGATATCCTTACCTCACGGCGTAAATTTGATATCCTCACAATACGCTTGTTCGGCGGAGCAAGAACCCGGATAACAATGATGATGATCATCAGTTATTTTGACGATCTGCCTGAAGGTAAAGATTTATAAGGAGTTGTCGATGTCAATACCTCTCATTCTGACCATCCTGGCGTGGGGAGcgcc                                           >  2:409602/1‑248 (MQ=255)
                                                                                                                                                                                                 aGCACGTTCATCGGACTACCGGGACCTAAAAACAATGCTGACATACGTGTTGTCGTCATGATGATATCCTTAACTAAAGGTGTCATCTTGATATTCTCAAAATACGCTTGATCGGCGGAGAAAGAACCCGATTAACAATGATGATGATCATCAGTTATTTTAACAATCTGCTCGAAGGTTAAGATTTATAAGGAGTTGTCGAGGTCCGTACCTCTCATTCTGACCATACTAGAGGGGGCactcaacatttt                                     >  2:318503/1‑240 (MQ=255)
                                                                                                                                                                                                       ttCATCGGACTACCGTGACCTAAAAACAATGCTGACATACGTGTTGAAGACATGATGATATCCTTAACTAAAGGTGTCATTTTGATATCCTCACAATACGCTTGTTCGGCGGAGTAAGAACCCGGATAACAATGATGATGATCATCAGTTATTTTGACGATCTGCCTGAAGGTGAAGATTTATAAGGAGTTGTCGATGTCAGTACCTCTCATTCTGACCATACTGGCGGGGGCGGCCACGTGTTTTGgcgc                               >  1:350741/1‑251 (MQ=255)
                                                                                                                                                                                                       ttCATAGGACTACCGTGACCTAAAAACAATGCTGACATACGTGTTGAAGACATGATGATATCCTTAACTAACGGTGTCATTTTGATATCCTCACAATACGCTTGTTCGGCGGAGTAAGAACCCGGCTAACAATGATGATGATCATCAGTTATTTTGAAGATCTGACTGAAGGTGAAGATTTATAAGGAGTTGTCGATGTCAGTACCTCTCATTCTGAACATACTGGCGGGGGGAGcaccggttattttggg                               >  2:377578/1‑236 (MQ=255)
                                                                                                                                                                                                            cGGACTACCGTGACCTAAACACAATGCTGACATACGCGTTGAAGACATGATGATATACTTAACTAAAGGTGTCATTTTGATATCCTCACAATACGCTTGTTCGGCGGAGTAAGAACCCGGATAACAATGATGATGATCATCAGTTATTTTGACGATCTGCCTGAAGGTGAAGATTTATAAGGAGTTGTCGATGTCAGTACCTCTCATTCTGACCATACTGGCgggggagcgcacggttttttgggcgattc                          >  2:225343/1‑227 (MQ=255)
                                                                                                                                                                                                                  aCCGTGACCTAAAAACAATGCTGACATACGTGTTGAAGACATGATGATATCCTTAACTAAAGGTGTCATTTTGATATCCTCACAATACGCTTGTTCGGCGGAGTAAGAACCCGGATAACAATGATGATGATCATCAGTTATTTTGACGATCTGCCTGAAGGTGAAGATTTATAAGGAGTTGTCGATGTCAGTACCTCTCATTCTGACCATACTGGCGGGGGCAGCCACGTTTATTGGCGCGTTTCTCGGCg                    <  1:235525/251‑1 (MQ=255)
                                                                                                                                                                                                                     gTGACCTAAAAACAATGCTGACATACGTGTTGAAGACATGATGATATCCTTAACTAAAGGTGTCATTTTGATATCCTCACAATACTCTTGTTCGGCGGAGTAAGAACCCGGATAACAATAATGATGATCATCAGTTATTTTGACGATCTGCCTGAAGGTGAAGATTTATAAGGAGTTGTCGATGTCAGTACCTCTCATTCTGACCATACTGGCGGGGGCAGCCACGTTTATTGGCGCGTTTCTCGGCGTtc                 <  2:21875/251‑1 (MQ=255)
                                                                                                                                                                                                                             aaaaCAATTCTGACATACGTGGTGAAGACATGATGATATCCTGAACTAAAGGTGTCATTTTGATATCCTCACAATACGCGTGTTCGGCGGAGTAAGAACCCGGATAACAATGATGATGATCATCAGTTATTTTGACGATCTGCCTGAAGGTGAAGAGGTATAAGTAGTTGTCGATGTCAGTACCTCTCATTCTTACCATACTGGCGGGGGCAGCCACGTGTATTGGCGCGTTTCTCGGCGTTCTCGGGCaa         <  2:356676/251‑1 (MQ=255)
                                                                                                                                                                                                                                 cAATGCTGACATACGTGTTGAAGACATGATGATATCCTTAACTAAAGGTGTCATTTTGATATCCTCACAATACGCTTGTTCGGCGGAGTAAGAACCCGGATAACAATGATGATGATCATCAGTTATTTTGACGATCTGCCTGAAGGTGAAGATTTATAAGGAGTTGTCGATGTCAGTACCTCTCATTCTGACCATACTGGCGGGGGCAGCCACGTTTATTGGCGCGTTTCTCGGCGTTCTCGGGCAAAAAc     <  1:433442/251‑1 (MQ=255)
                                                                                                                                                                                                                                    tGCTGACATACGTGTTGAAGACATGATGATATCCTTAACTAAAGGTGTCATTTTGATATCCTCACAATACGCTTGTTCGGCGGAGTAAGAACCCGGATAACAATGATGATGATCATCAGTTATTTTGACGATCTGCCTGAAGGTGAAGATTTATAAGGAGTTGTCGATGTCAGTACCTCTCATTCTGACCATACTGGCGGGGGCCGCCACGTTTATTGGCGCGTTTTTCGGCGTTCTCGGGCAAAAAcccc  >  1:314239/1‑250 (MQ=255)
                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                     
CCGGAGCAGGATAATGCGTATCGTACAGCGCCTGCGGGAAGCCAC‑CAAAGTCATGAATCGTGGGCGGCGTCTCCATCGCGGTCACTCCTGTTCCACGGGTAAACCAGTGAGCCGAAACCACCACAATCGCTTGCGGGCGTGGCAATGTCATCCCCAACTTCTGCCAGCTGCGG‑GTATACAAATTATCTTCCAGCACGTTCATCGGACTACCGTGACCTAAAAACAATGCTGGCATACGTGTTGAAGACATGATGATATCCTTAACTAAAGGTGTCATTTTGATATCCTCACAATACGCTTGTTCGGCGGAGTAAGAACCCGGATAACAATGATGATGATCATCAGTTATTTTGACGATCTGCCTGAAGGTGAAGATTTATAAGGAGTTGTCGATGTCAGTACCTCTCATTCTGACCATACTGGCGGGGGCAGCCACGTTTATTGGCGCGTTTCTCGGCGTTCTCGGGCAAAAACCCT  >  NC_000913/3182158‑3182634

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑