Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,205,789 A→T 6.5% pseudogene (59/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,7890AT 6.5% 314.9 / 8.8 93pseudogene (59/412 nt)beeEpseudogene, portal protein family, e14 prophage;Phage or Prophage Related
Reads supporting (aligned to +/- strand):  ref base A (40/46);  new base T (3/3);  total (44/49)
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.92e-01

TGTTGTTGCCGGGGCGCTGTGCCTGTTCTGGTCGTGGCTGGTGGCGCGATATCTCGACCGTACACAGTCGTCTGTCGGCGGAGGTAAATAGTGTTCTTTTCGGGATTATTTCAACGAAAAAGTGACGCACCGGTGACCACGCCAGCAGAGCTGGCGGATGCCATCGGGCTGTCGTATGACACCTATACCGGAAAGCAGATCAGCAGTCAGCGGGCTATGCGACTGACGGCGGTTTTTTCCTGCGTCAGAGTGCTGGCAGAGTCGGTCGGGATGTTGCCCTGCAATCTGTAT  >  NC_000913/1205641‑1205931
                                                                                                                                                    |                                                                                                                                              
tgttgttgCCGGGGCGCTGTGCCTGTTCTGGTCGTGGCTGGTGGCGCGATATCTCGACCGTACACAGTCGTCTGTCGGCGGAGGTAAATAGTGTTCTTTTCGGGATTATTTCAACGAAAAAGTGACGCACCGGTGACCACGCCAGCaga                                                                                                                                                <  3:670761/149‑1 (MQ=255)
   tgttgCCGGGGCGCTGTGCCTGTTCTGGTCGTGGCTGGTGGCGCGATCTCTCGACCGTACACAGTCGTCTGTCGGCGGAGGTAAATAGTGTTCTTTTCGGGATTATTTCAACGAAAAAGTGACGCACCGGTGACCACGCCAGCAGAGCt                                                                                                                                             <  3:807833/149‑1 (MQ=255)
        ccGGGGCGCTGTGCCTGTTCTGGTCGTGGCTGGTGGCGCGATATCTCGACCGTACACAGTCGTCTGTCGGCGGAGGTAAATAGTGTTCTTTTCGGGATTATTTCAACGAAAAAGTGACGCACCGGTGACCACGCCAGCTCTGCTGGCgt                                                                                                                                        <  4:739235/149‑2 (MQ=255)
         cGGGGCGCTGTGCCTGTTCTGGTCGTGGCTGGTGGCGCGATATCTCGACCGTACACAGTCGTCTGTCGGCGGAGGTAAATAGTGTTCTTTTCGGGATTATTTCAACGAAAAAGTGACGCACCGGTGACCACGCCAGCTCTGCTGGCgt                                                                                                                                        <  1:97967/148‑2 (MQ=255)
          ggggCGCTGTGCCTGTTCTGGTCGTGGCTGGTGGCGCGATATCTCGACCGTACACAGTCGTCTGTCGGCGGAGGTAAATAGTGTTCTTTTCGGGATTATTTCAACGAAAAAGTGACGCACCGGTGACCACGCCAGCAGAGCTGGCGGAt                                                                                                                                      >  4:843490/1‑149 (MQ=255)
          ggggCGCTGTGCCTGTTCTGGTCGTGGCTGGTGGCGCGATATCTCGACCGTACACAGTCGTCTGTCGGCGGAGGTAAATAGTGTTCTTTTCGGGATTATTTCAACGAAAAAGTGACGCACCGGTGACCACGCCAGCAGAGCTGGCGGAt                                                                                                                                      <  3:458817/149‑1 (MQ=255)
                  gtgCCTGTTCTGGTCGTGGCTGGTGGCGCGATATCTCGACCGTACACAGTCGTCTGTCGGCGGAGGTAAATAGTGTTCTTTTCGGGATTATTTCAACGAAAAAGTGACGCACCGGTGACCACGCCAGCAGAGCTGGCGGATGCCATCg                                                                                                                               <  2:633964/148‑1 (MQ=255)
                     ccTGTTCTGGTCGTGGCTGGTGGCGCGATATCTCGACCGTACACAGTCGTCTGTCGGCGGAGGTAAATAGTGTTCTTTTCGGGATTATTTCAACGAAAAAGTGACGCACCGGTGACCACGCCAGCAGAGCTGGCGGGTGCCATCGGGCt                                                                                                                           >  3:57198/1‑149 (MQ=255)
                      cTGTTCTGGTCGTGGCTGGTGGCGCGATATCTCGACCGTACACAGTCGTCTGTCGGCGGAGGTAAATAGTGTTCTTTTCGGGATTATTTCAACGAAAAAGTGACGCACCGGTGACCACGCCAGCAGAGCTGGCGGATGCCATCGGGCt                                                                                                                           >  2:774707/1‑148 (MQ=255)
                        gTTCTGGTCGTGGCTGGTGGCGCGATATCTCGACCGTACACAGTCGTCTGTCGGCGGAGGTAAATAGTGTTCTTTTCGGGATTATTTCAACGAAAAAGTGACGCACCGGTGACCACGCCAGCAGAGCTGGCGGATGCCATCGGGCTGt                                                                                                                         <  1:202620/148‑1 (MQ=255)
                        gTTCTGGTCGTGGCTGGTGGCGCGATATCTCGACCGTACACAGTCGTCTGTCGGCGGAGGTAAATAGTGTTCTTTTCGGGATTATTTCAACGAAAAAGTGACGCACCGGTGACCACGCCAGCAGAGCTGGCGGATGCCATCGGGCTGt                                                                                                                         >  2:96330/1‑148 (MQ=255)
                         ttCTGGTCGTGGCTGGTGGCGCGATATCTCGACCGTACACAGTCGTCTGTCGGCGGAGGTAAATAGTGTTCTTTTCGGGATTATTTCAACGAAAAAGTGACGCACCGGTGACCACGCCAGCAGAGCTGGCGGATGCCATCGGGCTGTc                                                                                                                        >  2:664618/1‑148 (MQ=255)
                            tGGTCGTGGCTGGTGGCGCGATATCTCGACCGTACACAGTCGTCTGTCGGCGGAGGTAAATAGTGTTCTTTTCGGGATTATTTCAACGAAAAAGTGACGCACCGGTGACCACGCCAGCAGAGCTGGCGGATGCCATCGGGCTGTCGTa                                                                                                                     <  2:190850/148‑1 (MQ=255)
                            tGGTCGTGGCTGGTGGCGCGATATCTCGACCGTACACAGTCGTCTGTCGGCGGAGGTAAATAGTGTTCTTTTCGGGATTATTTCAACGAAAAAGTGACGCACCGGTGACCACGCCAGCAGAGCTGGCGGATGCCATCGGGCTGTCGTa                                                                                                                     <  1:65118/148‑1 (MQ=255)
                            tGGTCGTGGCTGGTGGCGCGATATCTCGACCGTACACAGTCGTCTGTCGGCGGAGGTAAATAGTGTTCTTTTCGGGATTATTTCAACGAAAAAGTGACGCACCGGTGACCACGCCAGCAGAGCTGGCGGATGCCATCGGGCTGTCGTAt                                                                                                                    <  3:197971/149‑1 (MQ=255)
                            tGGTCGTGGCTGGTGGCGCGATATCTCGACCGTACACAGTCGTCTGTCGGCGGAGGTAAATAGTGTTCTTTTCGGGATTATTTCAACGAAAAAGTGACGCACCGGTGACCACGCCAGCAGAGCTGGCGGATGCCATCGGGCTGTCGTAt                                                                                                                    <  3:347859/149‑1 (MQ=255)
                                cGTGGCCGGTGGCGCGATATCTCGACCGCTAACAGTCGTCTGTCGGCGGAGCTAAATAGTGTTCTGTTCGGGATTATTTCAACGAAAAAGTGACGCACCAGTGACCACGCCAGCAGAGCTGGCGGATGCCATCGGGCTGTCGTAAGaca                                                                                                                >  4:80598/1‑149 (MQ=255)
                                  tGGCTGGTGGCGCGATATCTCGACCGTACACAGTCGTCTGTCGGCGGAGGTAAATAGTGTTCTTTTCGGGATTATTTCAACGAAAAAGTGACGCACCGGTGACCACGCCAGCAGAGCTGGCGGATGCCATCGGGCTGTCGTATGacac                                                                                                               >  1:91163/1‑148 (MQ=255)
                                    gCTGGTGGCGCGATATCTCGACCGTACACAGTCGTCTGTCGGCGGAGGTAAATAGTGTTCTTTTCGGGATTATTTCAACGAAAAAGTGACGCACCGGTGACCACGCCAGCAGAGCTGGCGGATGCCATCGGGCTGTCGTATGACACCt                                                                                                             >  2:163789/1‑148 (MQ=255)
                                        gtggCGCGATATCTCGACCGTACACAGTCGTCTGTCGGCGGAGGTAAATAGTGTTCTTTTCGGGATTATTTCAACGAAAAAGTGACGCACCGGTGACCACGCCAGCAGAGCTGGCGGATGCCATCGGGCTGTCGTATGACACCTATAcc                                                                                                        <  3:226042/149‑1 (MQ=255)
                                         tggCGCGATATCTCGACCGTACACAGTCGTCTGTCGGCGGAGGTAAATAGTGTTCTTTTCGGGATTATTTCAACGAAAAAGTGACGCACCGGTGACCACGCCAGCAGAGCTGGCGGATGCCATCGGGCTGTCGTATGACACCTATAcc                                                                                                        <  1:774707/148‑1 (MQ=255)
                                         tggCGCGATATCTCGACCGTACACAGTCGTCTGTCGGCGGAGGTAAATAGTGTTCTTTTCGGGATTATTTCAACGAAAAAGTGACGCACCGGTGACCACGCCAGCAGAGCTGGCGGATGCCATCGGGCTGTCGTATGACACCTATAcc                                                                                                        >  1:206734/1‑148 (MQ=255)
                                           gcgcGATATCTCGACCGTACACAGTCGTCTGTCGGCGGAGGTAAATAGTGTTCTTTTCGGGATTATTTCAACGAAAAAGTGACGCACCGGTGACCACGCCAGCAGAGCTGGCGGATGCCATCGGGCTGTCGTATGACACCTATACCGGa                                                                                                     <  3:291827/149‑1 (MQ=255)
                                              cgATATCTCGACCGTACACAGTCGTCTGTCGGCGGAGGTAAATAGTGTTCTTTTCGGGATTATTTCAACGAAAAAGTGACGCACCGGTGACCACGCCAGCAGAGCTGGCGGATGCCATCGGGCTGTCGTATGACACCTATACCGGaa                                                                                                    >  1:12880/1‑147 (MQ=255)
                                              cgATATCTCGACCGTACACAGTCGTCTGTCGGCGGAGGTAAATAGTGTTCTTTTCGGGATTATTTCAACGAAAAAGTGACGCACCGGTGACCACGCCAGCAGAGCTGGCGGATGCCATCGGGCTGTCGTATGACACCTATACCGGaa                                                                                                    <  2:12880/147‑1 (MQ=255)
                                               gATATCTCGACCGTACACAGTCGTCTGTCGGCGGAGGTAAATAGTGTTCTTTTCGGGATTATTTCAACGAAAAAGTGACGCACCGGTGACCACGCCAGCAGAGCTGGCGGATGCCATCGGGCTGTCGTATGACACCTATACCGGAAAg                                                                                                  >  2:782352/1‑148 (MQ=255)
                                                atatCTCGACCGTACACAGTCGTCTGTCGGCGGAGGTAAATAGTGTTCTTTTCGGGATTATTTCAACGAAAAAGTGACGCACCGGTGACCACGCCAGCAGAGCTGGCGGATGCCATCGGGCTGTCGTATGACACCTATACCGGAAAGCa                                                                                                <  3:398599/149‑1 (MQ=255)
                                                atatCTCGACCGTACACAGTCGTCTGTCGGCGGAGGTAAATAGTGTTCTTTTCGGGATTATTTCAACGAAAAAGTGACGCACCGGTGACCACGCCAGCAGAGCTGGCGGATGCCATCGGGCTGTCGTATGACACCTATACCGGAAAGCa                                                                                                <  3:843490/149‑1 (MQ=255)
                                                atatCTCGACCGTACACAGTCGTCTGTCGGCGGAGGTAAATAGTGTTCTTTTCGGGATTATTTCAACGAAAAAGTGACGCACCGGTGACCACGCCAGCAGAGCTGGCGGATGCCATCGGGCTGTCGTATGACACCTATACCGGAAAGCa                                                                                                >  4:949071/1‑149 (MQ=255)
                                                  atCTCGACCGTACACGATCGTCTGTCGGCGGCGGTAAATAGTGTTCTTATCAGGATTATTTCAACGAAAAAGTGACGCACCAGTGACTACGCCAGCAGGGCTGGCGGATGCCATCGGGCTGTCGTATGACACGTATACCGGAACGCAGa                                                                                              >  4:856878/1‑149 (MQ=255)
                                                   tctcGACCGTACACAGTCGTCTGTCGGCGGAGGTAAATAGTGTTCTTTTCGGGATTATTTCAACGAAAAAGTGACGCACCGGTGACCACGCCAGCAGAGCTGGCGGATGCCATCGGGCTGTCGTATGACACCTATACCGGAAAGCAGAt                                                                                             <  4:221659/149‑1 (MQ=255)
                                                    ctcGACCGTACACAGTCGTCTGTCGGCGGAGGTAAATAGTGTTCTTTTCGGGATTATTTCAACGAAAAAGTGACGCACCGGTGACCACGCCAGCAGAGCTGGCGGATGCCATCGGGCTGTCGTATGACACCTATACCGGAAAGCAGAt                                                                                             >  2:306854/1‑148 (MQ=255)
                                                     tcGACCGTACACAGTCGTCTGTCGGCGGAGGTAAATAGTGTTCTTTTCGGGATTATTTCAACGAAAAAGTGACGCACCGGTGACCACGCCAGCAGAGCTGGCGGATGCCATCGGGCTGTCGTATGACACCTATACCGGAAAGCAGATca                                                                                           >  3:339455/1‑149 (MQ=255)
                                                     tcGACCGTACACAGTCGTCTGTCGGCGGAGGTAAATAGTGTTCTTTTCGGGATTATTTCAACGAAAAAGTGACGCACCGGTGACCACGCCAGCAGAGCTGGCGGATGCCATCGGGCTGTCGTATGACACCTATACCGGAAAGCAGATca                                                                                           <  3:949071/149‑1 (MQ=255)
                                                      cGACCGTACACAGTCGTCTGTCGGCGGAGGTAAATAGTGTTCTTTTCGGGATTATTTCAACGAAAAAGTGACGCACCGGTGACCACGCCAGCAGAGCTGGCGGATGCCATCGGGCTGTCGTATGACACCTATACCGGAAAGCAGATgag                                                                                          <  3:9626/149‑1 (MQ=255)
                                                      cGACCGTACACAGTCGTCTGTCGGCGGAGGTAAATAGTGTTCTTTTCGGGATTATTTCAACGAAAAAGTGACGCACCGGTGACCACGCCAGCAGAGCTGGCGGATGCCATCGGGCTGTCGTATGACACCTATACCGGAAAGCAGATca                                                                                           >  2:317594/1‑148 (MQ=255)
                                                      cGACCGTACACAGTCGTCTGTCGGCGGAGGTAAATAGTGTTCTTTTCGGGATTATTTCAACGAAAAAGTGACGCACCGGTGACCACGCCAGCAGAGCTGGCGGATGCCATCGGGCTGTCGTATGACACCTATACCGGAAAGCAGATca                                                                                           >  2:327209/1‑148 (MQ=255)
                                                      cGACCGTACACAGTCGTCTGTCGGCGGAGGTAAATAGTGTTCTTTTCGGGATTATTTCAACGAAAAAGTGACGCACCGGTGACCACGCCAGCAGAGCTGGCGGATGCCATCGGGCTGTCGTATGACACCTATACCGGAAAGCAGATca                                                                                           >  2:695017/1‑148 (MQ=255)
                                                               acaGTCGTCTGTCGGCGGAGGTAAATAGTGTTCTTTTCGGGATTATTTCAACGAAAAAGTGACGCACCGGTGACCACGCCAGCAGAGCTGGCGGATGCCATCGGGCTGTCGTATGACACCTATACCGGAAAGCAGATcagca                                                                                        >  4:432946/1‑142 (MQ=255)
                                                               acaGTCGTCTGTCGGCGGAGGTAAATAGTGTTCTTTTCGGGATTATTTCAACGAAAAAGTGACGCACCGGTGACCACGCCAGCAGAGCTGGCGGATGCCATCGGGCTGTCGTATGACACCTATACCGGAAAGCAGATcagca                                                                                        <  3:432946/142‑1 (MQ=255)
                                                                caGTCGTCTGTCGGCGGAGGTAAATAGTGTTCTTTTCGGGATTATTTCAACGAAAAAGTGACGCACCGGTGACCACGCCAGCAGAGCTGGCGGATGCCATCGGGCTGTCGTGTGACACCTATACCGGAAAGCAGATCAGCAGTCAGCg                                                                                 <  1:795026/148‑1 (MQ=255)
                                                                 aGTCGTCTGTCGGCGGAGGTAAATAGTGTTCTTTTCGGGATTATTTCAACGAAAAACTGACGCACCGGTGACCAAGCCAGCAGAGCTGGCGGATGCCATCGGGCTGTCGTATGACACCTATACCGGAAAGCAGATCAGCACTCAGCggg                                                                               >  3:798817/1‑149 (MQ=255)
                                                                  gtcgtcTGTCGGCGGAGGTAAATAGTGTTCTTTTCGGGATTATTTCAACGAAAAAGTGACGCACCGGTGACCACGCCAGCAGAGCTGGCGGATGCCATCGGGCTGTCGTATGACACCTATACCGGAAAGCAGATCAGCAGTCAGCGGGc                                                                              >  3:229782/1‑149 (MQ=255)
                                                                       cTGTCGGCGGAGGTAAATAGTGTTCTTTTCGGGATTATTTCAACGAAAAAGTGACGCACCGGTGACCACGCCAGCAGAGCTGGCGGATGCCATCGGGCTGTCGTATGACACCTATACCGGAAAGCAGATCAGCAGTCAGCGGGCTATGc                                                                         <  3:8109/149‑1 (MQ=255)
                                                                           cggcggAGGTAAATAGTGTTCTTTTCGGGATTATTTCAACGAAAAAGTGACGCACCGGTGACCACGCCAGCAGAGCTGGCGGATGCCATCGGGCTGTCGTATGACACCTATACCGGAAAGCAGATCAGCAGTCAGCGGGCTATGCGAc                                                                      <  1:122439/148‑1 (MQ=255)
                                                                                gAGGTAAATAGTGTTCTTTTCGGGATTATTTCAACGAAAAAGTGACGCACCGGTGACCACGCCAGCAGAGCTGGCGGATGCCATCGGGCTGTCGTATGACACCTATACCGGAAAGCAGATCAGCAGTCAGCGGGCTATGCGACTGAcgg                                                                <  3:759228/149‑1 (MQ=255)
                                                                                 aGGTAAATAGTGTTCTTTTCGGGATTATTTCAACGAAAAAGTGACGCACCGGTGACCACGCCAGCAGAGCTGGCGGATGCCATCGGGCTGTCGTATGACACCTATACCGGAAAGCAGATCAGCAGTCAGCGGGCTATGCGACTGAcgg                                                                >  1:906484/1‑148 (MQ=255)
                                                                                 aGGTAAATAGTGTTCTTTGCGGGATTATTTCAACGAAAAAGTGACGCACCGGTGACCACGCCAGCAGAGCTGGCGGATGCCATCGGGCTGTCGTATGACACCTATACCGGAAAGCAGATCAGCAGTCAGCGGGCTATGCGACTGAcgg                                                                <  1:535910/148‑1 (MQ=255)
                                                                                   gTAAATAGTGTTCTTTTCGGGATTATTTCAACGAAAAAGTGACGCACCGGTGACCACGCCAGCAGAGCTGGCGGATGCCATCGGGCTGTCGTATGACACCTATACCGGAAAGCAGATCAGCAGTCAGCGGGCTATGCGACTGAcggcgg                                                             >  3:10080/1‑149 (MQ=255)
                                                                                   gTAAATAGTGTTCTTTTCGGGATTATTTCAACGAAAAAGTGACGCACCGGTGACCACGCCAGCAGAGCTGGCGGATGCCATCGGGCTGTCGTATGACACCTATACCGGAAAGCAGATCAGCAGTCAGCGGGCTATGCGACTGAcggcgg                                                             >  3:412819/1‑149 (MQ=255)
                                                                                    tAAATAGTGTTCTTTTCGGGATTATTTCAACGAAAAAGTGACGCACCGGTGACCACGCCAGCAGAGCTGGCGGATGCCATCGGGCTGTCGTATGACACCTATACCGGAAAGCAGATCAGCAGTCAGCGGGCTATGCGACTGAcggcgg                                                             <  1:695017/148‑1 (MQ=255)
                                                                                        tAGTGTTCTTTTCGGGATTATTTCAACGAAAAAGTGACGCACCGGTGACCACGCCAGCAGAGCTGGCGGATGCCATCGGGCTGTCGTATGACACCTATACCGGAAAGCAGATCAGCAGTCAGCGGGCTATGCGACTGACGGCGGttttt                                                        >  3:434986/1‑149 (MQ=255)
                                                                                        tAGTGTTCTTTTCGGGATTATTTCAACGAAAAAGTGACGCACCGGTGACCACGCCAGCAGAGCTGGCGGATGCCATCGGGCTGTCGTATGACACCTATACCGGAAAGCAGATCAGCAGTCAGCGGGCTATGCGACTGACGGCGGtttt                                                         <  1:468701/148‑1 (MQ=255)
                                                                                         aGTGTTCTTTTCGGGATTATTTCAACGAAAAAGTGACGCACCGGTGACCACGCCAGCAGAGCTGGCGGATGCCATCGGGCTGTCGTATGACACCTATACCGGAAAGCAGATCAGCAGTCAGCGGGCTATGCGACTGACGGCGGtttttt                                                       <  3:758862/149‑1 (MQ=255)
                                                                                           tgtTCTTTTCGGGATTATTTCAACGAAAAAGTGACGCACCGGTGACCACGCCAGCAGAGCTGGCGGATGCCATCGGGCTGTCGTATGACACCTATACCGGAAAGCa                                                                                                >  4:258342/1‑106 (MQ=255)
                                                                                           tgtTCTTTTCGGGATTATTTCAACGAAAAAGTGACGCACCGGTGACCACGCCAGCAGAGCTGGCGGATGCCATCGGGCTGTCGTATGACACCTATACCGGAAAGCa                                                                                                <  3:258342/106‑1 (MQ=255)
                                                                                            gtTCTTTTCGGGATTATTTCAACGAAAAAGTGACGCACCGGTGACCACGCCAGCAGAGCTGGCGGATGCCATCGGGCTGTCGTATGACACCTATACCGGAAAGCAGATCAGCAGTCAGCGGGCTATGCGACTGACGGCGGTTTTTTCCt                                                    <  4:478033/149‑1 (MQ=255)
                                                                                             ttCTTTTCGGGATTATTTCAACGAAAAAGTGACGCACCGGTGACCACGCCAGCAGAGCTGGCGGATGCCATCGGGCTGTCGTATGACACCTATACCGGAAAGCAGATCAGCAGTCAGCGGGCTATGCGa                                                                       <  3:46541/129‑1 (MQ=255)
                                                                                             ttCTTTTCGGGATTATTTCAACGAAAAAGTGACGCACCGGTGACCACGCCAGCAGAGCTGGCGGATGCCATCGGGCTGTCGTATGACACCTATACCGGAAAGCAGATCAGCAGTCAGCGGGCTATGCGa                                                                       >  4:46541/1‑129 (MQ=255)
                                                                                                ttttCGGGATTATTTCAACGAAAAAGTGACGCACCGGTGACCACGCCAGCAGAGCTGGCGGATGCCATCGGGCTGTCGTATGACACCTATACCGGAAAGCAGATCAGCAGTCAGCGGGCTATGCGACTGACGGCGGTTTTTTCCTGCg                                                 <  2:528090/148‑1 (MQ=255)
                                                                                                ttttCGGGATTATTTCAACGAAAAAGTGACGCACCGGTGACCACGCCAGCAGAGCTGGCGGATGCCATCGGGCTGTCGTATGACACCTATACCGGAAAGCAGATCAGCAGTCAGCGGGCTATGCGACTGACGGCGGTTTTTTCCTGCGt                                                >  4:253145/1‑149 (MQ=255)
                                                                                                      ggATTATTTCAACGAAAAAGTGACGCACCGGTGACCACGCCAGCAGAGCTGGCGGATGCCATCGGGCTGTCGTATGACACCTATACCGGAAAGCAGATCAGCAGTCAGCGGGCTATGCGACTGAcggc                                                               <  3:512933/128‑1 (MQ=255)
                                                                                                      ggATTATTTCAACGAAAAAGTGACGCACCGGTGACCACGCCAGCAGAGCTGGCGGATGCCATCGGGCTGTCGTATGACACCTATACCGGAAAGCAGATCAGCAGTCAGCGGGCTATGCGACTGAcggc                                                               >  4:512933/1‑128 (MQ=255)
                                                                                                      ggATTATTTCAACGAAAAAGTGACGCACCGGTGACCACGCCAGCAGAGCTGGCGGATGCCATCGGGCTGTCGTATGACACCTATACCGGAAAGCAGATCAGCAGTCAGCGGGCTATGCGACGGACGGCGGTTTTTTCCTGCGTAAGAGt                                          >  4:825986/1‑149 (MQ=255)
                                                                                                        attattTCAACGAAAAAGTGACGCACCGGTGACCACGCCAGCAGAGCTGGCGGATGCCATCGGGCTGTCGTATGACACCTATACCGGAAAGCAGATCAGCAGTCAGCGGGCTATGCGACTGACGGCGGTTTTTTCCTGCGTCAGAGTGc                                        >  4:700830/1‑149 (MQ=255)
                                                                                                         ttattTCAACGAAAAAGTGACGCACCGGTGACCACGCCAGCAGAGCTGGCGGATGCCATCGGGCTGTCGTATGACACCTATACCGGAAAGCAGATCAGCAGTCAGCGGGCTATGCGACTGACGGCGGTTTTTTCCTGCGTCAGAGTGc                                        >  1:851527/1‑148 (MQ=255)
                                                                                                         ttattTCAACGAAAAAGTGACGCACCGGTGACCACGCCAGCAGAGCTGGCGGATGCCATCGGGCTGTCGTATGACACCTATACCGGAAAGCAGATCAGCAGTCAGCGGGCTATGCGACTGACGGCGGTTTTTTCCTGCGTCAGAGTGCt                                       <  4:229782/149‑1 (MQ=255)
                                                                                                             ttCAACGAAAAAGTGACGCACCGGTGACCACGCCAGCAGAGCTGGCGGATGCCATCGGGCTGTCGTATGACACCTATACCGGAAAGCAGATCAGCAGTCAGCGGGCTATGCGACTGACGGCGGTTTTTTCCTGCGTCAGAGTGCTGGCa                                   <  3:704078/149‑1 (MQ=255)
                                                                                                               cAACGAAAAAGTGACGCACCGGTGACCACGCCAGCAGAGCTGGCGGATGCCATCGGGCTGTCGTATGACACCTATACCGGAAAGCAGATCAGCAGTCAGCGGGCTATGCGACTGACGGCGGTTTTTTCCTGCGTCAGAGTGCTGGCag                                  <  2:667666/148‑1 (MQ=255)
                                                                                                               cAACGAAAAAGTGACGCACCGGTGACCACGCCAGCAGAGCTGGCGGATGCCATCGGGCTGTCGTATGACACCTATACCGGAAAGCAGATCAGCAGTCAGCGGGCTATGCGACTGACGGCGGTTTTTTCCTGCGTCAGAGTGCTGGCag                                  <  1:782352/148‑1 (MQ=255)
                                                                                                                aaCGAAAAAGTGACGCACCGGTGACCACGCCAGCAGAGCTGGCGGATGCCATCGGGCTGTCGTATGACACCTATACCGGAAAGCAGATCAGCAGTCAGCGGGCTATGCGACTGACGGCGGTTTTTTCCTGCGTCAGAGTGCTGGCagag                                >  3:781475/1‑149 (MQ=255)
                                                                                                                 aCGAAAAAGTGACGCACCGGTGACCACGCCAGCAGAGCTGGCGGATGCCATCGGGCTGTCGTATGACACCTATACCGGAAAGCAGATCAGCAGTCAGCGGGCTATGCGACTGACGGCGGTTTTTTCCTGCGTCAGAGTGCTGGCAGAgt                               >  4:748464/1‑149 (MQ=255)
                                                                                                                  cGAAAAAGTGACGCACCGGTGACCACGCCAGCAGAGCTGGCGGATGCCATCGGGCTGTCGTATGACACCTATACCGGAAAGCAGATCAGCAGTCAGCGGGCTATGCGACTGACGGCGGTTTTTTCCTGCGTCAGAGTGCTGGCAGAgtc                              <  3:169489/149‑1 (MQ=255)
                                                                                                                      aaaGTGACGCACCGGTGACCACGCCAGCAGAGCTGGCGGATGCCATCGGGCTGTCGTATGACACCTATACCGGAAAGCAGATCAGCAGTCAGCGGGCTATGCGACTGACGGCGGTTTTTTCCTGCGTCAGAGTGCTGGCAGAgtcggt                           <  2:647272/148‑1 (MQ=255)
                                                                                                                          tGACGCACCGGTGACCACGCCAGCAGAGCTGGCGGATGCCATCGGGCTGTCGTATGACACCTATACCGGaaa                                                                                                   >  2:267625/1‑72 (MQ=255)
                                                                                                                          tGACGCACCGGTGACCACGCCAGCAGAGCTGGCGGATGCCATCGGGCTGTCGTATGACACCTATACCGGaaa                                                                                                   <  1:267625/72‑1 (MQ=255)
                                                                                                                          tGACGCACCGGTGACCACGCCAGCAGAGCTGGCGGATGCCATCGGGCTGTCGTATGACACCTATACCGGAAAGCAGATCAGCAGTCAGCGGGCTATGCGACTGACGGCGGTTTTTTCCTGCGTCAGAGTGCTGGCAGAGTCGGTCGGGa                      <  3:916779/149‑1 (MQ=255)
                                                                                                                          tGACGCACCGGTGACCACGCCAGCAGAGCTGGCGGATGCCATCGGGCTGTCGTATGACACCTATACCGGAAAGCAGATCAGCAGTCAGCGGGCTATGCGACTGACGGCGGTTTTTTCCTGCGTCAGAGTGCTGGCAGAGTCGGTCGGGa                      >  4:935270/1‑149 (MQ=255)
                                                                                                                                aCCGGTGACCACGCCAGCAGAGCTGGCGGATGCCATCGGGCTGTCGTATGACACCTATACCGGAAAGCAGa                                                                                              >  3:824200/1‑71 (MQ=255)
                                                                                                                                aCCGGTGACCACGCCAGCAGAGCTGGCGGATGCCATCGGGCTGTCGTATGACACCTATACCGGAAAGCAGa                                                                                              <  4:824200/71‑1 (MQ=255)
                                                                                                                                 ccGGTGACCACGCCAGCAGAGCTGGCGGATGCCATCGGGCTGTCGTATGACACCTATACCGGAAAGCAGATCAGCAGTCAGCGGGCTATGCGACTGACGGCGGTTTTTTCCTGCGTCAGAGTGCTGGCAGAGTCGGTCGGGATGTTGc                >  1:152303/1‑148 (MQ=255)
                                                                                                                                      gACCACGCCAGCAGAGCTGGCGGATGCCATCGGGCTGTCGTATGACACCTATACCGGAAAGCAGATCAGCAGTCAGCGGGCTATGCGACTGACGGCGGTTTTTTCCTGCGTCAGAGTGCTGGCAGAGTCGGTCGGGATGTTGCCCTGc           >  2:596556/1‑148 (MQ=255)
                                                                                                                                         tccGCCAGCTCTGCTGGCGGATGCCATCGGGCTGTCGTATGACACCTATACCGGAAAGCAGATCAGCAGTCAGCGGGCTATg                                                                          >  2:488837/3‑82 (MQ=255)
                                                                                                                                         tccGCCAGCTCTGCTGGCGGATGCCATCGGGCTGTCGTATGACACCTATACCGGAAAGCAGATCAGCAGTCAGCGGGCTATGCGACTGACGGCGGTTTTTTCCTGCGTCAGAGTGCTg                                      >  3:482775/3‑118 (MQ=255)
                                                                                                                                         tccGCCAGCTCTGCTGGCGGATGCCATCGGGCTGTCGTATGACACCTATACCGGAAAGCAGATCAGCAGTCAGCGGGCTATGCGACTGACGGCGGTTTTTTCCTGCGTCAGAGTGCTg                                      <  4:482775/116‑1 (MQ=255)
                                                                                                                                         tccGCCAGCTCTGCTGGCGGATGCCATCGGGCTGTCGTATGACACCTATACCGGAAAGCAGATCAGCAGTCAGCGGGCTATGCGACTGACGGCGGTTTTTTCCTGCGTCAGAGTGCTGGCAGAGGCGGTCGGGATGTTGCGCTGCAATc       >  4:930148/3‑149 (MQ=255)
                                                                                                                                         cACGCCAGGAGAGCTGGCGGATGCCATCGGGCTGTCGTATGACACCTATACCGGAAAGCAGATCAGCAGTCAGCGGGCTATGCGGCTGACGGCGGTTTTTTCCTGCGTCAGAGTGCTGGCAGAGTCGGTCGGGATGTTGCCCTGCAAt        <  2:91163/148‑1 (MQ=255)
                                                                                                                                          aCGCCAGCAGAGCTGGCGGATGCCATCGGGCTGTCGTATGACACCTATACCGGAAAGCAGATCAGCAGTCAGCGGGCTATGCGACTGACGGCGGTTTTTTCCTGCGTCAGAGTGCTGGCAGAGTCGGTCGGGATGTTGCCCTGCAATc       >  1:538097/1‑148 (MQ=255)
                                                                                                                                             ccAGCAGAGCTGGCGGATGCCATCGGGCTGTCGTATGACACCTATACCGGAAAGCAGATCAGCAGTCAGCGGGCTATGCGAc                                                                      <  1:587290/82‑1 (MQ=255)
                                                                                                                                             ccAGCAGAGCTGGCGGATGCCATCGGGCTGTCGTATGACACCTATACCGGAAAGCAGATCAGCAGTCAGCGGGCTATGCGAc                                                                      >  2:587290/1‑82 (MQ=255)
                                                                                                                                              cagcagAGCTGGCGGATGCCATCGGGCTGTCGTATGACACCTATACCGGAAAGCAGATCAGCAGTCAGCGGGCTATGCGACTGACGGCGGTTTTTTCCTGCGTCAGAGTGCTGGCAGAGTCGGTCGGGATGTTGCCCTGCAATCTGTa   <  2:152303/148‑1 (MQ=255)
                                                                                                                                              cagcagAGCTGGCGGATGCCATCGGGCTGTCGTATGACACCTATACCGGAAAGCAGATCAGCAGTCAGCGGGCTATGCGACTGACGGCGGTTTTTTCCTGCGTCAGAGTGCTGGCAGAGTCGGTCGGGATGTTGCCCTGCAATCTGTAt  <  4:412819/149‑1 (MQ=255)
                                                                                                                                              cagcagAGCTGGCGGATGCCATCGGGCTGTCGTATGACACCTATACCGGAAAGCAGATCAGCAGTCAGCGGGCTATGCGACTGACGGCGGTTTTTTCCTGCGTCAGAGTGCTGGCAGAGTCGGTCGGGATGTTGCCCTGCAATCTGTAt  <  4:186167/149‑1 (MQ=255)
                                                                                                                                              cagcagAGCTGGCGGATGCCATCGGGCTGTCGTATGACACCTATACCGGAAAGCAGATCAGCAGTCAGCGGGCTATGCGACTGACGGCGGTTTTTTCCTGCGTCAGAGTGCTGGCAGAGTCGGTCGGGATGTTGCCCTGCAATCTGTAt  <  4:10080/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                              
TGTTGTTGCCGGGGCGCTGTGCCTGTTCTGGTCGTGGCTGGTGGCGCGATATCTCGACCGTACACAGTCGTCTGTCGGCGGAGGTAAATAGTGTTCTTTTCGGGATTATTTCAACGAAAAAGTGACGCACCGGTGACCACGCCAGCAGAGCTGGCGGATGCCATCGGGCTGTCGTATGACACCTATACCGGAAAGCAGATCAGCAGTCAGCGGGCTATGCGACTGACGGCGGTTTTTTCCTGCGTCAGAGTGCTGGCAGAGTCGGTCGGGATGTTGCCCTGCAATCTGTAT  >  NC_000913/1205641‑1205931

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: