Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,777,571 A→G 5.9% D40G (GAC→GGC)  ypjF → CP4‑57 prophage; toxin of the YpjF‑YfjZ toxin‑antitoxin system

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,777,5710AG 5.9% 667.0 / 16.8 197D40G (GAC→GGC) ypjFCP4‑57 prophage; toxin of the YpjF‑YfjZ toxin‑antitoxin system
Reads supporting (aligned to +/- strand):  ref base A (88/97);  new base G (6/6);  total (94/103)
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

Only 100 of 212 total aligned reads displayed.

TCAACGCTAACTAATTTCACGAGAGCAAGCATGAACACTCTACCTGCTACAATTTCGCAGGCGGCGAAGCCCTGCCTGTCGCCAGTGGCTGTCTGGCAAATGCTACTGACACGCCTGCTGGAACAACACTATGGCCTGACACTGAACGACACGCCGTTCAGTGATGAAACTGTTATTAAGGAACATATCGATGCTGGTATCACTCTGGCCGATGCAGTCAATTTTCTGGTGGAAAAGTACGAACTGGTACGTATCGATCACAGAGGATTTTCGTGGCAACAACAGTCTCCA  >  NC_000913/2777423‑2777713
                                                                                                                                                    |                                                                                                                                              
tCAACGCTAACTAATTTCACGAGAGCAAGCATGAACACTCTACCTGCTACAATTTCGCAGGCGGCGAAGCCCTGCCTGTCGCCAGTGGCTGTCTGGCAAATGCTACTGACACGCCTGCTGGAACAACACTATGGCCTGACACTGAACGa                                                                                                                                                >  3:406114/1‑149 (MQ=255)
 cAACGCTAACTAATTTCACGAGAGCAAGCATGAACACTCTACCTGCTACAATTTCGCAGGCGGCGAAGCCCTGCCTGTCGCCAGTGGCTGTCTGGCAAATGCTACTGACACGCCTGCTGGAACAACACTATGGCCTGACACTGAACGac                                                                                                                                               <  3:295355/149‑1 (MQ=255)
  aaCGCTAACTAATTTCACGAGAGCAAGCATGAACACTCGACCTGCTACAATTTCGCAGGCGGCGAAGCCCTGCCTGTCGCCAGTGGCTGTCTGGCAAATGCTACTGACACGCCTGCTGGAACAACACTATGGCCTGACACTGAACGaca                                                                                                                                              >  4:1210765/1‑149 (MQ=255)
         actaaTTTCACGAGAGCAAGCATGAACACTCTACCTGCTACAATTTCGCAGGCGGCGAAGCCCTGCCTGTCGCCAGTGGCTGTCTGGCAAATGCTACTGACACGCCTGCTGGAACAACACTATGGCCTGACACTGACCGa                                                                                                                                                >  3:653834/1‑140 (MQ=255)
          ctaaTTTCACGAGAGCAAGCATGAACACTCTACCTGCTACAATTTCGCAGGCGGCGAAGCCCTGCCTGTCGCCAGTGGCTGTCTGGCAAATGCTACTGACACGCCTGCTGGAACAACACTATGGCCTGACACTGAACGACACGCCGTTc                                                                                                                                      >  1:565262/1‑149 (MQ=255)
           taaTTTCACGAGAGCAAGCATGAACACTCTACCTGCTACAATTTCGCAGGCGGCGAAGCCCTGCCTGTCGCCAGTGGCTGTCTGGCAAATGCTACTGACACGCCTGCTGGAACAACACTATGGCCTGACACTGAACGACACGCCGTTCa                                                                                                                                     <  1:480796/149‑1 (MQ=255)
             aTTTCACGAGAGCAAGCATGAACACTCTACCTGCTACAATTTCGCAGGCGGCGAAGCCCTGCCTGTCGCCAGTGGCTGTCTGGCAAATGCTACTGACACGCCTGCTGGAACAACACTATGGCCTGACACTGAACGACACGCCGTTCAGt                                                                                                                                   <  1:672037/149‑1 (MQ=255)
              tttCACGAGAGCAAGCATGAACACTCTACCTGCTACAATTTCGCAGGCGGCGAAGCCCTGCCTGTCGCCAGTGGCTGTCTGGCAAATGCTACTGACACGCCTGCTGGAACAACACTATGGCCTGACACTGAACGaca                                                                                                                                              >  1:795818/1‑137 (MQ=255)
              tttCACGAGAGCAAGCATGAACACTCTACCTGCTACAATTTCGCAGGCGGCGAAGCCCTGCCTGTCGCCAGTGGCTGTCTGGCAAATGCTACTGACACGCCTGCTGGAACAACACTATGGCCTGACACTGAACGACACGCCGTTCAGtg                                                                                                                                  >  1:1351971/1‑149 (MQ=255)
                tCACGAGAGCAAGCATGAACACTCTACCTGCTACAATTTCGCAGGCGGTGAAGCCCTGCCTGTCGCCAGTGGCTGTCTGGCAAATGCTACTGACACGCCTGCTGGAACAACACTATGGCCTGACACTGAACGACACGCCGTTCAGtgat                                                                                                                                <  3:1210765/149‑1 (MQ=255)
                tCACGAGAGCAAGCATGAACACTCTACCTGCTACAATTTCGCAGGCGGCGAAGCCCTGCCTGTCGCCAGAGGCTGTCTGGCAAATGCTACTGACACGCCTGCTGGAACAACACTATGGCCTGACACTGAACGACACGCCGTTCAGtgat                                                                                                                                >  2:1134359/1‑149 (MQ=255)
                  aCGAGAGCAAGCATGAACACTCTACCTGCTACAATTTCGCAGGCGGCGAAGCCCTGCCTGTCGCCAGTGGCTGTCTGGCAAATGCTACTGACACGCCTGCTGGAACAACACTATGGCCTGACACTGAACGACACGCCGTTCAGtgatga                                                                                                                              >  2:351360/1‑149 (MQ=255)
                    gagaGCAAGCATGAACACTCTACCTGCTACAATTTCGCAGGCGGCGAAGCCCTGCCTGTCGCCAGTGGCTGTCTGGCAAATGCTACTGACACGCCTGCTGGAACAACACTATGGCCTGACACTGAACGACACGCCGTTCAGTGATGaaa                                                                                                                            <  4:406114/149‑1 (MQ=255)
                    gagaGCAAGCATGAACACTCTACCTGCTACAATTTCGCAGGCGGCGAAGCCCTGCCTGTCGCCAGTGGCTGTCTGGCAAATGCTACTGACACGCCTGCTGGAACAACACTATGGCCTGACACTGAACGACACGCCGTTCAGTGATGaaa                                                                                                                            >  4:767560/1‑149 (MQ=255)
                        gcaagcaTGAACACTCTACCTGCTACAATTTCGCAGGCGGCGAAGCCCTGCCTGTCGCCAGTGGCTGTCTGGCAAATGCTACTGACACGCCTGCTGGAACAACACTATGGCCTGACACTGAACGACACGCCGTTCAGTGATGAAACTGt                                                                                                                        <  4:29436/149‑1 (MQ=255)
                         caagcaTGAACACTCTACCTGCTACAATTTCGCAGGCGGCGAAGCCCTGCCTGTCGCCAGTGGCTGTCTGGCAAATGCTACTGACACGCCTGCTGGAACAACACTATGGCCTGACACTGAACGACACGCCGTTCAGTGATGAAACTGtt                                                                                                                       <  2:1351971/149‑1 (MQ=255)
                          aagcaTGAACACTCTACCTGCTACAATTTCGCAGGCGGCGAAGCCCTGCCTGTCGCCAGTGGCTGTCTGGCAAATGCTACTGACACGCCTGCTGGAACAACACTATGGCCTGACACTGAACGACACGCCGTTCAGTGATGAAACTGtta                                                                                                                      <  1:1132788/149‑1 (MQ=255)
                            gcaTGAACACTCTACCTGCTACAATTTCGCAGGCGGCGAAGCCCTGCCTGTCGCCAGTGGCTGTCTGGCAAATGCTACTGACACGCCTGCTGGAACAACACTATGGCCTGACACTGAACGGCGTGTCGTTCAGt                                                                                                                                   >  1:849473/1‑134 (MQ=255)
                            gcaTGAACACTCTACCTGCTACAATTTCGCAGGCGGCGAAGCCCTGCCTGTCGCCAGTGGCTGTCTGGCAAATGCTACTGACACGCCTGCTGGAACAACACTATGGCCTGACACTGAACGACACGCCGTTc                                                                                                                                      <  1:399259/131‑1 (MQ=255)
                                 aaCACTCTACCTGCTACAATTTCGCAGGCGGCGAAGCCCTGCCTGTCGCCAGTGGCTGTCTGGCAAATGCTACTGACACGCCTGCTGGAACAACACTATGGCCTGACACTGAACGACACGCCGTTCAGTGATGAAACTGTTATTAAGGa                                                                                                               <  1:609799/149‑1 (MQ=255)
                                 aaCACTCTACCTGCTACAATTTCGCAGGCGGCGAAGCCCTGCCTGTCGCCAGTGGCTGTCTGGCAAATGCTACTGACACGCCTGCTGGAACAACACTATGGCCTGACACTGAACGACACGCCGTTCAGTGATGAAACTGTTATTAAGGa                                                                                                               <  2:78417/149‑1 (MQ=255)
                                  acacTCTACCTGCTACAATTTCGCAGGCGGCGAAGCCCTGCCTGTCGCCAGTGGCTGTCTGGCAAATGCTACTGACACGCCTGCTGGAACAACACTATGGCCTGACACTGAACGACACGCCGTTCAGTGATGAAACTGTTATTAAGGaa                                                                                                              <  3:52972/149‑1 (MQ=255)
                                  acacTCTACCTGCTACAATTTCGCAGGCGGCGAAGCCCTGCCTGTCGCCAGTGGCTGTCTGGCAAATGCTACTGACACGCCTGCTGGAACAACACTATGGCCTGACACTGAACGACACGCCGTTCAGTGATGAAACTGTTATTAAGGaa                                                                                                              >  4:545982/1‑149 (MQ=255)
                                     ctctACCTGCTACAATTTCGCAGGCGGCGAAGCCCTGCCTGTCGCCAGTGGCTGTCTGGCAAATGCTACTGACACGCCTGCTGGAACAACACTATGGCCTGACACTGAACGACACGCCGTTCAGTGATGAAACTGtta                                                                                                                      >  4:739608/1‑138 (MQ=255)
                                      tctACCTGCTACAATTTCGCAGGCGGCGAAGCCCTGCCTGTCGCCAGTGGCTGTCTGGCAAATGCTACTGACACGCCTGCTGGAACAACACTATGGCCTGACACTGAACGACACGCCGTTCAGTGATGAAACTg                                                                                                                         <  3:759355/134‑1 (MQ=255)
                                      tctACCTGCTACAATTTCGCAGGCGGCGAAGCCCTGCCTGTCGCCAGTGGCTGTCTGGCAAATGCTACTGACACGCCTGCTGGAACAACACTATGGCCTGACACTGAACGACACGCCGTTCAGTGATGAAACTGTTATTAAGGAACata                                                                                                          >  4:190716/1‑149 (MQ=255)
                                           cTGCTACAATTTCGCAGGCGGCGAAGCCCTGCCTGTCGCCAGTGGCTGTCTGGCAAATGCTACTGACACGCCTGCTGGAACAACACTATGGCCTGACACTGAACGACACGCCGTTCAGTGATGAAACTGTTATTAAGGAACATATCGAt                                                                                                     <  2:322619/149‑1 (MQ=255)
                                            tGCTACAATTTCGCAGGCGGCGAAGCCCTGCCTGTCGCCAGTGGCTGTCTGGCAAATGCTACTGACACGCCTGCTGGAACAACACTATGGCCTGACACTGAACGACACGCCGTTCAGtgatg                                                                                                                               <  2:260372/122‑1 (MQ=255)
                                               tACAATTTCGCAGGCGGCGAAGCCCTGCCTGTCGCCAGTGGCTGTCTGGCAAATGCTACTGACACGCCTGCTGGAACAACACTATGGCCTGACACTGAACGACACGCCGt                                                                                                                                        >  3:816985/1‑110 (MQ=255)
                                               tACAATTTCGCAGGCGGCGAAGCCCTGCCTGTCGCCAGTGGCTGTCTGGCAAATGCTACTGACACGCCTGCTGGAACAACACTATGGCCTGACACTGAACGACACGCCGTTCAGTGATGAAACTGTTATTAAGGAACATATCGATGCTg                                                                                                 <  4:678915/149‑1 (MQ=255)
                                                  aaTTTCGCAGGCGGCGAAGCCCTGCCTGTCGCCAGTGGCTGTCTGGCAAATGCTACTGACACGCCTGCTGGAACAACACTATGGCCTGACACTGAACGACACGCCGTTCAGTGATGAAACTGTTATTAAGGAACATATCGATGCTg                                                                                                 <  2:145958/146‑1 (MQ=255)
                                                  aaTTTCGCAGGCGGCGAAGCCCTGCCTGTCGCCAGTGGCTGTCTGGCAAATGCTACTGACACGCCTGCTGGAACAACACTATGGCCTGACACTGAACGACACGCCGTTCAGTGATGAAACTGTTATTAAGGAACATATCGATGCTg                                                                                                 <  4:1376067/146‑1 (MQ=255)
                                                  aaTTTCGCAGGCGGCGAAGCCCTGCCTGTCGCCAGTGGCTGTCTGGCAAATGCTACTGACACGCCTGCTGGAACAACACTATGGCCTGACACTGAACGACACGCCGTTCAGTGATGAAACTGTTATTAAGGAACATATCGATGCTGGTa                                                                                              >  3:317815/1‑149 (MQ=255)
                                                     ttCGCAGGCGGCGAAGCCCTGCCTGTCGCCAGTGGCTGTCTGGCAAATGCTACTGACACGCCTGCTGGAACAACACTATGGCCTGACACTGAACGACACGCCGTTCAGTGATGAAACTGTTATTAAGGAACATATCGATGCTGGTATCa                                                                                           <  2:935950/149‑1 (MQ=255)
                                                                 gAAGCCCTGCCTGTCGCCAGTGGCTGTCTGGCAAATGCTACTGACACGCCTGCTGGAACAACACTATGGCCTGACACTGAACGACACGCCGTTCAGTGATGAAACTGtta                                                                                                                      <  2:538585/110‑1 (MQ=255)
                                                                      cctgcctgTCGCCAGTGGCTGTCTGGCAAATGCTACTGACACGCCTGCTGGAACAACACTATGGCCTGACACTGAACGACACGCCGTTCAGTGATGAAACTGTTATTAAGGAACATATCGATGCTGGTATCACTCTGGCCGATGCAGTc                                                                          >  4:1203080/1‑149 (MQ=255)
                                                                       ctgcctgTCGCCAGTGGCTGTCTGGCAAATGCTACTGACACGCCTGCTGGAACAACACTATGGCCTGACACTGAACGACACGCCGtt                                                                                                                                       <  3:912442/87‑1 (MQ=255)
                                                                       ctgcctgTCGCCAGTGGCTGTCTGGCAAATGCTACTGACACGCCTGCTGGAACAACACTATGGCCTGACACTGAACGACACGCCGTTCAGTGATGAAACTGtt                                                                                                                       >  1:1059591/1‑103 (MQ=255)
                                                                       ctgcctgTCGCCAGTGGCTGTCTGGCAAATGCTACTGACACGCCTGCTGGAACAACACTATGGCCTGACACTGAACGACACGCCGTTCAGTGATGAAACTGtt                                                                                                                       <  2:1059591/103‑1 (MQ=255)
                                                                           ctgTCGCCAGTGGCTGTCTGGCAAATGCTACTGACACGCCTGCTGGAACAACACTATGGCCTGACACTGAACGACACGCCGTTCAGTGATGAAAc                                                                                                                           >  2:1089580/1‑95 (MQ=255)
                                                                           ctgTCGCCAGTGGCTGTCTGGCAAATGCTACTGACACGCCTGCTGGAACAACACTATGGCCTGACACTGAACGACACGCCGTTCAGTGATGAAAc                                                                                                                           <  1:1089580/95‑1 (MQ=255)
                                                                           ctgTCGCCAGTGGCTGTCTGGCAAATGCTACTGACACGCCTGCTGGAACAACACTATGGCCTGACACTGAACGACACGCCGTTCAGTGATGAAACTGTTATTAAGGAACATATCGATGCTGGTATCACTCTGGCCGATGCAGTCAAttt                                                                     >  3:221318/1‑149 (MQ=255)
                                                                            tgTCGCCAGTGGCTGTCTGGCAAATGCTACTGACACGCCTGCTGGAACAACACTATGGCCTGACACTGAACGACACGCCGTTCAGTGATGAAACTGTTATTAAGGAACATATCGATGCTGGTATCACTCTGGCCGATGCAGTCAAtttt                                                                    <  4:317815/149‑1 (MQ=255)
                                                                             gTCGCCAGTGGCTGTCTGGCAAATGCTACTGACACGCCTGCTGGAACAACACTATGGCCTGACACTGAACGACACGCCGTTCAGTGATGAAACTGTTATTAAGGAACATATCGATGCTg                                                                                                 >  2:852699/1‑119 (MQ=255)
                                                                             gTCGCCAGTGGCTGTCTGGCAAATGCTACTGACACGCCTGCTGGAACAACACTATGGCCTGACACTGAACGACACGCCGTTCAGTGATGAAACTGTTATTAAGGAACATATCGATGCTGGTAt                                                                                             <  4:1228576/123‑1 (MQ=255)
                                                                             gTCGCCAGTGGCTGTCTGGCAAATGCTACTGACACGCCTGCTGGAACAACACTATGGCCTGACACTGAACGACACGCCGTTCAGTGATGAAACTGTTATTAAGGAACATATCGATGCTGGTAt                                                                                             >  3:1228576/1‑123 (MQ=255)
                                                                             gTCGCCAGTGGCTGTCTGGCAAATGCTACTGACACGCCTGCTGGAACAACACTATGGCCTGACACTGAACGACACGCCGTTCAGTGATGAAACTGTTATTAAGGAACATATCGATGCTGGTATCACTCTGGCCGATGCAGTCAATTTTc                                                                   <  3:174762/149‑1 (MQ=255)
                                                                               cGCCAGTGGCTGTCTGGCAAATGCTACTGACACGCCTGCTGGAACAACACTATGGCCTGACACTGAACGACACGCCGTTCAGTGATGAAACTGTTATTAAGGAACATATCGATGCTGGTATCACTCTGGCCGATGCAGTCAATTTTCtg                                                                 <  3:1102990/149‑1 (MQ=255)
                                                                                  cAGTGGCTGTCTGGCAAATGCTACTGACACGCCTGCTGGAACAACACTATGGCCTGACACTGAACGACACGCCGTTCAGTGATGAAACTGTTATTAAGGAACata                                                                                                          >  3:775048/1‑105 (MQ=255)
                                                                                  cAGTGGCTGTCTGGCAAATGCTACTGACACGCCTGCTGGAACAACACTATGGCCTGACACTGAACGACACGCCGTTCAGTGATGAAACTGTTATTAAGGAACata                                                                                                          <  2:570667/105‑1 (MQ=255)
                                                                                       gCTGTCTGGCAAATGCTACTGACACGCCTGCTGGAACAACACTATGGCCTGACACTGAACGACACGCCGTTCAGTGATGAAACTGTTATTAAGGAACATATCGATGCTGGTATCACTCTGGCCGATGCAGTCAATTTTCTGGTGGaaaa                                                         >  2:415912/1‑149 (MQ=255)
                                                                                          gTCTGGCAAATGCTACTGACACGCCTGCTGGAACAACACTATGGCCTGACACTGAACGACACGCCGTTCAGTGATGAAACTGTTATTAAGGAACATATCGATGCTGGTATCACTCTGGCCGATGCAGTCAATTTTCTGGTGGAAAAGTa                                                      <  2:493434/149‑1 (MQ=255)
                                                                                           tCTGGCAAATGCTACTGACACGCCTGCTGGAACAACACTATGGCCTGACACTGAACGACACGCCGTTCAGTGATGAAACTGTTATTAAGGAACatat                                                                                                         <  4:58886/97‑1 (MQ=255)
                                                                                             tGGCAAATGCTACTGACACGCCTGCTGGAACAACACTATGGCCTGACACTGAACGACACGCCGTTCAGTGATGAAACTGTTATTAAGGAACATATCGATGCTGGTATCACTCTGGCCGATGCAGTCAATTTTCTGGTGGAAAAGTACGa                                                   >  4:709522/1‑149 (MQ=255)
                                                                                             tGGCAAATGCTACTGACACGCCTGCTGGAACAACACTATGGCCTGACACTGAACGACACGCCGTTCAGTGATGAAACTGTTATTAAGGAACATATCGATGCTGGTATCACTCTGGCCGATGCAGTCAATTTTCTGGTGGAAAAGTACGa                                                   >  3:699009/1‑149 (MQ=255)
                                                                                              ggCAAATGCTACTGACACGCCTGCTGGAACAACACTATGGCCTGACACTGAACGACACGCCGTTCAGTGATGAAACTGTTATTAAGGAACATATCGATGc                                                                                                   >  2:991095/1‑100 (MQ=255)
                                                                                               gCAAATGCTACTGACACGCCTGCTGGAACAACACTATGGCCTGACACTGAACGACACGCCGTTCAGTGATGAAACTGTTATTAAGGAACATATCGATGCTGGTATCACTCTGGCCGATGCAGTCAATTTTCTGGTGGAAAAGTACGAAc                                                 <  3:192571/149‑1 (MQ=255)
                                                                                                  aaTGCTACTGACACGCCTGCTGGAACAACACTATGGCCTGACACTGAACGACACGCCGTTCAGTGATGAAACTGTTATTAAGGAACATATCGATGCTGGTATCACTCTGGCCGATGCAGTCAATTTTCTGGTGGAAAAGTACGAACTgg                                              >  3:1367008/1‑149 (MQ=255)
                                                                                                    tGCTACTGACACGCCTGCTGGAACAACACTATGGCCTGACACTGAACGACACGCCGTTCAGTGATGAAACTGTTATTAAGGAACATATCGATGCTGGTATCACTCTGGCCGATGCAGTCAATTTTCTGGTGGAAAAGTACGAACTGGTa                                            <  2:611727/149‑1 (MQ=255)
                                                                                                       tACTGACACGCCTGCTGGAACAACACTATGGCCTGACACTGAACGACACGCCGTTCAGTGATGAAACTGTTATTAAGGAACATATCGATg                                                                                                    <  4:277377/90‑1 (MQ=255)
                                                                                                       tACTGACACGCCTGCTGGAACAACACTATGGCCTGACACTGAACGACACGCCGTTCAGTGATGAAACTGTTATTAAGGAACATATCGATGCTGGTATCACTCTGGCCGATGCAGTCAATTTTCTGGTGGAAAAGTACGAACTGGTa                                            >  4:216889/1‑146 (MQ=255)
                                                                                                         cTGACACGCCTGCTGGAACAACACTATGGCCTGACACTGAACGACACGCCGTTCAGTGATGAAACTGTTATTAAGGAACATATCGATGCTGGTATCACTCTGGCCGATGCAGTCAATTTTCTGGTGGAAAAGTACGAACTGGTACGTAt                                       >  4:834584/1‑149 (MQ=255)
                                                                                                          tGACACGCCTGCTGGAACAACACTATGGCCTGACACTGAACGACACGCCGTTCAGt                                                                                                                                   >  2:1404866/1‑56 (MQ=255)
                                                                                                          tGACACGCCTGCTGGAACAACACTATGGCCTGACACTGAACGACACGCCGTTCAGt                                                                                                                                   <  1:1404866/56‑1 (MQ=255)
                                                                                                          tGACACGCCTGCTGGAACAACACTATGGCCTGACACTGAACGACACGCCGTTCAGTGATGAAACTGTTATTAAGGAACATATCGATGCTGGTATCACTCTGGCCGATGCAGTCAATTTTCTGGTGGAAAAGTAc                                                     <  4:322723/134‑1 (MQ=255)
                                                                                                          tGACACGCCTGCTGGAACAACACTATGGCCTGACACTGAACGACACGCCGTTCAGTGATGAAACTGTTATTAAGGAACATATCGATGCTGGTATCACTCTGGCCGATGCAGTCAATTTTCTGGTGGAAAAGTAc                                                     >  3:322723/1‑134 (MQ=255)
                                                                                                            acacGCCTGCTGGAACAACACTATGGCCTGACACTGAACGACACGCCGTTCAGTGATGAAACTGTTATTAAGGAACATATCGATGCTGGTATCACTCTGGCCGATGCAGTCAATTTTCTGGTGGAAAAGTACGAACTGGTAc                                           <  1:124855/142‑1 (MQ=255)
                                                                                                            acacGCCTGCTGGAACAACACTATGGCCTGACACTGAACGACACGCCGTTCAGTGATGAAACTGTTATTAAGGAACATATCGATGCTGGTATCACTCTGGCCGATGCAGTCAATTTTCTGGTGGAAAAGTACGAACTGGTAc                                           >  2:124855/1‑142 (MQ=255)
                                                                                                              acGCCTGCTGGAACAACACTATGGCCTGACACTGAACGACACGCCGTTCAGTGATGAAACTGTTATTAAGGATCATATCGATGCTGGTATCACTCTGGCCGATGCAGTCAATTTTCTGGTGGaaaa                                                         >  4:80418/1‑126 (MQ=255)
                                                                                                                  ctgctgGAACAACACTATGGCCTGACACTGAACGACACGCCGTTCAGTGATGAAACTGTTATTAAGGAACATATCGATGCTGGTATCACTCTGGCCGATGCAGTCAATTTTCTGGTGGGAAAGTACGAAc                                                 <  1:767946/130‑1 (MQ=255)
                                                                                                                  ctgctgGAACAACACTATGGCCTGACACTGAACGACACGCCGTTCAGTGATGAAACTGTTATTAAGGAACATATCGATGCTGGTATCACTCTGGCCGATGCAGTCAATTTTCTGGTGGAAAAGTACGAAc                                                 >  2:767946/1‑130 (MQ=255)
                                                                                                                        gAACAACACTATGGCCTGACACTGAACGACACGCCGTTCAGTGATGAAACTGTTATTAAGGAACATATCGATGCTGGTATCACTCTGGCCGATGCAGTCAATTTTCtg                                                                 >  1:274924/1‑108 (MQ=255)
                                                                                                                           caacaCTATGGCCTGACACTGAACGACACGCCGTTCAGTGATGAAACTGTTATTAAGGAACATATCGATGCTGGTATCACTCTGGCCGATGCAGTCAATTTTCTGGTGGAAAAGTACGAACTGGTACGTATCGATCACAGAGGATTTTc                     >  1:1013905/1‑149 (MQ=255)
                                                                                                                           caacaCTATGGCCTGACACTGAACGACACGCCGTTCAGTGATGAAACTGTTATTAAGGAACATATCGATGCTGGTATCACTCTGGCCGATGCAGTCAATTTTCTGGTGGAAAAGTACGAACTGGTACGTATCGATCACAGAGGATTTTc                     <  2:1013905/149‑1 (MQ=255)
                                                                                                                            aacaCTATGGCCTGACACTGAACGACACGCCGTTCAGTGATGAAACTGTTATTAAGGAACATATCGATGCTGGTATCACTCTGGCCGATGCAGTCAATTTTCTGGTGGAAAAGTACGAACTGGTACGTATCGATCACAGAGGATTTTCg                    >  4:1082751/1‑149 (MQ=255)
                                                                                                                             acacTATGGCCTGACACTGAACGACACGCCGTTCAGTGATGAAACTGTTATTAAGGAACATATCGATGCTGGTATCACTCTGGCCGATGCAGTCAATTTTCTGGTGGAAAAGTACGAACTGGTACGTATCGATCACAGAGGATTTTCGt                   <  1:1654/149‑1 (MQ=255)
                                                                                                                              cacTATGGCCTGACACTGAACGACACGCCGTTCAGTGATGAAACTGTTATTAAGGAACATATCGATGCTGGTATCACTCTGGCCGATGCAGTCAATTTTCTGGTGGAAAAGTACGAACTGGTACGTATCGATCACAGAGGATTTTCGTg                  <  1:49928/149‑1 (MQ=255)
                                                                                                                               acTATGGCCTGACACTGAACGACACGCCGTTCAGTGATGAAACTGTTATTAAGGAACATATCGATGCTGGTATCACTCTGGCCGATGCAGTCAATTTTCTGGTGGAAAAGTACGAACTGGTa                                            >  3:278722/1‑122 (MQ=255)
                                                                                                                               acTATGGCCTGACACTGAACGACACGCCGTTCAGTGATGAAACTGTTATTAAGGAACATATCGATGCTGGTATCACTCTGGCCGATGCAGTCAATTTTCTGGTGGAAAAGTACGAACTGGTACGTATCg                                     >  1:978243/1‑129 (MQ=255)
                                                                                                                                cTATGGCCTGACACTGAACGACACGCCGTTCAGTGATGAAACTGTTATTAAGGAACATGTCGATGCTGGTATCACTCTGGCCGATGCAGTCAATTTTCTGGTGGAAAAGTACGAACTGGTACGTATCGATCACAGAGGATTTTCGTCGc                >  4:903547/1‑149 (MQ=255)
                                                                                                                                 tATGGCCTGACACTGAACGACACGCCGTTCAGTGATGAAACTGTTATTAAGGAACATATCGATGCTGGTATCACTCTGGCCGATGCAGTCAATTTTCTGGTGGaaaa                                                         <  1:1261613/107‑1 (MQ=255)
                                                                                                                                  aTGGCCTGACACTGAACGACACGCCGTTCAGTGATGAAACTGTTATTAAGGAACATATCGATGCTGGTATCACTCTGGCCGATGCAGTCAATTTTCTGGTGGAAAAGTACGAACTGGTACGTATCGATCACAGAGGATTTTCGTGGcaa              <  4:1118017/149‑1 (MQ=255)
                                                                                                                                    ggCCTGACACTGAACGACACGCCGTTCAGTGATGAAACTGTTATTAAGGAACATATCGATGCTGGTATCACTCTGGCCGATGCAGTCAATTTTCTGGTGGAAAAGTACGAACTGGTACGTATCGATCACAGAGGATTTTCGTGGcaa              >  2:721638/1‑147 (MQ=255)
                                                                                                                                    ggCCTGACACTGAACGACACGCCGTTCAGTGATGAAACTGTTATTAAGGAACATATCGATGCTGGTATCACTCTGGCCGATGCAGTCAATTTTCTGGTGGAAAAGTACGAACTGGTACGTATCGATCACAGAGGATTTTCGTGGcaa              <  1:721638/147‑1 (MQ=255)
                                                                                                                                      ccTGACACTGAACGACACGCCGTTCAGTGATGAAACTGTTATTAAGGAACATATCGATGCTGGTATCACTCTGGCCGATGCAGTCAATTTTCTGGTGGAAAAGTACGAACTGGTACGTATCGATCACAGAGGATTTTCGTGGCaacaac          <  1:925120/149‑1 (MQ=255)
                                                                                                                                        tGACACTGAACGACACGCCGTTCAGTGATGAAACTGTTATTAAGGAACATATCGATGCTGGTATCACTCTGGCCGATGCAGTCAATTTTCt                                                                  >  1:1042841/1‑91 (MQ=255)
                                                                                                                                        tGACACTGAACGACACGCCGTTCAGTGATGAAACTGTTATTAAGGAACATATCGATGCTGGTATCACTCTGGCCGATGCAGTCAATTTTCt                                                                  <  2:1042841/91‑1 (MQ=255)
                                                                                                                                        tGACACTGAACGACACGCCGTTCAGTGATGAAACTGTTATTAAGGAACATATCGATGCTGGTATCACTCTGGCCGATGCAGTCAATTTTCTGGCGGAAAAGTACGAACTGGTAc                                           >  3:1384235/1‑114 (MQ=255)
                                                                                                                                          acacTGAACGACACGCCGTTCAGTGATGAAACTGTTATTAAGGAACATATCGATGCTGGTATCACTCTGGCCGATGCAGTCAATTTTCTGGTGGAAAAGTACGAACTgg                                              >  1:1453343/1‑109 (MQ=255)
                                                                                                                                           cacTGAACGGCGTGTCGTTCAGTGATGAAACTGTTATTAAGGAACATATCGATGCTGGTATCACTCTGGCCGATGCAGTCAATTTTCTGGTGGAAAAGTACGAACTGGTACGTATCGATCACAGAGGa                          >  3:57301/1‑128 (MQ=255)
                                                                                                                                            acTGAACGGCGTGTCGTTCAGTGATGAAACTGTTATTAAGGAACATATCGATGCTGGTATCACTCTGGCCGATGCAGTCAATTTTCTGGTGGAAAAGTACGAACTGGTACGTATCGATCACAGAGGATTTTCGTGGCAACAACAGtctc    >  2:1298743/1‑149 (MQ=255)
                                                                                                                                            acTGAACGACACGCCGTTCAGTGATGAAACTGTTATTAAGGAACATATCGATGCTGGTATCACTCTGGCCGATGCAGTCAATTTTCTGGTGGAAAAGTACGAACTGGTACGTAt                                       <  3:874000/114‑1 (MQ=255)
                                                                                                                                            acTGAACGACACGCCGTTCAGTGATGAAACTGTTATTAAGGAACATATCGATGCTGGTATCACTCTGGCCGATGCAGTCAATTTTCTGGTGGAAAAGTACGAACTGGTACGTAt                                       >  4:874000/1‑114 (MQ=255)
                                                                                                                                            acTGAACGACACGCCGTTCAGTGATGAAACTGTTATTAAGGAACATATCGATGCTGGTATCACTCTGGCCGATGCAGTCAATTTTCTGGTGGAAAAGTACGAACTGGTACGTATCGa                                    <  4:243377/117‑1 (MQ=255)
                                                                                                                                            acTGAACGACACGCCGTTCAGTGATGAAACTGTTATTAAGGAACATATCGATGCTGGTATCACTCTGGCCGATGCAGTCAATTTTCTGGTGGAAAAGTACGAACTGGTACGTATCGa                                    >  3:243377/1‑117 (MQ=255)
                                                                                                                                            acTGAACGACACGCCGTTCAGTGATGAAACTGTTATTAAGGAACATATCGATGCTGGTATCACTCTGGCCGATGCAGTCAATTTTCTGGTGGAAAAGTACGAACTGGTACGTATCGATCACAGAGGATTTTCGTGGCAACAACAGtctc    <  4:1400834/149‑1 (MQ=255)
                                                                                                                                             cTGAACGGCGTGTCGTTCAGTGATGAAACTGTTATTAAGGAACATATCGATGCTGGTATCACTCTGGCCGATGCAGTCAATTTTCTGGTGGAAAAGTACGAACTGGTACGTATCGATCACAGAGGATTTTCGTGGcaa              <  4:1294396/138‑1 (MQ=255)
                                                                                                                                             cTGAACGACACGCCGTTCAGTGATGAAACTGTTATTAAGGAACATATCGATGCTGGTATCACTCTGGCCGATGCAGTCAATTTTCTGGTGGAAAAGTACGAACTGGTACGTATCGATCACAGAGGATTTTCGTGGCAACAACAGTCTcc   <  3:1203080/149‑1 (MQ=255)
                                                                                                                                              tGAACGACACGCCGTTCAGTGATGAAACTGTTATTAAGGAACATATCGATGCTGGTATCACTCTGGCCGATGCAGTCAATTTTCTGGTGGAAAAGTACGAACTgg                                              >  1:340017/1‑105 (MQ=255)
                                                                                                                                              tGAACGACACGCCGTTCAGTGATGAAACTGTTATTAAGGAACATATCGATGCTGGTATCACTCTGGCCGATGCAGTCAATTTTCTGGTGGAAAAGTACGAACTGGTACGTATCGATCACAGAGGATTTTCGTGGCAACAACAGTCTCCa  <  3:713342/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                              
TCAACGCTAACTAATTTCACGAGAGCAAGCATGAACACTCTACCTGCTACAATTTCGCAGGCGGCGAAGCCCTGCCTGTCGCCAGTGGCTGTCTGGCAAATGCTACTGACACGCCTGCTGGAACAACACTATGGCCTGACACTGAACGACACGCCGTTCAGTGATGAAACTGTTATTAAGGAACATATCGATGCTGGTATCACTCTGGCCGATGCAGTCAATTTTCTGGTGGAAAAGTACGAACTGGTACGTATCGATCACAGAGGATTTTCGTGGCAACAACAGTCTCCA  >  NC_000913/2777423‑2777713

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