Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,778,124 A→C 9.0% intergenic (+342/+22) ypjF → / ← ypjA CP4‑57 prophage; toxin of the YpjF‑YfjZ toxin‑antitoxin system/adhesin‑like autotransporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,778,1240AC 9.0% 243.6 / 17.0 78intergenic (+342/+22)ypjF/ypjACP4‑57 prophage; toxin of the YpjF‑YfjZ toxin‑antitoxin system/adhesin‑like autotransporter
Reads supporting (aligned to +/- strand):  ref base A (39/32);  new base C (3/4);  total (42/36)
Fisher's exact test for biased strand distribution p-value = 6.97e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.92e-01

AGCCCACCAATCATGATTGGACGGTGTAAGGACAACACCAACAAAAACAGGAAGTTAGAAGTCTCAGCAAAACACCGACCAGACGGTGAGGAGACATAAAAGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCAGTTCACACCAGCTACCGCGTTCCACGGGGATTCCACACCGGCACTATGGCTATACCCCACCCCAAGATGCCCGCTTAACGTACTGCTGAATGAGGCTTTAATACCTGCCTGGTAT  >  NC_000913/2777977‑2778272
                                                                                                                                                   |                                                                                                                                                    
aGCCCACCAATCATGATTGGACGGTGTAAGGACAACACCAACAAAAACAGGAAGTTAGAAGTCTCAGCAAAACACCGACCAGACGGTGAGGAGACATAAAAGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGAt                                                                                                                                                     <  2:339484/149‑1 (MQ=255)
   ccaccaATCATGATTGGACGGTGTAAGGACAACACCAACAAAAACAGGAAGTTAGAAGTCTCAGCAAAACACCGACCAGACGGTGAGGAGACATAAAAGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGAtgtg                                                                                                                                                  <  3:348862/149‑1 (MQ=255)
        aaTCATGATTGGACGGTGTAAGGACAACACCAACAAAAACAGGAAGTTAGAAGTCTCAGCAAAACACCGACCAGACGGTGAGGAGACATAAAAGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACATCTGTGGGCtt                                                                                                                                             <  1:238564/149‑1 (MQ=255)
          tCATGATTGGACGGTGTAAGGACAACACCAACAAAAACAGGAAGTTAGAAGTCTCAGCAAAACACCGACCAGACGGTGAGGAGACATAAAAGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCtttt                                                                                                                                           <  2:41605/149‑1 (MQ=255)
           cATGATTGGACGGTGTAAGGACAACACCAACAAAAACAGGAAGTTAGAAGTCTCAGCAAAACACCGACCAGACGGTGAGGAGACATAAAAGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCttttt                                                                                                                                          >  2:458000/1‑149 (MQ=255)
            aTGATTGGACGGTGTAAGGACAACACCAACAAAAACAGGAAGTTAGAAGTCTCAGCAAAACACCGACCAGACGGTGAGGAGACATAAAAGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACATCTGTGGGCtttcat                                                                                                                                         <  2:218837/149‑4 (MQ=255)
            aTGATTGGACGGTGTAAGGACAACACCAACAAAAACAGGAAGTTAGAAGTATCAGCAAAACACCGACCAGACGGTGAGGAGACATAAAAGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACATCTGTGGGCtttcat                                                                                                                                         <  4:321256/149‑4 (MQ=255)
              gATTGGACGGTGTAAGGACAACACCAACAAAAACAGGAAGTTAGAAGTCTCAGCAAAACACCGACCAGACGGTGAGGAGACATAAAAGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGt                                                                                                                                       <  2:217507/149‑1 (MQ=255)
                  ggACGGTGTAAGGACAACACCAACAAAAACAGGAAGTTAGAAGTCTCAGCAAAACACCGACCAGACGGTGAGGAGACATAAAAGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATTAAAGCCCACAGATGTGGGCTTTTTCGTTGAt                                                                                                                                   <  2:3197/149‑1 (MQ=255)
                   ggcGGTGGAAGGACAAGACCAACAAAAACAGAATGTTAGAAGCCTCGGAAAAACAACGACCACCCGGTGAGAAGAGATAAAAGGGTACGGAAAGGAGCCGCGGCGCCTGGTATCATGAAAGCCCACCGATGTGGGCTTTTTCGTTGATg                                                                                                                                  <  4:416596/147‑1 (MQ=255)
                   gACGGTGTAAGGACAACACCAACAAAAACAGGAAGTTAGAAGTCTCAGCAAAACACCGACCAGACGGTGAGGAGACATAAAAGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATg                                                                                                                                  <  3:208899/149‑1 (MQ=255)
                   gACGGTGTAAGGACAACACCAACAAAAACAGGAAGTTAGAAGTCTCAGCAAAACACCGACCAGACGGTGAGGAGACATAAAAGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATg                                                                                                                                  <  2:51143/149‑1 (MQ=255)
                                caacaCCAACAAAAACAGGAAGTTAGAAGTCTCAGCAAAACACCGACCAGACGGTGAGGACACATAAAAGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCa                                                                                                                     <  3:473026/149‑1 (MQ=255)
                                 aacaCCAACAAAAACAGGAAGTTAGAAGTCTCAGCAAAACACCGACCAGACGGTGAGGAGACATAAAAGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCAg                                                                                                                    >  3:315898/1‑149 (MQ=255)
                                      caacaaAAACAGGAAGTTAGAAGTCTCAGCAAAACACCGACCAGACGGTGAGGAGACATAAAAGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCAGTTcac                                                                                                               >  1:240574/1‑149 (MQ=255)
                                        acaacaacaGGAAGTTAGAAGTCTCAGCAAAACACCGACCAGACGGTAAGGAGACATAAAAGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAATGCCCACAGAt                                                                                                                                                     >  2:577960/1‑109 (MQ=255)
                                        acaaAAACAGGAAGTTAGAAGTCTCAGCAAAACACCGACCAGACGGTGAGGAGACATAAAAGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGAt                                                                                                                                                     <  1:577960/109‑1 (MQ=255)
                                        acaaAAACAGGAAGTTAGAAGTCTCAGCAAAACACCGACCAGACGGTGAGGAGACATAAAAGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCAGTTCacac                                                                                                             >  4:348264/1‑149 (MQ=255)
                                          aaaaaCAGGAAGTTAGAAGTCTCAGCAAAACACCGACCAGACGGTGAGGAGACATAAAAGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGaa                                                                                                                           <  1:276579/133‑1 (MQ=255)
                                          aaaaaCAGGAAGTTAGAAGTCTCAGCAAAACACCGACCAGACGGTGAGGAGACATAAAAGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGaa                                                                                                                           >  2:276579/1‑133 (MQ=255)
                                                    aGTTAGAAGTCTCAGCAAAACACCGACCAGACGGTGAGGAGACATAAAAGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCAGTTCACACCAGCTACCGCGt                                                                                                 <  3:384889/149‑1 (MQ=255)
                                                    aGTTAGAAGTCTCAGCAAAACACCGACCAGACGGTGAGGAGACATAAAAGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCAGTTCACACCAGCTACCGCGt                                                                                                 <  4:77998/149‑1 (MQ=255)
                                                     gTTAGAAGTCTCAGCAAAACACCGACCAGACGGTGAGGAGACATAAAAGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCt                                                                                                                                              >  1:24050/1‑103 (MQ=255)
                                                     gTTAGAAGTCTCAGCAAAACACCGACCAGACGGTGAGGAGACATAAAAGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCt                                                                                                                                              <  2:24050/103‑1 (MQ=255)
                                                     gTTAGAAGTCTCAGCAAAACACCGACCAGACGGTGAGGAGACATAAAAGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCAGTTCACACCAGCTACCGCGtt                                                                                                <  2:109631/149‑1 (MQ=255)
                                                     gTTAGAAGTCTCAGCAAAACACCGACCAGACGGTGAGGAGACATAAAAGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCAGTTCACACCAGCTACCGCGtt                                                                                                >  3:601912/1‑149 (MQ=255)
                                                            gTCTCAGCAAAACACCGACCAGACGGTGAGGAGACATAAAAGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCAGTTCACACCAGCTACCGCGTTCCACggg                                                                                         <  3:344816/149‑1 (MQ=255)
                                                             tctcAGCAAAACACCGACCAGACGGTGAGGAGACATAAAAGGATACGCAAAGGAGTCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCAGTTCACACCAGCTACCGCGTTCCACgggg                                                                                        >  4:629938/1‑149 (MQ=255)
                                                               tcAGCAAAACACCGACCAGACGGTGAGGAGACATAAAAGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCAGTTCACACCAGCTACCGCGTTCCACGGGGAt                                                                                      >  3:180665/1‑149 (MQ=255)
                                                               tcAGCAAAACACCGACCAGACGGTGAGGAGACATAAAAGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCAGTTCACACCAGCTACCGCGTTCCACGGGGAt                                                                                      >  2:325600/1‑149 (MQ=255)
                                                               tcAGCAAAACACCGACCAGACGGTGAGGAGACATAAAAGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCAGTTCACACCAGCTACCGCGTTCCACGGGGAt                                                                                      >  3:604606/1‑149 (MQ=255)
                                                               tcAGCAAAACACCGACCAGACGGTGAGGAGACATAAAAGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCAGTTCACACCAGCTACCGCGTTCCACGGGGAt                                                                                      >  2:7948/1‑149 (MQ=255)
                                                               tcAGCAAAACACCGACCAGACGGTGAGGAGACATAAAAGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCAGTTCACACCAGCTACCGCGTTCCACGGGGAt                                                                                      >  3:610278/1‑149 (MQ=255)
                                                                 aGCAAAACACCGACCAGACGGTGAGGAGACATAAAAGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCAGTTCACACCAGCTACCGCGTTCCACGGGGATTc                                                                                    <  1:356936/149‑1 (MQ=255)
                                                                   cAAAACACCGACCAGACGGTGAGGAGACATAAAAGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCAGTTCACACCAGCTACCGCGTTCCACGGGGATTCca                                                                                  <  2:568159/149‑1 (MQ=255)
                                                                    aaaaCACCGACCAGACGGTGAGGAGACATAAAAGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCAGTTCACACCAGCTACCGCGTTCCACGGGGATTCcac                                                                                 >  3:434054/1‑149 (MQ=255)
                                                                         acCGACCAGACGGTGAGGAGACATAAAAGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACATCTGTGGGCtttcat                                                                                                                                         >  2:374498/1‑85 (MQ=39)
                                                                         acCGACCAGACGGTGAGGAAACATAAAAGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACATCTGTGGGCtttcat                                                                                                                                         <  1:374498/88‑4 (MQ=38)
                                                                              ccAGACGGTGAGGAGACATAAAAGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAAc                                                                                                                          >  1:439364/1‑98 (MQ=255)
                                                                              ccAGACGGTGAGGAGACATAAAAGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAAc                                                                                                                          <  2:439364/98‑1 (MQ=255)
                                                                               cAGACGGTGAGGAGACATAAAAGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCAGTTCACACCAGCTACCGCGTTCCACGGGGATTCCACACCGGCACTAt                                                                      <  4:601912/149‑1 (MQ=255)
                                                                                aGACGGTGAGGAGACATAAAAGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCTCAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCAGTTCACACCAGCTACCGCGTTCCACGGGGATTCCACACCGGCACTATg                                                                     <  3:358116/149‑1 (MQ=255)
                                                                                aGACGGTGAGGAGACATAAAAGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCAGTTCACACCAGCTACCGCGTTCCACGGGGATTCCACACCGGCACTATg                                                                     >  1:284820/1‑149 (MQ=255)
                                                                                 gACGGTGAGGAGACATAAAAGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCAGTTCACACCAGCTACCGCGTTCCACGGGGATTCCACACCGGCACTATgg                                                                    <  3:119396/149‑1 (MQ=255)
                                                                                  aCGGTGAGGAGACATAAAAGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCAGTTCACACCAGCTACCGCGTTCCACGGGGATTCCACACCGGCACTATGGc                                                                   <  1:406318/149‑1 (MQ=255)
                                                                                      tGAGGAGACATAAAAGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCAGTTCACACCAGCTACCGCGTTCCACGGGGATTCCACACCGGCACTATGGCtata                                                               >  2:595879/1‑149 (MQ=255)
                                                                                         ggagACATAAAAGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCAGTTCACACCAGCTACCGCGTTCCACGGGGATTCCACACCGGCACTATGGCTATaccc                                                            >  4:34188/1‑149 (MQ=255)
                                                                                                 aaaaGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCAGTTCACACCAGCTACCGCGTTCCACGGGGATTCCACACCGGCACTATGGCTATACCCCACCCCaa                                                    >  2:311590/1‑149 (MQ=255)
                                                                                                 aaaaGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCAGTTCACACCAGCTACCGCGTTCCACGGGGATTCCACACCGGCACTATGGCTATACCCCACCCCaa                                                    >  2:503306/1‑149 (MQ=255)
                                                                                                  aaaGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCAGTTCACACCAGCTACCGCGTTCCACGGGGATTCCACACCGGCACTATGGCTATACCCCACCCCAAg                                                   >  1:699212/1‑149 (MQ=255)
                                                                                                   aaGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCAGTTCACACCAGCTACCGCGTTCCACGGGGATTCCACACCGGCACTATGGCTATACCCCACCCCAAGa                                                  >  2:279379/1‑149 (MQ=255)
                                                                                                       aTACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCAGTTCACACCAGCTACCGCGTTCCACGGGGATTCCACACCGGCACTATGGCTATACCCCACCCCAAGATGcc                                              >  4:663085/1‑149 (MQ=255)
                                                                                                        tACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCAGTTCACACCAGCTACCGCGTTCCACGGGGATTCCACACCGGCACTATGGCTATACCCCACCCCAAGATGccc                                             <  1:377734/149‑1 (MQ=255)
                                                                                                          ttaccaGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCAGTTCACACCAGCTACCGCGTTCCACGGGGATTCCACACCGGCACTATGGCTATACCCCACCCCAAGATGCCCGc                                           >  4:130088/6‑149 (MQ=255)
                                                                                                          cGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCAGTTCACACCAGCTACCGCGTTCCACGGGGATTCCACACCGGCACTATGGCTATACCCCACCCCAAGATGCCCGc                                           <  1:458000/149‑1 (MQ=255)
                                                                                                              caGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCAGTTCACACCAGCTACCGCGTTCCACGGGGATTCCACACCGGCACTATGGCTATACCCCACCCCAAGATGCCCGCTTaa                                       >  1:353356/2‑149 (MQ=255)
                                                                                                               aGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCAGTTCACACCAGCTACCGCGTTCCACGGGGATTCCACACCGGCACTATGGCTATACCCCACCCCAAGATGCCCGCTTAAc                                      >  3:69998/1‑149 (MQ=255)
                                                                                                               aGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCAGTTCACACCAGCTACCGCGTTCCACGGGGATTCCACACCGGCACTATGGCTATACCCCACCCCAAGATGCCCGCTTAAc                                      >  4:256703/1‑149 (MQ=255)
                                                                                                               aGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCAGTTCACACCAGCTACCGCGTTCCACGGGGATTCCACACCGGCACTATGGCTATACCCCACCCCAAGATGCCCGCTTAAc                                      <  3:193975/149‑1 (MQ=255)
                                                                                                               aGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCAGTTCACACCAGCTACCGCGTTCCACGGGGATTCCACACCGGCACTATGGCTATACCCCACCCCAAGATGCCCGCTTAAc                                      >  2:460567/1‑149 (MQ=255)
                                                                                                                      gcgGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCAGTTCACACCAGCTACCGCGTTCCACGGGGATTCCACACCGGCACTATGGCTATACCCCACCCCAAGATGCCCGCTTAACGTActgc                               >  1:628221/1‑149 (MQ=255)
                                                                                                                      gcgGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCAGTTCACACCAGCTACCGCGTTCCACGGGGATTCCACACCGGCACTATGGCTATACCCCACCCCAAGATGCCCGCTTAACGTActgc                               >  1:523479/1‑149 (MQ=255)
                                                                                                                               ggTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCAGTTCACACCAGCTACCGCGTTCCACGGGGATTCCACACCGGCACTATGGCTATACCCCACCCCAAGATGCCCGCTTAACGTACTGCTGAATGAgg                      >  2:274970/1‑149 (MQ=255)
                                                                                                                                   aCATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCAGTTCACACCAGCTACCGCGTTCCa                                                                                             <  1:666863/74‑1 (MQ=255)
                                                                                                                                   aCATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCAGTTCACACCAGCTACCGCGTTCCa                                                                                             >  2:666863/1‑74 (MQ=255)
                                                                                                                                     aTGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCAGTTCACACCAGCTACCGCGTTCCACGGGGATTCCACACCGGCACTATGGCTATACCCCACCCCAAGATGCCCGCTTAACGTACTGCTGAATGAGGCTTTaa                <  2:652158/149‑1 (MQ=255)
                                                                                                                                     aTGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCAGTTCACACCAGCTACCGCGTTCCACGGGGATTCCACACCGGCACTATGGCTATACCCCACCCCAAGATGCCCGCTTAACGTACTGCTGAATGAGGCTTTaa                >  1:261788/1‑149 (MQ=255)
                                                                                                                                       aaaaGCCCACATCTGTGGGCTTTTTCGTTGATGGTCAGAACGACCAGTTCACACCAGCTACCGCGTTCCACGGGGATTCCACACCGGCACTATGGCTATACCCCACCCCAAGATGCCCGCTTAACGTACTGCTGAATGAGGCTTTAATa              >  3:478416/2‑149 (MQ=255)
                                                                                                                                        aaaGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCAGTTCACACCAGCTACCGCGTTCCACGGGGATTCCACACCGGCACTATGGCTATACCCCACCCCAAGATGCCCGCTTAACGTACTGCTGAATGAGGCTTTAATAc             <  3:639968/149‑1 (MQ=255)
                                                                                                                                        aaaGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCAGTTCACACCAGCTACCGCGTTCCACGGGGATTCCACACCGGCACTATGGCTATACCCCACCCCAAGATGCCCGCTTAACGTACTGCTGAATGAGGCTTTAATAc             <  3:412566/149‑1 (MQ=255)
                                                                                                                                         aaGCCCACATCTGTGGGCTTTTTCGTTGATGGTCAGAACGACCAGTTCACACCAGCTACCGCGTTCCACGGGGATTCCACACCGGCACTATGGCTATACCCCACCCCAAGATGCCCGCTTAACGTACTGCTGAATAAGGCTTTAATAcc            >  2:691669/1‑149 (MQ=255)
                                                                                                                                         aaGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCAGTTCACACCAGCTACCGCGTTCCACGGGGATTCCACACCGGCACTATGGCTATACCCCACCCCAAGATGCCCGCTTAACGTACTGCTGAATGAGGCTTTAATAcc            >  3:178005/1‑149 (MQ=255)
                                                                                                                                           gCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCAGTTCACACCAGCTACCGCGTTCCACGGGGATTCCACACCGGCACTATGGCTATACCCCACCCCAAGATGCCCGCTTAACGTACTGCTGAATGAGGCTTTAATAcctg          >  3:32481/1‑149 (MQ=255)
                                                                                                                                            cccACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCAGTTCACACCAGCTACCGCGTTCCACGGGGATTCCACACCGGCACTATGGCTATACCCCACCCCAAGATGCCCGCTTAACGTACTGCTGAATGAGGCTTTAATAcctgc         >  1:264902/1‑149 (MQ=255)
                                                                                                                                              cacaGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCAGTTCACACCAGCTACCGCGTTCCACGGGGATTCCACACCGGCACTATGGCTATACCCCACCCCAAGATGCCCGCTTAACGTACTGCTGAATGAGGCTTTAATAcctgcct       >  4:214197/1‑149 (MQ=255)
                                                                                                                                              cacaGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCAGTTCACACCAGCTACCGCGTTCCACGGGGATTCCACACCGGCACTATGGCTATACCCCACCCCAAGATGCCCGCTTAACGTACTGCTGAATGAGGCTTTAATAcctgcct       >  2:589475/1‑149 (MQ=255)
                                                                                                                                              cacaGATGTAGGCTTTTTCGTTGATGGTCAGAACGACCAGTTCACACCAGCTACCGCGTTCCACGGGGATTCCACACGGGCACTATGGCTATACCCCACCCCAAGATGCCCGCTTAACGTACTGCTGAATGAGGCTTTAATAcctgcct       <  3:335239/149‑1 (MQ=255)
                                                                                                                                               acaGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCAGTTCACACCAGCTACCGCGTTCCACGGGGATTCCACACCGGCACTATGGCTATACCCCACCCCAAGATGCCCGCTTAACGTACTGCTGAATGAGGCTTTAATAcctgcctg      <  3:130088/149‑1 (MQ=255)
                                                                                                                                                   aTGTGGGCTTTTTCGTTGATGGTCAGAACGACCAGTTCACACCAGCTACCGCGTTCCACGGGGATTCCACACCGGCACTATGGCTATACCCCACCCCAAGATGCCCGCTTAACGTACTGCTGAATGAGGCTTTAATACCTGCCTGGtat  >  2:475957/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
AGCCCACCAATCATGATTGGACGGTGTAAGGACAACACCAACAAAAACAGGAAGTTAGAAGTCTCAGCAAAACACCGACCAGACGGTGAGGAGACATAAAAGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCAGTTCACACCAGCTACCGCGTTCCACGGGGATTCCACACCGGCACTATGGCTATACCCCACCCCAAGATGCCCGCTTAACGTACTGCTGAATGAGGCTTTAATACCTGCCTGGTAT  >  NC_000913/2777977‑2778272

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