Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,116,854 A→C 5.1% K16N (AAA→AAC yceA → putative rhodanese‑related sulfurtransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,116,8540AC 5.1% 539.3 / 15.5 155K16N (AAA→AACyceAputative rhodanese‑related sulfurtransferase
Reads supporting (aligned to +/- strand):  ref base A (82/65);  new base C (3/5);  total (85/70)
Fisher's exact test for biased strand distribution p-value = 4.69e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.79e-01

Only 100 of 175 total aligned reads displayed.

TTTTTCAACCTATTCGGGCAATTGTATGTATTGTCGCATTTTTTCGCCCGCAACCAAAATTTGTGGCTGAAGACTGGGCGAAATTGCCGCGCTTGTAAATAACAAATAATTTTTAATGCGCAAATGTAGCGTAAAATGTGTGGATGTTAATTATCGATAATTGCTATATCATGCCGCGGATTTTTACTTTCCCATCTCGCAGGAACCGTACACCATGCCAGTGTTACACAACCGCATTTCCAACGACGCGCTAAAAGCCAAAATGTTGGCTGAGAGCGAACCGCGAACCACCATTTCGTTTTACAAGTATTTCCACATCGCCGATCCTAAGGCGACCCGTGACGCTTTATATCAGCTGTTTACCGCGCTGAATGTTTTTGGGCGAGTGTATCTGGCGCATGAGGGCATTAACGCGCAAATCAGCGTACCTGCGAGCAATGTTGAAACATTTCGCGCGCAGCTCTATGCCTTCGACCCGGCTTTAGAGGGTTTACGCCTGAATATCGCGTTGGATGATGACGGGAAATC  >  NC_000913/1116593‑1117120
                                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                                          
          tATTCGGGCAATTGTATGTATTGTCGCATTTTTTCGCCCGCAACCAAAATTTGTGGCTGAAGACTGGGCGAAATTGCCGCGCTTGTAAATAACAAATAATTTTTACTGCGCAAATGTAGCGTAAAATGTGTGGATGTTAATTATCGATAATTGCTATATCATGCCGCGGATTTTTACTTTCCCATCTCGCAGGAACCNTACACCATCCCAGTGTTACACAACCTCATTTCCAACAACCCCCTAAAAGCCAAAACGTTAGCTAAGAgtgaacccca                                                                                                                                                                                                                                                     >  3:123725/1‑272 (MQ=255)
              cGGGCAATTGTATGTATTGTCGCATTTTTTCGCCCGCAACCAAAATTTGTGGCTGAAGAATGGGCGAAATTGCCGCGCTTGTAAATAACAAATAATTTTTAATGCGCAAATGTAGCGTAAAATGTGTGGATGTTAATTATCGATAATTGCTATATCATGCCGCGGATTTTTACTTTCCCATCTCGCAGGAACCGTACACCATGCCAGTGTTACACAACCGCATTTCCAACGACGCGCTAAAAGCCAAAATGTTGGCTGAGAgcgaaccgcgaac                                                                                                                                                                                                                                                  <  1:123725/274‑1 (MQ=255)
                                                                   tGAAGACTGGGCGAAATTGCCGCGATTGTAAATAACAAATAATTTTTAATGCGCAAATGTAGCGAAAAATGTGTGGATGTTAATTATCGAGAATTGCTATATCATGCCGCGGATTTTTACTTTCCCATCTCGCAGGAACCGTACACCATGCCAGTGTTACACAACCGCATTTCCAACGACGCGCTAAAAGCCAAAATGTTGGCTGAGAGCGAACCGCGAACCACCATTTCGTTTTACAAGTATTTCCACATCGCCGATCCTAAGGCGACCCGTGa                                                                                                                                                                                            <  2:186595/275‑1 (MQ=255)
                                                                   agaaagatGGGTGAAATTGCGGCATTGTTNNATAAATAATAATTTTTTATGCGACAAGGTGNCGTAAAATGTTTGGCTGTTAATTATGGATAAGTGCTATATCATGCCGCTGATTTTTACTTTCCCATCTCGCAGGNACCGTGCTCCCTGCCAGTGTTACACAACCGCATTTCCAACGACGCGCTAAAAGCCAAAATTTTGGCTGAGAGCGAACCGCGAACCACCATTTCGGTTTACAAGTATTTCCACATCGCCGATCCTAAGGCGACCCGTGa                                                                                                                                                                                            <  3:71801/268‑1 (MQ=255)
                                                                           gggCGAAATTGCCGCGCTTGTAAATAACAAATAATTTTTAATGCGCAAATGTAGCGTAAAATGTGTGGATGTTAATTATCGATAATTGCTATATCATGCCGCGGATTTTTACTTTCCCATCTCGCAGGAACCGTACACCATGCCAGTGTTACACAACCGCATTTCCAACGACGCGCTAAAAGCCAAAATGTTGGCTGAGAGCGAACCGCGAACCACCATTTCGTTTTAAAAGTATTTCCACATCCCCAATCCTAAGGCCAACCGGGAGCtttta                                                                                                                                                                                    >  4:277248/1‑272 (MQ=255)
                                                                            ggCGAAATTGCCGCGCTTGTAAATAACAAATAATTTTTAATGCGCAAATGTAGCGTAAAATGTGTGGATGTTAATTATCGATAATTGCTATATCATGCCGCGGATTTTTACTTTCCCATCTCGCAGGAACCGTACACCATGCCAGTGTTACACAACCGCATTTCCAACGACGCGCTAAAAGCAAAAATGTTGGCTGAGAGCGAACCGCGAACCACCATCCCGTTTTAAAAGTATTTCAACCTANNNCATCCTAAAGCCACCCGTGACGCTCta                                                                                                                                                                                     >  3:146567/1‑273 (MQ=255)
                                                                                                                                gCGTAAAATGTGTGGATGTTAATTATCGATAATTGCTATATCATGCCGCGGATTTTTACTTTCCCATCTCGCAGGAACCGTACACCATGCCAGTGTTACACAACCGCATTTCCAACGACGCGCTAAAAGCCAAAATGTTTGCTGAGAGCGAACCGCGCACCACCATTTCGTTTTAAAAGTATTTCCACATCGCCGATCCTAAGGCGACCCGTGACGCTTTTTATCAGCTGTTTACCGCGCTGAATGTTTTTGGGGTGGTGTATCTGGGGCAagag                                                                                                                               >  4:524483/1‑275 (MQ=255)
                                                                                                                                    aaaaTGGGTGGATGTTAATTATCGATAATTGCTATATCATGCCGCGGATTTTTACTTTCCCATCTCGCAGGAACCGTACACCATGCCAGTGTTACACAACCGCATTTCCAACGACGCGCTAAAAGCCAAAATGTTGGCTGAGAGCGAACCGCGAACCACCATTTCGTTTTACAAGTATTTCCACATCGCCGATCCTAAGGCGACCCGTGACGCTTTATATCAGCTGTTTACCGCGCTGAATGTTTTTGGGCGAGTGTATc                                                                                                                                          <  2:524484/260‑1 (MQ=255)
                                                                                                                                           gtgGATGTTAATTATCGATAATTGCTATATCATGCCGCGGATTTTTACTTTCCCATCTCGCAGGAACCGTACACCATGCCAGTGTTACACAACCGCATTTCCAACGACGCGCTAAAAGCCAAAATGTTGGCTGAGAGCGAACCGCGAACCACCATTTCGTTTTACAAGTATTTCCACATCGCCGATCCTAAGGCGACCCGTGACGCTTTATATCAGCTGTTTACCGCGCTGAATGTTTTTGGGCGAGTGTATCTGGCGCATGAGGGCATTAAcg                                                                                                                     >  1:76732/1‑274 (MQ=255)
                                                                                                                                           gtgGATGTTAATTATCGATAATTGCTATATCATGCCGCGGATTTTTACTTTCCCATCTCGCAGGAACCGTACACCATGCCAGTGTTACACAACCGCATTTCCAACGACGCGCTAAAAGCCAAAATGTTGGCTGAGAGCGAACCGCGAACCACCATTTCGTTTTACAAGTATTTCCACATCGCCGATCCTAAGGCGACCCGTGACGCTTTATATCAGCTGTTTACCGCGCTGAATGTTTTTGGGCGAGTGTATCTGGCGCATGAGGGCATTAAcg                                                                                                                     >  2:51865/1‑274 (MQ=255)
                                                                                                                                            tgGATGTTAATTATCGATAATTGCTATATCATGCCGCGGATTTTTACTTTCCCATCTCGCAGGAACCGTACACCATGCCAGTGTTACACAACCGCATTTCCAACGACGCGCTAAAAGCCAAAATGTTGGCTGAGAGCGAACCGCGAACCACCATTTCGTTTTACAAGTATTTCCACATCGCCGATCCTAAGGCGACCCGTGACGCttt                                                                                                                                                                                      >  3:88037/1‑208 (MQ=255)
                                                                                                                                               aTGTTAATTATCGATAATTGCTATATCATGCCGCGGATTTTTACTTTCCCATCTCGCAGGAACCGTACACCATGCCAGTGTTACACAACCGCATTTCCAACGACGCGCTAAAAGCCAAAATGTTGGCTGAGAgcgaaccg                                                                                                                                                                                                                                                       <  7:882934/140‑1 (MQ=255)
                                                                                                                                                tGTTAATTATCGATAATTGCTATATCATGCCGCGGATTTTTACTTTCCCATCTCGCAGGAACCGTACACCATGCCAGTGTTACACAACCGCATTTCCAACGACGCGCTAAAAGCCAAAATGTTGGCTGAGAGCGAACCGCGAACCACCATTTCGTTTTACAAGTATTTCCACATCGCCGATCCTAAGGCGACCCGTGACGCTTTatat                                                                                                                                                                                  <  1:88037/208‑1 (MQ=255)
                                                                                                                                                  ttAATTATCGATAATTGCTATATCATGCCGCGGATTTTTACTTTCCCATCTCGCAGGAACCGTACACCATGCCAGTGTTACACAACCGCATTTCCAACGACGCGCTAAAAGCCAAAATGTTGGCTGAGAGCGAACCGCGAACCACCAt                                                                                                                                                                                                                                            >  7:945188/1‑148 (MQ=255)
                                                                                                                                                   tAATTATCGATAATTGCTATATCATGCCGCGGATTTTTACTTTCCCATCTCGCAGGAACCGTACACCATGCCAGTGTTACACAACCGCATTTCCAACGACGCGCTAAAAGCCAACATTTTGGCt                                                                                                                                                                                                                                                                   >  5:81783/1‑124 (MQ=255)
                                                                                                                                                   tAATTATCGATAATTGCTATATCATGCCGCGGATTTTTACTTTCCCATCTCGCAGGAACCGTACACCATGCCAGTGTTACACAACCGCATTTCCAACGACGCGCTAAAAGCCAAAATGtt                                                                                                                                                                                                                                                                       <  8:23667/120‑1 (MQ=255)
                                                                                                                                                   tAATTATCGATAATTGCTATATCATGCCGCGGATTTTTACTTTCCCATCTCGCAGGAACCGTACACCATGCCAGTGTTACACAACCGCATTTCCAACGACGCGCTAAAAGCCAAAATGtt                                                                                                                                                                                                                                                                       >  7:23667/1‑120 (MQ=255)
                                                                                                                                                   tAATTATCCATAATTGCTATATCATGCCGCGGATTTTTACTTTCCCATCTCGCAGGAACCGTACACCATGCCAGTGTTACACAACCGCATTTCCAACGACGCGCTAAAAGCCAAAATGta                                                                                                                                                                                                                                                                       >  8:980716/1‑119 (MQ=255)
                                                                                                                                                    aaTTATCGATAATTGCTATATCATGCCGCGGATTTTTACTTTCCCATCTCGCAGGAACCGTACACCATGCCAGTGTTACACAACCGCATTTCCAACGACGCGCTAAAAGCCAAAATGTTGGCTgag                                                                                                                                                                                                                                                                <  6:467529/126‑1 (MQ=255)
                                                                                                                                                    aaTTATCGATAATTGCTATATCATGCCGCGGATTTTTACTTTCCCATCTCGCAGGAACCGTACACCATGCCAGTGTTACACAACCGCATTTCCAACGACGCGCTAAAAGCCAAAATGTTGGCTgag                                                                                                                                                                                                                                                                >  5:467529/1‑126 (MQ=255)
                                                                                                                                                     aTTATCGATAATTGCTATATCATGCCGCGGATTTTTACTTTCCCATCTCGCAGGAACCGTACACCATGCCAGTGTTACACAACCGCATTTCCAACGACGCGCTAAAAGCCAAAATGTTGGCTGAGAgcga                                                                                                                                                                                                                                                           <  7:401264/130‑1 (MQ=255)
                                                                                                                                                            aTAATTGCTATATCATGCCGCGGATTTTTACTTTCCCATCTCGCAGGAACCGTACACCATGCCAGTGTTACACAACCGCATTTCCAACGACGCGCTAAAAGCCAAAATGTTGGCTGAGAGCGAATCGCGAACCACCATTTCGTTTTACAATTATTTCCACATCGCCGATCCTACGGCGACCCGTGACGCTTTATATCAGCTTTTTACCGCGCTGAATGTTTTTGGGCGAGTGTATCTGGCGCATGTGGGCCTTACCGCGCAATTGAGCGTACCTg                                                                                                   >  4:166256/1‑275 (MQ=255)
                                                                                                                                                                  gCTATATCATGCCGCGGATTTTTACTTTCCCATCTCGCAGGAACCGTACACCATGCCAGTGTTACACAACCGCATTTCCAACGACGCGCTAAAAGCCAAAATGTTGGCTGAGAGCGAACCGCGAACCACCATTTCGTTTTACAAGTAt                                                                                                                                                                                                                            <  8:945188/148‑1 (MQ=255)
                                                                                                                                                                     atatCATGCCGCGGATTTTTACTTTCCCATCTCGCAGGAACCGTACACCATGCCAGTGTTACACAACCGCATTTCCAACGACGCGCTAAAAGCCAAAATGTTGGCTGAGAGCGAACCGCGAACCACCATTTCGTTTTACAAGTATTTCCACATCGCCGATCCTAAGGCGACCCGTGACGCTTTATATCAGCTGTTTACCGCGCTGAATGTTTTTGGGCGAGTGTATCTGGCGCATGAGGGCATTAACGCGCAAATCAGCGTACCTGCGAGCAATg                                                                                          >  2:256608/1‑275 (MQ=255)
                                                                                                                                                                           tGCCGCGGATTTTTACTTTCCCATCTCGCAGGAACCGTACACCATGCCAGTGTTACACAACCGCATTTCCAACGACGCGCTAAAAGCCaaa                                                                                                                                                                                                                                                                            >  5:365695/1‑91 (MQ=255)
                                                                                                                                                                           tGCCGCGGATTTTTACTTTCCCATCTCGCAGGAACCGTACACCATGCCAGTGTTACACAACCGCATTTCCAACGACGCGCTAAAAGCCAAAATGTTGGCTGAGAGCGAACCGCGAACCACCATTTCGTTTTACAAGTATTTCCACATCGCCGATCCTAAGGCGACCCGTGACGCTTTATATCAGCTGTTTACCGCGCTGAATGTTTTTGGGCGAGTGTATCTGGCGCATGAGGGCATTAACGCGCACATCAGAGTACCTGCGAGAAATGTTGaa                                                                                     >  4:390307/1‑274 (MQ=255)
                                                                                                                                                                           tGCCGCGGATTTTTACTTTCCCATCTCGCAGGAACCGTACACCATGCCAGTGTTACACAACCGCATTTCCAACGACGCGCTAAAAGCCAAAATGTTGGCTGAGAGCGAACCGCGAACCACCATTTCGTTTTACAAGTATTTCCACATCGCCGATCCTAAGGCGACCCGTGACGCTTTATATCAGCTGTTTACCGCGCTGAATGTTTTTGGGCGAGTGTATCTGGCGCATGAGGGCATTAACGCGCAAATCAGCGTACCTGCGAGCAATGTTGaa                                                                                     >  2:429821/1‑274 (MQ=255)
                                                                                                                                                                           tGCCGCGGATTTTTACTTTCCCATCTCGCAGGAACCGTACACCATGCCAGTGTTACACAACCGCATTTCCAACGACGCGCTAAAAGCCAAAATGTTGGCTGAGAGCGAACCGCGAACCACCATTTCGTTTTACAAGTATTTCCACATCGCCGACCCTAAGGCGACCCGTGCCGCTTTATATCAGCTTTTTACCGCGTTGAATGTTTTTGGGCGAGGTTATCCGGCTCATGCGGCCATTAACGCNNNAATAAGGGTACATGTGAGCAATGTTGaa                                                                                     >  3:145784/1‑274 (MQ=255)
                                                                                                                                                                            gCCGCGGATTTTTACTTTCCCATCTCGCAGGAACCGTACACCATGCCAGTGTTACACAACCGCATTTCCAACGACGCGCTAAAAGCCACAATGtt                                                                                                                                                                                                                                                                       >  8:1105040/1‑95 (MQ=255)
                                                                                                                                                                            gCCGCGGATTTTTACTTTCCCATCTCGCAGGAACCGTACACCATGCCAGTGTTACACAACCGCATTTCCAACGACGCGCTAAAAGCCAAAATGtt                                                                                                                                                                                                                                                                       <  7:1105040/95‑1 (MQ=255)
                                                                                                                                                                             ccGCGGATTTTTACTTTCCCATCTCGCAGGAACCGTACACCATGCCAGTGTTACACAACCGCATTTCCAACGACGCGCTAAAAGCCAAAATGTTGGCTGAGAGCGAACCGCGAACCACCATTTCGTTTTACAAGTATTTCCACATCGCCGATCCTAAGGCGACCCGTGACGCTTTATATCAGCTGTTTACCGCGCTGAATGTTTTTGGGCGAGTGTATCTGGCGCATGAGGGCATTAACGCGCAAATCAGCGTACCTGCGAGCAATGTTGAAAc                                                                                   >  2:123785/1‑274 (MQ=255)
                                                                                                                                                                               gcgGATTTTTACTTTCCCATCTCGCAGGAACCGTACACCATGCCAGTGTTACACAACCGCATTTCCAACGACGCGCTAAAAGCCAAAATGTTGGCTGAGAgcgaa                                                                                                                                                                                                                                                          <  8:837354/105‑1 (MQ=255)
                                                                                                                                                                                cgGATTTTTACTTTCCCATCTCGCAGGAACCGTACACCATGCCAGTGTTACACAACCGCATTTCCAACGACGCGCTAAAAGCCAAAATGTTGGCTGAGAGCGAACCGCGAAccacca                                                                                                                                                                                                                                             <  7:17022/117‑1 (MQ=255)
                                                                                                                                                                                cgGATTTTTACTTTCCCATCTCGCAGGAACCGTACACCATGCCAGTGTTACACAACCGCATTTCCAACGACGCGCTAAAAGCCAAAATGTTGGCTGAGAGCGAACCGCGAACCACCATTTCGTTTTACAAGTATTTCCACATCGCCGATCCTAAGGCGACCCGTGACGCTTTATATCAGCTGTTTACCGCGCTGAATGTTTTTGGGCGAGTGTATCTGGCCCATGCGGGCATTAACGCGCAAATCAGCGTACCTCCGAGCAATGTTGAAACAttt                                                                               >  2:131760/1‑275 (MQ=255)
                                                                                                                                                                                            ttCCCATCTCGCAGGAACCGTACACCATGCCAGTGTTACACAACCGCATTTCCAACGACGCGCTAAAAGCCAACAtt                                                                                                                                                                                                                                                                         >  8:535189/1‑76 (MQ=255)
                                                                                                                                                                                             tCCCATCTCGCAGGAACCGTACACCATGCCAGTGTTACACAACCGCATTTCCAACGACGCGCTAAAAGCCAAAATGTTTGCTGAGAGCGAACCGCGAACCACCATTTCGTTTTACAAGTATTTCCACATCGCCGATCCTAAGACGACCCGTGACGCTTTATATCAGCTGTTTACCGCGCTGTATGTTTTTGGGCGAGTGTATCTGGCCCATGAGGGCATTACCGCGCAAATAATCGTCACTGCGAGCAATTTTTAACCATTTCCCGCGCaacact                                                                  >  4:480333/1‑270 (MQ=255)
                                                                                                                                                                                                cATCTCGCAGGAACCGTACACCATGCCAGTGTTACACAACCGCATTTCCAACGACGCGCTAAAAGCCAAAATGt                                                                                                                                                                                                                                                                        >  6:81819/1‑74 (MQ=255)
                                                                                                                                                                                                cATCTCGCAGGAACCGTACACCATGCCAGTGTTACACAACCGCATTTCCAACGACGCGCTAAAAGCCAAAATGTTGGCTgag                                                                                                                                                                                                                                                                >  8:557708/1‑82 (MQ=255)
                                                                                                                                                                                                         gacgccGTAACCCCTGCCAGGTTGACACAACCGCATTCCCAACGACGCGCTAAAAGCCAAAATGTTGGCTGAGAGCGAACCGCGAACCACCATTTCGTTTTAAAAGTATTTCCACATCGCCGATCCTAAGGCGACCCGTGACGCTTTATATCAGCTGTTTACCGCGCTGAATGTTTTTGGGCGAGTGTATCTGGCGCATGAGGGCATTAACGCGCAAATCAGCGTACCTGCGAGCAATGTTGAAACATTTCGCGCGCAGCTCTATGCCTTCGAcc                                                      <  4:269186/271‑1 (MQ=255)
                                                                                                                                                                                                            aCCGTACACCATGCCAGTGTTACACAACCGCATTTCCAACGACGCGCTAAAAGCCAAAATGTTGGCTGAGAGCGAACCGCGAACCACCATTTCGTTTTACAAGTATTTc                                                                                                                                                                                                                         <  7:1047172/109‑1 (MQ=255)
                                                                                                                                                                                                                  caacATGCCAGGGTTTACAAACCGAATTTCCAACCGCGCGCTAAAAGCCAAAATGTTGGCTGAGAGCGAACCGCGAACCACCATTTCGTTTTACAAGTATTTCCACATCGCCGATCCTAAGGCGTCCCGGGACGCTTTATATCAGCTGTTTACCGCGCTGAATGTTTTTGGGCGAGTGTATCTGGCGCATGAGGGCATTAACGCGCAAATCAGCGTACCTGCGAGCAATGTTGAAACATTTCGCGCGCAGCTCTATGCCTTCGACCCGGCTTTag                                             <  4:212134/272‑1 (MQ=255)
                                                                                                                                                                                                                    ccATGCCAGTGTTACACAACCGCATTTCCAACGACGCGCTAAAAGCCAAAATGTTGGCTGAGAGCGAACCGCGAACCACCATTTCGt                                                                                                                                                                                                                                       >  7:991179/1‑87 (MQ=255)
                                                                                                                                                                                                                       tGCCAGTGTTACACAACCGCATTTCCAACGACGCGCTAAAAGCCAAAATGTTGGCTGAGGGCGAACCGCGAACCACCATTTCGTTTTACAAGTATTTCCACATCGCCGATCCTAAGGCGACCCGTGACGCTTTATATCAGCTGTTTAc                                                                                                                                                                       >  8:345331/1‑148 (MQ=255)
                                                                                                                                                                                                                       tGCCAGTGTTACACAACCGCATTTCCAACGACGCGCTAAAAGCCAAAATGTTGGCTGAGAGCGAACCGCGAACCACCATTTCGTTTTACAAGTAttt                                                                                                                                                                                                                          <  6:85734/97‑1 (MQ=255)
                                                                                                                                                                                                                       tGCCAGTGTTACACAACCGCATTTCCAACGACGCGCTAAAAGCCAAAATGTTGGCTGAGAGCGAACCGCGAACCACCATTTCGTTTTACAAGTAttt                                                                                                                                                                                                                          >  5:85734/1‑97 (MQ=255)
                                                                                                                                                                                                                       tGCCAGTGTTACACAACCGCATTTCCAACGACGCGCTAAAAGCCAAAATGTTGGCTGAGAGCGAACCGCGAACCACCATTTCGTTTTACAAGTATTTCCACATCGCCGATCCTAAGGCGACCCGTGACGCTTTATATCAGCTGTTTACCGCGCTGAATGTTTTTGGGCGAGTGTATCTGGCGCATGAGGGCATTAACGCGCAAATCAGCGTACCTGCGAGCAATGTTGAAACATTTCGCGCGCAGCTCTATGCCTTCGACCCGGCTTTAGAGGGt                                        >  2:4052/1‑275 (MQ=255)
                                                                                                                                                                                                                        gCCAGTGTTACACAACCGCATTTCCAACGACGCGCTAAAAGCCAAAATGTTGGCTGAGAGCGAACCGCGAACCACCATTTCGTTTTACAAGTATTTCCACATCGCCGATCCTaa                                                                                                                                                                                                        >  6:403634/1‑114 (MQ=255)
                                                                                                                                                                                                                        gCCAGTGTTACACAACCGCATTTCCAACGACGCGCTAAAAGCCAAAATGTTGGCTGAGAGCGAACCGCGAACCACCATTTCGTTTTACAAGTATTTCCACATCGCCGATCCTAAGGCGACCCGTGACGCTTTATATCAGCTGTTTACCGCGCTGAATGTTTTTGGGCGAGTGTATCTGGCGCATGAGGGCATTAACGCGCAAATCAGCGTACCTGCGa                                                                                                >  2:62461/1‑218 (MQ=255)
                                                                                                                                                                                                                        gCCAGTGTTACACAACCGCATTTCCAACGACGCGCTAAAAGCCAAAATGTTGGCTGAGAGCGAACCGCGAACCACCATTTCGTTTTACAAGTATTTCCACATCGCCGATCCTAAGGCGACCCGTGACGCTTTATATCAGCTGTTTACCGCGCTGAATGTTTTTGGGCGAGTGTATCTGGCGCATGAGGGCATTAACGCGCAAATCAGCGTACCTGCGAGCAATGTTGAAACATTTCGCGCGCAGCTCTATGCCTTCGACCCGGCTTTaaggctgt                                       >  2:269187/1‑268 (MQ=255)
                                                                                                                                                                                                                          cAGTGTTACACAACCGCATTTCCAACGACGCGCTAAAAGCCAAAATGTTGGCTGAGAGCGAACCGCGAACCACCATTTCGTTTTACAAGTATTTCCACATCGCCGATCCTAAGGCGACCCGTGACGCTTTATATCAGCTGTTTACCGCGCTGAATGTTTTTGGGCGAGTGTATCTGGCGCATGAGGGCATTAACGCGCAAATCAGCGTACCTGCGAGCAATGTTGAAACATTTCGCGCGCAGCTCTATGCCTTCGACCCGGCTTTAGAGGGttt                                      >  4:337193/1‑274 (MQ=255)
                                                                                                                                                                                                                           aGTGTTACACAACCGCATTTCCAACGACGCGCTAAAAGCCAAAATGTTGGCTGAGAGCGAACCGCGAACCACCATTTCGTTTTACGAGTATTTCCACATCGCCGATCCTAAGGCGACCCGTGACGCttt                                                                                                                                                                                      <  5:941925/129‑1 (MQ=255)
                                                                                                                                                                                                                           aGTGTTACACAACCGCATTTCCAACGACGCGCTAAAAGCCAAAATGTTGGCTGAGAGCGAACCGCGAACCACCATTTCGTTTTACAAGTATTTCCACATCGCCGATCCTAAGGCGACCCGTGACGCTTTATATCAGCTGTTTACcgcg                                                                                                                                                                   >  8:159602/1‑148 (MQ=255)
                                                                                                                                                                                                                            gtgtTACACAACCGCATTTCCAACGACGCGCTAAAAGCCAAAATGTTGGCTGAGAGCGAACCGCGAACCACCATTTCGTTTTACAAGTATTTCCACATCGCCGATCCTAAGGCGACCCGTGACGCTTTATATCAGCTGTTTACCGCGCTGAATGTTTTTGGGCGAGTGTATCTGGCGCATGAGGGCATTAACGCGCAAATCAGCGTACCTGCGAGCaa                                                                                            <  4:62460/218‑1 (MQ=255)
                                                                                                                                                                                                                              gtTACACAACCGCATTTCCAACGACGCGCTAAAAGCCAAAATGTTGGCTGAGAGCGAACCGCGAACCACCATTTCGTTTTACAAGTATTTCCACATCGCCGATCCTAAGGCGACCCGTGACGCTTTATATCAGCTGTTTACCGCGCTGAATGTTTTTGGGCGAGTGTATCTGGCGCATGAGGGCATTAACGCGCAAATCAGCGTACCTGCGAGCAATGTTGAAACATTTCGCGCGCAGCTCTATGCCTTCGACCCGGCTTTAGAGGGTTTACGc                                  <  2:337194/274‑1 (MQ=255)
                                                                                                                                                                                                                                 tACACAACCGCATTTCCAACGACGCGCTAAAAGCCAAAATGTTGGCTGAGAGCGAACCGCGAACCACCATTTCGTTTTACAAGTATTTCCACATCGCCGATCCTAAgg                                                                                                                                                                                                      <  6:879478/108‑1 (MQ=255)
                                                                                                                                                                                                                                 tACACAACCGCATTTCCAACGACGCGCTAAAAGCCAAAATGTTGGCTGAGAGCGAACCGCGAACCACCATTTCGTTTTACAAGTATTTCCACATCGCCGATCCTAAgg                                                                                                                                                                                                      >  5:879478/1‑108 (MQ=255)
                                                                                                                                                                                                                                  acacAACCGCATTTCCAACGACGCGCTAAAAGACAAAATGTTGGCTGAGAGCGAACCGCGAACCACCATTTCGTTTTACAAGTATTTCCACATCGCCGATCCTAAGGCGACCCGTGACGCTTTATATCAGCTGTTTACCGCGCTGAATGTTTTTGGGCGAGTGTATCTGGCGCATGAGGGCATTAACGCGCACATCAGCGTACCTGCGAGCAATGTTGAAACATTTCGCGCGCAGCTCTATGCCTTCGACCCGGCTTTAGAGGGTTTAAGACTGa                              >  2:327038/1‑275 (MQ=255)
                                                                                                                                                                                                                                   cacaACCGCATTTCCAACGACGCGCTAAAAGCCaaa                                                                                                                                                                                                                                                                            >  6:145721/1‑36 (MQ=255)
                                                                                                                                                                                                                                   cacaACCGCATTTCCAACGACGCGCTAAAAGCCaaa                                                                                                                                                                                                                                                                            <  5:145721/36‑1 (MQ=255)
                                                                                                                                                                                                                                   cacaACCGCATTTCCAACGACGCGCTAAAAGCCAAAATGTTGGCTGAGAGCGAACCGCGAACCACCATTTCGTTTTACAAGTATTTCCACATCGCCGATCCTAAGGCGACCCGTGACGCtt                                                                                                                                                                                       <  5:678285/121‑1 (MQ=255)
                                                                                                                                                                                                                                    acaACCGCATTTCCAACGACGCGCTAAAAGCCAAAATGTTg                                                                                                                                                                                                                                                                      >  7:764371/1‑41 (MQ=255)
                                                                                                                                                                                                                                    acaACCGCATTTCCAACGACGCGCTAAAAGCCAAAATGTTg                                                                                                                                                                                                                                                                      <  8:764371/41‑1 (MQ=255)
                                                                                                                                                                                                                                    acaACCGCATTTCCAACGACGCGCTAAAAGCCAAAATGTTGGCTGAGAGCGAACCGCGAACCACCATTTCGTTTTACAAGTATTTCCACATCGCCGATCCTAAGGCGACCCGTGACGCTTTATATCAGCTGTTTACCGCGCTGAATGTTTTTGGGCGAGTGTATCTGGCGCATGAGGGCATTAACGCGCAAATCAGCGTACCTGCGAGCAATGTTGAAACATTTCGCGCGCAGCTCTATGCCTTCGAcc                                                      <  1:38189/249‑1 (MQ=255)
                                                                                                                                                                                                                                     caACCGCATTTCCAACGACGCGCTAAAAGCCAAAATGTTGGCTGAGAGCGAACCGCGAACCACCATTTCGTTTTACAAGTATTTCCACATCGCCGATCCTAAGGCGACCCGTGACGCTTTATATCAGCTGTTTACCGCGCTGAATGtt                                                                                                                                                          <  5:202557/148‑1 (MQ=255)
                                                                                                                                                                                                                                       aCCGCATTTCCAACGACGCGCTAAAAGCCAAAATGTTGGCTGAGAGCGAACCGCGAACCACCATTTCGTTTTACAAGTATTTCCACATCGCCGATCCTAAGGCGACCCGTGACGCTTTATATc                                                                                                                                                                                 >  6:983776/1‑123 (MQ=255)
                                                                                                                                                                                                                                         cGCATTTCCAACGACGCGCTAAAAGCCAAAATGTTGGCTGAGAGCGAACCGCGAACCACCATTTCGTTTTACa                                                                                                                                                                                                                                 >  8:1063201/1‑73 (MQ=255)
                                                                                                                                                                                                                                         cGCATTTCCAACGACGCGCTAAAAGCCAAAATGTTGGCTGAGAGCGAACCGCGAACCACCATTTCGTTTTACa                                                                                                                                                                                                                                 <  7:1063201/73‑1 (MQ=255)
                                                                                                                                                                                                                                          gCATTTCCAACGACGCGCTAAAAGCCAAAATGTTGGCTGAGAGCGAACCGCGAACCACCATTTCGTTTTACAAGTATTTCCACATCGCCGATCCTAAGGCGACCCGTGACGCTTTATATCAGCTGTTTACCGCGCt                                                                                                                                                                 <  6:167222/136‑1 (MQ=255)
                                                                                                                                                                                                                                          gCATTTCCAACGACGCGCTAAAAGCCAAAATGTTGGCTGAGAGCGAACCGCGAACCACCATTTCGTTTTACAAGTATTTCCACATCGCCGATCCTAAGGCGACCCGTGACGCTTTATATCAGCTGTTTACCGCGCt                                                                                                                                                                 >  5:167222/1‑136 (MQ=255)
                                                                                                                                                                                                                                           cATTTCCAACGACGCGCTAAAAGCCAAAATGTTGGCTGAGAGCGAACCGCGAACCACCATTTCGTTTTACAAGTCTTTCCACATCGCCGATCCTAAGGCGACCCGTGACGCTTTATATCAGCTGTTTACCGCGCTGCATGTTTTTGGGCGAGGGTATCTGGCGCATGGGGGCATTAACGCGCAAATCAGCGTACCTGCGAGCAATGTTGAAACATTTCGCGCGCAGCGCTATGACTTCGACCCGGCTTTAGAGGGTTTAACCCGTAATATAGCGt                     >  4:148096/1‑275 (MQ=255)
                                                                                                                                                                                                                                             tttCCAACGACGCGCTAAAAGCCAAAATGTTGGCTGAGAGCGAACCGCGAACCACCATTTCGTTTTACAAGTATTTCCACATCGCCGAt                                                                                                                                                                                                             <  8:265183/89‑1 (MQ=255)
                                                                                                                                                                                                                                             tttCCAACGACGCGCTAAAAGCCAAAATGTTGGCTGAGAGCGAACCGCGAACCACCATTTCGTTTTACAAGTATTTCCACATCGCCGAt                                                                                                                                                                                                             >  7:265183/1‑89 (MQ=255)
                                                                                                                                                                                                                                              ttCCAACGACGCGCTAAAAGCCAAAATGTTGGCTGAGAGCGAACCGCGAACCACCATTTCGTTTTACAAGTATTTCCACATCGCCGATCCTAAGGCGACCCGTGACGCTTTATATCAGCTGTTTACCGCGCTGAATGTTTTTg                                                                                                                                                      <  5:457562/143‑1 (MQ=255)
                                                                                                                                                                                                                                                caaaCGACGCGCTAAAAGCCAAAATGTTGGCTGAGAGCGAACCGCGAACCACCATCTCGTTTTACAAGTATTTCCACATCGCCGATCCTAAGGCGACCCGTGACGCTTTATATCAGCTGTTTACCGCGCTGAATGTTTTTGGGCGAGTGTATCTGGCGCATGAGGGCATTAACGCGCAAATCAGCGTACCTGCGAGCAATGTTGAAACATTTCGCGCGCAGCTCTATGCCTTCGACCCGGCTTTAGAGGGTTTACGCCTGAATATCGCGTTGgat                <  1:146567/273‑1 (MQ=255)
                                                                                                                                                                                                                                                  aaCGACGCGCTAAAAGCCAAAATGTTGGCTGAGAGCGAACCGCGAACCACCATTTCGt                                                                                                                                                                                                                                       <  7:1127791/58‑1 (MQ=255)
                                                                                                                                                                                                                                                    cgacgCGCTAAAAGCCAAAATGTTGGCTGAGAGCGAACCGCGAACCACCATTTCGTTTTACAAGTATTTCCACATCGCCGATCCTAAg                                                                                                                                                                                                       >  7:540621/1‑88 (MQ=255)
                                                                                                                                                                                                                                                       cgcgCTAAAAGCCAAAATGTTGGCTGAGAGCGAACCGCGAACCACCATTTCGTTTTACAAGTATTTCCACATCGCCGATCCTAAGGCGACCCGTGACGCTTTATATCAGCTGTTTACCGCGCTGAATGTTTTTGGGCGAGTGTATCTGGCGCATGAGGGCATTAACGCGCAAATCAGCGTACCTGCGAGCAATGTTGAAACATTTCGCGCGCAGCTCTATGCCTTCGACCCGGCTTTAGAGGGTTTACGCCTGAATATCGCGTTGGATGATGACg         <  1:145784/275‑1 (MQ=255)
                                                                                                                                                                                                                                                           cTAAAAGCCAAAATGTTGGCTGAGAGCGAACCGCGAACCACCATTTCGTTTTACAAGTATTTCCACATCGCCGATCCTAAGGCGACCCGTGACGCTTTATATCAGCTGTTTACCGCGCTGAATGTTTTTGGGCGAGTGTATCTGgcgc                                                                                                                                    >  5:37516/1‑148 (MQ=255)
                                                                                                                                                                                                                                                            tAAAAGCCAAAATGTTGGCTGAGAGCGAACCGCGAACCACCATTTCGTTTTACAAGTAtt                                                                                                                                                                                                                           <  8:391912/60‑1 (MQ=255)
                                                                                                                                                                                                                                                            tAAAAGCCAAAATGTTGGCTGAGAGCGAACCGCGAACCACCATTTCGTTTTACAAGTATTTCCACATCGCCGATCCTAAGGCGACCCGTGACGCTTTATATCAGCTGTTTACCGCGCTGAATGTTTTTg                                                                                                                                                      >  7:825172/1‑129 (MQ=255)
                                                                                                                                                                                                                                                            tAAAAGCCAAAATGTTGGCTGAGAGCGAACCGCGAACCACCATTTCGTTTTACAAGTATTTCCACATCGCCGATCCTAAGGCGACCCGTGACGCTTTATATCAGCTGTTTACCGCGCTGAATGTTTTTg                                                                                                                                                      <  8:825172/129‑1 (MQ=255)
                                                                                                                                                                                                                                                            tAAAAGCCAAAATGTTGGCTGAGAGCGAACCGCGAACAACCATTTCGTTTTACAAGTATTTCCACATCGCCGATCCTAAGGCGACCCGTGACGCTTTATATCAGCTGTTTACCGCGCTGAATGTTTTTGGGCGAGTGTATCTGGCGCATGAGGGCATTAACGCGCAAATCAGCGTACCTGCGAGCAATGTTGAAACATTTCGCGCGCAGCTCTATGCCTTCGACCCGGCTTTAGAGGGTTTACGCCTGAATATCGCGTTGGATGATGACGGGaaa    <  1:71478/275‑1 (MQ=255)
                                                                                                                                                                                                                                                              aaaGCCAAAATGTTGGCGGAGGGCGAACCCGGAACCACCATTTCGTTTTACAAGTATTTCCACATCGCCGATCCTAAGGCGACCCGTGACGCTTTATATCAGCTGTTTACCGCGCTGAATGTTTTTGGGCGAGTGTATCTGGCGCATGAGGGCATTAACGCGCAAATCAGCGTACCTGCGAGCAATGTTGAAACATTTCGCGCGCAGCTCTATGCCTTCGACCCGGCTTTAGAGGGTTTACGCCTGAATATCGCGTTGGATGATGACGGGAAATc  <  4:51864/275‑1 (MQ=255)
                                                                                                                                                                                                                                                                 gCCAAAATGTTGGCTGAGAGCGAACCGCGAACCACCATTTCGTTTTACAAGTATTTCCACATCGCCGATCCTAAGGCGACCCGTGACGCTTTATATCAGCTGTTTACCGCGCTGAATGtttt                                                                                                                                                        <  6:296118/122‑1 (MQ=255)
                                                                                                                                                                                                                                                                  ccAAAATGTTGGCTGAGAGCGAACCGCGAACCACCATTTCGTTTTACAAGTATTTCCACATCGCCGATCCTAAGGCGACCCGTGACGCTTtat                                                                                                                                                                                    <  6:554228/93‑1 (MQ=255)
                                                                                                                                                                                                                                                                  ccAAAATGTTGGCTGAGAGCGAACCGCGAACCACCATTTCGTTTTACAAGTATTTCCACATCGCCGATCCTAAGGCGACCCGTGACGCTTtat                                                                                                                                                                                    <  6:121692/93‑1 (MQ=255)
                                                                                                                                                                                                                                                                  ccAAAATGTTGGCTGAGAGCGAACCGCGAACCACCATTTCGTTTTACAAGTATTTCCACATCGCCGATCCTAAGGCGACCCGTGACGCTTtat                                                                                                                                                                                    >  5:121692/1‑93 (MQ=255)
                                                                                                                                                                                                                                                                  ccAAAATGTTGGCTGAGAGCGAACCGCGAACCACCATTTCGTTTTACAAGTATTTCCACATCGCCGATCCTAAGGCGACCCGTGACGCTTtat                                                                                                                                                                                    >  5:554228/1‑93 (MQ=255)
                                                                                                                                                                                                                                                                  ccAAAATGTTGGCTGAGAGCGAACCGCGAACCACCATTTCGTTTTACAAGTATTTCCACATCGCCGATCCTAAGGCGACCCGTGACGCTTTATATCAGCTGTTTACCGCGCTGAATGTTTTTGGGCGAGTGTATCTGGCGCATGAggg                                                                                                                             >  5:281761/1‑148 (MQ=255)
                                                                                                                                                                                                                                                                   cAAAATGTTGGCTGAGAGCGAACCGCGAACCACCATTTCGTTTTACAAGTATTTCCACATCGCCGATCCTAAGGCGACCCGTGACGCt                                                                                                                                                                                        >  6:805076/1‑88 (MQ=255)
                                                                                                                                                                                                                                                                   cAAAATGTTGGCTGAGAGCGAACCGCGAACCACCATTTCGTTTTACAAGTATTTCCACATCGCCGATCCTAAGGCGACCCGTGACGCt                                                                                                                                                                                        >  6:459273/1‑88 (MQ=255)
                                                                                                                                                                                                                                                                    aaaaTGTTGGCTGAGAGCGAACCGCGAACCACCATTTCGTTTTACAAGTATTTCCACATCGCCGATCCTAAGGCGACCCGTGACGCTTTATATCAGCTg                                                                                                                                                                            <  7:62766/99‑1 (MQ=255)
                                                                                                                                                                                                                                                                    aaaaTGTTGGCTGAGAGCGAACCGCGAACCACCATTTCGTTTTACAAGTATTTCCACATCGCCGATCCTAAGGCGACCCGTGACGCTTTATATCAGCTg                                                                                                                                                                            <  5:984190/99‑1 (MQ=255)
                                                                                                                                                                                                                                                                     aaaTGTTGGCTGAGAGCGAACCGCGAACCACCATTTCGTTTTACAAGTATTTCCACATCGCCGATCCTAAGGCGACCCGTGACGCTTTATATCAGCTGTTTACCGCGCTGAATg                                                                                                                                                            >  8:468784/1‑114 (MQ=255)
                                                                                                                                                                                                                                                                     aaaTGTTGGCTGAGAGCGAACCGCGAACCACCATTTCGTTTTACAAGTATTTCCACATCGCCGATCCTAAGGCGACCCGTGACGCTTTATATCAGCTGTTTACCGCGCTGAATg                                                                                                                                                            <  7:468784/114‑1 (MQ=255)
                                                                                                                                                                                                                                                                     aaaTGTTGGCTGAGAGCGAACCGCGAACCACCATTTCGTTTTACAAGTATTTCCACATCGCCGATCCTAAGGCGACCCGTGACGCTTTATATCAGCTGTTTACCGCGCTGAATGTTTTTggg                                                                                                                                                    <  8:466355/122‑1 (MQ=255)
                                                                                                                                                                                                                                                                     aaaTGTTGGCTGAGAGCGAACCGCGAACCACCATTTCGTTTTACAAGTATTTCCACATCGCCGATCCTAAGGCGACCCGTGACGCTTTATATCAGCTGTTTACCGCGCTGAATGTTTTTggg                                                                                                                                                    >  7:466355/1‑122 (MQ=255)
                                                                                                                                                                                                                                                                     aaaTGTTGGCTGAGAGCGAACCGCGAACCACCATTTCGTTTTACAAGTATTTCCACATCGCCGATCCTAAGGCGACCCGTGACGCTTTATATCAGCTGTTTACCGCGCTGAATGTTTTTGGGCg                                                                                                                                                  <  8:1112989/124‑1 (MQ=255)
                                                                                                                                                                                                                                                                     aaaTGTTGGCTGAGAGCGAACCGCGAACCACCATTTCGTTTTACAAGTATTTCCACATCGCCGATCCTAAGGCGACCCGTGACGCTTTATATCAGCTGTTTACCGCGCTGAATGTTTTTGGGCGAGTGTATCTg                                                                                                                                        <  7:94743/134‑1 (MQ=255)
                                                                                                                                                                                                                                                                      aaTGTTGGCTGAGAGCGAACCGCGAACCACCATTTCGTTTTACAAGTATTTCCACATCGCCGATCCTAAGGCGACCCGTGACGCTTTATATCAGCTGTTTACCGCGCTGAATGTTTTTGGGCGAGTGTATCTGGCGCATg                                                                                                                                 <  7:380788/140‑1 (MQ=255)
                                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                                          
TTTTTCAACCTATTCGGGCAATTGTATGTATTGTCGCATTTTTTCGCCCGCAACCAAAATTTGTGGCTGAAGACTGGGCGAAATTGCCGCGCTTGTAAATAACAAATAATTTTTAATGCGCAAATGTAGCGTAAAATGTGTGGATGTTAATTATCGATAATTGCTATATCATGCCGCGGATTTTTACTTTCCCATCTCGCAGGAACCGTACACCATGCCAGTGTTACACAACCGCATTTCCAACGACGCGCTAAAAGCCAAAATGTTGGCTGAGAGCGAACCGCGAACCACCATTTCGTTTTACAAGTATTTCCACATCGCCGATCCTAAGGCGACCCGTGACGCTTTATATCAGCTGTTTACCGCGCTGAATGTTTTTGGGCGAGTGTATCTGGCGCATGAGGGCATTAACGCGCAAATCAGCGTACCTGCGAGCAATGTTGAAACATTTCGCGCGCAGCTCTATGCCTTCGACCCGGCTTTAGAGGGTTTACGCCTGAATATCGCGTTGGATGATGACGGGAAATC  >  NC_000913/1116593‑1117120

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