Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,098,519 A→G 6.4% intergenic (+116/+39) yggW → / ← yggM HemN family putative oxidoreductase/DUF1202 family putative secreted protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,098,5190AG 6.4% 413.5 / 12.0 125intergenic (+116/+39)yggW/yggMHemN family putative oxidoreductase/DUF1202 family putative secreted protein
Reads supporting (aligned to +/- strand):  ref base A (56/61);  new base G (3/5);  total (59/66)
Fisher's exact test for biased strand distribution p-value = 7.21e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.90e-01

Only 100 of 128 total aligned reads displayed.

ATTCGCTGCTGGAGCTTTTTCTGGCTGAGTAAACTTGTATTGCCGGATGCGGCGTGAACGCCTTATCCAGCCGACATGTGGCAGCGGTTGTAGGTCTGATAAGACGCGCAAGCGTCGCATCAGACGTTGATTGCCGGATGCGGGGTCAACGCCGCATCCTGCTACAAATCGTGCACTATATCAAACTTACTTAATCACCCCATTGCGCCGCAGTTGTTCGGTATCTTTTTTCCAGCGTTCATTAGCTTGGTTGATATCGTCAATTTTCATCAAAATCGCGACTTTTTTACTGTTTA  >  NC_000913/3098372‑3098667
                                                                                                                                                   |                                                                                                                                                    
 ttCGCTGCTGGAGCTTTTTCTGGCTGAGTAAACTTGTATTGCCGGATGCGGCGTGAACGCCTTATCCAGCCGACATGTGGCAGCGGTTGTAGGTCTGATAAGACGCGCAAGCGTCGCATCAGACGTTGATTGCCGGATGCGGGGTCAAc                                                                                                                                                    >  3:1006580/1‑149 (MQ=255)
    gctgctGGATCTTTTTCTGGCTGAGTAAACTTGTATTGCCGGATGCGGCGTGAACGCCTTATCCAGCCGACATGTGGCAGCGGTTGTAGGTCTGATAAGACGCGCAAGCGTCGCATCAGACGTTGATTGCCGGATGCGGGGTCAAcgcc                                                                                                                                                 <  2:1071788/149‑1 (MQ=255)
         tgGAGCTTTTTCTGGCTGAGTAAACTTGTATTGCCGGATGCGGCGTGAACGCCTTATCCAGCCGACATGTGGCAGCGGTTGTAGGTCTGATAAGACGCGCAAGCGTCGCATCAGACGTTGATTGCCGGATGCGGCGTTGACCCCGCATc                                                                                                                                            <  1:618024/149‑1 (MQ=255)
         tgGAGCTTTTTCTGGCTGAGTAAACTTGTATTGCCGGATGCGGCGTGAACGCCTTATCCAGCCGACATGTGGCAGCGGTTGTAGGTCTGATAAGACGCGCAAGCGTCGCATCAGACGTTGATTGCCGGATGCGGCGTTGACCCCGCATc                                                                                                                                            <  3:210043/149‑1 (MQ=255)
          ggAGCTTTTTCTGGCTGAGTAAACTTGTATTGCCGGATGCGGCGTGAACGCCTTATCCAGCCGACATGTGGCAGCGGTTGTAGGTCTGATAAGACGCGCAAGCGTCGCATCAGACGTTGATTGCCGGATGCGGGGTCAACGCCGCATcc                                                                                                                                           >  3:470389/1‑149 (MQ=255)
          ggAGCTTTTTCTGGCTGAGTAAACTTGTATTGCCGGATGCGGCGTGAACGCCTTATCCAGCCGACATGTGGCAGCGGTTGTAGGTCTGATAAGACGCGCAAGCGTCGCATCAGACGTTGATTGCCGGATGCGGCGTTGACCCCGCATcc                                                                                                                                           <  2:874433/149‑1 (MQ=255)
           gAGCTTTTTCTGGCTGAGTAAACTTGTATTGCCGGATGCGGCGTGAACGCCTTATCCAGCCGACATGTGGCAGCGGTTGTAGGTCTGATAAGACGCGCAAGCGTCGCATCAGACGTTGATTGCCGGATGCGGGGTCAACGCCGCATCCt                                                                                                                                          >  4:715253/1‑149 (MQ=255)
            aGCTTTTTCTGGCTGAGTAAACTTGTATTGCCGGATGCGGCGTGAACGCCTTATCCAGCCGACATGTGGCAGCGGTTGTAGGTCTGATAAGACGCGCAAGCGTCGCATCAGACGTTGATTGCCGGATGCGGGGTCAACGCCGCATCCTg                                                                                                                                         <  3:896587/149‑1 (MQ=255)
             gCTTTTTCTGGCTGAGTAAACTTGTATTGCCGGATGCGGCGTGAACGCCTTATCCAGCCGACATGTGGCAGCGGTTGTAGGTCTGATAAGACGCGCAAGCGTCGCATCAGACGTTGATTGCCGGATGCGGGGTCAACGCCGCATCCTGc                                                                                                                                        <  4:186601/149‑1 (MQ=255)
             gCTTTTTCTGGCTGAGTAAACTTGTATTGCCGGATGCGGCGTGAACGCCTTATCCAGCCGACATGTGGCAGCGGTTGTAGGTCTGATAAGACGCGCAAGCGTCGCATCAGACGTTGATTGCCGGATGCGGGGTCAACGCCGCATCCTGc                                                                                                                                        <  1:1431836/149‑1 (MQ=255)
               tttttCTGGCTGAGTAAACTTGTATTGCCGGATGCGGCGTGAACGCCTTATCCAGCCGACATGTGGCAGCGGTTGTAGGTCTGATAAGACGCGCAAGCGTCGCATCAGACGTTGATTGCCGGATGCGGCGTTGACCCCGCATc                                                                                                                                            >  2:383073/1‑143 (MQ=255)
               tttttCTGGCTGAGTAAACTTGTATTGCCGGATGCGGCGTGAACGCCTTATCCAGCCGACATGTGGCAGCGGTTGTAGGTCTGATAAGACGCGCAAGCGTCGCATCAGACGTTGATTGCCGGATGCGGCGTTGACCCCGCATc                                                                                                                                            <  1:383073/143‑1 (MQ=255)
                     tGGCTGAGTAAACTTGTATTGCCGGATGCGGCGTGAACGCCTTATCCAGCCGACATGTGGCAGCGGTTGTAGGTCTGATAAGACGCGCAAGCGTCGCATCAGACGTTGATTGCCGGATGCGGGGTCAACGCCGCATCCTGCTACAAATc                                                                                                                                <  2:240591/149‑1 (MQ=255)
                     tGGCTGAGTAAACTTGTATTGCCGGATGCGGCGTGAACGCCTTATCCAGCCGACATGTGGCAGCGGTTGTAGGTCTGATAAGACGCGCAAGCGTCGCATCAGACGTTGATTGCCGGATGCGGGGTCAACGCCGCATCCTGCTACAAATc                                                                                                                                <  2:483830/149‑1 (MQ=255)
                        cTGAGTAAACTTGTATTGCCGGATGCGGCGTGAACGCCTTATCCAGCCGACATGTGGCAGCGGTTGTAGGTCTGATAAGACGCGCAAGCGTCGCATCAGACGTTGATTGCCGGATGCGGGGTCAACGCCGCATCCTGCTACAAATCGTg                                                                                                                             <  4:829005/149‑1 (MQ=255)
                         tGAGTAAACTTGTATTGCCGGATGCGGCGTGAACGCCTTATCCAGCCGACATGTGGCAGCGGTTGTAGGTCTGATAAGACGCGCAAGCGTCGCATCAGACGTTGATTGCCGGATGCGGGGTCAACGCCGCATCCTGCTACAAATCGTGc                                                                                                                            >  3:149841/1‑149 (MQ=255)
                                 cTTGTATTGCCGGATGCGGCGTGAACGCCTTATCCAGCCGACATGTGGCAGCGGTTGTAGGTCTGATAAGACGCGCAAGCGTCGCATCAGACGTTGATTGCCGGATGCGGGGTCAACGCCGCATCCTGCTACAAATCGTGCACtata                                                                                                                      >  2:922586/1‑147 (MQ=255)
                                 cTTGTATTGCCGGATGCGGCGTGAACGCCTTATCCAGCCGACATGTGGCAGCGGTTGTAGGTCTGATAAGACGCGCAAGCGTCGCATCAGACGTTGATTGCCGGATGCGGGGTCAACGCCGCATCCTGCTACAAATCGTGCACtata                                                                                                                      <  1:922586/147‑1 (MQ=255)
                                 cTTGTATTGCCGGATGCGGCGTGAACGCCTTATCCAGCCGACATGTGGCAGCGGTTGTAGGTCTGATAAGACGCGCAAGCGTCGCATCAGACGTTGATTGCCGGATGCGGGGTCAACGCCGCATCCTGCTACAAATCGTGCACTATATc                                                                                                                    >  3:340231/1‑149 (MQ=255)
                                    gTATTGCCGGATGCGGCGTGAACGCCTTATCCAGCCGACATGTGGCAGCGGTTGTAGGTCTGATAAGACGCGCAAGCGTCGCATCAGACGTTGATTGCCGGATGCGGGGTCAACGCCGCATCCTGCTACAAATCGTGCACTATATCaaa                                                                                                                 <  2:1159434/149‑1 (MQ=255)
                                       ttGCCGGATGCGGCGTGAACGCCTTATCCAGCCGACATGTGGCAGCGGTTGTAGGTCTGATAAGACGCGCAAGCGTCGCATCAGACGTTGATTGCCGGATGCGGGGTCAACGCCGCATCCTGCTACAAATCGTGCACTATATCAAactt                                                                                                              >  4:635932/1‑149 (MQ=255)
                                        tGCCGGATGCGGCGTGAACGCCTTATCCAGCCGACATGTGGCAGCGGTTGTAGGTCTGATAAGACGCGCAAGCGTCGCATCAGACGTTGATTGCCGGATGCGGGGTCAACGCCGCATCCTGCTACAAATCGTGCACTATATCAAactta                                                                                                             >  2:598807/1‑149 (MQ=255)
                                           cGGATGCGGCGTGAACGCCTTATCCAGCCGACATGTGGCAGCGGTTGTAGGTCTGATAAGACGCGCAAGCGTCGCATCAGACGTTGATTGCCGGATGCGGGGTCAACGCCGCATCCTGCTACAAATCGTGCACTATATCAAacttactt                                                                                                          <  2:185041/149‑1 (MQ=255)
                                             gATGCGGCGTGAACGCCTTATCCAGCCGACATGTGGCAGCGGTTGTAGGTCTGATAAGACGCGCAAGCGTCGCATCAGACGTTGATTGCCGGATGCGGGGTCAACGCCGCATCCTGCTACAAATCGTGCACTATATCAAACTTACTTaa                                                                                                        <  1:1329794/149‑1 (MQ=255)
                                              aTGCGGCGTGAACGCCTTATCCAGCCGACATGTGGCAGCGGTTGTAGGTCTGATAAGACGCGCAAGCGTCGCATCAGACGTTGATTGCCGGATGCGGGGTCAACGCCGCATCCTGCTACAAATCGTGCACTATATCAAACTTACTTAAt                                                                                                       <  3:654221/149‑1 (MQ=255)
                                                 cggcgTGAACGCCTTATCCAGCCGACATGTGGCAGCGGTTGTAGGTCTGATAAGACGCGCAAGCGTCGCATCAGACGTTGATTGCCGGATGCGGGGTCAACGCCGCATCCTGCTACAAATCGTGCACTATATCAAACTTACTTAATCAc                                                                                                    <  3:859215/149‑1 (MQ=255)
                                                 cggcgTGAACGCCTTATCCAGCCGACATGTGGCAGCGGTTGTAGGTCTGATAAGACGCGCAAGCGTCGCATCAGACGTTGATTGCCGGATGCGGGGTCAACGCCGCATCCTGCTACAAATCGTGCACTATATCAAACTTACTTAATCAc                                                                                                    <  2:341414/149‑1 (MQ=255)
                                                  ggcgTGAACGCCTTATCCAGCCGACATGTGGCAGCGGTTGTAGGTCTGATAAGACGCGCAAGCGTCGCATCAGACGTTGATTGCCGGATGCGGGGTCAACGCCGCATCCTGCTACAAATCGTGCACTATATCAAACTTACTTAATCAcc                                                                                                   >  4:565011/1‑149 (MQ=255)
                                                      tGAACGCCTTATCCAGCCGACATGTGGCAGCGGTTGTAGGTCTGATAAGACGCGCAAGCGTCGCATCAGACGTTGATTGCCGGATGCGGGGTCAACGCCGCATCCTGCTACAAATCGTGCACTATATCAAACTTACTTAATCACCCCAt                                                                                               >  1:731607/1‑149 (MQ=255)
                                                       gAACGCCTTATCCAGCCGACATGTGGCAGCGGTTGTAGGTCTGATAAGACGCGCAAGCGTCGCATCAGACGTTGATTGCCGGATGCGGGGTCAACGCCGCATCCTGCTACAAATCGTGCACTATATCAAACTTACTTAATCACCCCAtt                                                                                              >  4:1427415/1‑149 (MQ=255)
                                                           gCCTTTTCCAGCCGACATGTGGCAGCGGTTGTAGGTCTGATAAGACGCGCAAGCGTCGCATAAGACGTTGATTGCCGGATGCGGGGTCAACGCGGCATCCTGCTACAAATCGTGCACTATATCAAACTTACTTAATAACCCCATTgcgc                                                                                          >  1:425440/1‑149 (MQ=255)
                                                           gCCTTATCCAGCCGACATGTGGCAGCGGTTGTAGGTCTGATAAGACGCGCAAGCGTCGCATCAGACGTTGATTGCCGGATGCGGGGTCAACGCCGCATCCTGCTACAAATCGTGCACTATATCAAACTTACTTAATCACCCCATTgcgc                                                                                          >  1:759575/1‑149 (MQ=255)
                                                            ccTTATCCAGCCGACATGTGGCAGCGGTTGTAGGTCTGATAAGACGCGCAAGCGTCGCATCAGACGTTGATTGCCGGATGCGGGGTCAACGCCGCATCCTGCTACAAATCGTGCACTATATCAAACTTACTTAATCACCCCAtt                                                                                              >  4:1202495/1‑144 (MQ=255)
                                                            ccTTATCCAGCCGACATGTGGCAGCGGTTGTAGGTCTGATAAGACGCGCAAGCGTCGCATCAGACGTTGATTGCCGGATGCGGGGTCAACGCCGCATCCTGCTACAAATCGTGCACTATATCAAACTTACTTAATCACCCCAtt                                                                                              <  3:1202495/144‑1 (MQ=255)
                                                               tATCCAGCCGACATGTGGCAGCGGTTGTAGGTCTGATAAGACGCGCAAGCGTCGCATCAGACGTTGATTGCCGGATGCGGGGTCAACGCCGCATCCTGCTACAAATCGTGCACTATATCAAACTTACTTAATCACCCCATTGCGCCGCa                                                                                      <  3:672394/149‑1 (MQ=255)
                                                                 tCCAGCCGACATGTGGCAGCGGTTGTAGGTCTGATAAGACGCGCAAGCGTCGCATCAGACGTTGATTGCCGGATGCGGGGTCAACGCCGCATCCTGCTACAAATCGTGCACTATATCAAACTTACTTAATCACCCCATTGCGCCGCAgt                                                                                    >  2:841733/1‑149 (MQ=255)
                                                                     gCCGACATGTGGCAGCGGTTGTAGGTCTGATAAGACGCGCAAGCGTCGCATCAGACGTTGATTGCCGGATGCGGGGTCAACGCCGCATCCTGCTACAAATCGTGCACTATATCAAACTTACTTAATCACCCCAtt                                                                                              >  1:580252/1‑135 (MQ=255)
                                                                     gCCGACATGTGGCAGCGGTTGTAGGTCTGATAAGACGCGCAAGCGTCGCATCAGACGTTGATTGCCGGATGCGGGGTCAACGCCGCATCCTGCTACAAATCGTGCACTATATCAAACTTACTTAATCACCCCAtt                                                                                              <  2:580252/135‑1 (MQ=255)
                                                                      ccGACATGTGGCAGCGGTTGTAGGTCTGATAAGACGCGCAAGCGTCGCATCAGACGTTGATTGCCGGATGCGGGGTCAACGCCGCATCCTGCTACAAATCGTGCACTATATCAAACTTACTTAATCACCCCATTGCGCCGCAGTTGTTc                                                                               <  1:63403/149‑1 (MQ=255)
                                                                      ccGACATGTGGCAGCGGTTGTAGGTCTGATAAGACGCGCAAGCGTCGCATCAGACGTTGATTGCCGGATGCGGGGTCAACGCCGCATCCTGCTACAAATCGTGCACTATATCAAACTTACTTAATCACCCCATTGCGCCGCAGTTGTTc                                                                               <  3:1427415/149‑1 (MQ=255)
                                                                      ccGACATGTGGCAGCGGTTGTAGGTCTGATAAGACGCGCAAGCGTCGCATCAGACGTTGATTGCCGGATGCGGGGTCAACGCCGCATCCTGCTACAAATCGTGCACTATATCAAACTTACTTAATCACCCCATTGCGCCGCAGTTGTTc                                                                               <  4:1053043/149‑1 (MQ=255)
                                                                      ccGACATGTGGCAGCGGTTGTAGGTCTGATAAGACGCGCAAGCGTCGCATCAGACGTTGATTGCCGGATGCGGGGTCAACGCCGCATCCTGCTACAAATCGTGCACTATATCAAACTTACTTAATCACCCCATTGCGCCGCAGTTGTTc                                                                               <  3:322593/149‑1 (MQ=255)
                                                                       cGACATGTGGCAGCGGTTGTAGGTCTGATAAGACGCGCAAGCGTCGCATCAGACGTTGATTGCCGGATGCGGGGTCAACGCCGCATCCTGCTACAAATCGTGCACTATATCAAACTTACTTAATCACCCCATTGCGCCGCAGTTGTTCg                                                                              <  1:1170119/149‑1 (MQ=255)
                                                                         aCATGTGGCAGCGGTTGTAGGTCTGATAAGACGCGCAAGCGTCGCATCAGACGTTGATTGCCGGATGCGGGGTCAACGCCGCATCCTGCTACAAATCGTGCACTATATCAAACTTACTTAATCACCCCATTGCGCCGCAGTTGTTCGGt                                                                            <  1:516447/149‑1 (MQ=255)
                                                                           aTGTGGCAGCGGTTGTAGGTCTGATAAGACGCGCAAGCGTCGCATCAGACGTTGATTGCCGGATGCGGGGTCAACGCCGCATCCTGCTACAAATCGTGCACTATATCAAACTTACTTAATCACCCCATTGCGCCGCAGTTGTTCGGTAt                                                                          >  3:443149/1‑149 (MQ=255)
                                                                           aTGTGGCAGCGGTTGTAGGTCTGATAAGACGCGCAAGCGTCGCATCAGACGTTGATTGCCGGATGCGGGGTCAACGCCGCATCCTGCTACAAATCGTGCACTATATCAAACTTACTTAATCACCCCATTGCGCCGCAGTTGTTCGGTAt                                                                          >  1:675987/1‑149 (MQ=255)
                                                                             gtgGCAGCGGTTGTAGGTCTGATAAGACGCGCAAGCGTCGCATCAGACGTTGATTGCCGGATGCGGGGTCAACGCCGCATCCTGCTACAAATCGTGCACTATATCAAACTTACTTAATCACCCCATTGCGCCGCAGTTGTTCGGTATCt                                                                        <  4:679334/149‑1 (MQ=255)
                                                                               ggCAGCGGTTGTAGGTCTGATAAGACGCGCAAGCGTCGCATCAGACGTTGATTGCCGGATGCGGGGTCAACGCCGCATCCTGCTACAAATCGTGCACTATATCAAACTTACTTAATCACCCCATTGCGCCGCAGTTGTTCGGTATCttt                                                                      >  2:76008/1‑149 (MQ=255)
                                                                                gCAGCGGTTGTAGGTCTGATAAGACGCGCAAGCGTCGCATCAGACGTTGATTGCCGGATGCGGGGTCAAcgcc                                                                                                                                                 <  2:766969/73‑1 (MQ=255)
                                                                                gCAGCGGTTGTAGGTCTGATAAGACGCGCAAGCGTCGCATCAGACGTTGATTGCCGGATGCGGGGTCAAcgcc                                                                                                                                                 >  1:766969/1‑73 (MQ=255)
                                                                                 cAGCGGTTGTAGGTCTGATAAGACGCGCAAGCGTCGCATCAGACGTTGATTGCCGGATGCGGGGTCAACGCCGCATCCTGCTACAAATCGTGCACTATATCAAACTTACTTAATCACCCCATTGCGCCGCAGTTGTTCGGTATCttttt                                                                    <  1:1458193/149‑1 (MQ=255)
                                                                                  aGCGGTTGTAGGTCTGATAAGACGCGCAAGCGTCGCATCAGACGTTGATTGACGGATGCGGGGTCAACGCCGCATCCTGCTACAAATCGTGCACTATATCAAACTTACTTAATCACCCCATTGCGCCGCAGTTGTTCGGTATCtttttt                                                                   <  2:1181840/149‑1 (MQ=255)
                                                                                   gCGGTTGTAGGTCTGATAAGACGCGCAAGCGTCGCATCAGACGTTGATTGCCGGATGCGGGGTCAACGCCGCATCCTGCTACAAATCGTGCACTATATCAAACTTACTTAATCa                                                                                                     <  4:1371099/114‑1 (MQ=255)
                                                                                               cTGATAAGACGCGCAAGCGTCGCATCAGACGTTGATTGCCGGATGCGGGGTCAACGCCGCATCCTGCTACAAATCGTGCACTATATCAAACTTACTTAATCACCCCATTGCGCCGCAGTTGTTCGGTATCTTTTTTCCAGCGTTCATTa                                                      <  2:754322/149‑1 (MQ=255)
                                                                                                      gACGCGCAAGCGTCGCATCAGACGTTGATTGCCGGATGCGGGGTCAACGCCGCATCCTGCTACAAATCGTGCACTATATCAAACTTACTTAATCACCCCATTGCGCCGCAGTTGTTCGGTATCTTTTTTCCAGCGTTCATTAGCTTGGt                                               >  1:288145/1‑149 (MQ=255)
                                                                                                      gACGCGCAAGCGTCGCATCAGACGTTGATTGCCGGATGCGGGGTCAACGCCGCATCCTGCTACAAATCGTGCACTATATCAAACTTACTTAATCACCCCATTGCGCCGCAGTTGTTCGGTATCTTTTTTCCAGCGTTCATTAGCTTGGt                                               <  4:1006580/149‑1 (MQ=255)
                                                                                                       aCGCGCAAGCGTCGCATCAGACGTTGATTGCCGGATGCGGGGTCAACGCCGCATCCTGCTACAAATCGTGCACTATATCAAACTTACTTAATCACCCCATTGCGCCGCAGTTGTTCGGTATCTTTTTTCCAGCGTTCATTAGCTTGGtt                                              <  3:608876/149‑1 (MQ=255)
                                                                                                       aCGCGCAAGCGTCGCATCAGACGTTGATTGCCGGATGCGGGGTCAACGCCGCATCCTGCTACAAATCGTGCACTATATCAAACTTACTTAATCACCCCATTGCGCCGCAGTTGTTCGGTATCTTTTTTCCAGCGTTCATTAGCTTGGtt                                              <  1:1430397/149‑1 (MQ=255)
                                                                                                       aCGCGCAAGCGTCGCATCAGACGTTGATTGCCGGATGCGGGGTCAACGCCGCATCCTGCTACAAATCGTGCACTATATCAAACTTACTTAATCACCCCATTGCGCCGCAGTTGTTCGGTATCTTTTTTCCAGCGTTCATTAGCTTGGt                                               <  2:376457/148‑1 (MQ=255)
                                                                                                        cgcgCAAGCGTCGCATCAGACGTTGATTGCCGGATGCGGGGTCAACGCCGCATCCTGCTACAAATCGTGCACTATATCAAACTTACTTAATCACCCCATTGCGCCGCAGTTGTTCGGTATCTTTTTTCCAGCGTTCATTAGCTTGGTTg                                             <  4:699686/149‑1 (MQ=255)
                                                                                                         gcgcAAGCGTCGCATCAGACGTTGATTGCCGGATGCGGGGTCAACGCCGCATCCTGCTACAAATCGTGCACTATATCAAACTTACTTAATCACCCCATTGCGCCGCAGTTGTTCGGTATCTTTTTTCCAGCGTTCATTAGCTTGGTTGa                                            <  4:564841/149‑1 (MQ=255)
                                                                                                               gCGTCGCATCAGACGTTGATTGCCGGATGCGGGGTCAACGCCGCATCCTGCTACAAATCGTGCACTATATCAAACTTACTTAATCACCCCATTGCGCCGCAGTTGTTCGGTATCTTTTTTCCAGCGTTCATTAGCTTGGTTGATATCGt                                      <  1:1091279/149‑1 (MQ=255)
                                                                                                                  tCGCATCAGACGTTGATTGCCGGATGCGGGGTCAACGCCGCATCCTGCTACAAATCGTGCACTATATCAAACTTACTTAATCACCCCATTGCGCCGCAGTTGTTCGGTATCTTTTTTCCAGCGTTCATTAGCTTGGt                                               <  4:935630/137‑1 (MQ=255)
                                                                                                                  tCGCATCAGACGTTGATTGCCGGATGCGGGGTCAACGCCGCATCCTGCTACAAATCGTGCACTATATCAAACTTACTTAATCACCCCATTGCGCCGCAGTTGTTCGGTATCTTTTTTCCAGCGTTCATTAGCTTGGTTGATATCGTCaa                                   <  2:884240/149‑1 (MQ=255)
                                                                                                                  tCGCATCAGACGTTGATTGCCGGATGCGGGGTCAACGCCGCATCCTGCTACAAATCGTGCACTATATCAAACTTACTTAATCACCCCATTGCGCCGCAGTTGTTCGGTATCTTTTTTCCAGCGTTCATTAGCTTGGTTGATATCGTCaa                                   <  2:533394/149‑1 (MQ=255)
                                                                                                                  tCGCATCAGACGTTGATTGCCGGATGCGGGGTCAACGCCGCATCCTGCTACAAATCGTGCACTATATCAAACTTACTTAATCACCCCATTGCGCCGCAGTTGTTCGGTATCTTTTTTCCAGCGTTCATTAGCTTGGTTGATATCGTCaa                                   <  3:772569/149‑1 (MQ=255)
                                                                                                                       tCAGACGTTGATTGCCGGATGCGGGGTCAACGCCGCATCCTGCTACAAATCGTGCACTATATCAAACTTACTTAATCACCCCATTGCGCCGCAGTTGTTCGGTATCTTTTTTCCAGCGTTCATTAGCTTGGTTGATATCGTCAATTTtc                              >  4:876605/1‑149 (MQ=255)
                                                                                                                         aGACGTTGATTGCCGGATGCGGGGTCAACGCCGCATCCTGCTACAAATCGTGCACTATATCAAACTTACTTAATCACCCCATTGCGCCGCAGTTGTTCGGTATCTTTTTTCCAGCGTTCATTAGCTTGGTTGATATCGTCAATTTtcat                            <  4:1370054/149‑1 (MQ=255)
                                                                                                                          gACGTTGATTGCCGGATGCGGGGTCAACGCCGCATCCTGCTACAAATCGTGCACTATATCAAacttactt                                                                                                          >  2:32571/1‑70 (MQ=255)
                                                                                                                          gACGTTGATTGCCGGATGCGGGGTCAACGCCGCATCCTGCTACAAATCGTGCACTATATCAAacttactt                                                                                                          <  1:32571/70‑1 (MQ=255)
                                                                                                                          gACGTTGATTGCCGGATGCGGGGTCAACGCCGCATCCTGCTACAAATCGTGCACTATATCAAACTTACTTAATCACCCCATTGCGCCGCAGTTGTTCGGTATCTTTTTTCCAGCGTTCATTAGCTTGGTTGATATCGTCAATTTtcatc                           >  2:1120600/1‑149 (MQ=255)
                                                                                                                          gACGTTGATTGCCGGATGCGGGGTCAACGCCGCATCCTGCTACAAATCGTGCACTATATCAAACTTACTTAATCACCCCATTGCGCCGCAGTTGTTCGGTATCTTTTTTCCAGCGTTCATTAGCTTGGTTGATATCGTCAATTTtcatc                           <  2:84808/149‑1 (MQ=255)
                                                                                                                           aCGTTGATTGCCGGATGCGGGGTCAACGCCGCATCCTGCTACAAATCGTGCACTATATCAAACTTACTTAATCACCCCATTGCGCCGCAGTTGTTCGGTATCTTTTTTCCAGCGTTCATTAGCTTGGTTGatat                                         >  3:58163/1‑134 (MQ=255)
                                                                                                                           aCGTTGATTGCCGGATGCGGGGTCAACGCCGCATCCTGCTACAAATCGTGCACTATATCAAACTTACTTAATCACCCCATTGCGCCGCAGTTGTTCGGTATCTTTTTTCCAGCGTTCATTAGCTTGGTTGatat                                         <  4:58163/134‑1 (MQ=255)
                                                                                                                             gTTGATTGCCGGATGCGGGGTCAACGCCGCATCCTGCTACAAATCGTGCACTATATCAAACTTACTTAATCACCCCATTGCGCCGCAGTTGTTCGGTATCTTTTTTCCAGCGTTCATTAGCTTGGTTGATATCGTCAATTTTCATCaaa                        >  3:35060/1‑149 (MQ=255)
                                                                                                                                gATTGCCGGATGCGGGGTCAACGCCGCATCCTGCTACAAATCGTGCACTATATCAAACTTACTTAATCACCCCATTGCGCCGCAGTTGTTCGGTATCTTTTTTCCAGCGTTCATTAGCTTGGTTGATATCGTCAATTTTCATCAAAATc                     >  3:835918/1‑149 (MQ=255)
                                                                                                                                gATTGCCGGATGCGGGGTCAACGCCGCATCCTGCTACAAATCGTGCACTATATCAAACTTACTTAATCACCCCATTGCGCCGCAGTTGTTCGGTATCTTTTTTCCAGCGTTCATTAGCTTGGTTGATATCGTCAATTTTCATCAAAATc                     >  3:301990/1‑149 (MQ=255)
                                                                                                                                  ttGCCGGATGCGGGGTCAACGCCGCATCCTGCTACAAATCGTGCACTATATCAAACTTACTTAATCACCCCATTGCGCCGCAGTTGTTCGGTATCTTTTTTCCAGCGTTCATTAGCTTGGTTGATATCGTCAATTTTCATCAAAATcgc                   <  4:470389/149‑1 (MQ=255)
                                                                                                                                  ttGCCGGATGCGGGGTCAACGCCGCATCCTGCTACAAATCGTGCACTATATCAAACTTACTTAATCACCCCATTGCGCCGCAGTTGTTCGGTATCTTTTTTCCAGCGTTCATTAGCTTGGTTGATATCGTCAATTTTCATCAAAATcgc                   <  1:1155139/149‑1 (MQ=255)
                                                                                                                                     ccGGATGCGGGGTCAACGCCGCATCCTGCTACAAATCGTGCACTATATCAAACTTACTTAATCACCCCATTGCGCCGCAGTTGTTCGGTATCTTTTTTCCAGCGTTCATTAGCTTGGTTGATATCGTCAATTTTCATCAAAATCGCGAc                >  3:623707/1‑149 (MQ=255)
                                                                                                                                      cGGATGCGGGGTCAACGCCGCATCCTGCTACAAATCGTGCACTATATCAAACTTACTTAATCACCCCATTGCGCCGCAGTTGTTCGGTATCTTTTTTCCAGCGTTCATTAGCTTGGt                                               <  4:1317063/117‑1 (MQ=255)
                                                                                                                                      cGGATGCGGGGTCAACGCCGCATCCTGCTACAAATCGTGCACTATATCAAACTTACTTAATCACCCCATTGCGCCGCAGTTGTTCGGTATCTTTTTTCCAGCGTTCATTAGCTTGGt                                               >  3:1317063/1‑117 (MQ=255)
                                                                                                                                      cGGATGCGGGGTCAACGCCGCATCCTGCTACAAATCGTGCACTATATCAAACTTACTTAATCACCCCATTGCGCCGCAGTTGTTCGGTATCTTTTTTCCAGCGTTCATTAGCTTGGTTGATATCGTCAATTTTCATCAAAATCGCGACt               >  1:204205/1‑149 (MQ=255)
                                                                                                                                      cGGATGCGGGGTCAACGCCGCATCCTGCTACAAATCGTGCACTATATCAAACTTACTTAATCACCCCATTGCGCCGCAGTTGTTCGGTATCTTTTTTCCAGCGTTCATTAGCTTGGTTGATATCGTCAATTTTCATCAAAATCGCGACt               >  1:311568/1‑149 (MQ=255)
                                                                                                                                      cGGATGCGGGGGCAACGCCGCATCCTGCTACAAATCGTGCACTATATCAAACTTACTTAATCACCCCATTGCGCCGCAGTTGTTCGGTATCTTTTTTCCAGCGTTCATTAGCTTGGTTGATATCGTCAATTTTCATCAAAATCGCGACt               >  3:1145597/1‑149 (MQ=255)
                                                                                                                                      aggATGCGGCGTTGACCCCGCATCCTGCTACAAATCGTGCACTATATCAAACTTACTTAATCACCCCATTGCGCCGCAGTTGTTCGGTATCTTTTTTCCAGCGTTCATTAGCTTGGTTGATATCGTCAATTTTCATCAAAATCGCGACt               >  2:608688/2‑149 (MQ=255)
                                                                                                                                        gATGCGGGGTCAACGCCGCATCCTGCTACAAATCGTGCACTATATCAAACTTACTTAATCACCCCATTGCGCCGCAGTTGTTCGGTATCTTTTTTCCAGCGTTCATTAGCTTGGTTGATATCGTCAATTTTCATCAAAATCGCGACttt             >  1:708425/1‑149 (MQ=255)
                                                                                                                                        gATGCGGGGTCAACGCCGCATCCTGCTACAAATCGTGCACTATATCAAACTTACTTAATCACCCCATTGCGCCGCAGTTGTTCGGTATCTTTTTTCCAGCGTTCATTAGCTTGGTTGATATCGTCAATTTTCATCAAAATCGCGACttt             >  3:1188297/1‑149 (MQ=255)
                                                                                                                                        gATGCGGCGTTGACCCCGCATCCTGCTACAAATCGTGCACTATATCAAACTTACTTAATCACCCCATTGCGCCGCAGTTGTTCGGTATCTTTTTTCCAGCGTTCATTAGCTTGGTTGATATc                                        >  4:1115117/1‑122 (MQ=255)
                                                                                                                                          tGCGGGGTCAACGCCGCATCCTGCTACAAATCGTGCACTATATCAAACTTACTTAATCACCCCATTGCGCCGCAGTTGTTCGGTATCTTTTTTCCAGCGTTCATTAGCTTGGTTGATATCGTCAATTTTCATCAAAATCGCGACttttt           <  4:425801/149‑1 (MQ=255)
                                                                                                                                             ggggTCAACGCCGCATCCTGCTACAAATCGTGCACTATATCAAactta                                                                                                             <  4:505000/48‑1 (MQ=255)
                                                                                                                                             ggggTCAACGCCGCATCCTGCTACAAATCGTGCACTATATCAAactta                                                                                                             >  3:505000/1‑48 (MQ=255)
                                                                                                                                             ggggTCAACGCCGCATCCTGCTACAAATCGTGCACTATATCAAACTTACTTAATCACCCCATTGCGCCGCAGTTGTTCGGTATCTTTTTTCCAGCGTTCATTAGCTTGGTTGATATCGTCAATTTTCATCAAAATCGCGACTTTTTTAc        >  2:1134801/1‑149 (MQ=255)
                                                                                                                                              gggTCAACGCCGCATCCTGCTACAAATCGTGCACTATATCAAACTTACTTAATCACCCCATTGCGCCGCAGTTGTTCGGTATCTTTTTTCCAGCGTTCATTAGCTTGGTTGATATCGTCAATTTTCATCAAAATCGCGACTTTTTTACt       >  4:1448588/1‑149 (MQ=255)
                                                                                                                                                gTCAACGCCGCATCCTGCTACAAATCGTGCACTATATCAAACTTACTTAATCACCCCATTGCGCCGCAGTTGTTCGGTATCTTTTTTCCAGCGTTCATTAGCTTGGTTGATATCGTCAATTTTCATCAAAATCGCGACTTTTTTACTGt     <  4:314948/149‑1 (MQ=255)
                                                                                                                                                 tCAACGCCGCATCCTGCTACAAATCGTGCACTATATCAAACTTACTTAATCACCCCATTGCGCCGCAGTTGTTCGGTATCTTTTTTCCAGCGTTCATTAGCTTGGTTGATATCGTCAATTTTCATCAAAATCGCGACTTTTTTACTGtt    >  4:354401/1‑149 (MQ=255)
                                                                                                                                                  cAACGCCGCATCCTGCTACAAATCGTGCACTATATCAAACTTACTTAATCACCCCATTGCGCCGCAGTTGTTc                                                                               >  2:937657/1‑73 (MQ=255)
                                                                                                                                                  cAACGCCGCATCCTGCTACAAATCGTGCACTATATCAAACTTACTTAATCACCCCATTGCGCCGCAGTTGTTc                                                                               <  1:937657/73‑1 (MQ=255)
                                                                                                                                                   aaCGCCGCATCCTGCTACAAATCGTGCACTATATCAAACTTACTTAATCACCCCATTGCGCCGCAGTTGTTCGGTATCTTTTTTCCAGCGTTCATTAGCTTGGTTGATATCGTCAATTTTCATCAAAATCGCGACTTTTTTACTGTTTa  >  3:457669/1‑149 (MQ=255)
                                                                                                                                                   aaCGCCGCATCCTGCTACAAATCGTGCACTATATCAAACTTACTTAATCACCCCATTGCGCCGCAGTTGTTCGGTATCTTTTTTCCAGCGTTCATTAGCTTGGTTGATATCGTCAATTTTCATCAAAATCGCGACTTTTTTACTGTTTa  >  3:861705/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
ATTCGCTGCTGGAGCTTTTTCTGGCTGAGTAAACTTGTATTGCCGGATGCGGCGTGAACGCCTTATCCAGCCGACATGTGGCAGCGGTTGTAGGTCTGATAAGACGCGCAAGCGTCGCATCAGACGTTGATTGCCGGATGCGGGGTCAACGCCGCATCCTGCTACAAATCGTGCACTATATCAAACTTACTTAATCACCCCATTGCGCCGCAGTTGTTCGGTATCTTTTTTCCAGCGTTCATTAGCTTGGTTGATATCGTCAATTTTCATCAAAATCGCGACTTTTTTACTGTTTA  >  NC_000913/3098372‑3098667

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