Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,368,982 G→C 6.0% intergenic (+28/+32) yhcG → / ← yhcH DUF1016 family protein in the PD‑(D/E)XK nuclease superfamily/DUF386 family protein, cupin superfamily

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,368,9820GC 6.0% 646.7 / 23.8 183intergenic (+28/+32)yhcG/yhcHDUF1016 family protein in the PD‑(D/E)XK nuclease superfamily/DUF386 family protein, cupin superfamily
Reads supporting (aligned to +/- strand):  ref base G (86/86);  new base C (5/6);  total (91/92)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.32e-01

Only 100 of 188 total aligned reads displayed.

GTTTGCCTAACACCGTTCTGGCAAGCGAATATAAGATGCAACTACCTGATGAGAAACGACTCGCAGATGAACTCGTTCGAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTAACCTTCACCACAACCTTTTTAATCTCTCCAGGCTCGCCGACAACGCATCCTGGTTTATGCGGTTCACCTGGCATAAACACGGCGAACATTCCCGG  >  NC_000913/3368836‑3369130
                                                                                                                                                  |                                                                                                                                                    
gTTTGCCTAACACCGTTCTGGCAAGCGAATATAAGATGCAACTACCTGATGAGAAACGACTCGCAGATGAACTCGTTCGAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAgtg                                                                                                                                                    >  1:1201424/1‑149 (MQ=255)
 tttGCCTAACACCGTTCTGGCAAGCGAATATAAGATGCAACTACCTGATGAGAAACGACTCGCAGATGAACTCGTTCGAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAgtgt                                                                                                                                                   <  1:1348127/149‑1 (MQ=255)
 tttGCCTAACACCGTTCTGGCAAGCGAATATAAGATGCAACTACCTGATGAGAAACGACTCGCAGATGAACTCGTTCGAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACActgt                                                                                                                                                   <  1:1180455/149‑1 (MQ=255)
  ttGCCTAACACCGTTCTGGCAAGCGAATATAAGATGCAACTACCTGATGAGAAACGACTCGCAGATGAACTCGTTCGAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTc                                                                                                                                                  <  3:1424955/149‑1 (MQ=255)
   tGCCTAACACCGTTCTGGCAAGCGAATATAAGATGCAACTACCTGATGAGAAACGACTCGCAGATGAACTCGTTCGAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCt                                                                                                                                                 <  3:5051/149‑1 (MQ=255)
   tGCCTAACACCGTTCTGGCAAGCGAATATAAGATGCAACTACCTGATGAGAAACGACTCGCAGATGAACTCGTTCGAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCt                                                                                                                                                 <  2:858749/149‑1 (MQ=255)
    gCCTAACACCGTTCTGGCAAGCGAATATAAGATGCAACTACCTGATGAGAAACGACTCGCAGATGAACTCGTTCGAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTg                                                                                                                                                >  1:535166/1‑149 (MQ=255)
      cTAACACCGTTCTGGCAAGCGAATATACGATGCAACTACCTGGTGAGAAACGACTCGCAGTTGAACTCGTTCGAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCa                                                                                                                                              <  1:1561659/149‑1 (MQ=255)
      cTAACACCGTTCTGGCAAGCGAATATAAGATGCAACTACCTGATGAGAAACGACTCGCAGATGAACTCGTTCGAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACACTGTCTGCa                                                                                                                                              <  4:717890/149‑1 (MQ=255)
         acacCGTTCTGGCAAGCGAATATAAGATGCAACTACCTGATGAGAAACGACTCGCAGATGAACTCGTTCGAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACACTGTCTGCAGaa                                                                                                                                           <  3:1528338/149‑1 (MQ=255)
          cacCGTTCTGGCAAGCGAATATAAGATGCAACTACCTGATGAGAAACGACTCGCAGATGAACTCGTTCGAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTc                                                                                                                                                  <  2:1049592/141‑1 (MQ=255)
          cacCGTTCTGGCAAGCGAATATAAGATGCAACTACCTGATGAGAAACGACTCGCAGATGAACTCGTTCGAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACACTGTCTGCAGAAt                                                                                                                                          <  3:135361/149‑1 (MQ=255)
           acCGTTCTGGCAAGCGAATATAAGATGCAACTACCTGATGAGAAACGACTCGCAGATGAACTCGTTCGAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAATc                                                                                                                                         >  3:105155/1‑149 (MQ=255)
           acCGTTCTGGCAAGCGAATATAAGATGCAACTACCTGATGAGAAACGACTCGCAGATGAACTCGTTCGAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAATc                                                                                                                                         >  2:1514620/1‑149 (MQ=255)
               ttCTGGCAAGCGAATATAAGATGCAACTACCTGATGAGAAACGACTCGCAGATGAACTCGTTCGAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAATCTTTa                                                                                                                                     <  2:621909/149‑1 (MQ=255)
                   ggCAAGCGAATATAAGATGCAATTACCTGATGAGAAACGACTCGCAGATGAACTCGTTCGAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAATCTTTAAAAc                                                                                                                                 >  1:1042967/1‑149 (MQ=255)
                     cAAGCGAATATAAGATGCAACTACCTGATGAGAAACGACTCGCAGATGAACTCGTTCGAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGa                                                                                                                               <  3:625302/149‑1 (MQ=255)
                      aaGCGAATATAATATGCAACTACCTGATGAGAAACGACTCGCAGATGAACGCGTTCGAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAGCGGTTACGAATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGaa                                                                                                                              <  2:969169/149‑1 (MQ=255)
                      aaGCGAATATAAGATGCAACTACCTGATGAGAAACGACTCGCAGATGAACTCGTTCGAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGaa                                                                                                                              <  2:832750/149‑1 (MQ=255)
                      aaGCGAATATAAGATGCAACTACCTGATGAGAAACGACTCGCAGATGAACTCGTTCGAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGaa                                                                                                                              <  3:1371397/149‑1 (MQ=255)
                                aaGATGCAACTACCTGATGAGAAACGACTCGCAGATGAACTCGTTCGAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTc                                                                                                                                                  >  2:490020/1‑119 (MQ=255)
                                aaGATGCAACTACCTGATGAGAAACGACTCGCAGATGAACTCGTTCGAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTc                                                                                                                                                  <  1:490020/119‑1 (MQ=255)
                                 aGATGCAACTACCTGATGAGAAACGACTCGCAGATGAACTCGTTCGAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTaa                                                                                                                   <  3:718539/149‑1 (MQ=255)
                                 aGATGCAACTACCTGATGAGAAACGACTCGCAGATGAACTCGTTCGAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTaa                                                                                                                   <  2:582389/149‑1 (MQ=255)
                                 aGATGCAACTACCTGATGAGAAACGACTCGCAGATGAACTCGTTCGAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTaa                                                                                                                   <  1:1514620/149‑1 (MQ=255)
                                 aGATGCAACTACCTGATGAGAAACGACTCGCAGATGAACTCGTTCGAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTaa                                                                                                                   <  1:793341/149‑1 (MQ=255)
                                 aGATGCAACTACCTGATGAGAAACGACTCGCAGATGAACTCGTTCGAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTaa                                                                                                                   <  3:1496736/149‑1 (MQ=255)
                                  gATGCAACTACCTGATGAGAAACGACTCGCAGATGAACTCGTTCGAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAg                                                                                                                  >  3:646132/1‑149 (MQ=255)
                                       aaCTACCTGATGAGAAACGACTCGCAGATGAACTCGTTCGAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCAtt                                                                                                             >  1:776512/1‑149 (MQ=255)
                                         cTACCTGATGAGAAACGACTCGCAGATGAACTCGTTCGAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTaa                                                                                                           <  4:1334310/149‑1 (MQ=255)
                                          tACCTGATGAGAAACGACTCGCAGATGAACTCGTTCGAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTaaa                                                                                                          >  2:21458/1‑149 (MQ=255)
                                           aCCTGATGAGAAACGACTCGCAGATGAACTCGTTCGAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAAt                                                                                                         >  3:292418/1‑149 (MQ=255)
                                               gatgaGAAACGACTCGCAGATGAACTCGTTCGAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGc                                                                                                     >  2:724655/1‑149 (MQ=255)
                                               gatgaGAAACGACTCGCAGATGAACTCGTTCGAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGc                                                                                                     >  1:1609246/1‑149 (MQ=255)
                                                   agaAACGACTCGCAGATGAACTCGTTCGAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTa                                                                                                 <  4:541771/149‑1 (MQ=255)
                                                   agaAACGACTCGCAGATGAACTCGTTCGAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTa                                                                                                 >  4:772818/1‑149 (MQ=255)
                                                    gaAACGACTCGCAGATGAACTCGTTCGAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTaa                                                                                                <  2:898748/149‑1 (MQ=255)
                                                     aaaCGACTCGCAGATGAACTCGTTCGAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTAAc                                                                                               <  2:710369/149‑1 (MQ=255)
                                                     aaaCGACTCGCAGATGAACTCGTTCGAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTAAc                                                                                               <  3:979815/149‑1 (MQ=255)
                                                     aaaCGACTCGCAGATGAACTCGTTCGAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTAAc                                                                                               >  1:607946/1‑149 (MQ=255)
                                                      aaCGACTCGCAGATGAACTCGTTCGAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATc                                                                                                        >  4:825738/1‑139 (MQ=255)
                                                           cTCGCAGATGAACTCGTTCGAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTAACCTTcac                                                                                         <  1:1066213/149‑1 (MQ=255)
                                                            tCGCAGATGAACTCGTTCGAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTAACCTTcacc                                                                                        >  4:21550/1‑149 (MQ=255)
                                                            tCGCAGATGAACTCGTTCGAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTAACCTTcacc                                                                                        >  1:1154261/1‑149 (MQ=255)
                                                            tCGCAGATGAACTCGTTCGAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTAACCTTcacc                                                                                        >  3:1256408/1‑149 (MQ=255)
                                                            tCGCAGATGAACTCGTTCGAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTAACCTTcacc                                                                                        >  3:1190419/1‑149 (MQ=255)
                                                                aGATGAACTCGTTCGAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACACTGTCTGCAGaa                                                                                                                                           <  3:312712/94‑1 (MQ=255)
                                                                aGATGAACTCGTTCGAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACACTGTCTGCAGaa                                                                                                                                           >  4:312712/1‑94 (MQ=255)
                                                                  aTGAACTCGTTCGAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTaa                                                                                                                   <  1:186357/116‑1 (MQ=255)
                                                                   tGAACTCGTTCGAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTAACCTTCACCACAACCt                                                                                 <  2:535166/149‑1 (MQ=255)
                                                                    gAACTCGTTCGAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAACCATTAAATCAGCCTTAACCTTCACCACAACCtt                                                                                >  1:894421/1‑149 (MQ=255)
                                                                     aaCTCGTTCGAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTAACCTTCACCACAACCttt                                                                               <  1:1034349/149‑1 (MQ=255)
                                                                      aCTCGTTCGAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTAACCTTCACCACAACCtttt                                                                              <  4:292418/149‑1 (MQ=255)
                                                                      aCTCGTTCGAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTAACCTTCACCACAACCtttt                                                                              <  1:1457477/149‑1 (MQ=255)
                                                                        tCGTTCGAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAATCtt                                                                                                                                       >  2:394336/1‑90 (MQ=255)
                                                                        tCGTTCGAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTAACCTTCACCACAACCTTTTTa                                                                            >  2:1139201/1‑149 (MQ=255)
                                                                        tCGTTCGAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTAACCTTCACCACAACCTTTTTa                                                                            >  3:243105/1‑149 (MQ=255)
                                                                          gTTCGAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTAACCTTCACCACAACCTTTTTAAt                                                                          <  1:1139201/149‑1 (MQ=255)
                                                                               aaCACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTAACCTTCACCACAACCTTTTTAAtctctc                                                                     >  3:1025721/1‑149 (MQ=255)
                                                                                  acaGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTAACCTTCACCACAACCTTTTTAATCTCTCCAg                                                                  >  4:1230245/1‑149 (MQ=255)
                                                                                      gCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAAt                                                                                                                                          >  1:909212/1‑73 (MQ=255)
                                                                                        ggTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTAACCTTCACCACAACCTTTTTAATCTCTCCAGGCTCGc                                                            <  4:256484/149‑1 (MQ=255)
                                                                                            cTAGAGGAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTAACCTTCACCACAACCTTTTTAATCTCTCCAGGCTCGCCGAc                                                        >  3:753027/1‑149 (MQ=255)
                                                                                              agagGAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTAACCTTCACCACAACCTTTTTAATCTCTCCAGGCTCGCCGACaa                                                      <  4:631563/149‑1 (MQ=255)
                                                                                               gagGAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAg                                                                                                                  >  1:270092/1‑88 (MQ=255)
                                                                                                  gAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTAACCTTCACCACAACCTTTTTAAt                                                                          >  4:449130/1‑125 (MQ=255)
                                                                                                     ggCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTAACCTTCACCACAACCTTTTTAATCTCTCCAGGCTCGCCGACAACGCATcc                                               <  3:1304220/149‑1 (MQ=255)
                                                                                                     ggCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTAACCTTCACCACAACCTTTTTAATCTCTCCAGGCTCGCCGACAACGCATcc                                               <  4:1524144/149‑1 (MQ=255)
                                                                                                      gCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTAACCTTCACCACAACCTTTTTAATCTCTCCAGGCTCGCCGACAACGCATCCt                                              <  4:1025721/149‑1 (MQ=255)
                                                                                                        tataGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTAACCTTcacc                                                                                        <  1:1501404/105‑1 (MQ=255)
                                                                                                        tataGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTAACCTTCACCACAACCTTTTTAATCTCTCCAGGCTCGCCGACAACGCATCCTgg                                            >  2:1360215/1‑149 (MQ=255)
                                                                                                         ataGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTAACCTTCACCACAACCTTTTTAATCTCTCCAGGCTCGCCGACAACGCATCCTGGt                                           >  2:1705212/1‑149 (MQ=255)
                                                                                                         ataGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTAACCTTCACCACAACCTTTTTAATCTCTCCAGGCTCGCCGACAACGCATCCTGGt                                           <  2:1575988/149‑1 (MQ=255)
                                                                                                           aGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAATCTGTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTAACCTTCACCACAACCTTTTTAATCTCTCCAGGCTCGCCGACAACGCATCCTGGttt                                         >  2:481559/1‑149 (MQ=255)
                                                                                                               gCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTATGCCATTAAATCAGCCTTAACCTTCACCACAACCTTTTTAATCTCTCCAGGCTCGCCGACAACGCATCCTGGTTTATGc                                     >  2:387323/1‑149 (MQ=255)
                                                                                                                 cGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTAACCTTCACCACAACCTTTTTAATCTCTCCAGGCTCGCCGACAACGCATCCTGGTTTATGCgg                                   <  3:21550/149‑1 (MQ=255)
                                                                                                                  gTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTAACCTTCACCACAACCTTTTTAATCTCTCCAGGCTCGCCGACAAc                                                     <  1:389041/130‑1 (MQ=255)
                                                                                                                    tAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTAACCTTCACCACAACCTTTTTAATCTCTCCAGGCTCGCCGACAACGCATCCTGGTTTATGCGGTTc                                <  1:189114/149‑1 (MQ=255)
                                                                                                                     aaTTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTAACCTTCACCACAACCTTTTTAATCTCTCCAGGCTCGCCGACAACGCATCCTGGTTTATGCGGTTCa                               >  3:1336059/1‑149 (MQ=255)
                                                                                                                         tCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTAACCTTCACCACAACCTTTTTAATCTCTCCAGGCTCGCCGACAACGCATCCTGGTTTATGCGGTTCACCTg                           <  4:753027/149‑1 (MQ=255)
                                                                                                                          cAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTAACCTTCACCACAACCTTTTTAATCTCTCCAGGCTCGCCGACAACGCATCCTGGTTTATGCGGTTCACCTgg                          >  1:1688597/1‑149 (MQ=255)
                                                                                                                                tttCGAATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTAACCTTCACCACAACCTTTTTAATCTCTCCAGGCTCGCCGACAACGCATCCTGGTTTATGCGGTTCACCTGGCATaaa                    >  1:292171/1‑149 (MQ=255)
                                                                                                                                tttCGAATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTAACCTTCACCACAACCTTTTTAATCTCTCCAGGCTCGCCGACAACGCATCCTGGTTTATGCGGTTCACCTGGCATaaa                    >  1:641245/1‑149 (MQ=255)
                                                                                                                                tttCGAATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTAACCTTCACCACAACCTTTTTAATCTCTCCAGGCTCGCCGACAACGCATCCTGGTTTATGCGGTTCACCTGGCATaaa                    >  4:280929/1‑149 (MQ=255)
                                                                                                                                  tCGAATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTAACCTTCACCACAACCTTTTTAATCTCTCCAGGCTCGCCGACAACGCATCCTGGTTTATGCGGTTCACCTGGCATAAaca                  >  3:865449/1‑149 (MQ=255)
                                                                                                                                   cGAATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTAACCTTCACCACAACCTTTTTAATCTCTCCAGGCTCGCCGACAACGCATCCTGGTTTATGCGGTTCACCTGGCATAAacac                 <  2:1042967/149‑1 (MQ=255)
                                                                                                                                    agaTTCTGCAGACACTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTAACCTTCACCACAACCTTTTTAATCTCTCCAGGCTCGCCGACAACGCATCCTGGTTTATGCGGTTCACCTGGCATAAACACg                >  4:1320758/3‑149 (MQ=255)
                                                                                                                                     aaTTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTAACCTTCACCACAACCTTTTTAATCTCTCCAGGCTCGCCGACAACGCATCCTGGTTTATGCGGTTCACCTGGCATAAACACgg               >  1:397089/1‑149 (MQ=255)
                                                                                                                                       ttCTGCAGACACTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTAACCTTCACCACAACCTTTTTAATCTCTCCAGGCTCGCCGACAACGCATCCTGGTTTATGCGGTTCACCTg                           >  3:467672/1‑135 (MQ=255)
                                                                                                                                       ttCTGCAGACACTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTAACCTTCACCACAACCTTTTTAATCTCTCCAGGCTCGCCGACAACGCATCCTGGTTTATGCGGTTCACCTg                           <  4:467672/135‑1 (MQ=255)
                                                                                                                                        tCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTAACCTTCACCACAACCTTTTTAATCTCTCCAGGCTCGCCGACAACGCATCCTGGTTTATGCGGTTCACCTGGCATAAACACGGCGa            <  2:1201424/149‑1 (MQ=255)
                                                                                                                                         cTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTAACCTTCACCACAACCTTTTTAATCTCTCCAGGCTCGCCGACAACGCATCCTGGTTTATGCGGTTCACCTGGCATAAACACGGCGaa           <  4:584707/149‑1 (MQ=255)
                                                                                                                                          tGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTAACCTTCACCACAACCTTTTTAATCTCTCCAGGCTCGCCGACAACGCATCCTGGTTTATGCGGtt                                 >  4:372432/1‑126 (MQ=255)
                                                                                                                                          tGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTAACCTTCACCACAACCTTTTTAATCTCTCCAGGCTCGCCGACAACGCATCCTGGTTTATGCGGtt                                 <  3:372432/126‑1 (MQ=255)
                                                                                                                                          tGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTAACCTTCACCACAACCTTTTTAATCTCTCCAGGCTCGCCGACAACGCATCCTGGTTTATGCGGTTCACCTGGCATAAACACGGCGAAc          >  3:296929/1‑149 (MQ=255)
                                                                                                                                          tGCAGACACTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTAACCTTCACCACAACCTTTTTAATCTCTCCAGGCTCGCCGTCAACGCATCCTGGTTTATGCGGTTCACCTGGCATAAACACGGCGAAc          >  1:340239/1‑149 (MQ=255)
                                                                                                                                            cAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTAACCTTCACCACAACCTTTTTAATCTCTCCAGGCTCGCCGACAACGCATCCTGGTTTATGCGGTTCACCTGGCATAAACACGGCGAACAt        >  4:800823/1‑149 (MQ=255)
                                                                                                                                               aCAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTAACCTTCACCACAACCTTTTTAATCTCTCCAGGCTCGCCGACAACGCATCCTGGGTTATGCGGTTCACCTGGCATAAACACGGCGAACATTcc     >  2:1020148/1‑149 (MQ=255)
                                                                                                                                                cAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTAACCTTCACCACAACCTTTTTAATCTCTCCAGGCTCGCCGACAACGCATCCTGGTTTATGCGGTTCACCTGGCATAAACACGGCGAACATTccc    <  2:728648/149‑1 (MQ=255)
                                                                                                                                                  gtgtCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTAACCTTCACCACAACCTTTTTAATCTCTCCAGGCTCGCCGACAACGCATCCTGGTTTATGCGGTTCACCTGGCATAAACACGGCGAACATTCCCgg  <  3:586464/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
GTTTGCCTAACACCGTTCTGGCAAGCGAATATAAGATGCAACTACCTGATGAGAAACGACTCGCAGATGAACTCGTTCGAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTAACCTTCACCACAACCTTTTTAATCTCTCCAGGCTCGCCGACAACGCATCCTGGTTTATGCGGTTCACCTGGCATAAACACGGCGAACATTCCCGG  >  NC_000913/3368836‑3369130

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: