Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,205,787 A→T 7.5% A19A (GCA→GCT beeE → e14 prophage; putative protein BeeE

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,205,7870AT 7.5% 232.5 / 5.9 75A19A (GCA→GCTbeeEe14 prophage; putative protein BeeE
Reads supporting (aligned to +/- strand):  ref base A (34/35);  new base T (3/3);  total (37/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGTTGTTGCCGGGGCGCTGTGCCTGTTCTGGTCGTGGCTGGTGGCGCGATATCTCGACCGTACACAGTCGTCTGTCGGCGGAGGTAAATAGTGTTCTTTTCGGGATTATTTCAACGAAAAAGTGACGCACCGGTGACCACGCCAGCAGAGCTGGCGGATGCCATCGGGCTGTCGTATGACACCTATACCGGAAAGCAGATCAGCAGTCAGCGGGCTATGCGACTGACGGCGGTTTTTTCCTGCGTCAGAGTGCTGGCAGAGTCGGTCGGGATGTTGCCCTGCAATCTGTATCACC  >  NC_000913/1205641‑1205935
                                                                                                                                                  |                                                                                                                                                    
tgttgttgCCGGGGCGCTGTGCCTGTTCTGGTCGTGGCTGGTGGCGCGATATCTCGACCGTACACAGTCGTCTGTCGGCGGAGGTAAATAGTGTTCTTTTCGGGATTATTTCAACGAAAAAGTGACGCACCGGTGACCACGCCAGCagag                                                                                                                                                   >  1:756239/1‑150 (MQ=255)
      tgCCGGGGCGCTGTGCCTGTTCTGGTCGTGGCTGGTGGCGCGATATCTCGACCGTACACAGTCGTCTGTCGGCGGAGGTAAATAGTGTTCTTTTCGGGATTATTTCAACGAAAAAGTGACGCACCGGTGACCACGCCAGCAGAGCTGGCg                                                                                                                                             <  2:80952/150‑1 (MQ=255)
      tgCCGGGGCGCTGTGCCTGTTCTGGTCGTGGCTGGTGGCGCGATATCTCGACCGTACACAGTCGTCTGTCGGCGGAGGTAAATAGTGTTCTTTTCGGGATTATTTCAACGAAAAAGTGACGCACCGGTGACCACGCCAGCAGAGCTGGCg                                                                                                                                             <  1:1559366/150‑1 (MQ=255)
        ccGGGGCGCTGTGCCTGTTCTGGTCGTGGCTGGTGGCGCGATATCTCGACCGTACACAGTCGTCTGTCGGCGGAGGTAAATAGTGTTCTTTTCGGGATTATTTCAACGAAAAAGTGACGCACCGGTGACCACGCCAGCAGAGCTGGCGGa                                                                                                                                           <  2:628963/150‑1 (MQ=255)
         cGGGGCGCTGTGCCTGTTCTGGTCGTGGCTGGTGGCGCGATATCTCGACCGTACACAGTCGTCTGTCGGCGGAGGTAAATAGTGTTCTTTTCGGGATTATTTCAACTAAAAAGTGACGCACCTGTGACCACGCCAGCTCTGCTGGCgtgg                                                                                                                                          <  2:60156/150‑4 (MQ=255)
           gggCGCTGTGCCTGTTCTGGTCGTGGCTGGTGGCGCGATATCTCGACCGTACACAGTCGTCTGTCGGCGGAGGTAAATAGTGTTCTTTTCGGGATTATTTCAACGAAAAAGTGACGCACCGGTGACCACGCCAGCTCTGCTGGCgtggtc                                                                                                                                        <  1:300893/150‑6 (MQ=255)
                   tgCCTGTTCTGGTCGTGGCTGGTGGCGCGATATCTCGACCGTACCCAGTCGTCTGTCGGCGGAGGTAAATAGTGTTCTTTTCGGGATTATTTCAACGAAAAAGTGACGCACCGGTGACCACGCCAGCagag                                                                                                                                                   >  1:768444/1‑131 (MQ=255)
                   tgCCTGTTCTGGTCGTGGCTGGTGGCGCGATATCTCGACCGTACCCAGTCGTCTGTCGGCGGAGGTAAATAGTGTTCTTTTCGGGATTATTTCAACGAAAAAGTGACGCACCGGTGACCACGCCAGCagag                                                                                                                                                   <  2:768444/131‑1 (MQ=255)
                   tgCCTGTTCTGGTCGTGGCTGGTGGCGCGATATCTCGACCGTACACAGTCGTCTGTCGGCGGAGGTAAATAGTGTTCTTTTCGGGATTATTTCAACGAAAAAGTGACGCACCGGTGACCACGCCAGCAGAGCTGGCGGATGCCATCGGGc                                                                                                                                <  1:1490138/150‑1 (MQ=255)
                    gCCTGTTCTGGTCGTGGCTGGTGGCGCGATATCTCGACCGTACACAGTCGTCTGTCGGCGGAGGTAAATAGTGTTCTTTTCGGGATTATTTCAACGAAAAAGTGACGCACCGGTGACCACGCCAGCAGAGCTGGCGGATGCCATCGGGCt                                                                                                                               >  1:1016365/1‑150 (MQ=255)
                     ccTGTTCTGGTCGTGGCTGGTGGCGCGATATCTCGACCGTACACAGTCGTCTGTCGGCGGAGGTAAATAGTGTTCTTTTCGGGATTATTTCAACGAAAAAGTGACGCACCGGTGACCACGCCAGCAGAGCTGGCGGATGCCATCGGGCTg                                                                                                                              <  1:1470984/150‑1 (MQ=255)
                          tCTGGTCGTGGCTGGTGGCGCGATATCTCGACCGTACACAGTCGTCTGTCGGCGGAGGTAAATAGTGTTCTTTTCGGGATTATTTCAACGAAAAAGTGACGCACCGGTGACCACGCCAGCAGAGCTGGCGGATGCCATCGGGCTGTCGTa                                                                                                                         >  1:868579/1‑150 (MQ=255)
                          tCTGGTCGTGGCTGGTGGCGCGATATCTCGACCGTACACAGTCGTCTGTCGGCGGAGGTAAATAGTGTTCTTTTCGGGATTATTTCAACGAAAAAGTGACGCACCGGTGACCACGCCAGCAGAGCTGGCGGATGCCATCGGGCTGTCGTa                                                                                                                         <  1:380016/150‑1 (MQ=255)
                          tCTGGGCGTGGCTGGTGGCGCGATATCTCGACCGTACACAGTCGTCTGTCGGCGGAGGTAAATAGTGTTCTTTTCGGGATTATTTCAACGAAAAAGTGACGCACCGGTGACCACGCCAGCAGAGCTGGCGGATGCCATCGGGCTGTCGTa                                                                                                                         <  1:1324692/150‑1 (MQ=255)
                            tGGTCGTGGCTGGTGGCGCGATATCTCGACCGTACACAGTCGTCTGTCGGCGGAGGTAAATAGTGTTCTTTTCGGGATTATTTCAACGAAAAAGTGACGCACCGGTGACCACGCCAGCAGAGCTGGCGGATGCCATCGGGCTGTCGTATg                                                                                                                       <  1:588380/150‑1 (MQ=255)
                            tGGTCGTGGCTGGTGGCGCGATATCTCGACCGTACACAGTCGTCTGTCGGCGGAGGTAAATAGTGTTCTTTTCGGGATTATTTCAACGAAAAAGTGACGCACCGGTGACCACGCCAGCAGAGCTGGCGGATGCCATCGGGCTGTCGTATg                                                                                                                       <  1:353144/150‑1 (MQ=255)
                             ggTCGTGGCTGGTGGCGCGATATCTCGACCGTACACAGTCGTCTGTTGGCGGATGTAAATAGTGTTCTTTTCGGGTTTATTTCAACGAAATTTTGTCGCAGCGGTGACCTCGCGAGCTGTGCTGGCGGATGCCATCGGGCTGTCGTATGa                                                                                                                      <  2:671371/150‑1 (MQ=255)
                              gTCGTGGCTGGTGGCGCGATATCTCGACCGTACACAGTCGTCTGTCGGCGGAGGTAAATAGTGTTCTTTTCGGGATTATTTCAACGAAAAAGTGACGCACCGGTGACCACGCCAGCAGAGCTGGCGGATGCCATCGGGCTGTCGTATGac                                                                                                                     >  1:863945/1‑150 (MQ=255)
                                cGTGGCTGGTGGCGCGATATCTCGACCGTACACAGTCGTCTGTCGGCGGAGGTAAATAGTGTTCTTTTCGGGATTATTTCAACGAAAAAGTGACGCACCGGTGACCACGCCAGCAGAGCTGGCGGATGCCATCGGGCTGTCGTATGacac                                                                                                                   >  2:1061731/1‑150 (MQ=255)
                                 gTGGCTGGTGGCGCGATATCTCGACCGTACACAGTCGTCTGTCGGCGGAGGTAAATAGTGTTCTTTTCGGGATTATTTCAACGAAAAAGTGACGCACCGGTGACCACGCCAGCAGAGCTGGCGGATGCCATCGGGCTGTCGTATGACAcc                                                                                                                  >  2:807413/1‑150 (MQ=255)
                                  tGGCTGGTGGCGCGATATCTCGACCGTACACAGTCGTCTGTCGGCGGAGGTAAATAGTGTTCTTTTCGGGATTATTTCAACGAAAAAGTGACGCACCGGTGACCACGCCAGCAGAGCTGGCGGATGCCATCGGGCTGTCGTATGACACCt                                                                                                                 <  2:1585235/150‑1 (MQ=255)
                                   ggCTGGTGGCGCGATATCTCGACCGTACACAGTCGTCTGTCGGCGGAGGTAAATAGTGTTCTTTTCGGGATTATTTCAACGAAAAAGTGACGCACCGGTGACCACGCCAGCAGAGCTGGCGGATGCCATCGGGCTGTCGTATGACACCta                                                                                                                <  1:763939/150‑1 (MQ=255)
                                    gCTGGTGGCGCGATATCTCGACCGTACACAGTCGTCTGTCGGCGGAGGTAAATAGTGTTCTTTTCGGGATTATTTCAACGAAAAAGTGACGCACCGGTGACCACGCCAGCAGAGCTGGCGGATGCCATCGGGCTGTCGTATGACACCtat                                                                                                               >  1:1567286/1‑150 (MQ=255)
                                       ggtggCGCGATATCTCGACCGTACACAGTCGTCTGTCGGCGGAGGTAAATAGTGTTCTTTTCGGGATTATTTCAACGAAAAAGTGACGCACCGGTGACCACGCCAGCAGAGCTGGCGGATGCCATCGGGCTGTCGTATGACACCTATAcc                                                                                                            <  1:4687/150‑1 (MQ=255)
                                            cgcgATATCTCGACCGTACACAGTCGTCTGTCGGCGGAGGTAAATAGTGTTCTTTTCGGGATTATTTCAACGAAAAAGTGACGCACCGGTGACCACGCCAGCAGAGCTGGCGGATGCCATCGGGCTGTCGTATGACACCTATACCGGaaa                                                                                                       >  1:1003981/1‑150 (MQ=255)
                                            cgcgATATCTCGACCGTACACAGTCGTCTGTCGGCGGAGGTAAATAGTGTTCTTTTCGGGATTATTTCAACGAAAAAGTGACGCACCGGTGACCACGCCAGCAGAGCTGGCGGATGCCATCGGGCTGTCGTATGACACCTATACCGGaaa                                                                                                       >  1:149794/1‑150 (MQ=255)
                                            cgcgATATCTCGACCGTACACAGTCGTCTGTCAGCGGAGGTAAATAGTGTTCTTTTCGGGATTATTTCAACGAAAAAGTGACGCACCGGTGACCACGCCAGCAGAGCTGGCGGATGCCATCGGGCTGTCGTATGACACCTATACCGGaaa                                                                                                       >  2:226660/1‑150 (MQ=255)
                                               gATATCTCGACCGTACACAGTCGTCTGTCGGCGGAGGTAAATAGTGTTCTTTTCGGGATTATTTCAACGAAAAAGTGACGCACCGGTGACCACGCCAGCAGAGCTGGCGGATGCCATCGGGCTGTCGTATGACACCTATACCGGAAAGCa                                                                                                    >  2:665492/1‑150 (MQ=255)
                                                   tctcGACCGTACACAGTCGTCTGGCGGCGGAGGTAAATACTGTTCTTTTCGGGATTATTTCAACGAAAAAGTGACGCACCGGTGACCACGCCAGCAGAGCTGGCGGATGCCATCGGGCTGTCGTATGACACCTATACCGGAAAGCAGATc                                                                                                <  2:169350/150‑1 (MQ=255)
                                                    ctcGACCGTACACAGTCGTCTGTCGGCGGAGGTAAATAGTGTTCTTTTCGGGATTATTTCAACGAAAAAGTGACGCACCGGTGACCACGCCAGCAGAGCTGGCGGATGCCATCGGGCTGTCGTATg                                                                                                                       >  2:1381012/1‑126 (MQ=255)
                                                    ctcGACCGTACACAGTCGTCTGTCGGCGGAGGTAAATAGTGTTCTTTTCGGGATTATTTCAACGAAAAAGTGACGCACCGGTGACCACGCCAGCAGAGCTGGCGGATGCCATCGGGCTGTCGTATg                                                                                                                       <  1:1381012/126‑1 (MQ=255)
                                                      cGACCGTACACAGTCGTCTGTCGGCGGAGGTAAATAGTGTTCTTTTCGGGATTATTTCAACGAAAAAGTGACGCACCGGTGACCACGCCAGCTCTGCTGGCGGATGCCATCGGGCTGTCGTATGACACCTATACCGGAAAGCAGATcagc                                                                                             >  1:892222/1‑150 (MQ=255)
                                                        aCCGTACACAGTCGTCTGTCGGCGGAGGTAAATAGTGTTCTTTTCGGGATTATTTCAACGAAAAAGTGACGCACCGGTGACCACGCCAGCAGAGCTGGCGGATGCCATCGGGCTGTCGTATGACACCTATACCGGAAAGCAGATcagcag                                                                                           <  1:936326/150‑1 (MQ=255)
                                                         ccGTACACAGTCGTCTGTCGGCGGAGGTAAATAGTGTTCTTTTCGGGATTATTTCAACGAAAAAGTGACGCACCGGTGACCACGCCAGCAGAGCTGGCGGATGCCATCGGGCTGTCGTATGACACCTATACCGGAAAGCAGATCAGCAGt                                                                                          >  1:476747/1‑150 (MQ=255)
                                                           gTACACAGTCGTCTGTCGGCGGAGGTAAATAGTGTTCTTTTCGGGATTATTTCAACGAAAAAGTGACGCACCGGTGACCACGCCAGCAGAGCTGGCGGATGCCATCGGGCTGTCGTATGACACCTATACCGGAAAGCAGATCAGCAGTCa                                                                                        >  2:534531/1‑150 (MQ=255)
                                                              cacaGTCGTCTGTCGGCGGAGGTAAATAGTGTTCTTTTCGGGATTATTTCAACGAAAAAGTGACGCACCGGTGACCACGCCAGCAGAGCTGGCGGATGCCATCGGGCTGTCGTATGACACCTATACCGGAAAGCAGATCAGCAGTCAGCg                                                                                     >  1:683699/1‑150 (MQ=255)
                                                               acaGTCGTCTGTCGGCGGAGGTAAATAGTGTTCTTTTCGGGATTATTTCAACGAAAAAGTGACGCACCGGTGACCACGCCAGCAGAGCTGGCGGATGCCATCGGGCTGTCGTATGACACCTATACCGGAAAGCAGATCAGCAGTCAGCgg                                                                                    <  1:433806/150‑1 (MQ=255)
                                                                 aGTCGTCTGTCGGCGGAGGTAAATAGTGTTCTTTTCGGGATTATTTCAACGAAAAAGTGACGCACCGGTGACCACGCCAGCAGAGCTGGCGGATGCCATCGGGCTGTCGTATGACACCTATACCGGAAAGCAGATCAGCAGTCAGCGGGc                                                                                  <  2:1029720/150‑1 (MQ=255)
                                                                 aGTCGTCTGTCGGCGGAGGTAAATAGTGTTCTTTTCGGGATTATTTCAACGAAAAAGTGACGCACCGGTGACCACGCCAGCAGAGCTGGCGGATGCCATCGGGCTGTCGTATGACACCTATACCGGAAAGCAGATCAGCAGTCAGCGGGc                                                                                  <  2:149794/150‑1 (MQ=255)
                                                                    cgtcTGTCGGCGGAGGTAAATAGTGTTCTTTTCGGGATTATTTCAACGAAAAAGTGACGCACCGGTGACCACGCCAGCAGAGCTGGCGGATGCCATCGGGCTGTCGTATGACACCTATACCGGAAAGCAGATCAGCAGTCAGCGGGCTAt                                                                               >  2:27711/1‑150 (MQ=255)
                                                                     gtcTGTCGGCGGATGTAAATAGTGTTCTTTTCGGGATTATTTCAACGAAAAAGTGACGCACCGGTGACCACGCCAGCAGAGCTGGCGGATGCCATCGGGCTGTCGTATGACACCTATACCGGAAAGCAGATCAGCAGGCAGCGGGCTATg                                                                              <  1:745814/150‑1 (MQ=255)
                                                                     gtcTGTCGGCGGAGGTAAATAGTGTTCTTTTCGGGATTATTTCAACGAAAAAGTGACGCACCGGTGACCACGCCAGCAGAGCTGGCGGATGCCATCGGGCTGTCGTATGACACCTATACCGGAAAGCAGATCAGCAGTCAGCGGGCTATg                                                                              >  1:820838/1‑150 (MQ=255)
                                                                       cTGTCGGCGGAGGTAAATAGTGTTCTTTTCGGGATTATTTCAACGAAAAAGTGACGCACCGGTGACCACGCCAGCAGAGCTGGCGGATGCCATCGGGCTGTCGTATGACACCTATACCGGAAAGCAGATCAGCAGTCAGCGGGCTATGCg                                                                            >  1:15407/1‑150 (MQ=255)
                                                                         gTCGGCGGAGGTAAATAGTGTTCTTTTCGGGATTATTTCAACGAAAAAGTGACGCACCGGTGACCACGCCAGCAGAGCTGGCGGATGCCATCGGGCTGTCGTATGACACCTATACCGGAAAGCAGATCAGCAGTCAGCGGGCTATGCGAc                                                                          >  1:182656/1‑150 (MQ=255)
                                                                         gTCGGCGGAGGTAAATAGTGTTCTTTTCGGGATTATTTCAACGAAAAAGTGACGCACCGGTGACCACGCCAGCAGAGCTGGCGGATGCCATCGGGCTGTCGTATGACACCTATACCGGAAAGCAGATCAGCAGTCAGCGGGCTATGCGAc                                                                          <  1:303441/150‑1 (MQ=255)
                                                                            ggcggAGGTAAATAGTGTTATTTTCGGGATTATTTCAACGAAAAAGTGACGCACCGGTGACCACGCCAGCAGAGCTTGCGGATGCCATCGGGCTGTCGTATGACAACTATACCGGAAAGCAGATCAGCACTCAGCGGGCTATGCGACTGa                                                                       >  2:1543258/1‑150 (MQ=255)
                                                                             gcggAGGTAAATAGTGTTCTTTTCGGGATTATTTCAACGAAAAAGTGACGCACCGGTGACCACGCCAGCAGAGCTGGCGGATGCCATCGGGCTGTCGTATGACACCTATACCGGAAAGCAGATCAGCAGTCAGCGGGCTATGCGACTGAc                                                                      >  2:1431288/1‑150 (MQ=255)
                                                                                    tAAATAGTGTTCTTTTCGGGATTATTTCAACGAAAAAGTGACGCACCGGTGACCACGCCAGCAGAGCTGGCGGATGCCATCGGGCTGTCGTATGACACCTATACCGGAAAGCAGATCAGCAGTCAGCGGGCTATGCGACTGAcggc                                                                   >  2:976648/1‑146 (MQ=255)
                                                                                    tAAATAGTGTTCTTTTCGGGATTATTTCAACGAAAAAGTGACGCACCGGTGACCACGCCAGCAGAGCTGGCGGATGCCATCGGGCTGTCGTATGACACCTATACCGGAAAGCAGATCAGCAGTCAGCGGGCTATGCGACTGAcggc                                                                   <  1:976648/146‑1 (MQ=255)
                                                                                        tAGTGTTCTTTTCGGGATTATTTCAACGAAAAAGTGACGCACCGGTGACCACGCCAGCAGAGCTGGCGGATGCCATCGGGCTGTCGTATGACACCTATACCGGAAAGCAGATCAGCAGTCAGCGGGCTATGCGACTGACGGCGGtttttt                                                           <  2:820838/150‑1 (MQ=255)
                                                                                           tgtTCTTTTCGGGATTATTTCAACGAAAAAGTGACGCACCGGTGACCACGCCAGCAGAGCTGGCGGATGCCATCGGGCTGTCGTATGACACCTATACCGGAAAGCAGATCAGCAGTCAGCGGGCTATGCGACTGACGGCGGTTTTTTCCt                                                        <  2:949455/150‑1 (MQ=255)
                                                                                           tgtTCTTTTCGGGATTATTTCAACGAAAAAGTGACGCACCGGTGACCACGCCAGCAGAGCTGGCGGATGCCATCGGGCTGTCGTATGACACCTATACCGGAAAGCAGATCAGCAGTCAGCGGGCTATGCGACTGACGGCGGTTTTTTCCt                                                        >  2:1511572/1‑150 (MQ=255)
                                                                                                 tttCGGGATTATTTCAACGAAAAAGTGACGCACCGGTGACCACGCCAGCAGAGCTGGCGGATGCCATCGGGCTGTCGTATGACACCTATACCGGAAAGCAGATCAGCAGTCAGCGGGCTATGCGACTGACGGCGGTTTTTTCCTGCGTCa                                                  >  2:767816/1‑150 (MQ=255)
                                                                                                  ttCGGGATTATTTCAACGAAAAAGTGACGCACCGGTGACCACGCCAGCAGAGCTGGCGGATGCCATCGGGCTGTCGTATGACACCTATACCGGAAAGCAGATCAGCAGTCAGCGGGCTATGCGACTGACGGCGGTTTTTTCCTGCGTCag                                                 <  2:182656/150‑1 (MQ=255)
                                                                                                           attTCAACGAAAAAGTGACGCACCGGTGACCACGCCAGCAGAGCTGGCGGATGCCATCGGGCTGTCGTATGACACCTATACCGGAAAGCAGATCAGCAGTCAGCGGGCTATGCGACTGACGGCGGTTTTTTCCTGCGTCAGAGTGCTGGc                                        <  1:665492/150‑1 (MQ=255)
                                                                                                             ttCAACGAAAAAGTGACGCACCGGTGACCACGCCAGCAGAGCTGGCGGATGCCATCGGGCTGTCGTATGACACCTATACCGGAAAGCAGATCAGCAGTCAGCGGGCTATGCGACTGACGGCGGTTTTTTCCTGCGTCAGAGTGCTGGCag                                      <  2:1003981/150‑1 (MQ=255)
                                                                                                             ttCAACGAAAAAGTGACGCACCGGTGACCACGCCAGCAGAGCTGGCGGATGCCATCGGGCTGTCGTATGACACCTATACCGGAAAGCAGATCAGCAGTCAGCGGGCTATGCGACTGACGGCGGTTTTTTCCTGCGTCAGAGTGCTGGCag                                      <  2:1567286/150‑1 (MQ=255)
                                                                                                              tCAACGAAAAAGTGACGCACCGGTGACCACGCCAGCAGAGCTGGCGGATGCCATCGGGCTGTCGTATGACACCTATACCGGAAAGCAGATCAGCAGTCAGCGGGCTATGCGACTGACGGCGGTTTTTTCCTGCGTCAGAGTGCTGGCaga                                     <  2:1368776/150‑1 (MQ=255)
                                                                                                               cAACGAAAAAGTGACGCACCGGTGACCACGCCAGCAGAGCTGGCGGATGCCATCGGGCTGTCGTATGACACCTATACCGGAAAGCAGATCAGCAGTCAGCGGGCTATGCGACTGACGGCGGTTTTTTCCTGCGTCAGAGTGCTGGCagag                                    >  2:973064/1‑150 (MQ=255)
                                                                                                               cAACGAAAAAGTGACGCACCGGTGACCACGCCAGCAGAGCTGGCGGATGCCATCGGGCTGTCGTATGACACCTATACCGGAAAGCAGATCAGCAGTCAGCGGGCTATGCGACTGACGGCGGTTTTTTCCTGCGTCAGAGTGCTGGCagag                                    >  2:440096/1‑150 (MQ=255)
                                                                                                                 aCGAAAAAGTGACGCACCGGTGACCACGCCAGCAGAGCTGGCGGATGCCATCGGGCTGTCGTATGACACCTATACCGGAAAGCAGATCAGCAGTCAGCGGGCTATGCGACTGACGGCGGTTTTTTCCTGCGTCAGAGTGCTGGCAGAgtc                                  <  1:1061731/150‑1 (MQ=255)
                                                                                                                   gAAAAAGTGACGCACCGGTGACCACGCCAGCAGAGCTGGCGGATGCCATCGGGCTGTCGTATGACACCTATACCGGAAAGCAGATCAGCAGTCAGCGGGCTATGCGACTGACGGCGGTTTTTTCCTGCGTCAGAGTGCTGGCAGAgtcgg                                >  1:568591/1‑150 (MQ=255)
                                                                                                                       aaGTGACGCACCGGTGACCACGCCAGCAGAGCTGGCGGATGCCATCGGGCTGTCGTATGACACCTATACCGGAAAGCAGATCAGCAGTCAGCGGGCTATGCGACTGACGGCGGTTTTTTCCTGCGTCAGAGTGCTGGCAGAGtcggtcgg                            <  1:584303/150‑1 (MQ=255)
                                                                                                                        aGTGACGCACCGGTGACCACGCCAGCAGAGCTGGCGGATGCCATCGGGCTGTCGTATGACACCTATAGCGGAAAGCAGATCAGCAGTCAGCCGGCTATGCGACTGAAGGCGGTTTTTTCCTGCGTCAGAGTGCTGGCAGAGTCGGTCggg                           >  2:823736/1‑150 (MQ=255)
                                                                                                                        aGTGACGCACCGGTGACCACGCCAGCAGAGCTGGCGGATGCCATCGGGCTGTCGTATGACACCTATACCGGAAAGCAGATCAGCAGTCAGCGGGCTATGCGACTGACGGCGGTTTTTTCCTGCGTCAGAGTGCTGGCAGAGTCGGTCggg                           <  2:418546/150‑1 (MQ=255)
                                                                                                                             cGCACCGGTGACCACGCCAGCAGAGCTGGCGGATGCCATCGGGCTGTCGTATGACACCTATACCGGAAAGCAGATCAGCAGTCAGCGGGCTATGCGACTGACGGCGGTTTTTTCCTGCGTCAGAGTGCTGGCAGAGTCGGTCGGGATGtt                      >  2:866026/1‑150 (MQ=255)
                                                                                                                             cGCACCGGTGACCACGCCAGCAGAGCTGGCGGATGCCATCGGGCTGTCGTATGACACCTATACCGGAAAGCAGATCAGCAGTCAGCGGGCTATGCGACTGACGGCGGTTTTTTCCTGCGTCAGAGTGCTGGCAGAGTCGGTCGGGATGtt                      >  2:860393/1‑150 (MQ=255)
                                                                                                                               cACCGGTGACCACGCCAGCAGAGCTGGCGGATGCCATCGGGCTGTCGTATGACACCTATACCGGAAAGCAGATCAGCAGTCAGCGGGCTATGCGACTGACGGCGGTTTTTTCCTGCGTCAGAGTGCTGGCAGAGTCGGTCGGGATGTTGc                    <  1:1458352/150‑1 (MQ=255)
                                                                                                                                aCCGGTGACCACGCCAGCAGAGCTGGCGGATGCCATCGGGCTGTCGTATGACACCTATACCGGAAAGCAGATCAGCAGTCAGCGGGCTATGTGACTGACGGCGGTTTTTTCCTGCGTCAGAGTGCTGGCAGAGTCGGTCGGGATGTTGcc                   >  2:1133604/1‑150 (MQ=255)
                                                                                                                                 ccGGTGACCACGCCAGCAGAGCTGGCGGATGCCATCGGGCTGTCGTATGACACCTATACCGGAAAGCAGATCAGCAGTCAGCGGGCTATGCGACTGACGGCGGTTTTTTCCTGCGTCAGAGTGCTGGCAGAGTCGGTCGGGATGTTGccc                  <  1:1295102/150‑1 (MQ=255)
                                                                                                                                 ccGGTGACCACGCCAGCAGAGCTGGCGGATGCCATCGGGCTGTCGTATGACACCTATACCGGAAAGCAGATCAGCAGTCAGCGGGCTATGCGACTGACGGCGGTTTTTTCCTGCGTCAGAGTGCTGGCAGAGTCGGTCGGGATGTTGccc                  >  2:436550/1‑150 (MQ=255)
                                                                                                                                 ccGGTGACCACGCCAGCAGAGCTGGCGGATGCCATCGGGCTGTCGTATGACACCTATACCGGAAAGCAGATCAGCAGTCAGCGGGCTACGCGACTGACGGCGGTTTTTTCCTGCGTCAGAGTGCTGGCAGAGTCGGTCGGGATGTTGccc                  <  1:760525/150‑1 (MQ=255)
                                                                                                                                   ggTGACCACGCCAGCAGAGCTGGCGGATGCCATCGGGCTGTCGTATGACACCTATACCGGAAAGCAGATCAGCAGTCAGCGGGCTATGCGACTGACGGCGGTTTTTTCCTGCGTCAGAGTGCTGGCAGAGTCGGTCGGGATGTTGCCCTg                >  2:1480101/1‑150 (MQ=255)
                                                                                                                                   ggTGACCACGCCAGCAGAGCTGGCGGATGCCATCGGGCTGTCGTATGACACCTATACCGGAAAGCAGATCAGCAGTCAGCGGGCTATGCGACTGACAGCGGTTTTTTCCTGCGTCAGAGTGCTGGCAGAGTCGGTCGGGATGTTGCCCTg                >  1:1193661/1‑150 (MQ=255)
                                                                                                                                      gcatccGCCAGCTCTGCTGGCGGATGCCATCGGGCTGTCGTATGACACCTATACCGGAAAGCAGCTCAGCAGTCAGCGGGCTATGCGACTGACGGCGGTTTTTTCCTGCGTCAGAGTGCTGGCAGAGTCGGTCGGGATGTTGCCCTGCaa             >  1:751898/6‑150 (MQ=255)
                                                                                                                                      gcatccGCCAGCTCTGCTGGCGGATGCCATCGGGCTGTCGTATGACACCTATACCGGAAAGCAGATCAGCAGTCAGCGGGCTATGCGACTGACGGCGGTTTTTTCCTGCGTCAGAGTGCTGGCTGAGTCGGTCGGGATGTTGCCCTGCaa             >  1:1075543/6‑150 (MQ=255)
                                                                                                                                      gACCACGCCAGCAGAGCTGGCGGATGCCATCGGGCTGTCGTATGACACCTATACCGGAAAGCAGATCAGCAGTCAGCGGGCTATGCGACTGACGGCGGTTTTTTCCTGCGTCAGAGTGCTGGCAGAGTCGGTCGGGATGTTGCCCTGCaa             <  2:863945/150‑1 (MQ=255)
                                                                                                                                        ccACGCCAGCAGAGCTGGCGGATGCCATCGGGCTGTCGTATGACACCTATACCGGAAAGCAGATCAGCAGTCAGCGGGCTATGCGACTGACGGCGGTTTTTTCCTGCGTCAGAGTGCTGGCAGAGTCGGTCGGGATGTTGCCCTGCAATc           >  2:449310/1‑150 (MQ=255)
                                                                                                                                              cagcagAGCTGGCGGATGCCATCGGGCTGTCGTATGACACCTATACCGGAAAGCAGATCAGCAGTCAGCGGGCTATGCGACTGACGGCGGTTTTTTCCTGCGTCAGAGTGCTGGCAGAGTCGGTCGGGATGTTGCCCTGCAATCTGTATc     <  1:408478/150‑1 (MQ=255)
                                                                                                                                               agcagAGCTGGCGGATGCCATCGGGCTGTCGTATGACACCTATACCGGAAAGCAGATCAGCAGTCAGCGGGCTATGCGACTGACGGCGGTTTTTTCCTGCGTCAGAGTGCTGGCAGAGTCGGTCGGGATGTTGCCCTGCAATCTGTATcc    >  1:1436269/1‑149 (MQ=255)
                                                                                                                                                 cagAGCTGGCGGATGCCATCGGGCTGTCGTATGACACCTATACCGGAAAGCAGATCAGCAGTCAGCGGGCTATGCGACTGACGGCGGTTTTTTCCTGCGTCAGAGTGCTGGCAGAGTCGGTCGGGATGTTGCCCTGCAATCTGTATCAcc  >  1:27486/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
TGTTGTTGCCGGGGCGCTGTGCCTGTTCTGGTCGTGGCTGGTGGCGCGATATCTCGACCGTACACAGTCGTCTGTCGGCGGAGGTAAATAGTGTTCTTTTCGGGATTATTTCAACGAAAAAGTGACGCACCGGTGACCACGCCAGCAGAGCTGGCGGATGCCATCGGGCTGTCGTATGACACCTATACCGGAAAGCAGATCAGCAGTCAGCGGGCTATGCGACTGACGGCGGTTTTTTCCTGCGTCAGAGTGCTGGCAGAGTCGGTCGGGATGTTGCCCTGCAATCTGTATCACC  >  NC_000913/1205641‑1205935

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: