Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 770,805 G→A 11.1% intergenic (+194/‑653) mngB → / → cydA alpha‑mannosidase/cytochrome d terminal oxidase, subunit I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913770,8050GA11.1% 148.6 / 9.0 54intergenic (+194/‑653)mngB/cydAalpha‑mannosidase/cytochrome d terminal oxidase, subunit I
Reads supporting (aligned to +/- strand):  ref base G (19/29);  new base A (4/2);  total (23/31)
Fisher's exact test for biased strand distribution p-value = 3.84e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.06e-01

AGGAAAATCAAGGTTCAAATCTATCGGGGCCTTTTTTACCCGGCCAGTCACGGACGTTCAGTTACCGGCTTGCCTGAATAGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAGCGCATCAGGTAATTTTGCGTTTATCTTCACTCTCAAGCCACGTATATGTGGCTTTATTTTTAACAAAATAATAACCCTGGGTGAGTTAATTATAATATAATTATAAGTTAACTAAATGTTAATATTGGCGGGGTGGATTTATGCCTTTATTAGTAATCCTGAAACTCTGCGTCGTATTAGCCAG  >  NC_000913/770535‑770947
                                                                                                                                                                                                                                                                              |                                                                                                                                              
agaaaaaTAAAGGTTTAAAACTATCGGGTCCTTGTATACCCGGACAGTCACGGACGTTCAGTTAGCGGCTTGCCGTAATAGAAATCAAAACGAAGCCACATATGCGCGGCCAGAGTGTGGACAAAGGGCGCTTTGTTCATGCCGGATACGGCGTGAACGGTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAGCGCATCAGGTAATTTTGCGTTTATCTTCACTCTCATGCCACGtat                                                                                                                                             <  4:132215/271‑1 (MQ=255)
aGGAAAATCAAGGTTCAAATCTATCGGGGCCTTTTTTACCCGGCCAGTCACGGACGTTCAGTTACCGGCTTGCCTGAATAGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAGCGCATCAGGTAATTTTGCGTTTATCTTCACTCTCAAGCCACGtata                                                                                                                                            >  2:593038/1‑275 (MQ=255)
aGGAAAATCAAGGTTCAAATCTATCGGGGCCTTTTTTACCCGGCCAGTCACGGACGTTCAGTTACCGGCTTGCCTGAATAGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAGCGCATCAGGTAATTTTGCGTTTATCTTCACTCTCAAGCCACGtata                                                                                                                                            >  2:602782/1‑275 (MQ=255)
  gAAAATCAAGGTTCAAAGCTATCGGGGCCTTTTTTACCCGGCCAGTCACGGACGTTCAGTTACCGGCTTGCCTGAATAGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAGCGCATCAGGTAATTTTGCGTTTATCTTCACTCTCAAGCCACGTATAtg                                                                                                                                          <  2:217383/275‑1 (MQ=255)
       tctaGGTTCAAATCTATCGGTGCCTTTTTTACCCGGCCCGTACCGGACGTTCAGTTACCTGCTTACCTTAATCGCAATCAAACCGAATCCACCTATGCCCGGCGTGAATTTCGCCAAAGGGCTCTTTGTGAAGGCCGGATACGGCATGAACTCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATAATGTAGACACTGGCAAGCGTAGCGCATCAGGTAATTTTGCGTTTATATTCACTCTCAAGCCACGTATATGTGGCt                                                                                                                                     <  4:127692/272‑1 (MQ=255)
                    cTAGCGGGGCCTGTTTTAACCGGCCAGTCACGGACGTTACGTTACCGGCTTGCCTGAATAGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAGCGCATCAGGTAATTTTGCGTTTATCTTCACTCTCAAGCCACGTATATGTGGCTTTATTTTTAACaa                                                                                                                        <  4:681465/275‑1 (MQ=255)
                                  tttACCCGGCCAGGGACGGACGTTCAGTTACCGGCTGGCCTGAATAGCAAACGAACCGAAGCCACAAATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTGCATGCCGGATACGGCATGAACGATTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAGCGCATCAGGTAATTTTGCGTTTATCTTCACTCTCAAGCCACGTATATGTGGCTTTATTTTTAACAAAATAATAACCCTgg                                                                                                          <  4:633782/275‑1 (MQ=255)
                                    tACCCGGCCAGTACCGGATGTTCCGTTACCGGCTTGCCTGAATAGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCCATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAGCGCATCAGGTAATTTTGCGTTTATCTTCACTCTCAAGCCACGTATATGTGGCTTTATTTTTAACAAAATAATAACCCTGGGt                                                                                                        <  4:602780/275‑1 (MQ=255)
                                          gCCAGTCACGGACGTTCAGTTACCGGCTTGCCTGAATAGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATAAATTGCAGAGATGTTTTAGACACTTGCAAGCGTAGCGCATCAGGTAATTTTGCGTTTATCCTCACTCTCAAGCCACGTATATGTGGCTTTATTTTTAACCAAATAATAAcc                                                                                                              >  2:633784/1‑263 (MQ=255)
                                                                    ctcgCCTGAATAGAAATCAAACCTAAGCCACATAAGCGCGGTCAGATTTTTGACAAAGGGTGCTTTGTTCATGCCGGATACGGCATGAACGCTTGATTCTGTCTACAAAAGCAGGCACATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAGCGCATCAGGTAATTTTGCGTTTATCTTCACTCTCAAGCCACGTATATGTGGCTTTATTTTTAACAAAATAATAACCCTGGGTGAGTTAATTATGATATAATTATAAGTTAACTa                                                                        <  4:342013/272‑1 (MQ=255)
                                                                                       tccnnagtccACATATTGTCTGCCAGATTTTGGACAAAGGGCGCTTTGTACATNCCAGATACGGCATGTACGATTTATTCGGTCTACAAAAGAAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAGCGCATCAGGTAATTTTGCGTTTATCTTCACTCTCAAGCCACGTAGATGTGGCTTTATTTTTAAc                                                                                                                          <  4:726081/198‑1 (MQ=255)
                                                                                                cacaTATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAGCGCATCAGGTAATTTTGCGTTTATCTTCACTCTCAAGCCACGTATATGTGGCTTTATTTTTAACAAAATAATAACCCTGGGTGAGTTAATTATAATATAATTATAAGTTAACTAAATGTTAATATTGGCGGGGTGGATTTAt                                            <  4:593036/275‑1 (MQ=255)
                                                                                                                       gACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAGCGCATCAGGTAATTTTGCGTTTATCTTCACTCTCAAGCCACGTATATGTGGCTTTATTTTTAACAAaataat                                                                                                                  >  4:389406/1‑182 (MQ=255)
                                                                                                                           aaGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAGCGCATCAGGTAATTTTGCGTTTATCTTCACTCTCAAGCCACGTATATGTGGCTTTATTTTTAACAAAATAATAAcc                                                                                                              <  2:389408/182‑1 (MQ=255)
                                                                                                                                   tttGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAGCGCATCAGGTAATTTTGCGTTTATCTTCACTCTCAAGCCACGTATATGTgg                                                                                                                                       <  8:664059/149‑1 (MQ=255)
                                                                                                                                    ttGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAGCGCATCAGGTAATTTTGCGTTTATCTTCACTCTCAAGCCACGTATATGTgg                                                                                                                                       >  6:800513/1‑148 (MQ=255)
                                                                                                                                         cATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAGCGCATCAGGTAATTTTGCGTTTATCTTCACTCTCAAGCCACGTATATGTGGCTTTAt                                                                                                                                 <  7:50849/149‑1 (MQ=255)
                                                                                                                                           tGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAGCGCATCAGGTAATTTTGCGTTTATCTTCACTCTCAAGCCACGTATATGTGGCTTTAtt                                                                                                                                >  5:171512/1‑148 (MQ=255)
                                                                                                                                           tGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAGCGCATCAGGTAATTTTGCGTTTATCTTCACTCTCAAGCCACGTATATGTGGCTTTAtt                                                                                                                                <  6:532805/148‑1 (MQ=255)
                                                                                                                                            gCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAGCGCATCAGGTAATTTTGCGTTTATCTTCACTCTCAAGCCACGTATATGTGGCTTTAttt                                                                                                                               >  5:415028/1‑148 (MQ=255)
                                                                                                                                               ggATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAGCGCATCAGGTAATTTTGCGTTTATCTTCACTCTCAAGCCACGTATATGTGGCTTTATTTTTa                                                                                                                            >  6:770549/1‑148 (MQ=255)
                                                                                                                                                 aTACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAGCGCATCAGGTAATTTTGCGTTTATCTTCACTCTCAAGCCACGTATATGTGGCTTTATTTTTAAc                                                                                                                          >  5:126611/1‑148 (MQ=255)
                                                                                                                                                               cTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAGCGCATCAGGTAATTTTGCGTTTATCTTCACTCTCAAGCCACGTATAtgtg                                                                                                                                        >  5:455772/1‑120 (MQ=255)
                                                                                                                                                               cTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAGCGCATCAGGTAATTTTGCGTTTATCTTCACTCTCAAGCCACGTATAtgtg                                                                                                                                        <  6:455772/120‑1 (MQ=255)
                                                                                                                                                               cTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAGCGCATCAGGTAATTTTGCGTTTATCTTCACTCTCAAGCCACGTATATGTGGCTTTATTTTTAACAAAATAATAACCCt                                                                                                            <  6:171512/148‑1 (MQ=255)
                                                                                                                                                               cTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATCGCAGAGATGATGTAGACACTGGCAAGCGTAGCGCATCAGGTAATTTTGCGTTTATCTTCACTCTCAAGCCACATATACGTGGCt                                                                                                                                     >  6:304502/1‑123 (MQ=255)
                                                                                                                                                               cTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATCGCAGAGATGATGTAGACACTGGCAAGCGTAGCGCATCAGGTAATTTTGCGTTTATCTTCACTCTCAAGCCACATATACGTGGCt                                                                                                                                     <  5:304502/123‑1 (MQ=255)
                                                                                                                                                                tttATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAGCGCATCAGGTAATTTTGCGTTTATCTTCACTCTCAAGCCACGTATATGTGGCTTTATTTTTAACaaa                                                                                                                       >  1:184810/1‑136 (MQ=255)
                                                                                                                                                                    ttCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAGCGCATCAGGTAATTTTGCGTTTATCTTCACTCTCAAGCCACGTATATGTGGCTTTATTTTTAACAAaataa                                                                                                                   <  3:184809/136‑1 (MQ=255)
                                                                                                                                                                    ttCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAGCGCATCAGGTAATTTTGCGTTTATCTTCACTCTCAAGCCACGTATATGTGGCTTTATTTTTAACAAAATAATAACCCTGGGTGa                                                                                                      <  8:749540/149‑1 (MQ=255)
                                                                                                                                                                      cGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAGCGCATCAGGTAATTTTGCGTTTATCTTCACTCTCAAGCCACGTATATGTGGCTTTATTTTTAACAAAATAATAACCCTGGGTGAg                                                                                                     <  5:800513/148‑1 (MQ=255)
                                                                                                                                                                         tCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAGCGCATCAGGTAATTTTGCGTTTATCTTCACTCTCAAGCCACGTATATGTGGCTTTATTTTTAACAAAATAATAACCCTGGGTGAGTTa                                                                                                  <  6:393726/148‑1 (MQ=255)
                                                                                                                                                                          cTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAGCGCATCAGGTAATTTTGCGTTTATCTTCACTCTCAAGCCACGTATATGTGGCTTTATTTTTAACAAAATAATAACCCt                                                                                                            >  2:215796/1‑137 (MQ=255)
                                                                                                                                                                           tACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAGCGCATCAGGTAATTTTGCGTTTATCTTCACTCTCAAGCCACATATACGTGGCttga                                                                                                                                  >  7:124040/1‑112 (MQ=255)
                                                                                                                                                                           tACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAGCGCATCAGGTAATTTTGCGTTTATCTTCACTCTCAAGCCACATATACGTGGCttga                                                                                                                                  <  8:124040/114‑3 (MQ=255)
                                                                                                                                                                            aCAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAGCGCATCAGGTAATTTTGCGTTTATCTTCACTCTCAAGCCACGTATATGTGGCTTTATTTTTAACAAAATAATAACCCTGGGTGAGTTAAtt                                                                                               <  5:770549/148‑1 (MQ=255)
                                                                                                                                                                             cAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAGCGCATCAGGTAATTTTGCGTTTATCTTCACTCTCAAGCCACGTATATGTGGCTTTATTTTTAACAAAAtaataa                                                                                                                >  5:313375/1‑130 (MQ=255)
                                                                                                                                                                             cAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAGCGCATCAGGTAATTTTGCGTTTATCTTCACTCTCAAGCCACGTATATGTGGCTTTATTTTTAACAAAAtaataa                                                                                                                <  6:313375/130‑1 (MQ=255)
                                                                                                                                                                              aaaaGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAGCGCATCAGGTAATTTTGCGTTTATCTTCACTCTCAAGCCACGTATATGTGGCTTTATTTTTAACAAAATAATAACCCTGGGt                                                                                                        <  4:215795/137‑1 (MQ=255)
                                                                                                                                                                                     ggCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAGCGCATCAGGTAATTTTGCGTTTATCTTCACTCTCAAGCCACGTATATGTGGCTTTATTTTTAACAAAATAATAACCCTGGGTGAg                                                                                                     >  4:783336/1‑133 (MQ=255)
                                                                                                                                                                                       cAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAGCGCATCAGGTAATTTTGCGTTTATCTTCACTCTCAAGCCACGTATATGTGGCTTTATTTTTAACAAAATAATAACCCTGGGTGAGTTAATTATAATATAAtt                                                                                    >  6:7013/1‑148 (MQ=255)
                                                                                                                                                                                         aaTTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAGCGCATCAGGTAATTTTGCGTTTATCTTCACTCTCAAGCCACGTATATGTGGCTTTATTTTTAACAAAATAATAACCCTGGGTGAGTTaa                                                                                                 <  2:783338/133‑1 (MQ=255)
                                                                                                                                                                                                     tGCAGAGATGATGTAGACACTGGCAAGCGTAGCGCATCAGGTAATTTTGCGTTTATCTTCACTCTCAAGCCACGTATATGTGGCTTTATTTTTa                                                                                                                            >  4:354746/1‑94 (MQ=255)
                                                                                                                                                                                                     tGCAGAGATGATGTAGACACTGGCAAGCGTAGCGCATCAGGTAATTTTGCGTTTATCTTCACTCTCAAGCCACGTATATGTGGCTTTATTTTTAACaaaa                                                                                                                      >  4:20699/1‑100 (MQ=255)
                                                                                                                                                                                                         gagaTGATGTAGACACTGGCAAGCGTAGCGCATCAGGTAATTTTGCGTTTATCTTCACTCTCAAGCCACGTATATGTGGCTTTATTTTTAACaa                                                                                                                        <  2:354748/94‑1 (MQ=255)
                                                                                                                                                                                                         gagaTGATGTAGACACTGGCAAGCGTAGCGCATCAGGTAATTTTGCGTTTATCTTCACTCTCAAGCCACGTATATGTGGCTTTATTTTTAACAAaataat                                                                                                                  <  2:20699/100‑1 (MQ=255)
                                                                                                                                                                                                           gatgatGTAGACACTGGCAAGAGTAGCGCATCAGGTAATTTTGCGTTTATCTTCACTCTCAAGCCACGTATATGTGGCTTTATTTTTAACAAAATAATAACCCTGGGTGAGTTAATTATAATATAATTATAAGTTAACTAAATGTTaa                                                                <  6:415028/148‑1 (MQ=255)
                                                                                                                                                                                                                      cacTGGCAAGCGTAGCGCATCAGGTAATTTTGCGTTTATCTTCACTCTCAAGCCACGTATATGTGGCTTTATTTTTAACAAAATAATAACCCTGGGTGAGTTAATTATAATATAATTATAAGTTAACTAAATGTTAATATTGGCgggg                                                     <  6:126611/148‑1 (MQ=255)
                                                                                                                                                                                                                      cacTGGCAAGCGTAGCGCATCAGGTAATTTTGCGTTTATCTTCACTCTCAAGCCACGTATATGTGGCTTTATTTTTAACAAAATAATAACCCTGGGTGAGTTAATTATAATATAATTATAAGTTAACTAAATGTTAATATTGGCGGGGt                                                    >  7:663513/1‑149 (MQ=255)
                                                                                                                                                                                                                           gCAAGCGTAGCGCATCAGGTAATTTTGCGTTTATCTTCACTCTCAAGCCACGTATATGTGGCTTTATTTTTAACAAAATAATAACCCTGGGTGAGTTAATTATAATATAATTATAAGTTAACTAAATGTTAATATTGGCGGGGTGGAtt                                               <  8:663513/149‑1 (MQ=255)
                                                                                                                                                                                                                                                          tATCTTCACTCTCAAGCCACGTATATGTGGCTTTATTTTTAACAAAATAATAACCCTGGGTGAGTTAATTATAATATAATTATAAGTTAACTAAATGTTAATATTGGCGGGGTGGATTTATGCCTTTATTAGTAATCCTGAAACTCTg                 <  5:7013/148‑1 (MQ=255)
                                                                                                                                                                                                                                                                    ctcAAGCCACGTATATGTGGCTTTATTTTTAACAAAATAATAACCCTGGGTGAGTTAATTATAATATAATTATAAGTTAACTAAATGTTAATATTGGCGGGGTGGATTTATGCCTTTATTAGTAATCCTGAAACTCTGCGTCGTATTa       >  6:598050/1‑148 (MQ=255)
                                                                                                                                                                                                                                                                     taaaGCCACATATACGTGGCTTTATTTTTAACAAAATAATAACCCTGGGTGAGTTAATTATAATATAATTATAAGTTAACTAAATGTTAATATTGGCGGGGTGGATTTATGCCTTTATTAGTAATCCTTAAACTCTGCGTCGTATTAGc     >  8:238444/3‑149 (MQ=255)
                                                                                                                                                                                                                                                                         gCCACATATACGTGGCTTTATTTTTAACAAAATAATAACCCTGGGTGAGTTAATTATAATATAATTATAAGTTAACTAAATGTTAATATTGGCGGGGTGGATTTATGCCTTTATTAGTAATCCTGAAACTCTGCGTCGTATTAGCCAg  >  5:848146/1‑148 (MQ=255)
                                                                                                                                                                                                                                                                              |                                                                                                                                              
AGGAAAATCAAGGTTCAAATCTATCGGGGCCTTTTTTACCCGGCCAGTCACGGACGTTCAGTTACCGGCTTGCCTGAATAGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAGCGCATCAGGTAATTTTGCGTTTATCTTCACTCTCAAGCCACGTATATGTGGCTTTATTTTTAACAAAATAATAACCCTGGGTGAGTTAATTATAATATAATTATAAGTTAACTAAATGTTAATATTGGCGGGGTGGATTTATGCCTTTATTAGTAATCCTGAAACTCTGCGTCGTATTAGCCAG  >  NC_000913/770535‑770947

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