Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,205,173 A→G 5.5% pseudogene (770/1134 nt) aaaE → pseudogene, phage terminase protein A 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,1730AG 5.5% 202.2 / 3.3 74pseudogene (770/1134 nt)aaaEpseudogene, phage terminase protein A family, e14 prophage;Phage or Prophage Related
Reads supporting (aligned to +/- strand):  ref base A (36/33);  new base G (2/2);  total (39/35)
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.14e-01

AGAAAAATGAAGATCGCCGGACAGCCGAACGCTTTCAGAAATGGGTGGAAATGGGCGTTCTGACCGTTACCGATGGTGCGGAGGTGGATTATCGCTACATCCTCGAAGAGGCCAAAGCGGCGAACAAAATCAGCCCGGTCAGTGAGTCACCCATCGACCCCTTCGGGGCGACCGGGCTGTCACATGACCTTGCTGATGAAGACCTGAACCCCGTCACCATCATTCAGAACTACACCAACATGTCCGATCCGATGAAAGAGCTGGAAGCGGCGATTGAATCGGGGCGCTTTCA  >  NC_000913/1205029‑1205320
                                                                                                                                                |                                                                                                                                                    
agaAAAATGAAGATCGCCGGACAGCCGAACGCTTTCAGAAATGGGTGGAAATGGGCGTTCTGACCGTTACCGATGGTGCGGAGGTGGATTATCGCTACATCCTCGAAGAGGCCAAAGCGGCGAACAAAATCAGCCCGGTCAGTGAGTCa                                                                                                                                                  >  2:474709/1‑149 (MQ=255)
 gaAAAATGAAGATCGCCGGACAGCCGAACGCTTTCAGAAATGGGTGGAAATGGGCGTTCTGACCGTTACCGATGGTGCGGAGGTGGATTATCGCTACATCCTCGAAGAGGCCAAAGCGGCGAACAAAATCAGCCCGGTCAGTGAGTCAc                                                                                                                                                 <  2:419906/149‑1 (MQ=255)
 gaAAAATGAAGATCGCCGGACAGCCGAACGCTTTCAGAAATGGGTGGAAATGGGCGTTCTGACCGTTACCGATGGTGCGGAGGTGGATTATCGCTACATCCTCGAAGAGGCCAAAGCGGCGAACAAAATCAGCCCGGTCAGTGAGTCAc                                                                                                                                                 <  1:930996/149‑1 (MQ=255)
 gaAAAAGGAAGATCGCCGGACAGCCGAACGCTTTCAGAAATGGGTGGAAATGGGCGTTCTGACCGTTACCGATGGTGCGGAGGTGGATTATCGCTACATCCTCGAAGAGGCCAAAGCGGCGAACAAAATCAGCCCGGTCAGTGAGTCAc                                                                                                                                                 <  2:827981/149‑1 (MQ=255)
     aaTGAAGATCGCCGGACAGCCGAACGCTTTCAGAAATGGGTGGGAGTGGGCGTTCTGAGCGTTACCGATGGTGCGGAGGTGGATTATCGGTACATCCTCGAAGAGGCCACAGCGGCGAACAATATCAGCCCGGTCCGTGAGTCAcccct                                                                                                                                             >  2:760364/1‑147 (MQ=255)
        gatgATCGCCGGACAGCCGAACGCTTTCAGAAATGGGTGGAAATGGGCGTTCTGACCGTTACCGATGGTGCGGAGGTGGATTATCGCTACATCCTCGAAGAGGCCAAAGCGGCGAACAAAATCAGCCCGGTCAGTGAGTCACCCATCGa                                                                                                                                          <  1:757203/146‑1 (MQ=255)
        gAAGATCGCCGGACAGCCGAACGCTTTCAGAAATGGGTGGAAATGGGCGTTCTGACCGTTACCGATGGTGCGGAGGTGGATTATCGCTACATCCTCGAAGAGGCCAAAGCGGCGAACAAAATCAGCCCGGTCAGTGAGTCACCCATCGa                                                                                                                                          >  1:40464/1‑149 (MQ=255)
         aaGATCGCCGGACAGCCGAACGCTTTCAGAAATGGGTGGAAATGGGCGTTCTGACCGTTACCGATGGTGCGGAGGTGGATTATCGCTACATCCTCGAAGAGGCCAAAGCGGCGAACAAAATCAGCCCGGTCAGTGAGTCACCCATCGAc                                                                                                                                         <  2:452487/149‑1 (MQ=255)
          aGATCGCCGGACAGCCGAACGCTTTCAGAAATGGGTGGAAATGGGCGTTCTGACCGTTACCGATGGTGCGGAGGTGGATTATCGCTACATCCTCGAAGAGGCCAAAGCGGCGAACAAAATCAGCCGGGTCAGTGAGTCACCCATCGAcc                                                                                                                                        <  1:806971/149‑1 (MQ=255)
            agcGCCGGACAGCCGAACGCTTTCAGAAATGGGTGGAAATGGGCGTTCTGACCGTTACCGATGGTGCGGAGGTGGATTATCGCTACATCCTCGAAGAGGCCAAAGCGGCGAACAAAATCAGCCCGGTCAGTGAGTCACCCATCGAcccc                                                                                                                                      <  1:261687/147‑1 (MQ=255)
            aTCGCCGGACAGCCGAACGCTTTCAGAAATGGGTGGAAATGGGCGTTCTGACCGTTACCGATGGTGCGGAGGTGGATTATCGCTACATCCTCGAAGAGGCCAAAGCGGCGGACAAAATCAGCCCGGGCAGTGAGTCACCCGTCGAcccc                                                                                                                                      >  2:842953/1‑149 (MQ=255)
            aTCGCCGGACAGCCGAACGCTTTCAGAAATGGGTGGAAATGGGCGTTCTGACCGTTACCGATGGTGCGGAGGTGGATTATCGCTACATCCTCGAAGAGGCCAAAGCGGCGAACAAAATCAGCCCGGTCAGTGAGTCACCCATCGAcccc                                                                                                                                      >  2:842995/1‑149 (MQ=255)
            aTCGCCGGACAGCCGAACGCTTTCAGAAATGGGTGGAAATGGGCGTTCTGACCGTTACCGATGGTGCGGAGGTGGATTATCGCTACATCCTCGAAGAGGCCAAAGCGGCGAACAAAATCAGCCCGGTCAGTGAGTCACCCATCGAcccc                                                                                                                                      >  1:1138918/1‑149 (MQ=255)
              cGCCGGACAGCCGAACGCTTTCAGAAATGGGTGGAAATGGGCGTTCTGACCGTTACAGATGGTGCGGAGGTGGATTATCGCTACATCCTCGAAGAGGCCAAAGCGGCGAACAAAATCAGCCCGGACAGTGAGTCACCCATCGACCCCtt                                                                                                                                    >  2:1294743/1‑149 (MQ=255)
               gCCGGACAGCCGAACGCGTTCAGAAATGGGTGGAAATGGGCGTTCTGACCGTTACCGATGGTGCGGAGGTGGATTATCGCTACATCCTCGAAGAGGCCAAAGCGGCGAACAAAATCAGCCCGGTCAGTGAGTCACCCATCGACCCCTTc                                                                                                                                   <  1:608906/149‑1 (MQ=255)
                 cGGACAGCCGAACGCCTTCAGAAATGGGTGGAAATGGGCGTTCTGACCGTTACCGATGGTGCGGAGGTGGATTATCGCTACATCCTCGAAGAGGCCAAAGCGGCGAACAAAATCAGCCCGGTCAGTGAGTCACCCATCGACCCCTTCgg                                                                                                                                 <  1:145734/149‑1 (MQ=255)
                 cGGACAGCCGAACGCCTTCAGAAATGGGTGGAAATGGGCGTTCTGACCGTTACCGATGGTGCGGAGGTGGATTATCGCTACATCCTCGAAGAGGCCAAAGCGGCGAACAAAATCAGCCCGGTCAGTGAGTCACCCATCGACCCCTTCgg                                                                                                                                 <  2:26838/149‑1 (MQ=255)
                  ggACAGCCGAACGTTTTCAGAAATGGGTGGAAAGGGGCGTTCTGACCGTTACCGATGGTGCGGAGGTGGATTATCGCTACATCCTCGAAGAGGCCAAAGCGGCGAACAAAATCAGCCCGGTCAGTGAGTCACCCATCGACCCCTTCggg                                                                                                                                <  1:228104/149‑1 (MQ=255)
                  ggACAGCCGAACGCTTTCAGAAATGGGTGTAAATGGGCGTTCTGACCGTTACCGATGGTGCGGAGGTGGATTATCGCTACATCCTCGAAGAGGCCAAAGCGGCGAACAAAATCAGCCCGGTCAGTGAGTCACCCATCGACCCCTTCggg                                                                                                                                <  1:470459/149‑1 (MQ=255)
                  ggACAGCCGAACGCTTGCAGAAATGGGTGTAAATGGGCGTTCTGACCGTTACCGATGGTGCGGAGGTGGATTATCGCTACATCCTCGAAGAGGCCAAAGCGGCGAACAAAATCAGCCCGGTCAGTGAGTCACCCATCGACCCCTTCggg                                                                                                                                <  2:1198594/149‑1 (MQ=255)
                   gACAGCCGAACGCTTTCAGAAATGGGTGGAAATGGGCGTTCTGACCGTTACCGATGGTGCGGAGGTGGATTATCGCTACATCCTCGAAGAGGCCAAAGCGGCGAACAAAATCAGCCCGGTCAGTGGGTCACCCATCGACCCCTTCgggg                                                                                                                               >  2:1304436/1‑149 (MQ=255)
                   gACAGCCGAACGCTTTCAGAAATGGGAGGAAACGAGCGTACTCACCGTTACCGATGGTGCGGAGGTGGCTTATCGCTACATCCTCGAAGAGGCCAAAGCGGGGCACAAAATCAGCCCCGCCAGGCCGTACCCCACCGACCCCTTCcggg                                                                                                                               >  2:249264/1‑149 (MQ=255)
                   gACAGACGAACGCTTTCAGAAATGGGTGGAAATGGGCGTTCTGACCGTTACCGATGGTGCGGAGGTGGATTATCGCTACATCCTCGAAGAGGCCAAAGCGGCGAACAAAATCAGCCCGGTCAGTGGGTCACCCATCGACCCCTTCgggg                                                                                                                               >  2:1020517/1‑149 (MQ=255)
                      aGCCGAACGCTTTCAGAAATGGGTGGAAATGGGCGTTCTGACCGTTACCGATGGTGCGGAGGTGGATTATCGCTACATCCTCGAAGAGGCCAAAGCGGCGAACAAAATCAGCCCGGTCAGTGAGTCACCCATCGACCCCTTCGGGGCGa                                                                                                                            >  1:1039522/1‑149 (MQ=255)
                        ccGAACGCTTTCAGAAATGGGTGGAAATGGGCGTTCTGACCGTTACCGATGGTGCGGAGGTGGATTATCGCTACATCCTCGAAGAGGCCAAAGCGGCGAACAAAATCAGCCCGGTCAGTGAGTCACCCATCGACCCCTTCGGGGCgctg                                                                                                                          >  1:611612/1‑146 (MQ=255)
                        ccGAACGCTTTCAGAAATGGGTGGAAATGGGCGTTCTGACCGTTACCGATGGTGCGGAGGTGGATTATCGCTACATCCTCGAAGAGGCCAAAGCGGCGAACAAAATCAGCCCGGTCAGTGAGTCACCCATCGACCCCTTCGGGGCGAcc                                                                                                                          >  2:933826/1‑149 (MQ=255)
                               cTTTCAGAAATGGGTGGAAATGGGCGTTCTGACCGTTACCGATGGTGCGGAGGTGGATTATCGCTACATCCTCGAAGAGGCCAAAGCGGCGAACAAAATCAGCCCGGTCAGTGAGTCACCCATCGACCCCTTCGGGGCGACCGGGCTGt                                                                                                                   <  2:1164613/149‑1 (MQ=255)
                                tttCAGAAATGGGTGGAAATGGGCGTTCTGACCGTTACCGATGGTGCGGAGGTGGATTATCGCTACATCCTCGAAGAGGCCAAAGCGGCGAACAAAATCAGCCCGGTCAGTGAGTCACCCATCGACCCCTTCGGGGCGACCGGGCTGTc                                                                                                                  >  1:1098655/1‑149 (MQ=255)
                                 ttCAGAAATGGGTGGAAATGGGCGTTCTGACCGTTACCGATGGTGCGGAGGTGGATTATCGCTACATCCTCGAAGAGGCCAAAGCGGCGAACAAAATCAACCCGGTCAGTGAGTCACCCATCGACCCCTTCGGGGCGACCGGGCTGTca                                                                                                                 >  2:1093497/1‑149 (MQ=255)
                                    aGAAATGGGTGGAAATGGGCGTTCTGACCGTTACCGATGGTGCGGAGGTGGATTATCGCTACATCCTCGAAGAGGCCAAAGCGGCGAACAAAATCAGCCCGGTCAGTGAGTCACCCATCGACCCCTTCGGGGCGACCGGGCTGTCACAt                                                                                                              >  1:1262523/1‑149 (MQ=255)
                                      aaaTGGGTGGAAATGGGCGTTCTGACCGTTACCGATGGTGCGGAGGTGGATTATCGCTACATCCTCGAAGAGGCCAAAGCGGCGAACAAAATCAGCCCGGTCAGTGAGTCACCCATCGACCCCTTCGGGGCGACCGGGCTGTCACATGa                                                                                                            <  1:1096795/149‑1 (MQ=255)
                                      aaaTGGGTGGAAATGGGCGTTCTGACCGTTACCGATGGTGCGGAGGTGGATTATCGCTACATCCTCGAAGAGGCCAAAGCGGCGAACAAAATCAGCCCGGTCAGTGAGTCACCCATCGACCCCTTCGGGGCGACCGGGCTGTCACATGa                                                                                                            >  2:139745/1‑149 (MQ=255)
                                       aaggggTGGAAATGGGCGTTCCGACCGGTACCGATGGTGCGGAGGTGGATTATCGCTACATCCTCGAAGAGGCCAAAGCGGCGAACAAAATCAGCCCGGTCAGTGAGTCACCCATCGACCCCTTCGGGGCGACCGGGCTGTCACATGAc                                                                                                           <  1:386459/146‑1 (MQ=255)
                                        aTGGGTGGAAATGGGCGTTCTGACCGTTACCGATGGTGCGGAGGTGGATTATCGCTACATCCTCGAAGAGGCCAAAGCGGCGAACAAAATCAGCCCGGTCAGTGAGTCACCCATCGACCCCTTCGGGGCGACCGGGCTGTCACATGAcc                                                                                                          >  1:1223348/1‑149 (MQ=255)
                                         tGGGTGGAAATGGGCGTTCTGACCGTTACCGATGGTGCGGAGGTGGATTATCGCTACATCCTCGAAGAGGCCAAAGCGGCGAACAAAATCAGCCCGGTCAGTGAGTCACCCATCGACCCCTTCGGGGCGACCGGGCTGTCACATGACCt                                                                                                         <  1:1201554/149‑1 (MQ=255)
                                           ggTGGAAATGGGCGTTCTGACCGTTACCGATGGTGCGGAGGTGGATTATCGCTACATCCTCGAAGAGGCCAAAGCGGCGAACAAAATCAGCCCGGTCAGTGAGTCACCCATCGACCCCTTCGGGGCGACCGGGCTGTCACATGACCTTg                                                                                                       <  1:707814/149‑1 (MQ=255)
                                             tGGAAATGGGCGTTCTGACCGTTACCGATGGTGCGGAGGTGGATTATCGCTACATCCTCGAAGAGGCCAAAGCGGCGAACAAAATCAGCCCGGTCAGTGAGTCACCCATCGACCCCTTCGGGGCGACCGGGCTGGCACCTGACCTTGCt                                                                                                     >  1:267404/1‑149 (MQ=255)
                                              ggAAATGGGCGTTCTGACCGTTACCGATGGTGCGGAGGTGGATTATCGCTACATCCTCGAAGAGGCCAAAGCGGCGAACAAAATCAGCCCGGTCAGTGAGTCACCCATCGACCCCTTCGGGGCGACCGGGCTGTCACATGACCTTGCtg                                                                                                    <  2:82651/149‑1 (MQ=255)
                                                    gggCGTTCTGACCGTTACCGATGGTGCGGAGGTGGATTATCGCTACATCCTCGAAGAGGCCAAAGCGGCGAACAAAATCAGCCCGGTCAGTGAGTCACCCATCGACCCCTTCGGGGCGACCGGGCTGTCACATGACCTTGCTGATGAAg                                                                                              <  2:1162158/149‑1 (MQ=255)
                                                      gCGTTCTGACCGTTACCGATGGTGCGGAGGTGGATTATCGCTACATCCTCGAAGAGGCCAAAGCGGCGAACAAAATCAGCCCGGTCAGTGAGTCACCCATCGACCCCTTCGGGGCGACCGGGCTGTCACATGACCTTGCTGATGAAGAc                                                                                            >  1:473350/1‑149 (MQ=255)
                                                        gTTCTGACCGTTACCGATGGTGCGGAGGTGGATTATCGCTACATCCTCGAAGAGGCCAAAGCGGCGAACAAAATCAGCCCGGTCAGTGAGTCACCCATCGACCCCTTCGGGGCGACCGGGCTGGCACATGACCTTGCTGATGAAGACCt                                                                                          >  2:1195435/1‑149 (MQ=255)
                                                        gTTCTAACCGTTACCGATGGTGCGGAGGTGGATTATCGCTACATCCTCGAAGAGGCCAAAGCGGCGAACAAAATCAGCCCGGTCAGTGAGTCACCCATCGACCCCTTCGGGGCGACCGGGCTGTCACATGACCTTGCTGCTGAAGAGCt                                                                                          >  2:1020570/1‑149 (MQ=255)
                                                         ttCTGACCGTTACCGATGGTGCGGAGGTGGATTATCGCTACATCCTCGAAGAGGCCAAAGCGGCGAACAAAATCAGCCCGGTCAGTGAGTCACCCATCGACCCCTTCGGGGCGACCGGGCTGTCACATGACCTTGCTGATGAAGACCTg                                                                                         >  2:36330/1‑149 (MQ=255)
                                                                      cGATGGTGCGGAGGTGGATTATCGCTACATCCTCGAAGAGGCCAAAGCGGCGAACAAAATCAGCCCGGTCAGTGAGTCACCCATCGACCCCTTCGGGGCGACCGGGCTGTCACATGACCTTGCTGATGAAGACCTGAACCCCGTCACca                                                                            >  2:663400/1‑149 (MQ=255)
                                                                            tGCGGAGGTGGATTATCGCTACATCCTCGAAGAGGCCAAAGCGGCGAACAAAATCAGCCCGGTCAGTGAGTCACCCATCGACCCCTTCGGGGCGACCGGGCTGTCACATGACCTTGCTGATGAAGACCTGAACCCCGTCACCATCATTc                                                                      <  2:40464/149‑1 (MQ=255)
                                                                            tGCGGAGGTGGATTATCGCTACATCCTCGAAGAGGCAAAAGCGGCGAACAAAATCAGCCCGGTCAGTGAGTCAACCATCGAGTCCTTCGGGGCGAGCGGGCTGTCAAATGACCTTGCTGATGAAGACCTGAACCCCGTCACCATCATTc                                                                      >  1:1134754/1‑149 (MQ=255)
                                                                              cGGAGGTGGATTATCGCTACATCCTCGAAGAGGCCAAAGCGGCGAACAAAATCAGCCCGGTCAGTGAGTCACCCATCGACCCCTTCGGGGCGACCGGGCTGTCACATGACCTTGCTGATGAAGACCTGAACCCCGTCACCATCATTCAg                                                                    <  1:196602/149‑1 (MQ=255)
                                                                                gAGGTGGATTATCGCTACATCCTCGAAGAGGCCAAAGCGGCGAACAAAATCAGCCCGGTCAGTGAGTCACCCATCGACCCCTTCGGGGCGACCGGGCTGTCACATGACCTTGCTGCTGAAGACCTGAACCCCGTCACCATCATTCAGaa                                                                  <  1:1064789/149‑1 (MQ=255)
                                                                                gAGGTGGATTATCGCTACATCCTCGAAGAGGCCAAAGCGGCGAACAAAATCAGCCCGGTCAGTGAGTCACCCATCGACCCCTTCGGGGCGACCGGGCTGTCACATGACCTTGCTGATGAAGACCTGAACCCCGTCACCATCATTCAGaa                                                                  >  2:399898/1‑149 (MQ=255)
                                                                                gAGGTGGATTATCGCTACATCCTCGAAGAGGCCAAAGCGGCGAACAAAATCAGCCCGGTCAGTGAGTCACCCATCGACCCCTTCGGGGCGACCGGGCTGTCACATGACCTTGCTGATGAAGACCTGAACCCCGTCACCATCATTCAGaa                                                                  <  1:474709/149‑1 (MQ=255)
                                                                                    tGGATTATCGCTACATCCTCGAAGAGGCCAAAGCGGCGAACAAAATCAGCCCGGTCAGTGAGTCACCCATCGGCCCCTTCGGGGCGACCGGGCTGTCACATGACGTTGCTGATGAAGACCTGAACCCCGTCACCAGCATTCAGAACtta                                                              >  1:78207/1‑147 (MQ=255)
                                                                                            cGCTACATCCTCGAAGAGGCCAAAGCGGCGAACAAAATCAGCCCGGTCAGTGGGTCACCCATCGACCCCTTCGGGGCGACCGGGCTGTCACATGACCTTGCTGATGAAGACCTGAACCCCGTCACCATCATTCAGAACTACACCAACAt                                                      <  1:1304436/149‑1 (MQ=255)
                                                                                            cGCTACATCCTCGAAGAGGCCAAAGCGGCGAACAAAATCAGCCCGGTCAGTGGGTCACCCATCGACCCCTTCGGGGCGACCGGGCTGTCACATGACCTTGCTGATGAAGACCTGAACCCCGTCACCATCATTCAGAACTACACCAACAt                                                      <  1:1020517/149‑1 (MQ=255)
                                                                                             gCTACATCCTCGAAGAGGCCAAAGCGGCGAACAAAATCAGCCCGGTCAGTGAGTCACCCATCGACCCCTTCGGGGCGACCGGGCTGTCACATGACCTTGCTGATGAAGACCCGAACCCCGTCACCATCATTCAGAACTACACCAACATg                                                     >  1:438616/1‑149 (MQ=255)
                                                                                              cTACATCCTCGAAGAGGCCAAAGCGGCGAACAAAATCAGCCCGGTCAGTGAGTCACCCATCGACCCCTTCGGGGCGACCGGGCTGTCACATGACCTTGCTGATGAAGACCTGAACCCCGTCACCATCATTCAGAACTACACCACCATGt                                                   >  1:982700/1‑149 (MQ=255)
                                                                                               tACATCCTCGAAGACGCCAAAGCGGCGAACAAAATCAGCCCGGTCAGTGAGTCACCCATCGACCCCTTCGGGGCGACCGGGCTGTCACATGACCTTGCTGATGAAGACCTGAACCCCGTCACCATCATTCAGAACTACACCAACATGtc                                                  <  1:1225427/149‑1 (MQ=255)
                                                                                               cacacCCTCCAAGAGGCCAAAGTGGCGAACAAAATCAGCCCGCTCAGTGAGTCACCCATCGACCCCTTCGGGCCGACCGGGCTGTCACATGACCTTGCTGATGAAGACCTGAACCCCGTCACCATCATTCAGAACTACACCAACATGtc                                                  <  1:555496/148‑1 (MQ=255)
                                                                                                 cATCCTCGAAGAGGCCAAAGCGGCGAACAAAATCAGCCCGGTCAGTGAGTCACCCATCGACCCCTTCGGGGCGACCGGGCTGTCACATGACCTTGCTGATGAAGACCTGAACCCCGTCACCATCATTCAGAACTACACCAACATGtccg                                                >  1:1453135/1‑149 (MQ=255)
                                                                                                    ccTCGAAGAGGCCAAAGCGGCGAACAAAATCAGCCCGGACAGTGAGGCACCCATCGACCCCTTCGGGGCGACCGGGCTGTCACATGACCTTGCTGATGAAGACCTGAACCCCGTCCCCATCATTCAGAACTACACCAACCTGtccgatc                                             >  1:350137/1‑149 (MQ=255)
                                                                                                      tCGAAGAGGCCAAAGCGGCGAACAAAATCAGCCCGGTCAGTGAGTCACCCATCGACCCCTTCGGGGCGACCGGGCTGTCACATGACCTTGCTGATGAAGACCTGAACCCCGTCACCATCATTCAGAACTACACCAACATGtccgatccg                                           >  1:677164/1‑149 (MQ=255)
                                                                                                          agagGCCCAAGCGGCGAACAAAATCAGCCCGGCCAGTGAGTCACCCATCGACCCCTTCGGGGCGACCGGGCTGTCACATGACCTTGCTGATGAAGACCTGAACCCCGTCACCATCATTCAGAACTACACCAACATGTCCGATCCGATGa                                       <  2:240460/149‑1 (MQ=255)
                                                                                                              gCCAAAGCGGCGAACAAAATCAGCCCGGTCAGTGAGTCACCCATCGACCCCTTCGGGGCGACCGGGCTGTCACATGACCTTGCTGATGAAGACCTGAACCCCGTCACCATCATTCAGAACTACACCAACATGTCCGATCCGATGAAaga                                   >  1:107856/1‑149 (MQ=255)
                                                                                                                 aaaGCGGCGAACAAAATCAGCCCGGTCAGTGAGTCACCCATCGACCCCTTCGGGGCGACCGGGCTGTCACATGACCTTGCTGATGAAGACCTGAACCCCGTCACCATCATTCAGAACTACACCAACATGTCCGATCCGATGAAAGAGCt                                <  1:1240864/149‑1 (MQ=255)
                                                                                                                  aaGCGGCGAACAAAATCGGCCCGGTCAGTGAGTCACCCATCGACCCCTTCGGGGCGACCGGGCTGTCACATGACCTTGCTGATGAAGACCTGAACCCCGTCACCATCATTCAGAACTACACCAACATGTCCGATCCGATGAAAGAGCTg                               <  2:1449924/149‑1 (MQ=255)
                                                                                                                    gcggcgAACAAAGTCAGCCCGGTCAGTGAGTCACCCATCGACCCCTTCGGGGCGACCGGGCTGTCACATGACCTTGCTGATGAAGACCTGAACCCCGTCACCATCATTCAGAACTACACCAACATGTCCGATCCGATGAAAGAGCTGGa                             <  2:1051748/149‑1 (MQ=255)
                                                                                                                          aaCAAAATCAGCCCGGTCAGTGAGTCACCCATCGATCCCTTCGGGGCGACCGGGCTGTCACATGACCTTGCTGATGAAGACCTGAACCCCGTCACCATCATTCAGAACTACACCAACATGTCCGATCCGCTGAAAGAGCTGGAAgcggc                       >  2:905481/1‑149 (MQ=255)
                                                                                                                          aaCAAAATCAGCCCGGTCAGTGAGTCACCCATCGACCCCTTCGGGGCGACCGGGCTGTCACATGACCTTGCTGATGAAGACCTGAACCCCGTCACCATCATTCAGAACTACACCAACATGTCCGATCCGATGAAAGAGCTGGAAgcggc                       <  1:1341002/149‑1 (MQ=255)
                                                                                                                          aaCAAAATCAGCCCGGTCAGTAAGTCACCCAGCGACCCCTTCGGGGCGACCGGGCTGTCACATGACCTTGCTGATGAAGACCTGAACCCCGTCACCATCATTCAGAACTACACCAACATGTCCGATCCGATGAAAGAGCTGGGAgcggg                       >  2:1103157/1‑148 (MQ=255)
                                                                                                                              aaaTCAGCCCGGTCAGTGAGTCACCCATCGACCCCTTCGGGGCGACCGGGCTGTCACATGACCTTGCTGATGAAGACCTGAACCCCGTCACCATCATTCAGAACTACACCAACATGTCCGATCCGATGAAAGAGCTGGAAGCGGCGAtt                   <  1:663400/149‑1 (MQ=255)
                                                                                                                              aaaTCAGCCCGGTCAGTGAGTCACCCATCGACCCCTTCGGGGCGACCGGGCTGTCACATGACCTTGCTGATGAAGACCTGAACCCCGTCACCATCATTCAGAACTAACAATTTTGTTCCGATCCGATGAAAGAGCTGGAAGCGGCGAtt                   <  1:710943/149‑1 (MQ=255)
                                                                                                                                     cccGGTCAGTGAGTCACCCATCGACCCCTTCGGGGCGACCGGGCTGTCACATGACCTTGCTGATGAAGACCTGAACCCCGTCACCATCATTCAGAACTACACCAACATGTCCGATCCGATGAAAGAGCTGGAAGCGGCGATTGAATCgg            >  2:617353/1‑149 (MQ=255)
                                                                                                                                            aGTGAGTCACCCATCGACCCCTTCGGGGCGACCGGGCTGTCACATGACCTTGCTGATGAAGACCTGAACCCCGTCACCATCATTCAGAACTACACCAACATGTCCGATCCGATGAAAGAGCTGGACGCGGCGGTTGAATCGGGGCGCtt     >  1:281682/1‑149 (MQ=255)
                                                                                                                                             gTGAGTCACCCATCGACCCCTTCGGGGCGACCGGGCTGTCACATGACCTTGCTGATGAAGACCTGAACCCCGTCACCATCATTCAGAACTACACCAACATGTCCGATCCGATGAAAGAGCTGGAAGCGGCGATTGAGTCGGGGCGCttt    >  1:622640/1‑149 (MQ=255)
                                                                                                                                               gAGTCACCCATCGACCCCTTCGGGGCGACCGGGCTGTCACATGACCTTGCTGATGAAGACCTGAACCCCGTCACCATCATTCAGAACTACACCAACATGTCCGATCCGATGAAAGAGCTGAACGCGGCGATTGAACCGCGGCGCTTgca  >  1:780277/1‑149 (MQ=255)
                                                                                                                                                |                                                                                                                                                    
AGAAAAATGAAGATCGCCGGACAGCCGAACGCTTTCAGAAATGGGTGGAAATGGGCGTTCTGACCGTTACCGATGGTGCGGAGGTGGATTATCGCTACATCCTCGAAGAGGCCAAAGCGGCGAACAAAATCAGCCCGGTCAGTGAGTCACCCATCGACCCCTTCGGGGCGACCGGGCTGTCACATGACCTTGCTGATGAAGACCTGAACCCCGTCACCATCATTCAGAACTACACCAACATGTCCGATCCGATGAAAGAGCTGGAAGCGGCGATTGAATCGGGGCGCTTTCA  >  NC_000913/1205029‑1205320

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