Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 644,030 A→G 5.8% intergenic (‑63/+167) rnk ← / ← rna regulator of nucleoside diphosphate kinase/ribonuclease I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913644,0300AG 5.8% 490.8 / 14.3 156intergenic (‑63/+167)rnk/rnaregulator of nucleoside diphosphate kinase/ribonuclease I
Reads supporting (aligned to +/- strand):  ref base A (82/64);  new base G (3/6);  total (85/71)
Fisher's exact test for biased strand distribution p-value = 3.01e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.78e-01

Only 100 of 157 total aligned reads displayed.

GCAAACCAGCATAGGCGGGTTGCTCCAGCAGAATATCGATGCGTTCGGCATCCAGGTCGTTAATGATGATAGTTGGTCTGGACATTTTTACTCCATGTCGTCGGTGCTGCGAGTGTCGCAGATAAACATACCCAAAAGAAAACCCTCACCGTCAGGCGGCGAGGGTTTAACTCACATGATGATACTGACTGTTGCTCACTCTTTGAAGTGATTTGCGTCACATTCAGGGAATTCCTCAATGCAGCGCATTATGTATAAATCTTAATCGCCTTGGTTTATGGAAGACGAATTAGCGT  >  NC_000913/643883‑644178
                                                                                                                                                   |                                                                                                                                                    
gCAAACCAGCATAGGCGGGTTGCTCCAGCAGAATATCGATGCGTTCGGCATCCAGGTCGTTAATGATGATAGTTGGTCTGGACATTTTTACTCCATGTCGTCGGTGCTGCGAGTGTCGCAGATAAACATACCCAAAAGAAAACCCTCAc                                                                                                                                                     <  4:395744/149‑1 (MQ=255)
gCAAACCAGCATAGGCGGGTTGCTCCAGCAGAATATCGATGCGTTCGGCATCCAGGTCGTTAATGATGATAGTTGGTCTGGACATTTTTACTCCATGTCGTCGGTGCTGCGAGTGTCGCAGATAAACATACCCAAAAGAAAACCCTCAc                                                                                                                                                     >  3:79314/1‑149 (MQ=255)
  aaaCCAGCATAGGCGGGTTGCTCCAGCAGAATATCGATGCGTTCGGCATCCAGGTCGTTAATGATGATAGTTGGTCTGGACATTTTTACTCCATGTCGTCGGTGCTGCGAGTGTCGCAGATAAACATACCCAAAAGAAAACCCTCACCg                                                                                                                                                   >  1:957836/1‑149 (MQ=255)
  aaaCCAGCATAGGCGGGTTGCTCCAGCAGAATATCGATGCGTTCGGCATCCAGGTCGTTAATGATGATAGTTGGTCTGGACATTTTTACTCCATGTCGTCGGTGCTGCGAGTGTCGCAGATAAACATACCCAAAAGAAAACCCTCACCg                                                                                                                                                   >  1:338350/1‑149 (MQ=255)
  aaaCCAGCATAGGCGGGTTGCTCCAGCAGAATATCGATGCGTTCGGCATCCAGGTCGTTAATGATGATAGTTGGTCTGGACATTTTTACTCCATGTCGTCGGTGCTGCGAGTGGCGCAGATAAACATACCCAAAAGAAAACCCTCACCg                                                                                                                                                   >  3:739911/1‑149 (MQ=255)
   aaCCAGCATAGGCGGGTTGCTCCAGCAGAATATCGATGCGTTCGGCATCCAGGTCGTTAATGATGATAGTTGGTCTGGACATTTTTACTCCATGTCGTCGGTGCTGCGAGTGTCGCAGATAAACATACCCAAAAGAAAACCCTCACCGt                                                                                                                                                  >  3:302711/1‑149 (MQ=255)
    aCCAGCATAGGCGGGTTGCTCCAGCAGAATATCGATGCGTTCGGCATCCAGGTCGTTAATGATGATAGTTGGTCTGGACATTTTTACTCCATGTCGTCGGTGCTGCGAGTGTCGCAGATAAACATACCCAAAAGAAAACCCTCACCGTc                                                                                                                                                 >  3:393763/1‑149 (MQ=255)
     ccAGCATAGGCGGGTTGCTCCAGCAGAATATCGATGCGTTCGGCATCCAGGTCGTTAATGATGATAGTTGGTCTGGACATTTTTACTCCATGTCGTCGGTGCTGCGAGTGTCGCAGATAAACATACCCAAAAGAAAACCCTCACCGTCa                                                                                                                                                <  3:6302/149‑1 (MQ=255)
        gCATAGGCGGGTTGCTCCAGCAGAATATCGATGCGTTCGGCATCCAGGTCGTTAATGATGATAGTTGGTCTGGACATTTTTACTCCATGTCGTCGGTGCTGCGAGTGTCGCAGATAAACATACCCAAAAGAAAACCCTCACCGTCAggc                                                                                                                                             <  4:378176/149‑1 (MQ=255)
         cATAGGCGGGTTGCTCCAGCAGAATATCGATGCGTTCGGCATCCAGGTCGTTAATGATGATAGTTGGTCTGGACATTTTTACTCCATGTCGTCGGTGCTGCGAGTGTCGCAGATAAACATACCCAAAAGAAAACCCTCACCGTCAggcg                                                                                                                                            <  2:950842/149‑1 (MQ=255)
               cGGGTTGCTCCAGCAGAATATCGATGCGTTCGGCATCCAGGTCGTTAATGATGATAGTTGGTCTGGACATTTTTACTCCATGTCGTCGGTGCTGCGAGTGTCGCAGATAAACATACCCAAAAGAAAACCCTCGCCGCCTGACGGTGAgg                                                                                                                                      <  4:71689/149‑1 (MQ=255)
               cGGGTTGCTCCAGCAGAATATCGATGCGTTCGGCATCCAGGTCGTTAATGATGATAGTTGGTCTGGACATTTTTACTCCATGTCGTCGGTGCTGCGAGTGTCGCAGATAAACATACCCAAAAGAAAACCCTCACCGTCAGGCGGCGAgg                                                                                                                                      <  4:1073563/149‑1 (MQ=255)
               cGGGTTGCTCCAGCAGAATATCGATGCGTTCGGCATCCAGGTCGTTAATGATGATAGTTGGTCTGGACATTTTTACTCCATGTCGTCGGTGCTGCGAGTGTCGCAGATAAACATACCCAAAAGAAAACCCTCACCGTCAGGCGGCGAgg                                                                                                                                      <  2:404014/149‑1 (MQ=255)
                gggTTGCTCCAGCAGAATATCGATGCGTTCGGCATCCAGGTCGTTAATGATGATAGTTGGTCTGGACATTTTTACTCCATGTCGTCGGTGCTGCGAGTGTCGCAGATAAACATACCCAAAAGAAAACCCTCGCCGCCTGACGGTGAggg                                                                                                                                     <  3:641461/149‑1 (MQ=255)
                   ttGCTCCAGCAGAATATCGATGCGTTCGGCATCCAGGTCGTTAATGATGATAGTTGGTCTGGACATTTTTACTCCATGTCGTCGGTGCTGCGAGTGTCGCAGATAAACATACCCAAAAGAAAACCCTCGCCGCCTGACGGTGAggg                                                                                                                                     <  2:587837/146‑1 (MQ=255)
                   ttGCTCCAGCAGAATATCGATGCGTTCGGCATCCAGGTCGTTAATGATGATAGTTGGTCTGGACATTTTTACTCCATGTCGTCGGTGCTGCGAGTGTCGCAGATAAACATACCCAAAAGAAAACCCTCGCCGCCTGACGGTGAggg                                                                                                                                     >  3:203299/1‑146 (MQ=255)
                   ttGCTCCAGCAGAATATCGATGCGTTCGGCATCCAGGTCGTTAATGATGATAGTTGGTCTGGACATTTTTACTCCATGTCGTCGGTGCTGCGAGTGTCGCAGATAAACATACCCAAAAGAAAACCCTCGCCGCCTGACGGTGAggg                                                                                                                                     >  1:587837/1‑146 (MQ=255)
                   ttGCTCCAGCAGAATATCGATGCGTTCGGCATCCAGGTCGTTAATGATGATAGTTGGTCTGGACATTTTTACTCCATGTCGTCGGTGCTGCGAGTGTCGCAGATAAACATACCCAAAAGAAAACCCTCGCCGCCTGACGGTGAggg                                                                                                                                     <  4:203299/146‑1 (MQ=255)
                       tCCAGCAGAATATCGATGCGTTCGGCATCCAGGTCGTTAATGATGATAGTTGGTCTGGACATTTTTACTCCATGTCGTCGGTGCTGCGAGTGTCGCAGATAAACATACCCAAAAGAAAACCCTCACCGTCAGGCGGCGAGGGTTTAACt                                                                                                                              >  2:828824/1‑149 (MQ=255)
                       tCCAGCAGAATATCGATGCGTTCGGCATCCAGGTCGTTAATGATGATAGTTGGTCTGGACATTTTTACTCCATGTCGTCGGTGCTGCGAGTGTCGCAGATAAACATACCCAAAAGAAAACCCTCACCGTCAGGCGGCGAGGGTGTAACt                                                                                                                              >  1:1008244/1‑149 (MQ=255)
                          agcagAATATCGATGCGTTCGGCATCCAGGTCGTTAATGATGATAGTTGGTCTGGACATTTTTACTCCATGTCGTCGGTGCTGCGAGTGTCGCAGATAAACATACCCAAAAGAAAACCCTCACCGTCAGGCGGCGAGGGTTTAACTcac                                                                                                                           >  3:26226/1‑149 (MQ=255)
                            caggatatCGATGCGTTCGGCATCCAGGTCGTTAATGATGATAGTTGGTCTGGACATTTTTACTCCATGTCGTCGGTGCTGCGAGTGTCGCAGATAAACATACCCAAAAGAAAACCCTCACCGTCAGGCGGCGAGGGTTTAACTCACat                                                                                                                         <  1:154861/145‑1 (MQ=255)
                            cagAATATCGATGCGTTCGGCATCCAGGTCGTTAATGATGATAGTTGGTCTGGACATTTTTACTCCATGTCGTCGGTGCTGCGAGTGTCGCAGATAAACATACCCAAAAGAAAACCCTCACCGTCAGGCGGCGAGGGTTTAACTCACat                                                                                                                         >  2:1015678/1‑149 (MQ=255)
                            cagAATATCGATGCGTTCGGCATCCAGGTCGTTAATGATGATAGTTGGTCTGGACATTTTTACTCCATGTCGTCGGTGCTGCGAGTGTCGCAGATAAACATACCCAAAAGAAAACCCTCACCGTCAGGCGGCGAGGGTTTAACTCACat                                                                                                                         <  1:289283/149‑1 (MQ=255)
                              gAATATCGATGCGTTCGGCATCCAGGTCGTTAATGATGATAGTTGGTCTGGACATTTTTACTCCATGTCGTCGGTGCTGCGAGTGTCGCAGATAAACATACCCAAAAGAAAACCCTCACCGTCAGGCGGCGAGGGTTTAACTCACatga                                                                                                                       >  4:420616/1‑149 (MQ=255)
                               aaTATCGATGCGTTCGGCATCCAGGTCGTTAATGATGATAGTTGGTCTGGACATTTTTACTCCATGTCGTCGGTGCTGCGAGTGTCGCAGATAAACATACCCAAAAGAAAACCCTCACCGTCAGGCGGCGAGGGTTTAACTCACatgat                                                                                                                      <  4:3983/149‑1 (MQ=255)
                                 tatCGATGCGTTCGGCATCCAGGTCGTTAATGATGATAGTTGGTCTGGACATTTTTACTCCATGTCGTCGGTGCTGCGAGTGTCGCAGATAAACATACCCAAAAGAAAACCCTCACCGTCAGGCGGCGAGGGTTTAACTCACAtgatga                                                                                                                    >  4:444274/1‑149 (MQ=255)
                                  atCGATGCGTTCGGCATCCAGGTCGTTAATGATGATAGTTGGTCTGGACATTTTTACTCCATGTCGTCGGTGCTGCGAGTGTCGCAGATAAACATACCCAAAAGAAAACCCTCACCGTCAGGCGGCGAGGGTTTAACTCACATgatgat                                                                                                                   >  1:868772/1‑149 (MQ=255)
                                    cGATGCGTTCGGCATCCAGGTCGTTAATGATGATAGTTGGTCTGGACATTTTTACTCCATGTCGTCGGTGCTGCGAGTGTCGCAGATAAACATACCCAAAAGAAAACCCTCACCGTCAGGCGGCGAGGGTTTAACTCACATGATGATac                                                                                                                 <  4:26226/149‑1 (MQ=255)
                                     gATGCGTTCGGCATCCAGGTCGTTAATGATGATAGTTGGTCTGGACATTTTTACTCCATGTCGTCGGTGCTGCGAGTGTCGCAGATAAACATACCCAAAAGAAAACCCTCACCGTCAGGCGGCGAGGGTTTAACTCACATGATGATact                                                                                                                >  4:757615/1‑149 (MQ=255)
                                     gATGCGTTCGGCATCCAGGTCGTTAATGATGATAGTTGGTCTGGACATTTTTACTCCATGTCGTCGGTGCTGCGAGTGTCGCAGATAAACATACCCAAAAGAAAACCCTCACCGTCAGGCGGCGAGGGTTTAACTCACATGATGATact                                                                                                                >  2:342569/1‑149 (MQ=255)
                                      aTGCGTTCGGCATCCAGGTCGTTAATGATGATAGTTGGTCTGGACATTTTTACTCCATGTCGTCGGTGCTGCGAGTGTCGCAGATAAACATACCCAAAAGAAAACCCTCACCGTCAGGCGGCGAGGGTTTAACTCACATGATGATactg                                                                                                               >  2:318888/1‑149 (MQ=255)
                                       tGCGTTCGGCATCCAGGTCGTTAATGATGATAGTTGGTCTGGACATTTTTACTCCATGTCGTCGGTGCTGCGAGTGTCGCAGATAAACATACCCAAAAGAAAACCCTCACCGTCAGGCGGCGAGGGTTTAACTCACATGATGATactga                                                                                                              >  1:423661/1‑149 (MQ=255)
                                           ttCGGCATCCAGGTCGTTAATGATGATAGTTGGTCTGGACATTTTTACTCCATGTCGTCGGTGCTGCGAGTGTCGCAGATAAACATACCCAAAAGAAAACCCTCACCGTCAGGCGGCGAGGGTTTAACTCACATGATGATACTGACTGt                                                                                                          >  2:165644/1‑149 (MQ=255)
                                               gCATCCAGGTCGTTAATGATGATAGTTGGTCTGGACATTTTTACTCCATGTCGTCGGTGCTGCGAGTGTCGCAGATAAACATACCCAAAAGAAAACCCTCACCGTCAGGCGGCGAGGGTTTAACTCACATGATGATACTGACTGTTGCt                                                                                                      >  3:981964/1‑149 (MQ=255)
                                                 aTCCAGGTCGTTAATGATGATAGTTGGTCTGGACATTTTTACTCCATGTCGTCGGTGCTGCGAGTGTCGCAGATAAACATACCCAAAAGAAAACCCTCACCGTCAGGCGGCGAGGGTTTAACTCACATGATGATACTGACTGTTGCTCa                                                                                                    >  1:222930/1‑149 (MQ=255)
                                                   ccAGGTCGTTAATGATGATAGTTGGTCTGGACATTTTTACTCCATGTCGTCGGTGCTGCGAGTGTCGCAGATAAACATACCCAAAAGAAAACCCTCACCGTCAGGCGGCGAGGGTTTAACTCACATGATGATACTGACTGTTGCTCAct                                                                                                  <  1:504202/149‑1 (MQ=255)
                                                         cGTTAATGATGATAGTTGGTCTGGACATTTTTACTCCATGTCGTCGGTGCTGCGAGTGTCGCAGATAAACATACCCAAAAGAAAACCCTCACCGTCAGGCGGCGAGGGTTTAACTCACATGATGATACTGACTGTTGCTCACTCTTTGa                                                                                            <  2:798446/149‑1 (MQ=255)
                                                          gTTAATGATGATAGTTGGTCTGGACATTTTTACTCCATGTCGTCGGTGCTGCGAGTGTCGCAGATAAACATACCCAAAAGAAAACCCTCACCGTCAGGCGGCGAGGGTTTAACTCACATGATGATACTGACTGTTGCTCACTCTTTGaa                                                                                           <  4:335481/149‑1 (MQ=255)
                                                                gatgatAGTTGGTCTGGACATTTTTACTCCATGTCGTCGGTGCTGCGAGTGTCGCAGATAGACATACCCAAAAGAAAACCCTCACCGTCAGGCGGCGAGGGTTTAACTCACATGATGATACTGACTGTTGCTCACTCTTTGAAGTGAtt                                                                                     <  1:1054692/149‑1 (MQ=255)
                                                                 atgatAGTTGGTCTGGACATTTTTACTCCATGTCGTCGGTGCTGCGAGTGTCGCAGATAAACATACCCAAAAGAAAACCCTCACCGTCAGGCGGCGAGGGTTTAACTCACATGATGATACTGACTGTTGCTCACTCTTTGAAGTGAttt                                                                                    <  3:477678/149‑1 (MQ=255)
                                                                   gatAGTTGGTCTGGACATTTTTACTCCATGTCGTCGGTGCTGCGAGTGTCGCAGATAAACATACCCAAAAGAAAACCCTCACCGTCAGGCGGCGAGGGTTTAACTCACATGATGATACTGACTGTTGCTCACTCTTTGAAGTGATTTGc                                                                                  >  3:1071703/1‑149 (MQ=255)
                                                                       gTTGGTCTGGACATTTTTACTCCATGTCGTCGGTGCTGCGAGTGTCGCAGATAAACATACCCAAAAGAAAACCCTCACCGTCAGGCGGCGAGGGTTTAACTCACATGATGATACTGACTGTTGCTCACTCTTTGAAGTGATTTGCGTca                                                                              >  3:366495/1‑149 (MQ=255)
                                                                        ttGGTCTGGACATTTTTACTCCATGTCGTCGGTGCTGCGAGTGTCGCAGATAAACATACCCAAAAGAAAACCCTCACCGTCAGGCGGCGAGGGTTTAACTCACATGATGATACTGACTGTTGCTCACTCTTTGAAGTGATTTGCGTcac                                                                             >  3:1003623/1‑149 (MQ=255)
                                                                          ggTCTGGACATTTTTACTCCATGTCGTCGGTGCTGCGAGTGTCGCAGATAAACATACCCAAAAGAAAACCCTCACCGTCAGGCGGCGAGGGTTTAACTCACATGATGATACTGACTGTTGCTCACTCTTTGAAGTGATTTGCGTCACAt                                                                           <  2:1022841/149‑1 (MQ=255)
                                                                            tCTGGACATTTTTACTCCATGTCGTCGGTGCTGCGAGTGTCGCAGATAAACATACACAAAAGAAAACCCTCACCGTCAGGCGGCGAGGGGTTAACTCACATGATGATACTGACTGTTGCTCACTCTTTGAAGTGATTTGCGTCACATTc                                                                         >  2:295998/1‑149 (MQ=255)
                                                                             cTGGACATTTTTACTCCATGTCGTCGGTGCTGCGAGTGTCGCAGATAAACATACCCAAAAGAAAACCCTCACCGTCAGGCGGCGAGGGTTTAACTCACATGATGATACTGACTGTTGCTCACTCTTTGAAGTGATTTGCGTCACATTCa                                                                        >  1:377767/1‑149 (MQ=255)
                                                                              tGGACATTTTTACTCCATGTCGTCGGTGCTGCGAGTGTCGCAGATAAACATACCCAAAAGAAAACCCTCACCGTCAGGCGGCGAGGGTTTAACTCACATGATGATACTGACTGTTGCTCACTCTTTGAAGTGATTTGCGTCACATTCAg                                                                       <  4:788870/149‑1 (MQ=255)
                                                                              tGGACATTTTTACTCCATGTCGTCGGTGCTGCGAGTGTCGCAGATAAACATACCCAAAAGAAAACCCTCACCGTCAGGCGGCGAGGGTTTAACTCACATGATGATACTGACTGTTGCTCACTCTTTGAAGTGATTTGCGTCACATTCAg                                                                       >  1:1021145/1‑149 (MQ=255)
                                                                                  cATTTTTACTCCATGTCGTCGGTGCTGCGAGTGTCGCAGATAAACATACCCAAAAGAAAACCCTCACCGTCAGGCGGCGAGGGTTTAACTCACATGATGATACTGACTGTTGCTCACTCTTTGAAGTGATTTGCGTCACATTCAGGGaa                                                                   <  2:475494/149‑1 (MQ=255)
                                                                                         aCTCCATGTCGTCGGTGCTGCGAGTGTCGCAGATAAACATACCCAAAAGAAAACCCTCACCGTCAGGCGGCGAGGGTTTAACTCACATGATGATACTGACTGTTGCTCACTCTTTGAAGTGATTTGCGTCACATTCAGGGAATTCCTCa                                                            >  2:730317/1‑149 (MQ=255)
                                                                                           tCCATGTCGTCGGTGCTGCGAGTGTCGCAGATAAACATACCCAAAAGAAAACCCTCACCGTCAGGCGGCGAGGGTTTAACTCACATGATGATac                                                                                                                 <  3:302850/94‑1 (MQ=255)
                                                                                           tCCATGTCGTCGGTGCTGCGAGTGTCGCAGATAAACATACCCAAAAGAAAACCCTCACCGTCAGGCGGCGAGGGTTTAACTCACATGATGATac                                                                                                                 >  4:302850/1‑94 (MQ=255)
                                                                                           tCCATGTCGTCGGTGCTGCGAGTGTCGCAGATAAACATACCCAAAAGAAAACCCTCACCGTCAGGCGGCGAGGGTTTAACTCACATGATGATACTGACTGTTGCTCACTCTTTGAAGTGATTTGCGTCACATTCAGGGAATTCCTCAAt                                                          >  3:222486/1‑149 (MQ=255)
                                                                                           tCCATGTCGTCGGTGCTGCGAGTGTCGCAGATAAACATACCCAAAAGAAAACCCTCACCGTCAGGCGGCGAGGGTTTAACTCACATGATGATACTGACTGTTGCTCACTCTTTGAAGTGATTTGCGTCACATTCAGGGAATTCCTCAAt                                                          >  2:454111/1‑149 (MQ=255)
                                                                                              aTGTCGTCGGTGCTGCGAGTGTCGCAGATAAACATACCCAAAAGAAAACCCTCACCGTCAGGCGGCGAGGGTTTAACTCACATGATGATACTGACTGTTGCTCACTCTTTGAAGTGATTTGCGTCACATTCAGGGAATTCCTCAATGCa                                                       >  3:1056715/1‑149 (MQ=255)
                                                                                              aTGTCGTCGGTGCTGCGAGTGTCGCAGATAAACATACCCAAAAGAAAACCCTCACCGTCAGGCGGCGAGGGTTTAACTCACATGATGATACTGACTGTTGCTCACTCTTTGAAGTGATTTGCGTCACATTCAGGGAATTCCTCAATGCa                                                       <  4:603584/149‑1 (MQ=255)
                                                                                               tGTCGTCGGTGCTGCGAGTGTCGCAGATAAACATACCCAAAAGAAAACCCTCACCGTCAGGCGGCGAGGGTTTAACTCACATGATGATACTGACTGTTGCTCACTCTTTGAAGTGATTTGCGTCACATTCAGGGAATTCCTCAATGCAg                                                      >  1:966374/1‑149 (MQ=255)
                                                                                                   gtcgGTGCTGCGAGTGTCGCAGATAAACATACCCAAAAGAAAACCCTCACCGTCAGGCGGCGAGGGTTTAACTCACATGATGATACTGACTGTTGCTCACTCTTTGAAGTGAt                                                                                      >  4:858638/1‑113 (MQ=255)
                                                                                                   gtcgGTGCTGCGAGTGTCGCAGATAAACATACCCAAAAGAAAACCCTCACCGTCAGGCGGCGAGGGTTTAACTCACATGATGATACTGACTGTTGCTCACTCTTTGAAGTGAt                                                                                      <  3:858638/113‑1 (MQ=255)
                                                                                                   gtcgGTGCTGCGAGTGTCGCAGATAAACATACCCAAAAGAAAACCCTCACCGTCAGGCGGCGAGGGTTTAACTCACATGATGATACTGACTGTTGCTCACTCTTTGAAGTGATTTGCGTCACATTCAGGGAATTCCTCAATGCAGCGCa                                                  <  4:294975/149‑1 (MQ=255)
                                                                                                    tcgGTGCTGCGAGTGTCGCAGATAAACATACCCAAAAGAAAACCCTCACCGTCAGGCGGCGAGGGTTTAACTCACATGATGATACTGACTGTTGCTCACTCTTTGAAGTGATTTGCGTCACATTCAGGGAATTCCTCAATGCAGCGCAt                                                 >  2:957830/1‑149 (MQ=255)
                                                                                                      ggTGCTGCGAGTGTCGCAGATAAACATACCCAAAAGAAAACCCTCACCGTCAGGCGGCGAGGGTTTAACTCACATGATGATACTGACTGTTGCTCACTCTTTGAAGTGATTTGCGTCACATTCAGGGAATTCCTCAATGCAGCGCATTa                                               <  2:341550/149‑1 (MQ=255)
                                                                                                         gctgcGAGTGTCGCAGATAAACATACCCAAAAGAAAACCCTCACCGTCAGGCGGCGAGGGTTTAACTCACATGATGATACTGACTGTTGCTCACTCTTTGAAGTGATTTGCGTCACATTCAGGGAATTCCTCAATGCAGCGCATTATGt                                            >  3:506543/1‑149 (MQ=255)
                                                                                                          ctgcGAGTGTCGCAGATAAACATACCCAAAAGAAAACCCTCACCGTCAGGCGGCGAGGGTTTAACTCACATGATGATACTGACTGTTGCTCACTCTTTGAAGTGATTTGCGTCACATTCAGGGAATTCCTCAATGCAGCGCATTATGta                                           <  3:420616/149‑1 (MQ=255)
                                                                                                                gtgtCGCAGATAAACATACCCAAAAGAAAACCCTCACCGTCAGGCGGCGAGGGTTTAACTCACATGATGATACTGACTGTTGCTCACTCTTTGAAGTGATTTGCGTCACATTCAGGGAATTCCTCAATGCAGCGCATTATGTATAAATc                                     <  3:491208/149‑1 (MQ=255)
                                                                                                                   tCGCAGATAAACATACCCAAAAGAAAACCCTCACCGTCAGGCGGCGAGGGTTTAACTCACATGATGATACTGACTGTTGCTCACTCTTTGAAGTGATTTGCGTCACATTCAGGGAATTCCTCAATGCAGCGCATTATGTATAAATCTTa                                  >  2:759795/1‑149 (MQ=255)
                                                                                                                    cGCAGATAAACATACCCAAAAGAAAACCCTCACCGTCAGGCGGCGAGGGTTTAACTCACATGATGATACTGACTGTTGCTCACTCTTTGAAGTGATTTGCGTCACATTCAGGGAATTCCTCAATGCAGCGCATTATGTATAAATCTTaa                                 <  1:1053957/149‑1 (MQ=255)
                                                                                                                     gCAGATAAACATACCCAAAAGAAAACCCTCACCGTCAGGCGGCGAGGGTTTAACTCACATGATGATACTGACTGTTGCTCACTCTTTGAAGTGATTTGCGTCACATTCAGGGAATTCCTCAATGCAGCGCATTATGTATAAATCTTAAt                                <  3:444274/149‑1 (MQ=255)
                                                                                                                             aCATACCCAAAAGAAAACCCTCACCGTCAGGCGGCGAGGGTTTAACTCACATGATGATACTGACTGTTGCTCACTCTTTGAAGTGATTTGCGTCACATTCAGGGAATTCCTCAATGCAGCGCATTATGTATAAATCTTAATCGCCTTgg                        <  4:506543/149‑1 (MQ=255)
                                                                                                                             aCATACCCAAAAGAAAACCCTCACCGTCAGGCGGCGAGGGTTTAACTCACATGATGATACTGACTGTTGCTCACTCTTTGAAGTGATTTGCGTCACATTCAGGGAATTCCTCAATGCAGCGCATTATGTATAAATCTTAATCGCCTTgg                        >  4:79271/1‑149 (MQ=255)
                                                                                                                               aTACCCAAAAGAAAACCCTCACCGTCAGGCGGCGAGGGTTTAACTCACATGATGATACTGACTGTTGCTCACTCttt                                                                                              >  1:240204/1‑77 (MQ=255)
                                                                                                                               aTACCCAAAAGAAAACCCTCACCGTCAGGCGGCGAGGGTTTAACTCACATGATGATACTGACTGTTGCTCACTCttt                                                                                              <  2:240204/77‑1 (MQ=255)
                                                                                                                                  cccAAAAGAAAACCCTCACCGTCAGGCGGCGAGGGTTTAACTCACATGATGATACTGACTGTTGCTCACTCTTTGAAGTGATTTGCGTCACATTCAGGGAATTCCTCAATGCAGCGCATTATGTATAAATCTTAATCGCCTTGGTTTAt                   <  2:63798/149‑1 (MQ=255)
                                                                                                                                   ccAAAAGAAAACCCTCACCGTCAGGCGGCGAGGGTTTAACTCACATGATGATACTGACTGTTGCTCACTCTTTGAAGTGATTTGCGTCACATTCAGGGAATTCCTCAATGCAGCGCATTATGTATAAATCTTAATCGCCTTGGTTTATg                  >  4:946484/1‑149 (MQ=255)
                                                                                                                                   ccAAAAGAAAACCCTCACCGTCAGGCGGCGAGGGTTTAACTCACATGATGATACTGACTGTTGCTCACTCTTTGAAGTGATTTGCGTCACATTCAGGGAATTCCTCAATGCAGCGCATTATGTATAAATCTTAATCGCCTTGGTTTATg                  <  3:757615/149‑1 (MQ=255)
                                                                                                                                    cAAAAGAAAACCCTCACCGTCAGGCGGCGAGGGTTTAACTCACATGATGATACTGACTGTTGCTCACTCTTTGAAGTGATTTGCGTCACATTCAGGGAATTCCt                                                              <  2:203246/104‑1 (MQ=255)
                                                                                                                                    cAAAAGAAAACCCTCACCGTCAGGCGGCGAGGGTTTAACTCACATGATGATACTGACTGTTGCTCACTCTTTGAAGTGATTTGCGTCACATTCAGGGAATTCCTCAATGCAGCGCATTATGTATAAATCTTAATCGCCTTGGt                       <  1:419146/143‑1 (MQ=255)
                                                                                                                                    cAAAAGAAAACCCTCACCGTCAGGCGGCGAGGGTTTAACTCACATGATGATACTGACTGTTGCTCACTCTTTGAAGTGATTTGCGTCACATTCAGGGAATTCCTCAATGCAGCGCATTATGTATAAATCTTAATCGCCTTGGt                       >  2:419146/1‑143 (MQ=255)
                                                                                                                                    cAAAAGAAAACCCTCACCGTCAGGCGGCGAGGGTTTAACTCACATGATGATACTGACTGTTGCTCACTCTTTGAAGTGATTTGCGTCACATTCAGGGAATTCCTCAATGCAGCGCATTATGTATAAATCTTAATCGCCTTGGTTTATgg                 >  4:895366/1‑149 (MQ=255)
                                                                                                                                    cAAAAGAAAACCCTCACCGTCAGGCGGCGAGGGTTTAACTCACATGATGATACTGACTGTTGCTCACTCTTTGAAGTGATTTGCGTCACATTCAGGGAATTCCTCAATGCAGCGCATTATGTATAAATCTTAATCGCCTTGGTTTATgg                 >  3:296169/1‑149 (MQ=255)
                                                                                                                                    cAAAAGAAAACCCTCACCGTCAGGCGGCGAGGGTTTAACTCACATGATGATACTGACTGTTGCTCACTCTTTGAAGTGATTTGCGTCACATTCAGGGAATTCCTCAATGCAGCGCATTATGTATAAATCTTAATCGCCTTGGTTTATgg                 >  1:87384/1‑149 (MQ=255)
                                                                                                                                     aaaaGAAAACCCTCACCGTCAGGCGGCGAGGGTTTAACTCACATGATGATACTGACTGTTGCTCACTCTTTGAAGTGATTTGCGTCACATTCAGGGAATTCCTCAATGCAGCGCATTATGTATAAATCTTAATCGCCTTGGTTTATGGa                >  1:731893/1‑149 (MQ=255)
                                                                                                                                         gAAAACCCTCACCGTCAGGCGGCGAGGGTTTAACTCACATGATGATACTGACTGTTGCTCACTCTTTGAAGTGATTTGCGTCACATTCAGGGAATTCCTCAATGCAGCGCATTATGTATAAATCTTAATCGCCTTGGTTTATGGAAGAc            >  3:726758/1‑149 (MQ=255)
                                                                                                                                         gAAAACCCTCACCGTCAGGCGGCGAGGGTTTAACTCACATGATGATACTGACTGTTGCTCACTCTTTGAAGTGATTTGCGTCACATTCAGGGAATTCCTCAATGCAGCGCATTATGTATAAATCTTAATCGCCTTGGCTTATGGCAGAc            >  3:299772/1‑149 (MQ=255)
                                                                                                                                          aaaaCCCTCACCGTCAGGCGGCGAGGGTTTAACTCACATGATGATACTGACTGTTGCTCACTCTTTGAAGTGATTTGCGTCACATTCAg                                                                       >  4:480276/1‑89 (MQ=255)
                                                                                                                                            aaCCCTCACCGTCAGGCGGCGAGGGTTTAACTCACATGATGATACTGACTGTTGCTCACTCTTTGAAGTGATTTGCGTCACATTCAGGGAATTCCTCAATGCAGCGCATTATGTATAAATCTTAATCGCCTTGGTTTATGGAAGACGaa         >  2:360728/1‑149 (MQ=255)
                                                                                                                                            aaCCCTCACCGTCAGGCGGCGAGGGTTTAACTCACATGATGATACTGACTGTTGCTCACTCTTTGAAGTGATTTGCGTCACATTCAGGGAATTCCTCAATGCAGCGCATTATGTATAAATCTTAATCGCCTTGGTTTATGGAAGACGaa         >  2:450991/1‑149 (MQ=255)
                                                                                                                                             aCCCTCACCGTCAGGCGGCGAGGGTTTAACTCACATGATGATACTGACTGTTGCTCACTCTTTGAAGTGATTTGCGTCACATTCAGGGAATTCCTCAATGCAGCGCATTATGTATAAATCTTAATCGCCTTGGTTTATGGAAGACGAAt        <  1:759795/149‑1 (MQ=255)
                                                                                                                                             aCCCTCACCGTCAGGCGGCGAGGGTTTAACTCACATGATGATACTGACTGTTGCTCACTCTTTGAAGTGATTTGCGTCACATTCAGGGAATTCCTCAATGCAGCACATTATGTATAAATCTTAATCGCCTTGGTTTATGGAAGACGa          <  1:175310/147‑1 (MQ=255)
                                                                                                                                              cccTCGCCGCCTGACGGTGAGGGTTTAACTCACATGATGATACTGACTGTTGCTCACTCTTTGAAGTGATTTGCGTCACATTCAGGGAATTCCTCAATGCAGCGCATTATGTATaaa                                       <  3:667027/117‑1 (MQ=255)
                                                                                                                                              cccTCGCCGCCTGACGGTGAGGGTTTAACTCACATGATGATACTGACTGTTGCTCACTCTTTGAAGTGATTTGCGTCACATTCAGGGAATTCCTCAATGCAGCGCATTATGTATAAATCTTaa                                 <  3:137464/123‑1 (MQ=255)
                                                                                                                                              cccTCACCGTCAGGCGGCGAGGGTTTAACTCACATGATGATACTGACTGTTGCTCACTCTTTGAAGTGATTTGCGTCACATTCAGGGAATTCCTCAATGCAGCGCATTATGTATAAATCTTAATCGCCTTGGTTTATGGAAGACGAAtt       >  1:250464/1‑149 (MQ=255)
                                                                                                                                               ccTCACCGTCAGGCGGCGAGGGTTTAACTCACATGATGATACTGACTGTTGCTCACTCTTTGAAGTGATTTGCGTCACATTCAGGGAATTCCTCAATGCAGCGCATTATGTATAAATCTTAATCGCCTTGGTTTATGGAAGACGAATTa      >  3:882575/1‑149 (MQ=255)
                                                                                                                                               ccTCACCGTCAGGCGGCGAGGGTTTAACTCACATGATGATACTGACTGTTGCTCACTCTTTGAAGTGATTTGCGTCACATTCAGGGAATTCCTCAATGCAGCGCATTATGTATAAATCTTAATCGCCTTGGTTTATGGAAGACGAATTa      <  1:318888/149‑1 (MQ=255)
                                                                                                                                               ccTCACCGTCAGGCGGCGAGGGTTTAACTCACATGATGATACTGACTGTTGCTCACTCTTTGAAGTGATTTGCGTCACATTCAGGGAATTCCTCAATGCAGCGCATTATGTATAAATCTTAATCGCCTTGGTTTATGGAAGACGAATTa      >  1:552095/1‑149 (MQ=255)
                                                                                                                                                cTCACCGTCAGGCGGCGAGGGTTTAACTCACATGATGATACTGACTGTTGCTCACTCTTTGAAGTGATTTGCGTCACATTCAGGGAATTCCTCAATGCAGCGCATTATGTATAAATCTTAATCGCCTTGGTTTATGGAAGACGAATTAg     >  1:654054/1‑149 (MQ=255)
                                                                                                                                                cTCACCGTCAGGCGGCGAGGGTTTAACTCACATGATGATACTGACTGTTGCTCACTCTTTGAAGTGATTTGCGTCACATTCAGGGAATTCCTCAATGCAGCGCATTATGTATAAATCTTAATCGCCTTGGTTTATGGAAGACGAATTAg     >  1:662628/1‑149 (MQ=255)
                                                                                                                                                 tCACCGTCAGGCGGCGAGGGTTTAACTCACATGATGATACTGACTGTTGCTCACTCTTTGAAGTGATTTGCGTCACATTCAGGGAATTCCTCAATGCAGCGCATTATGTATAAATCTTAATCGCCTTGGTTTATGGAAGACGAATTAGc    <  3:111463/149‑1 (MQ=255)
                                                                                                                                                   aCCGTCAGGCGGCGAGGGTTTAACTCACATGATGATACTGACTGTTGCTCACTCTTTGAAGTGATTTGCGTCACATTCAGGGAATTCCTCAATGCAGCGCATTATGTATAAATCTTAATCGCCTTGGTTTATGGAAGACGAATTAGCgt  <  1:684297/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
GCAAACCAGCATAGGCGGGTTGCTCCAGCAGAATATCGATGCGTTCGGCATCCAGGTCGTTAATGATGATAGTTGGTCTGGACATTTTTACTCCATGTCGTCGGTGCTGCGAGTGTCGCAGATAAACATACCCAAAAGAAAACCCTCACCGTCAGGCGGCGAGGGTTTAACTCACATGATGATACTGACTGTTGCTCACTCTTTGAAGTGATTTGCGTCACATTCAGGGAATTCCTCAATGCAGCGCATTATGTATAAATCTTAATCGCCTTGGTTTATGGAAGACGAATTAGCGT  >  NC_000913/643883‑644178

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