Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,205,788 G→C 5.3% pseudogene (58/412 nt) beeE → pseudogene, portal protein family, e14 prophage;Phage or Prophage Related

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,205,7880GC 5.3% 334.4 / 7.9 95pseudogene (58/412 nt)beeEpseudogene, portal protein family, e14 prophage;Phage or Prophage Related
Reads supporting (aligned to +/- strand):  ref base G (43/47);  new base C (2/3);  total (45/50)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.21e-01

TGTTGCCGGGGCGCTGTGCCTGTTCTGGTCGTGGCTGGTGGCGCGATATCTCGACCGTACACAGTCGTCTGTCGGCGGAGGTAAATAGTGTTCTTTTCGGGATTATTTCAACGAAAAAGTGACGCACCGGTGACCACGCCAGCAGAGCTGGCGGATGCCATCGGGCTGTCGTATGACACCTATACCGGAAAGCAGATCAGCAGTCAGCGGGCTATGCGACTGACGGCGGTTTTTTCCTGCGTCAGAGTGCTGGCAGAGTCGGTCGGGATGTTGCCCTGCAATCTGT  >  NC_000913/1205644‑1205929
                                                                                                                                                |                                                                                                                                             
tgttgCCGGGGCGCTGTGCCTGTTCTGGTCGTGGCTGGTGGCGCGATATCTCGACCGTACACAGTCGTCTGTCGGCGGAGGTAAATAGTGTTCTTTTCGGGATTATTTCAACGAAAAAGTGACGCACCGGTGACCACGCCAGCAGAGCt                                                                                                                                           <  3:156225/149‑1 (MQ=255)
  ttgCCGGGGCGCTGTGCCTGTTCTGGTCGTGGCTGGTGGCGCGATATCTCGACCGTACACAGTCGTCTGTCGGCGGAGGTAAATAGTGTTCTTTTCGGGATTATTTCAACGAAAAAGTGATGCACCGGTGACCACGCCAGCAGAGCTgg                                                                                                                                         <  1:994103/149‑1 (MQ=255)
      cGGGGCGCTGTGCCTGTTCTGGTCGTGGCTGGTGGCGCGATATCTCGACCGTACACAGTCGTCTGTCGGCGGAGGTAAATAGTGTTCTTTTCGGGATTATTTCAACGAAAAAGTGACGCACCGGTGACCACGCCAGCTCTGCTGGCgtg                                                                                                                                     <  2:655633/149‑3 (MQ=255)
       ggggCGCTGTGCCTGTTCTGGTCGTGGCTGGTGGCGCGATATCTCGACCGTACACAGTCGTCTGTCGGCGGAGGTAAATAGTGTTCTTTTCGGGATTATTTCAACGAAAAAGTGACGCACCGGTGACCACGCCAGCAGAGCTGGCGGAt                                                                                                                                    <  1:812725/149‑1 (MQ=255)
       ggggCGCTGTGCCTGTTCTGGTCGTGGCTGGTGGCGCGATATCTCGACCGTACACAGTCGTCTGTCGGCGGAGGTAAATAGTGTTCTTTTCGGGATTATTTCAACGAAAAAGTGACGCACCGGTGACCACGCCAGCAGAGCTGGCGGAt                                                                                                                                    <  3:276541/149‑1 (MQ=255)
       ggggCGCTGTGCCTGTTCTGGTCGTGGCTGGTGGCGCGATATCTCGACCGTACACAGTCGTCTGTCGGCGGAGGTAAATAGTGTTCTTTTCGGGATTATTTCAACGAAAAAGTGACGCACCGGTGACCACGCCAGCAGAGCTGGCGGAt                                                                                                                                    <  3:511754/149‑1 (MQ=255)
            gcTGTGCCTGTTCTGGTCGTGGCTGGTGGCGCGATATCTCGACCGTACACAGTCGTCTGTCGGCGGAGGTAAATAGTGTTCTTTTCGGGATTATTTCAACGAAAAAGTGACGCACCGGTGACCACGCCAGCAGAGCTGGCGGATGCCAt                                                                                                                               >  1:482850/1‑149 (MQ=255)
             cTGTGCCTGTTCTGGTCGTGGCTGGTGGCGCGATATCTCGACCGTACACAGTCGTCTGTCGGCGGAGGTAAATAGTGTTCTTTTCGGGATTATTTCAACGAAAAAGTGACGCACCGGTGACCACGCCAGCAGAGCTGGCGGATGCCATc                                                                                                                              >  1:117359/1‑149 (MQ=255)
              tgtgCCTGTTCTGGTCGTGGCTGGTGGCGCGATATCTCGACCGTACACAGTCGTCTGTCGGCGGAGGTAAATAGTGTTCTTTTCGGGATTATTTCAACGAAAAAGTGACGCACCGGTGACCACGCCAGCAGAGCTGGCGGATGCCATCg                                                                                                                             <  1:153159/149‑1 (MQ=255)
               gtgCCTGTTCTGGTCGTGGCTGGTGGCGCGATATCTCGACCGTACACAGTCGTCTGTCGGCGGAGGTAAATAGTGTTCTTTTCGGGATTATTTCAACGAAAAAGTGACGCACCGGTGACCACGCCAGCAGAGCTGGCGGATGCCATCgg                                                                                                                            <  3:922972/149‑1 (MQ=255)
                 gCCTGTTCTGGTCGTGGCTGGTGGCGCGATATCTCGACCGTACACAGTCGTCTGTCGGCGGAGGTAAATAGTGTTCTTTTCGGGATTATTTCAACGAAAAAGTGACGCACCGGTGACCACGCCAGCAGAGCTGGCGGATGCCATCGGGc                                                                                                                          >  3:873315/1‑149 (MQ=255)
                 gCCTGTTCTGGTCGTGGCTGGTGGCGCGATATCTCGACCGTACACAGTCGTCTGTCGGCGGAGGTAAATAGTGTTCTTTTCGGGATTATTTCAACGAAAAAGTGACGCACCGGTGACCACGCCAGCAGAGCTGGCGGATGCCATCGGGc                                                                                                                          >  2:829230/1‑149 (MQ=255)
                  ccTGTTCTGGTCGTGGCTGGTGGCGCGATATCTCGACCGTACACAGTCGTCTGTCGGCGGAGGTAAATAGTGTTCTTTTCGGGATTATTTCAACGAAAAAGTGACGCACCGGTGACCACGCCAGCAGAGCTGGCGGATGCCATCGGGCt                                                                                                                         >  1:867519/1‑149 (MQ=255)
                   cTGTTCTGGTCGTGGCTGGTGGCGCGATATCTCGACCGTACACAGTCGTCTGTCGGCGGAGGTAAATAGTGTTCTTTTCGGGATTATTTCAACGAAAAAGTGACGCACCGGTGACCACGCCAGCAGAGCTGGCGGATGCCATCGGGCTg                                                                                                                        >  3:106098/1‑149 (MQ=255)
                    tGTTCTGGTCGTGGCTGGTGGCGCGATATCTCGACCGTACACAGTCGTCTGTCGGCGGAGGTAAATAGTGTTCTTTTCGGGATTATTTCAACGAAAAAGTGACGCACCGGTGACCACGCCAGCAGAGCTGGCGGATGCCATCGGGCTgc                                                                                                                       >  4:170953/1‑148 (MQ=255)
                      ttCTGGTCGTGGCTGGTGGCGCGATATCTCGACCGTACACAGTCGTCTGTCGGCGGAGGTAAATAGTGTTCTTTTCGGGATTATTTCAACGAAAAAGTGACGCACCGGTGACCACGCCAGCAGAGCTGGCGGATGCCATCGGGCTGCCg                                                                                                                     >  2:809713/1‑149 (MQ=255)
                        cTGGTCGTGGCTGGTGGCGCGATATCTCGACCGTACACAGTCGTCTGTCGGCGGAGGTAAATAGTGTTCTTTTCGGGATTATTTCAACGAAAAAGTGACGCACCGGTGACCACGCCAGCAGAGCTGGCGGATGCCATCGGGCTGTCGTa                                                                                                                   <  1:292564/149‑1 (MQ=255)
                        cTGGTCGTGGCTGGTGGCGCGATATCTCGACCGTACACAGTCGTCTGTCGGCGGAGGTAAATAGTGTTCTTTTCGGGATTATTTCAACGAAAAAGTGACGCACCGGTGACCACGCCAGCAGAGCTGGCGGATGCCATCGGGCTGTCGTa                                                                                                                   <  4:873315/149‑1 (MQ=255)
                         tGGTCGTGGCTGGTGGCGCGATATCTCGACCGTACACAGTCGTCTGTCGGCGGAGGTAAATAGTGTTCTTTTCGGGATTATTTCAACGAAAAAGTGACGCACCGGTGACCACGCCAGCAGAGCTGGCGGATGCCATCGGGCTGTCGTAt                                                                                                                  <  4:114205/149‑1 (MQ=255)
                         tGGTCGTGGCTGGTGGCGCGATATCTCGACCGTACACAGTCGTCTGTCGGCGGAGGTAAATAGTATTCTCTTCGGGATTATTTCAACGAAAAAGTGACGCACCGGTGACCACGCCAGCAGAGCTGGCGGATGCCATCGGGCTGTCGTAt                                                                                                                  <  3:761348/149‑1 (MQ=255)
                         tGGTCGTGGCTGGTGGCGCGATATCTCGACCGTACACAGTCGTCGGTCGGCGGAGGTAAATAGTGTTCTTTTCGGGATTATTTCAACGAAAAAGTGACGCACCGGTGACCACGCCAGCAGAGCTGGCGGATGCCATCGGGCTGTCGTAt                                                                                                                  >  4:503302/1‑149 (MQ=255)
                               tGGCTGGTGGCGCGATATCTCGACCGTACACAGTCGTCTGTCGGCGGAGGTAAATAGTGTTCTTTTCGGGATTATTTCAACGAAAAAGTGACGCACCGGTGACCACGCCAGCAGAGCTGGCGGATGCCATCGGGCTGTCGTATGACAcc                                                                                                            >  4:357741/1‑149 (MQ=255)
                               tGGCTGGTGGCGCGATATCTCGACCGTACACAGTCGTCTGTCGGCGGAGGTAAATAGTGTTCTTTTCGGGATTATTTCAACGAAAAAGTGACGCACCGGTGACCACGCCAGCAGAGCTGGCGGATGCCATCGGGCTGTCGTATGACAcc                                                                                                            >  4:700208/1‑149 (MQ=255)
                                          gcgATATCTCGACCGTACACAGTCGTCTGTCGGCGGAGGTAAATAGTGTTCTTTTCGGGATTATTTCAACGAAAAAGTGACGCACCGGTGACCACGCCAGCAGAGCTGGCGGATGCCATCGGGCTGTCGTATGACACCTATACCGGaaa                                                                                                 <  1:975064/149‑1 (MQ=255)
                                          gcgATATCTCGACCGTACACAGTCGTCTGTCGGCGGAGGTAAATAGTGTTCTTTTCGGGATTATTTCAACGAAAAAGTGACGCACCGGTGACCACGCCAGCAGAGCTGGCGGATGCCATCGGGCTGTCGTATGACACCTATACCGGaaa                                                                                                 <  1:727730/149‑1 (MQ=255)
                                            gATATCTCGACCGTACACAGTCGTCTGTCGGCGGAGGTAAATAGTGTTCTTTTCGGGATTATTTCAACGAAAAAGTGACGCACCGGTGACCACGCCAGCAGAGCTGGCGGATGCCATCGGGCTGTCGTATGACACCTATACCGGAAAGc                                                                                               <  4:299750/149‑1 (MQ=255)
                                             atatCTCGACCGTACACAGTCGTCTGTCGGCGGAGGTAAATAGTGTTCTTTTCGGGATTATTTCAACGAAAAAGTGACGCACCGGTGACCACGCCAGCAGAGCTGGCGGATGCCATCGGGCTGTCGTATGACACCTATACCGGAAAGCa                                                                                              <  4:980998/149‑1 (MQ=255)
                                                tctcGACCGTACACAGTCGTCTGTCGGCGGAGGTAAATAGTGTTCTTTTCGGGATTATTTCAACGAAAAAGTGACGCACCGGTGACCACGCCAGCAGAGCTGGCGGATGCCATCGGGCTGTCGTATGACACCTATACCGGAAAGCAGAt                                                                                           <  2:567808/149‑1 (MQ=255)
                                                tctcGACCGTACACAGTCGTCTGTCGGCGGAGGTAAATAGTGTTCTTTTCGGGATTATTTCAACGAAAAAGTGACGCACCGGTGACCACGCCAGCAGAGCTGCGGATGCCATCGGGCTGTCGTATGACACCTATACCGGAAAGCAGATc                                                                                          >  4:312946/1‑149 (MQ=255)
                                                  tcGACCGTACACAGTCGTCTGTCGGCGGAGGTATATAGTGTTCTTTTCGGGATTATTTCAACGAAAAAGTGACGCACCGGTGACCACGCCAGCAGAGCTGGCGGATGCCATCGGGCTGACGTATGACACCTATACAGGAAAGCAGATca                                                                                         >  4:781928/1‑149 (MQ=255)
                                                  tcGACCGTACACAGTCGTCTGTCGGCGGAGGTAAATAGTGTTCTTTTCGGGATTATTTCAACGAAAAAGTGACGCACCGGTGACCACGCCAGCTCTGCTGGCgtg                                                                                                                                     >  3:662729/1‑103 (MQ=255)
                                                  tcGACCGTACACAGTCGTCTGTCGGCGGAGGTAAATAGTGTTCTTTTCGGGATTATTTCAACGAAAAAGTGACGCACCGGTGACCACGCCAGCTCTGCTGGCgtg                                                                                                                                     <  4:662729/105‑3 (MQ=255)
                                                  tcGACCGTACACAGTCGTCTGTCGGCGGAGGTAAATAGTGTTCTTTTCGGGATTATTTCAACGAAAAAGTGACGCACCGGTGACCACGCCAGCAGAGCTGGCGGATGCCATCGGGCTGTCGTATGACACCTATACCGGAAAGCa                                                                                              <  3:17075/144‑1 (MQ=255)
                                                  tcGACCGTACACAGTCGTCTGTCGGCGGAGGTAAATAGTGTTCTTTTCGGGATTATTTCAACGAAAAAGTGACGCACCGGTGACCACGCCAGCAGAGCTGGCGGATGCCATCGGGCTGTCGTATGACACCTATACCGGAAAGCa                                                                                              >  4:17075/1‑144 (MQ=255)
                                                  tcGACCGTACACAGTCGTCTGTCGGCGGAGGTAAATAGTGTTCTTTTCGGGATTATTTCAACGAAAAAGTGACGCACCGGTGACCACGCCAGCAGAGCTGGCGGATGCCATCGGGCTGTCGTATGACACCTATACCGGAAAGCAGATca                                                                                         >  1:1039431/1‑149 (MQ=255)
                                                  tcGACCGTACACAGTCGTCTGTCGGCGGAGGTAAATAGTGTTCTTTTCGGGATTATTTCAACGAAAAAGTGACGCACCGGTGACCACGCCAGCAGAGCTGGCGGATGCCATCGGGCTGTCGTATGACACCTATACCGGAAAGCAGATca                                                                                         <  1:85356/149‑1 (MQ=255)
                                                  tcGACCGTACACAGTCGTCTGTCGGCGGAGGTAAATAGTGTTCTTTTCGGGATTATTTCAACGAAAAAGTGACGCACCGGTGACCACGCCAGCAGAGCTGGCGGATGCCATCGGGCTGTCGTATGACACCTATACCGGAAAGCAGATca                                                                                         >  4:214522/1‑149 (MQ=255)
                                                  tcGACCGTACACAGTCGTCTGTCGGCGGAGGTAAATAGTGTTCTTTTCGGGATTATTTCAACGAAAAAGTGACGCACCGGTGACCACGCCAGCAGAGCTGGCGGATGCCATCGGGCTGTCGTATGACACCTATACCGGAAAGCAGATca                                                                                         >  3:759496/1‑149 (MQ=255)
                                                   cGACCGTACACAGTCGTCTGTCGGCGGAGGTAAATAGTGTTCTTTTCGGGATTATTTCAACGAAAAAGTGACGCACCGGTGACCACGCCAGCAGAGCTGGCGGATGCCATCGGGCTGTCGTATGACACCTATACCGGAAAGCAGATcag                                                                                        <  1:329993/149‑1 (MQ=255)
                                                      ccGTACACAGTCGTCTGTCGGCGGAGGTAAATAGTGTTCTTTTCGGGATTATTTCAACGAAAAAGTGACGCACCGGTGACCACGCCAGCAGAGCTGGCGGATGCCATCGGGCTGTCGTATGACACCTATACCGGAAAGCAGATcagcag                                                                                     >  1:360201/1‑149 (MQ=255)
                                                        gTACACAGTCGTCTGTCGGCGGAGGTAAATAGTGTTCTTTTCGGGATTATTTCAACGAAAAAGTGACGCACCGGTGACCACGCCAGCAGAGCTGGTGGATGCCATCGGGCTGTCGTATGACACCTATACCGGAAAGCAGATCAGCAGTc                                                                                   <  3:355109/149‑1 (MQ=255)
                                                        gTACACAGTCGTCTGTCGGCGGAGGTAAATAGTGTTCTTTTCGGGATTATTTCAACGAAAAAGTGACGCACCGGTGACCACGCCAGCAGAGCTGGTGGATGCCATCGGGCTGTCGTATGACACCTATACCGGAAAGCAGATCAGCAGTc                                                                                   >  4:355109/1‑149 (MQ=255)
                                                            acaGTCGTCTGTCGGCGGAGGTAAATAGTGTTCTTTTCGGGATTATTTCAACGAAAAAGTGACGCACCGGTGACCACGCCAGCTCTGCTGGCgtg                                                                                                                                     >  3:466902/1‑93 (MQ=255)
                                                            acaGTCGTCTGTCGGCGGAGGTAAATAGTGTTCTTTTCGGGATTATTTCAACGAAAAAGTGACGCACCGGTGACCACGCCAGCAGAGCTGGCGGATGCCATCGGGCTGTCGTATGACACCTATACCGGAAAGCAGATCAGCAGTCAGCg                                                                               >  4:241267/1‑149 (MQ=255)
                                                            acaGTCGTCTGTCGGCGGAGGTAAATAGTGTTCTTTTCGGGATTATTTCAACGAAAAAGTGACGCACCGGTGACCACGCCAGCAGAGCTGGCGGATGCCATCGGGCTGTCGTATGACACCTATACCGGAAAGCAGATCAGCAGTCAGCg                                                                               >  2:259243/1‑149 (MQ=255)
                                                               gtcgtcTGTCGGCGGAGGTAAATAGTGTTCTTTTCGGGATTATTTCAACGAAAAAGTGACGCACCGGTGACCACGCCAGCAGAGCTGGCGGATGCCATCGGGCTGTCGTATGACACCTATACCGGAAAGCAGATCAGCAGTCAGCGGGc                                                                            >  3:205409/1‑149 (MQ=255)
                                                                tcgtcTGTCGGCGGAGGTAAATAGTGTTCTTTTCGGGATTATTTCAACGAAAAAGTGACGCACCGGTGACCACGCCAGCAGAGCTGGCGGATGCCATCGGGCTGTCGTATGACACCTATACCGGAAAGCAGATCAGCAGTCAGCGGGCt                                                                           >  4:406698/1‑149 (MQ=255)
                                                                    cTGTCGGCGGAGGTAAATAGTGTTCTTTTCGGGATTATTTCAACGAAAAAGTGACGCACCGGTGACCACGCCAGCAGAGCTGGCGGATGCCATCGGGCTGTCGTATGACACCTATACCGGAAAGCAGATCAGCAGTCAGCGGGCTATGc                                                                       <  3:428861/149‑1 (MQ=255)
                                                                    cTGTCGGCGGAGGTAAATAGTGTTCTTTTCGGGATTATTTCAACGAAAAAGTGACGCACCGGTGACCACGCCAGCAGAGCTGGCGGATGCCATCGGGCTGTCGTATGACACCTATACCGGAAAGCAGATCAGCAGTCAGCGGGCTATGc                                                                       <  3:8727/149‑1 (MQ=255)
                                                                      gTCGGCGGAGGTAAATAGTGTTCTTTTCGGGATTATTTCAACGAAAAAGTGACGCACCGGTGACCACGCCAGCAGAGCTGGCGGATGCCATCGGGCTGTCGTATGACACCTATACCGGAAAGCAGATCAGCAGTCAGCGGGCTATGCGa                                                                     <  2:321708/149‑1 (MQ=255)
                                                                      gTCGGCGGAGGTAAATAGTGTTCTTTTCGGGATTATTTCAACGAAAAAGTGACGCACCGGTGACCACGCCAGCAGAGCTGGCGGATGCCATCGGGCTGTCGTATGACACCTATACCGGAAAGCAGATCAGCAGTCAGCGGGCTATGCGa                                                                     <  4:529926/149‑1 (MQ=255)
                                                                                gTAAATAGTGTTCTTTTCGGGATTATTTCAACGAAAAAGTGACGCACCGGTGACCACGCCAGCAGAGCTGGCGGATGCCATCGGGCTGTCGTATGACACCTATACCGGAAAGCAGATCAGCAGTCAGCGGGCTATGCGACTGAcggcgg                                                           >  1:614137/1‑149 (MQ=255)
                                                                                 tAAATAGTGTTCTTTTCGGGATTATTTCAACGAAAAAGTGACGCACCGGTGACCACGCCAGCAGAGCTGGCGGATGCCATCGGGCTGTCGTATGACACCTATACCGGAAAGCAGATCAGCAGTCAGCGGGCTATGCGACTGACGGCGGt                                                          >  3:996858/1‑149 (MQ=255)
                                                                                  aaaTAGTGTTCTTTTCGGGATTATTTCAACGAAAAAGTGACGCACCGGTGACCACGCCAGCAGAGCTGGCGGATGCCATCGGGCTGTCGTATGACACCTATACCGGAAAGCAGATCAGCAGTCAGCGGGCTATGCGACTGACGGCGGtt                                                         <  2:117359/149‑1 (MQ=255)
                                                                                     tAGTGTTCTTTTCGGGATTATTTCAACGAAAAAGTGACGCACCGGTGACCACGCCAGCAGAGCTGGCGGATGCCATCGGGCTGTCGTATGACACCTATACCGGAAAGCAGATCATCAGTCAGCGGGCTATGCGACTGACGGCGGttttt                                                      <  3:700208/149‑1 (MQ=255)
                                                                                     tAGTGTTCTTTTCGGGATTATTTCAACGAAAAAGTGACGCACCGGTGACCACGCCAGCAGAGCTGGCGGATGCCATCGGGCTGTCGTATGACACCTATACCGGAAAGCAGATCAGCAGTCAGCGGGCTATGCGACTGACGGCGGttttt                                                      <  3:357741/149‑1 (MQ=255)
                                                                                     tAGTGTTCTTTTCGGGATTATTTCAACGAAAAAGTGACGCACCGGTGACCACGCCAGCAGAGCTGGCGGATGCCATCGGGCTGTCGTATGACACCTATACCGGAAAGCAGATCAGCAGTCAGCGGGCTATGCGACTGACGGCGGttttt                                                      <  3:274583/149‑1 (MQ=255)
                                                                                     tAGTGTTCTTTTCGGGATTATTTCAACGAAAAAGTGACGCACCGGTGACCACGCCAGCAGAGCTGGCGGATGCCATCGGGCTGTCGTATGACACCTATACCGGAAAGCAGATCAGCAGTCAGCGGGCTATGCGACTGACGGCGGttttt                                                      <  2:37900/149‑1 (MQ=255)
                                                                                      aGTGTTCTTTTCGGGATTATTTCAACGAAAAAGTGACGCACCGGTGACCACGCCAGCAGAGCTGGCGGATGCCATCGGGCTGTCGTATGACACCTATACCGGAAAGCAGATca                                                                                         >  1:522989/1‑113 (MQ=255)
                                                                                      aGTGTTCTTTTCGGGATTATTTCAACGAAAAAGTGACGCACCGGTGACCACGCCAGCAGAGCTGGCGGATGCCATCGGGCTGTCGTATGACACCTATACCGGAAAGCAGATca                                                                                         <  2:522989/113‑1 (MQ=255)
                                                                                       gtgtTCTTTTCGGGATTATTTCAACGAAAAAGTGACGCACCGGTGACCACGCCAGCAGAGCTGGCGGATGCCATCGGGCTGTCGTATGACACCTATACCGGAAAGCAGATCAGCAGTCAGCGGGCTATGCGACTGACGGCGGTTTTTTc                                                    <  1:719622/149‑1 (MQ=255)
                                                                                        tgtTCTTTTCGGGATTATTTCAACGAAAAAGTGACGCACCGGTGACCACGCCAGCAGAGCTGGCGGATGCCATCGGGCTGTCGTATGACACCTATACCGGAAAGCAGATCAGCAGTCAGCgg                                                                              >  2:523636/1‑122 (MQ=255)
                                                                                        tgtTCTTTTCGGGATTATTTCAACGAAAAAGTGACGCACCGGTGACCACGCCAGCAGAGCTGGCGGATGCCATCGGGCTGTCGTATGACACCTATACCGGAAAGCAGATCAGCAGTCAGCgg                                                                              <  1:523636/122‑1 (MQ=255)
                                                                                        tgtTCTTTTCGGGATTATTTCAACGAAAAAGTGACGCACCGGTGACCACGCCAGCAGAGCTGGCGGATGCCATCGGGCTGTCGTATGACACCTATACCGGAAAGCAGATCAGCAGTCAGCGGGCTATGCGACTGACGGCGGTTTTTTcc                                                   >  2:232164/1‑149 (MQ=255)
                                                                                         gtTCTTTTCGGCATTATTTCAACAAAAAAGTGACGCACCGGTGACCACGCCAGCAGAGCTGGCGGATGCCATCGGGCTGTCGTATGACACCTATACCGGAAAGCAGAGCAGCATTCAGCGGGCTATGCGACTGACGGCGGTTTTTTGCt                                                  <  4:926319/149‑1 (MQ=255)
                                                                                          ttCTTTTCGGGATTATTTCAACGAAAAAGTGACGCACCGGTGACCACGCCAGCAGAGCTGGCGGATGCCATCGGGCTGTCGTATGACACCTATACCGGAAAGCAGATCAGCAGTCAGCGGGCTATGCGACTGACGGCGGTTTTTTCCTg                                                 >  3:302671/1‑149 (MQ=255)
                                                                                          ttCTTTTCGGGATTATTTCAACGAAAAAGTGACGCACCGGTGACCACGCCAGCAGAGCTGGCGGATGCCATCGGGCTGTCGTATGACACCTATACCGGAAAGCAGATCAGCAGTCAGCGGGCTATGCGACTGACGGCGGTTTTTTCCTg                                                 >  3:287976/1‑149 (MQ=255)
                                                                                          ttCTTTTCGGGATTATTTCAACGAAAAAGTGACGCACCGGTGACCACGCCAGCAGAGCTGGCGGATGCCATCGGGCTGTCGTATGACACCTATACCGGAAAGCAGATCAGCAGTCAGCGGGCTATGCGACTGACGGCGGTTTTTTCCTg                                                 >  3:884357/1‑149 (MQ=255)
                                                                                             ttttCGGGATTATTTCAACGAAAAAGTGACGCACCGGTGACCACGCCAGCAGAGCTGGCGGATACCATCGGGCTGTCGTATGACACCTATACCGGAAAGCAGATCAGCAGTCAGCGGGCTATGCGACTGACGGCGGTTTTTTCCTGCGt                                              >  4:262848/1‑149 (MQ=255)
                                                                                                 cGGGATTATTTCAACGAAAAAGTGACGCACCGGTGACCACGCCAGCAGAGCTGGCGGATGCCATCGGGCTGTCGTATGACACCTATACCGGAAAGCAGATCAGCAGTCAGCGGGCTATGAGACTGACGGCGGTTTTTTCCTGCGTCata                                          <  3:671871/149‑3 (MQ=255)
                                                                                                      ttattTCAACGAAAAAGTGACGCACCGGTGACCACGCCAGCAGAGCTGGCGGATGCCATCGGGCTGTCGTATGACACCTATACCGGAAAGCAGATCAGCAGTCAGCGGGCTATGCGACTGACGGCGGTTTTTTCCTGCGTCAGAGTGCt                                     >  4:631810/1‑149 (MQ=255)
                                                                                                       tattTCAACGAAAAAGTGACGCACCGGTGACCACGCCAGCAGAGCTGGCGGACGCCATCGGGCTGTCGTATGACACCTATACCGGAAAGCAGATCATCAGTCAGCGGGCTATGCGACTGACGGCGGTTTTTTCCTGCGTCAGAGTGCTg                                    <  4:888619/149‑1 (MQ=255)
                                                                                                         tttCAACGAAAAAGTGACGCACCGGTGACCACGCCAGCAGAGCTGGCGGATGCCATCGGGCTGTCGTATGACACCTATACCGGAAAGCAGATCAGCAGTCAGCGGGCTATGCGACTGACGGCGGTTTTTTCCTGCGTCAGAGTGCTGGc                                  >  2:974275/1‑149 (MQ=255)
                                                                                                            cAACGAAAAAGTGACGCACCGGTGACCACGCCAGCAGAGCTGGCGGATGCCATCGGGCTGTCGTATGACACCTATACCGGAAAGCAGATCAGCAGTCAGCGGGCTATGCGACTGACGGCGGTTTTTTCCTGCGTCAGAGTGCTGGCaga                               <  3:414320/149‑1 (MQ=255)
                                                                                                              aCGAAAAAGTGACGCACCGGTGACCACGCCAGCAGAGCTGGCGGATGCCATCGGGCTGTCGTATGACACCTATACCGGAAAGCAGATCAGCAGTCAGCGGGCTATGCGACTGACGGCGGTTTTTTCCTGCGTCAGAGTGCTGGCAGAgt                             >  4:311255/1‑149 (MQ=255)
                                                                                                                gAAAAAGTGACGCACCGGTGACCACGCCAGCAGAGCTGGCGGATGCCATCGGGCTGTCGTATGACACCTATACCGGAAAGCAGATCAGCAGTCAGCGGGCTATGCGACTGACGGCGGTTTTTTCCTGCGTCAGAGTGCTGGCAGAgtcg                           >  3:760244/1‑149 (MQ=255)
                                                                                                                gAAAAAGTGACGCACCGGTGACCACGCCAGCAGAGCTGGCGGATGCCATCGGGCTGTCGTATGACACCTATACCGGAAAGCAGATCAGCAGTCAGCGGGCTATGCGACTGACGGCGGTTTTTTCCTGCGTCAGAGTGCTGGCAGAgtcg                           >  3:190563/1‑149 (MQ=255)
                                                                                                                   aaaGTGACGCACCGGTGACCACGCCAGCAGAGCTGGCGGATGCCATCGGGCTGTCGTATGACACCTATACCGGAAAGCAGATCAGCAGTCAGCGGGCTATGCGACTGACGGCGGTTTTTTCCTGCGTCAGAGTGCTGGCAGAgtcggtc                        <  1:970463/149‑1 (MQ=255)
                                                                                                                      gTGACGCACCGGTGACCACGGCAGCAGAGCTGGCGGATGCCATCGGGCTGTCGTATGACACCTATACCGGAAAGCAGATCAGCAGTCAGCGGGCTATGCGACTGACGGCGGTTTTTTCCTGCGTCAGAGTGCTGGCAGAGTCGGTCggg                     <  3:745545/149‑1 (MQ=255)
                                                                                                                      gTGACGCACCGGTGACCACGCCAGCAGAGCTGGCGGATGCCATCGGGCTGTCGTATGACACCTATACCGGAAAGCAGATCAGCAGTCAGCGGGCTATGCGACTGACGGCGGTTTTTTCCTGCGTCAGAGTGCTGGCAGAGTCGGTCggg                     <  1:979934/149‑1 (MQ=255)
                                                                                                                      gTGACGCACCGGTGACCACGCCAGCAGAGCTGGCGGATGCCATCGGGCTGTCGTATGACACCTATACCGGAAAGCAGATCAGCAGTCAGCGGGCTATGCGACTGACGGCGGTTTTTTCCTGCGTCAGAGTGCTGGCAGAGTCGGTCggg                     <  3:322348/149‑1 (MQ=255)
                                                                                                                         aCGCACCGGTGACCACGCCAGCAGAGCTGGCGGATGCCATCGGGCTGTCGTATGACACCTATACCGGAAAGCAGATCAGCAGTCAGCGGGCTATGCGACTGACGGCGGTTTTTTCCTGCGTCAGAGTGCTGGCAGAGTCGGTCGGGATg                  >  2:647000/1‑149 (MQ=255)
                                                                                                                            cACCGGTGACCACGCCAGCAGAGCTGGCGGATGCCATCGGGCTGTCGTATGACACCTATACCGGAAAGCAGATCAGCAGTCAGCGGGCTATGCGACTGACGGCGGTTTTTTCCTGCGTCAGAGTGCTGGCAGAGTCGGTCGGGATGTTg               >  2:277695/1‑149 (MQ=255)
                                                                                                                            cACCGGTGACCACGCCAGCAGAGCTGGCGGATGCCATCGGGCTGTCGTATGACACCTATACCGGAAAGCAGATCAGCAGTCAGCGGGCTATGCGACTGACGGCGGTTTTTTCCTGCGTCAGAGTGCTGGCAGAGTCGGTCGGGATGTTg               <  1:499349/149‑1 (MQ=255)
                                                                                                                            cACCGGTGACCACGCCAGCAGAGCTGGCGGATGCCATCGGGCTGTCGTATGACACCTATACCGGAAAGCAGATCAGCAGTCAGCGGGCTATGCGACTGACGGCGGTTTTTTCCTGCGTCAGAGAGCTGGCAGAGTCGGTCGGGATGTTg               >  3:986160/1‑149 (MQ=255)
                                                                                                                             aCCGGTGACCACGCCAGCAGAGCTGGCGGATGCCATCGGGCTGTCGTATGACACCTATACCGGAAAGCAGATCAGCAGTCAGCGGGCTATGCGACTGACGGCGGTTTTTTCCTGCGTCAGAGTGCTGGCAGAGTCGGTCGGGATGTTGc              <  3:123011/149‑1 (MQ=255)
                                                                                                                             aCCGGTGACCACGCCAGCAGAGCTGGCGGATGCCATCGGGCTGTCGTATGACACCTATACCGGAAAGCAGATCAGCAGTCAGCGGGCTATGCGACTGACGGCGGTTTTTTCCTGCGTCAGAGTGCTGGCAGAGTCGGTCGGGATGTTGc              >  3:984379/1‑149 (MQ=255)
                                                                                                                              ccGGTGACCACGCCAGCAGAGCTGGCGGATGCCATCGGGCTGTCGTATGACACCTATACCGGAAAGCAGATCAGCAGTCAGCGGGCTATGCGACTGACGGCGGTTTTTTCCTGCGTCAGAGTGCTGGCAGAGTCGGTCGGGATGTTGcc             <  1:35464/149‑1 (MQ=255)
                                                                                                                                ggTGACCACGCCAGCAGAGCTGGCGGATGCCATCGGGCTGTCGTATGACACCTATACCGGAAAGCAGATCAGCAGTCAGCGGGCTATGCGACTGACGGCGGTTTTTTCCTGCGTCAGAGTGCTGGCAGAGTCGGTCGGGATGTTGCCCt           <  3:542993/149‑1 (MQ=255)
                                                                                                                                ggTGACCACGCCAGCAGAGCTGGCGGATGCCATCGGGCTGTCGTATGACACCTATACCGGAAAGCAGATCAGCAGTCAGCGGGCTATGCGACTGACGGCGGTTTTTTCCTGCGTCAGAGTGCTGGCAGAGTCGGTCGGGATGTTGCCCt           <  2:482850/149‑1 (MQ=255)
                                                                                                                                    aCCACGCCAGCAGAGCTGGCGGATGCCATCGGGCTGTCGTATGACACCTATACCGGAAAGCAGATCAGCAGTCAGCGGGCTATGCGACTGACGGCGGTTTTTTCCTGCGTCAGAGTGCTGGCAGAGTCGGTCGGGATGTTGCCCTGCaa       >  1:198627/1‑149 (MQ=255)
                                                                                                                                     atccGCCAGCTCTGCTGGCGGATGCCATCGGGCTGTCGTATGACACCTATACCGGAAAGCAGATCAGCAGTCAGCGGGCTATGCGACTGACGGCGGTTTTTTCCTGCGt                                              <  3:415583/106‑1 (MQ=255)
                                                                                                                                      cACGCCAGCAGAGCTGGCGGATGCCATCGGGCTGTCGTATGACACCTATACCGGAAAGCAGATCAGCAGTCAGCGGGCTATGCGACTGACGGCGGTTTTTTCCTGCGTCAGAGTGCTGGCAGAGTCGGTCGGGATGTTGCCCTGCAATc     <  2:498631/149‑1 (MQ=255)
                                                                                                                                       aCGCCAGCAGAGCTGGCGGATGCCATCGGGCTGTCGTATGACACCTATACCGGAAAGCAGATCAGCAGTCAGCGGGCTATGCGACTGACGGCGGTTTTTTCCTGCGTCAGAGTGCTGGCAGAGTCGGTCGGGATGTTGCCCTGCAATCt    >  1:567928/1‑149 (MQ=255)
                                                                                                                                         gCCAGCAGAGCTGGCGGATGCCATCGGGCTGTCGTATGACACCTATACCGGAAAGCAGATCAGCAGTCAGCGGGCTATGCGACTGACGGCGGTTTTTTCCTGCGTCAGAGTGCTGGCAGAGTCGGTCGGGATGTTGCCCTGCAATCTGt  <  1:578485/149‑1 (MQ=255)
                                                                                                                                            agcagAGCTGGCGGATGCCATCGGGCTGTCGTATGACACCTATACCGGAAAGCAGATCAGCAGTCAGCGGGCTATGCGACTGACGGCGGTTTTTTCCTGCGTCAGAGTGCTGGCAGAGTCGGTCGGGa                    <  3:1035338/128‑1 (MQ=255)
                                                                                                                                            agcagAGCTGGCGGATGCCATCGGGCTGTCGTATGACACCTATACCGGAAAGCAGATCAGCAGTCAGCGGGCTATGCGACTGACGGCGGTTTTTTCCTGCGTCAGAGTGCTGGCAGAGTCGGTCGGGa                    >  4:1035338/1‑128 (MQ=255)
                                                                                                                                            agcagAGCTGGCGGATGCCATCGGGCTGTCGTATGACACCTATACCGGAAAGCAGATCAGCAGTCAGCGGGCTATGCGACTGACGGCGGTTTTTTCCTGCGTCAGAGTGCTGGCAGAGTCGGTCGGGATATTg               >  2:568752/1‑133 (MQ=255)
                                                                                                                                            agcagAGCTGGCGGATGCCATCGGGCTGTCGTATGACACCTATACCGGAAAGCAGATCAGCAGTCAGCGGGCTATGCGACTGACGGCGGTTTTTTCCTGCGTCAGAGTGCTGGCAGAGTCGGTCGGGATATTg               <  1:568752/133‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                             
TGTTGCCGGGGCGCTGTGCCTGTTCTGGTCGTGGCTGGTGGCGCGATATCTCGACCGTACACAGTCGTCTGTCGGCGGAGGTAAATAGTGTTCTTTTCGGGATTATTTCAACGAAAAAGTGACGCACCGGTGACCACGCCAGCAGAGCTGGCGGATGCCATCGGGCTGTCGTATGACACCTATACCGGAAAGCAGATCAGCAGTCAGCGGGCTATGCGACTGACGGCGGTTTTTTCCTGCGTCAGAGTGCTGGCAGAGTCGGTCGGGATGTTGCCCTGCAATCTGT  >  NC_000913/1205644‑1205929

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